17 #ifndef __deal2__shifted_matrix_h
18 #define __deal2__shifted_matrix_h
21 #include <deal.II/base/config.h>
22 #include <deal.II/base/smartpointer.h>
39 template<
class MATRIX>
52 void shift (
const double sigma);
57 double shift ()
const;
62 template <
class VECTOR>
63 void vmult (VECTOR &dst,
const VECTOR &src)
const;
68 template <
class VECTOR>
69 double residual (VECTOR &dst,
const VECTOR &src,
const VECTOR &rhs)
const;
98 template<
class MATRIX,
class MASSMATRIX,
class VECTOR>
113 void shift (
const double sigma);
118 double shift ()
const;
123 void vmult (VECTOR &dst,
const VECTOR &src)
const;
128 double residual (VECTOR &dst,
const VECTOR &src,
const VECTOR &rhs)
const;
155 template <
class MATRIX>
164 template <
class MATRIX>
172 template <
class MATRIX>
181 template <
class MATRIX>
182 template <
class VECTOR>
192 template <
class MATRIX>
193 template <
class VECTOR>
197 const VECTOR &rhs)
const
202 dst.sadd(-1.,1.,rhs);
203 return dst.l2_norm ();
208 template <
class MATRIX,
class MASSMATRIX,
class VECTOR>
215 A(&A), M(&M), sigma(sigma)
219 template <
class MATRIX,
class MASSMATRIX,
class VECTOR>
226 template <
class MATRIX,
class MASSMATRIX,
class VECTOR>
234 template <
class MATRIX,
class MASSMATRIX,
class VECTOR>
237 const VECTOR &src)
const
249 template <
class MATRIX,
class MASSMATRIX,
class VECTOR>
253 const VECTOR &rhs)
const
262 dst.sadd(-1.,1.,rhs);
263 return dst.l2_norm ();
267 DEAL_II_NAMESPACE_CLOSE
SmartPointer< const MATRIX, ShiftedMatrixGeneralized< MATRIX, MASSMATRIX, VECTOR > > A
ShiftedMatrixGeneralized(const MATRIX &A, const MASSMATRIX &M, const double sigma)
void vmult(VECTOR &dst, const VECTOR &src) const
double residual(VECTOR &dst, const VECTOR &src, const VECTOR &rhs) const
SmartPointer< const MATRIX, ShiftedMatrix< MATRIX > > A
ShiftedMatrix(const MATRIX &A, const double sigma)
void vmult(VECTOR &dst, const VECTOR &src) const
double residual(VECTOR &dst, const VECTOR &src, const VECTOR &rhs) const
SmartPointer< const MASSMATRIX, ShiftedMatrixGeneralized< MATRIX, MASSMATRIX, VECTOR > > M