SparseMatrix<Number, MemorySpace>::matrix_norm_square(
const Vector<Number, MemorySpace> &v) const
{
+ AssertDimension(m(), v.size());
Assert(matrix->isFillComplete(), ExcMatrixNotCompressed());
Assert(matrix->getRowMap()->isSameAs(*matrix->getDomainMap()),
ExcNotQuadratic());
const Vector<Number, MemorySpace> &u,
const Vector<Number, MemorySpace> &v) const
{
+ AssertDimension(m(), u.size());
+ AssertDimension(m(), v.size());
Assert(matrix->isFillComplete(), ExcMatrixNotCompressed());
Assert(matrix->getRowMap()->isSameAs(*matrix->getDomainMap()),
ExcNotQuadratic());
PetscScalar
MatrixBase::matrix_norm_square(const VectorBase &v) const
{
+ AssertDimension(m(), v.size());
+
VectorBase tmp(v);
vmult(tmp, v);
return tmp * v;
MatrixBase::matrix_scalar_product(const VectorBase &u,
const VectorBase &v) const
{
+ AssertDimension(m(), u.size());
+ AssertDimension(m(), v.size());
+
VectorBase tmp(u);
vmult(tmp, v);
return u * tmp;
TrilinosScalar
SparseMatrix::matrix_norm_square(const MPI::Vector &v) const
{
+ AssertDimension(m(), v.size());
Assert(matrix->RowMap().SameAs(matrix->DomainMap()), ExcNotQuadratic());
MPI::Vector temp_vector;
SparseMatrix::matrix_scalar_product(const MPI::Vector &u,
const MPI::Vector &v) const
{
+ AssertDimension(m(), u.size());
+ AssertDimension(m(), v.size());
Assert(matrix->RowMap().SameAs(matrix->DomainMap()), ExcNotQuadratic());
MPI::Vector temp_vector;