mult() [1/10]

template<typename RT>

void MatrixOps::mult	(	const DMat< RT > &	A,
		const DMat< RT > &	B,
		DMat< RT > &	result_product )

Definition at line 437 of file mat-linalg.hpp.

{
  // printf("entering dmat mult\n");
  assert(A.numColumns() == B.numRows());
  assert(A.numRows() == result_product.numRows());
  assert(B.numColumns() == result_product.numColumns());
 
  typename RT::Element tmp(A.ring());
  for (size_t i = 0; i < A.numRows(); i++)
    for (size_t j = 0; j < B.numColumns(); j++)
      {
        auto& val = result_product.entry(i,j);
        for (size_t k = 0; k < A.numColumns(); ++k)
          {
            A.ring().mult(tmp, A.entry(i,k), B.entry(k,j));
            A.ring().add(val, val, tmp);
          }
      }
}

References DMat< ACoeffRing >::entry(), DMat< ACoeffRing >::numColumns(), DMat< ACoeffRing >::numRows(), and DMat< ACoeffRing >::ring().