.\"
.de Id
..
.de Sp
.if n .sp
.if t .sp 0.4
..
.TH solver_abtb 2rheolef "rheolef-6.5" "rheolef-6.5" "rheolef-6.5"
.\" label: /*Class:solver_abtb
.SH NAME
\fBsolver_abtb\fP -- direct or iterative solver iterface for mixed linear systems
.\" skip: @cindex mixed linear problem
.\" skip: @cindex conjugate gradien algorithm
.\" skip: @cindex finite element method
.\" skip: @cindex stabilized mixed finite element method
.\" skip: @cindex Stokes problem
.\" skip: @cindex incompresible elasticity
.SH SYNOPSIS

.\" begin_example
.Sp
.nf
    solver_abtb stokes     (a,b,mp);
    solver_abtb elasticity (a,b,c,mp);
.Sp
.fi
.\" end_example
.SH DESCRIPTION

The \fBsolver_abtb\fP class provides direct or iterative algorithms for some mixed problem:
.\" begin_example
.Sp
.nf
       [ A  B^T ] [ u ]    [ Mf ]
       [        ] [   ]  = [    ]
       [ B  -C  ] [ p ]    [ Mg ]
.Sp
.fi
.\" end_example
where A is symmetric positive definite and C is symmetric positive.
By default, iterative algorithms are considered for tridimensional problems
and direct methods otherwise.
Such mixed linear problems appears for instance with the discretization
of Stokes problems.
The C matrix can be zero and then the corresponding argument can be omitted
when invoking the constructor.
Non-zero C matrix appears for of Stokes problems with stabilized P1-P1 element,
or for nearly incompressible elasticity problems.
.SH DIRECT ALGORITHM

When the kernel of \fBB^T\fP is not reduced to zero, then the pressure p is defined up to a constant
and the system is singular. In the case of iterative methods, this is not a problem.
But when using direct method, the system is then completed to impose a constraint on the pressure term
and the whole matrix is factored one time for all.
.SH ITERATIVE ALGORITHM

The preconditionned conjugate gradient algorithm is used, where the \fBmp\fP matrix
is used as preconditionner. See see mixed_solver(4).
.SH EXAMPLES

See the user's manual for practical examples for the nearly incompressible
elasticity, the Stokes and the Navier-Stokes problems.
.\" skip start:AUTHOR:

.\" skip start:DATE:
.\" skip start:METHODS:
.\" END
.SH IMPLEMENTATION
.\" begin_example
.Sp
.nf
template <class T, class M = rheo_default_memory_model>
class solver_abtb_basic {
public:

// typedefs:

  typedef typename csr<T,M>::size_type size_type;

// allocators:

  solver_abtb_basic ();
  solver_abtb_basic (const csr<T,M>& a, const csr<T,M>& b, const csr<T,M>& mp,
       const solver_option_type& opt = solver_option_type());
  solver_abtb_basic (const csr<T,M>& a, const csr<T,M>& b, const csr<T,M>& c, const csr<T,M>& mp,
       const solver_option_type& opt = solver_option_type());

// accessors:

  void solve (const vec<T,M>& f, const vec<T,M>& g, vec<T,M>& u, vec<T,M>& p) const;

protected:
// internal
  void init();
// data:
  mutable solver_option_type _opt;
  csr<T,M>          _a;
  csr<T,M>          _b;
  csr<T,M>          _c;
  csr<T,M>          _mp;
  solver_basic<T,M> _sA;
  solver_basic<T,M> _sa;
  solver_basic<T,M> _smp;
  bool              _need_constraint;
};
typedef solver_abtb_basic<Float,rheo_default_memory_model> solver_abtb;
.Sp
.fi
.\" end_example

.\" LENGTH = 1
.SH SEE ALSO
mixed_solver(4)