Detailed Description

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

G+Smo/Eigen iterative solver for linear systems inside the incompressible flow solvers (classes derived from gsFlowSolverBase).

Template Parameters

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

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

Collaboration diagram for gsFlowLinSystSolver_iter< 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_iter (const gsFlowSolverParams< T > &params)
	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

template<class T , int MatOrder, class SolverType >

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

virtual

Solve the linear system.

Parameters

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

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