dotNxM()

TDmatrix linear::op::dotNxM	(	const TDmatrix::Tcoc &	,
		const TDmatrix::Tcoc &
	)

dot product for NxM matrices.

We cannot simply overload the operator*() function. The general libaff array underlying TDmatrix is an NxMxKxL matrix, even if used for an NxM matrix. In most cases this doesn't matter. But with matrix dot products we have to distinguish.

Definition at line 45 of file dot.cc.

References checkNxM(), and LINEAR_assert.

Referenced by main(), and test().

    {
      const int dim1=0;
      const int dim2=1;
      checkNxM(A);
      checkNxM(B);
      LINEAR_assert(A.first(dim2)==B.first(dim1),
                    "size of input matrices does not match!");
      LINEAR_assert(A.last(dim2)==B.last(dim1),
                    "size of input matrices does not match!");
      TDmatrix C(aff::Shaper(A.f(dim1),A.l(dim1))(B.f(dim2),B.l(dim2)));
      for (int i=C.f(dim1); i<= C.l(dim1); ++i)
      {
        for (int j=C.f(dim2); j<= C.l(dim2); ++j)
        {
          C(i,j)=0.;
          for (int k=A.f(dim2); k<=A.l(dim2); ++k)
          {
            C(i,j)+=A(i,k)*B(k,j);
          }
        }
      }
      return(C);
    }

Here is the call graph for this function:

Here is the caller graph for this function: