public override QsValue SubtractOperation(QsValue vl)
{
QsValue value;
if (vl is QsReference)
{
value = ((QsReference)vl).ContentValue;
}
else
{
value = vl;
}
if (value is QsScalar)
{
var scalar = value as QsScalar;
QsTensor NewTensor = new QsTensor();
if (this.Order > 3)
{
foreach (var iTensor in this.InnerTensors)
{
NewTensor.SubtractOperation((QsTensor)iTensor.AddOperation(value));
}
}
else
{
for (int il = 0; il < this.MatrixLayers.Count; il++)
{
QsMatrix ResultMatrix = (QsMatrix)this.MatrixLayers[il].SubtractOperation(scalar);
NewTensor.AddMatrix(ResultMatrix);
}
}
return(NewTensor);
}
if (value is QsTensor)
{
var tensor = value as QsTensor;
if (this.Order != (tensor.Order))
{
throw new QsException("Adding two different ranked tensors are not supported");
}
if (this.Order > 3)
{
QsTensor NewTensor = new QsTensor();
for (int i = 0; i < this.InnerTensors.Count(); i++)
{
var iTensor = this.InnerTensors[i];
NewTensor.AddInnerTensor((QsTensor)
iTensor.SubtractOperation(tensor.InnerTensors[i]));
}
return(NewTensor);
}
else
{
if (tensor.MatrixLayers.Count == this.MatrixLayers.Count)
{
// the operation will succeed
if (tensor.FaceRowsCount == this.FaceRowsCount)
{
if (tensor.FaceColumnsCount == this.FaceColumnsCount)
{
QsTensor NewTensor = new QsTensor();
for (int il = 0; il < this.MatrixLayers.Count; il++)
{
QsMatrix ResultMatrix = (QsMatrix)this.MatrixLayers[il].SubtractOperation(tensor.MatrixLayers[il]);
NewTensor.AddMatrix(ResultMatrix);
}
return(NewTensor);
}
}
}
}
}
throw new QsException("Tensor Subtract Operation with " + value.GetType().Name + " Failed");
}
