gsFlowLinSystSolver_iterSP< T, MatOrder, SolverType > Class Template Reference

Detailed Description

template<class T, int MatOrder, class SolverType>
class gismo::gsFlowLinSystSolver_iterSP< T, MatOrder, SolverType >

G+Smo/Eigen iterative solver for saddle-point linear systems inside the incompressible flow solvers (classes derived from gsFlowSolverBase) with block preconditioners.

Template Parameters

T	coefficient type
MatOrder	sparse matrix storage order (ColMajor/RowMajor)
SolverType	the G+Smp/Eigen iterative solver type

Inheritance diagram for gsFlowLinSystSolver_iterSP< T, MatOrder, SolverType >:

Collaboration diagram for gsFlowLinSystSolver_iterSP< T, MatOrder, SolverType >:

Public Member Functions
virtual void	applySolver (const gsSparseMatrix< T, MatOrder > &mat, const gsMatrix< T > &rhs, gsMatrix< T > &solution)
	Solve the linear system.
virtual void	applySolver (const gsSparseMatrix< T, MatOrder > &mat, const gsMatrix< T > &rhs, gsMatrix< T > &solution, real_t alpha_u, real_t alpha_p, index_t usize, index_t pdofs)
	Solve the Navier–Stokes linear system with underrelaxation.
T	getAvgLinIterations () const
	Returns the average iteration count of the linear solver per applySolver() call.
std::vector< index_t >	getLinIterVector () const
	Returns vector of iteration counts of the linear solver for each call of applySolver().
virtual const T	getSolverSetupTime () const
	Returns the total time spent on linear solver setup.
virtual const T	getSolveTime () const
	Returns the total time spent on solving of the linear systems.
	gsFlowLinSystSolver_iterSP (const gsFlowSolverParams< T > &params, const gsINSAssembler< T, MatOrder > *assemblerPtr)
	Constructor.
virtual void	setupPreconditioner (const gsSparseMatrix< T, MatOrder > &mat)
	Setup the preconditioner for a given matrix.
virtual void	setupSolver (const gsSparseMatrix< T, MatOrder > &mat)
	Setup the linear solver for a given matrix (nothing happens here).

Protected Member Functions
T	stopwatchStart ()
	Start measuring time (decides whether to use gsStopwatch or MPI_Wtime)
T	stopwatchStop ()
	Stop measuring time (decides whether to use gsStopwatch or MPI_Wtime)

Member Function Documentation

applySolver() [1/2]

template<class T , int MatOrder, class SolverType >

void applySolver ( const gsSparseMatrix< T, MatOrder > &  mat, const gsMatrix< T > &  rhs, gsMatrix< T > &  solution  )

virtualinherited

Solve the linear system.

Parameters

[out]

solution

a reference to the vector, where the computed solution will be stored

Reimplemented from gsFlowLinSystSolver< T, MatOrder >.

applySolver() [2/2]

template<class T , int MatOrder>

void applySolver ( const gsSparseMatrix< T, MatOrder > &  mat, const gsMatrix< T > &  rhs, gsMatrix< T > &  solution, real_t  alpha_u, real_t  alpha_p, index_t  usize, index_t  pdofs  )

virtualinherited

Solve the Navier–Stokes linear system with underrelaxation.

Parameters

[in]	mat	a const reference to the system matrix
[in]	rhs	a const reference to the system right-hand side
[out]	solution	a reference to the vector, where the computed solution will be stored
[in]	alpha_u	velocity relaxation parameter
[in]	alpha_p	pressure relaxation parameter
[in]	usize	size of the velocity part of the system
[in]	pdofs	number of pressure DOFs