subtractMultipleTo() [1/10]

template<typename RT>

void MatrixOps::subtractMultipleTo	(	DMat< RT > &	C,
		const DMat< RT > &	A,
		const DMat< RT > &	B )

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

{
  assert(A.numColumns() == B.numRows());
  assert(A.numRows() == C.numRows());
  assert(B.numColumns() == C.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 = C.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().subtract(val, val, tmp);
          }
      }
}

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