// Multiply by vector
NEWMAT::ReturnMatrix FullBFMatrix::MulByVec(const NEWMAT::ColumnVector& invec) const
{
if (invec.Nrows() != int(Ncols())) {throw BFMatrixException("FullBFMatrix::MulByVec: Matrix-vector size mismatch");}
NEWMAT::ColumnVector ret;
ret = (*mp)*invec;
ret.Release();
return(ret);
}
// Given A*x=b, solve for x
NEWMAT::ReturnMatrix FullBFMatrix::SolveForx(const NEWMAT::ColumnVector& b, // Ignoring all parameters except b
MISCMATHS::MatrixType type,
double tol,
int miter) const
{
if (int(Nrows()) != b.Nrows()) {throw BFMatrixException("FullBFMatrix::SolveForx: Matrix-vector size mismatch");}
NEWMAT::ColumnVector ret;
ret = mp->i()*b;
ret.Release();
return(ret);
}
