Calculate the cross between two matrixs as row wise or column wise
int Cross( matrixbase & mbSource, int nDim = 1 )
Returns 0 on success and a non-zero error code on failure.
EX1
void matrixbase_Cross_ex1() { int rc; matrix<double> mat1 = { {1, 2, 3} }; matrix<double> mat2 = { {4, 5, 6} }; rc = mat1.Cross(mat2); // Multiply mat2 to mat1 (row-wise) if(rc!=0) printf("Error: Cross failed. rc=%d\n"); else{ printf("The matrix is:\n"); for(int ii=0; ii< mat1.GetNumRows(); ii++){ for(int jj=0; jj< mat1.GetNumCols(); jj++) printf("%g ", mat1[ii][jj]); printf("\r\n"); } } }
EX2
// cross product of two (row) vectors void matrixbase_Cross_ex2() { int rc; matrix<double> mat1 = { {1, 2, 3} }; matrix<double> mat2 = { {4, 5, 6} }; // mat1 X mat2 (row-wise) = // {2*6-3*5, 3*4-1*6, 1*5-2*4} // {-3, 6, -3} // MatrixPage MatPg1; MatPg1.Create("Origin"); MatrixLayer MatLy1 = MatPg1.Layers(0); Matrix Mat1(MatLy1); Mat1 = mat1; MatrixPage MatPg2; MatPg2.Create("Origin"); MatrixLayer MatLy2 = MatPg2.Layers(0); Matrix Mat2(MatLy2); Mat2 = mat1; // Mat2 contents are initially identical to the original mat1 MatrixPage MatPg3; MatPg3.Create("Origin"); MatrixLayer MatLy3 = MatPg3.Layers(0); Matrix Mat3(MatLy3); Mat2 = mat2; // This is used to cross-product on Mat2 rc=mat1.Cross(mat2); // Multiply mat2 to mat1 (row-wise) if(rc!=0) printf(" Error: Cross on %s by %s failed. rc=%d\n", Mat1.GetName(),Mat2.GetName(),rc); else { Mat3 = mat1; printf(" %s X %s = %s (row wise product)\n", Mat1.GetName(), Mat2.GetName(),Mat3.GetName()); } }
EX3
// cross product of two matrices (column-wise product) void matrixbase_Cross_ex3() { int rc; matrix<double> mat1 = { {1, 1, 2, 3}, {2, 1, 2, 3}, {3, 1, 2, 3} }; matrix<double> mat2 = { {4, 4, 5, 6}, {5, 4, 5, 6}, {6, 4, 5, 6} }; // mat1 X mat2 (column-wise) = // {-3, 0 0 0} // { 6, 0 0 0} // {-3, 0 0 0} // // = transpose( {1 2 3}X(4 5 6) / {1 1 1}X{4 4 4} / {2 2 2}X{5 5 5} / {3 3 3}X{6 6 6} ) MatrixPage MatPg1; MatPg1.Create("Origin"); MatrixLayer MatLy1 = MatPg1.Layers(0); Matrix Mat1(MatLy1); Mat1 = mat1; MatrixPage MatPg2; MatPg2.Create("Origin"); MatrixLayer MatLy2 = MatPg2.Layers(0); Matrix Mat2(MatLy2); Mat2 = mat1; // Mat2 contents are initially identical to the original mat1 MatrixPage MatPg3; MatPg3.Create("Origin"); MatrixLayer MatLy3 = MatPg3.Layers(0); Matrix Mat3(MatLy3); Mat2 = mat2; // This is used to cross-product on Mat2 rc=mat1.Cross(mat2, 2); // Multiply mat2 to mat1 column-wise if(rc!=0) printf(" Error: Cross on %s by %s failed. rc=%d\n", Mat1.GetName(),Mat2.GetName(),rc); else { Mat3 = mat1; printf(" %s X %s = %s (column wise product)\n", Mat1.GetName(), Mat2.GetName(),Mat3.GetName()); } }
EX4
#define MAKEMAT(X) MatrixLayer MatLy##X;MatLy##X.Create();Matrix Mat##X(MatLy##X) void matrixbase_Cross_ex4() { matrix mat1 = { {1, 2, 3}, {1, 2, 3}, {1, 2, 3} }; matrix mat2 = { {2, 3, 4}, {3, 4, 5}, {4, 5, 6} }; MAKEMAT(1); MAKEMAT(2); MAKEMAT(3); MAKEMAT(4); MAKEMAT(5); MAKEMAT(6); MAKEMAT(7); Mat1=mat1; Mat2=mat2; Mat3=mat1*mat2; printf("%s = %s * %s <== proper multiplication of 2 matrices\n",Mat3.GetName(),Mat1.GetName(),Mat2.GetName()); Mat4=mat1*10; printf("%s = 10 * %s <== scalr multiplication of a matrix\n",Mat4.GetName(),Mat1.GetName()); Mat5=mat1; Mat5.DotMultiply(mat2); printf("%s = DotMultiply(%s,%s) <== element-wise multiplication of 2 matrices\n",Mat5.GetName(),Mat1.GetName(),Mat2.GetName()); Mat6=mat1; Mat6.Cross(mat2); printf("%s = Cross(%s,%s) <== cross product of 2 matrices\n",Mat6.GetName(),Mat1.GetName(),Mat2.GetName()); mat1.CumulativeProduct(Mat7); printf("%s = CumulativeProduct(%s) <== cumulative product of a matrix\n",Mat7.GetName(),Mat1.GetName()); }
Cross product of this matrix and the source matrix along the dimension identified by nDim. The nDim dimension of both matrices must be three and the underlying base type of the source matrix must be less than or equal to the underlying base type of this matrix.
Note
The following is a summary of various matrix operations:
1) Multiplication:
1-1) * ... Proper multiplication of 2 matrices (mat1*mat2)
1-2) * ... Scaler multiplication of a matrix (mat1*A or A*mat1)
1-3) DotMultiply ... Element-wise multiplication of 2 matrices
1-4) Cross ... Cross product of 2 matrices
1-5) CumulativeProduct ... cumulative product of a matrix
2) Division:
2-1) / ... Not defined for 2 matrices. For proper division, multiply Inverse(mat1)
2-2) / ... Scaler division of a matrix (mat1/A)
2-3) DotDivide ... Element-wise division of 2 matrices
3) Addition:
3-1) + ... Element-wise addition of 2 matrices
3-2) + ... Scaler addition of a matrix (mat1+A or A+mat1)
3-3) SumColumns ... Summation of each column in a matrix
3-4) CumulativeSum ... Cumulative product of a matrix
4) Subtraction:
4-1) - ... Element-wise subtraction of 2 matrices
4-2) - ... Scaler subtraction of a matrix (mat1-A or A-mat1)
4-3) Difference ... Difference of adjacent 2 rows in a matrix (1st order differential)
5) Power:
5-1) ^ (or pow) ... Not defined as element-wise power of 2 matrices
5-2) ^ (or pow) ... Not defined as scaler power of a matrix
5-3) DotPower ... Element-wise power of 2 matrices
To depict the differences of above multiplications, try the EX3.
matrixbase::DotMultiply, matrixbase::DotDivide, matrixbase::DotPower
origin.h