[sharonyue]

Hi guys,

I dont know how should I call this inner product in OpenFOAM syntax. you know that if I wanna operates inner product on a tensor(T) and a vector(v). I can use:

T & v

but this time I wanna operates inner product on a n×n matrix with a n×1 matrix.

First I define this n×n matrix as a scalarSquareMatrix and define this n×1 matrix as scalarDiagonalMatrix.

But I dont know which sign can be used as inner product on this two matrices?



BTW, I use for loop to calculate but failed:

for (int c = 0; c < nWeights_*2; c++)
    {
	for ( int r = 0; r < nWeights_*2; r++)  
	{
	      mT[c] = invA[c][r] * moments_[r];
	}
    }

Only mT[3] is rite. The others of it is null.