11/*
2- Copyright (C) 2011 - 2024 by the authors of the ASPECT code.
2+ Copyright (C) 2011 - 2025 by the authors of the ASPECT code.
33
44 This file is part of ASPECT.
55
@@ -27,9 +27,8 @@ namespace aspect
2727 {
2828 /* *
2929 * Implement the block Schur preconditioner
30- *(A B^T
31- * 0 S)^{-1}.
32- **/
30+ * (A B^T; 0 S)^{-1}.
31+ */
3332 template <class AInvOperator , class SInvOperator , class BTOperator , class VectorType >
3433 class BlockSchurPreconditioner : public Subscriptor
3534 {
@@ -39,7 +38,7 @@ namespace aspect
3938 * @param A_inverse_operator Approximation of the inverse of the velocity block.
4039 * @param S_inverse_operator Approximation for the inverse Schur complement.
4140 * @param BToperator Operator for the B^T block of the Stokes system.
42- ** /
41+ */
4342 BlockSchurPreconditioner (
4443 const AInvOperator &A_inverse_operator,
4544 const SInvOperator &S_inverse_operator,
@@ -48,15 +47,14 @@ namespace aspect
4847 /* *
4948 * Matrix vector product with this preconditioner object.
5049 */
51- void vmult (VectorType &dst,
50+ void vmult (VectorType &dst,
5251 const VectorType &src) const ;
5352
5453 private:
5554 /* *
5655 * References to the various operators this preconditioner works with.
5756 */
5857
59-
6058 const AInvOperator &A_inverse_operator;
6159 const SInvOperator &S_inverse_operator;
6260 const BTOperator &BT_operator;
@@ -70,18 +68,17 @@ namespace aspect
7068 const SInvOperator &S_inverse_operator,
7169 const BTOperator &BT_operator)
7270 :
73-
7471 A_inverse_operator (A_inverse_operator),
7572 S_inverse_operator (S_inverse_operator),
76- BT_operator (BT_operator)
73+ BT_operator (BT_operator)
7774 {}
7875
7976
8077
8178 template <class AInvOperator , class SInvOperator , class BTOperator , class VectorType >
8279 void
83- BlockSchurPreconditioner<AInvOperator, SInvOperator, BTOperator, VectorType>::
84- vmult (VectorType &dst,
80+ BlockSchurPreconditioner<AInvOperator, SInvOperator, BTOperator, VectorType>::
81+ vmult (VectorType &dst,
8582 const VectorType &src) const
8683 {
8784 typename VectorType::BlockType utmp (src.block (0 ));
@@ -103,4 +100,5 @@ namespace aspect
103100 }
104101 }
105102}
103+
106104#endif
0 commit comments