public IVariable Subtract(IVariable x, GekkoTime t)
{
//clone this Add() method and do with O.SubtractMatrixScalar, also for Matrix and matrix
EVariableType type = x.Type();
if (type == EVariableType.Matrix)
{
double[,] a = this.data;
double[,] b = ((Matrix)x).data;
int m = a.GetLength(0);
int k = a.GetLength(1);
if (b.GetLength(0) != m || b.GetLength(1) != k)
{
G.Writeln2("*** ERROR: The two matrices are not compatible for subtraction");
G.Writeln2(" " + m + "x" + k + " and " + b.GetLength(0) + "x" + b.GetLength(1) + " do not match");
throw new GekkoException();
}
double[,] c = O.SubtractMatrixMatrix(a, b, m, k);
Matrix z = new Matrix();
z.data = c;
return(z);
}
else
{
//subtraction of a scalar is not legal, this is like AREMOS
G.Writeln2("*** ERROR: You are trying to add a MATRIX and a " + type.ToString().ToUpper());
throw new GekkoException();
}
}
/// <summary>
/// 解析變數。
/// </summary>
/// <param name="value">解析值。</param>
/// <remarks>
/// 只對單一欄位解析，如:[@@Session]、[@FieldName]、[@Table.FieldName] 與常數。
/// </remarks>
public void Parse(string value)
{
Regex oRegex;
MatchCollection oMatches;
GroupCollection oGroups;
string sPattern;
string sWord;
string[] oArray;
sPattern = @"\[(?<word>@+[\w]*.?[\w]*)\]";
oRegex = new Regex(sPattern, RegexOptions.IgnoreCase);
oMatches = oRegex.Matches(value);
foreach (Match m in oMatches)
{
oGroups = m.Groups;
sWord = oGroups["word"].Value;
if (sWord.Contains("@@"))
{
// 系統變數
_Value = StrFunc.StrSubstring(sWord, 2);
_VariableType = EVariableType.Variable;
}
else if (sWord.Contains("@"))
{
// 欄位變數
_Value = StrFunc.StrSubstring(sWord, 1);
if (sWord.Contains("."))
{
// 資料欄位，格式為 [@Table.FieldName]，但其中 [@DB.NULL] 或 [@DB.NOTNULL] 為常數
oArray = StrFunc.StrSplit(this.Value, ".");
if (oArray.Length > 0)
{
if (StrFunc.SameText(StrFunc.StrUpper(oArray[0]), "DB"))
{
_VariableType = EVariableType.Constant; // DB.NULL 或 DB.NOTNULL 為常數
}
else
{
_VariableType = EVariableType.TableField;
}
}
}
else
{
// 欄位變數，格式為 [@FieldName]
_VariableType = EVariableType.Field; // 格式為 [@Table.FieldName]
}
}
return;
}
_Value = value;
_VariableType = EVariableType.Constant;
}
public EVariableDef(EVariableType t)
{
type = t;
}
public IVariable Multiply(IVariable x, GekkoTime t)
{
EVariableType type = x.Type();
if (type == EVariableType.Matrix)
{
double[,] a = this.data;
double[,] b = ((Matrix)x).data;
// m x k * p x n
int m = a.GetLength(0);
int k = a.GetLength(1);
int p = b.GetLength(0);
int n = b.GetLength(1);
if (p == 1 && n == 1)
{
//Special case, #a * #onebyone
double[,] c = O.MultiplyMatrixScalar(a, b[0, 0], m, k);
Matrix z = new Matrix();
z.data = c;
return(z);
}
else if (m == 1 && k == 1)
{
//Special case, #onebyone * #a
double[,] c = O.MultiplyMatrixScalar(b, a[0, 0], p, n);
Matrix z = new Matrix();
z.data = c;
return(z);
}
else
{
if (k != p)
{
G.Writeln2("*** ERROR: The two matrices are not compatible for multiplication");
G.Writeln2(" " + m + "x" + k + " and " + b.GetLength(0) + "x" + b.GetLength(1) + " do not match");
}
double[,] c = null;
if (false)
{
c = new double[m, n];
alglib.rmatrixgemm(m, n, k, 1, a, 0, 0, 0, b, 0, 0, 0, 0, ref c, 0, 0);
}
else
{
c = Program.MultiplyMatrices(a, b);
}
Matrix z = new Matrix();
z.data = c;
return(z);
//c = {{74, 80, 86, 92}, {173, 188, 203, 218}}
}
}
else if (type == EVariableType.Val)
{
//This is allowed in AREMOS, too
double[,] a = this.data;
double b = O.GetVal(x, t);
int m = a.GetLength(0);
int k = a.GetLength(1);
double[,] c = O.MultiplyMatrixScalar(a, b, m, k);
Matrix z = new Matrix();
z.data = c;
return(z);
}
else
{
G.Writeln2("*** ERROR: You are trying to multiply a MATRIX and a " + type.ToString().ToUpper());
throw new GekkoException();
}
}
