LibMesh::LinearSolver<T> & PetscLinearSolver<T> 类接口

This class provides a uniform interface for linear solvers, this base class is overloaded to provide linear solvers from different packages like PETSc / LASPACK

LinearSolver类继承关系：



共有成员函数：

LinearSovler(const libMesh::Communicator)

virtual ~LinearSolver()

clear(),

init(),

solver_type()

set_solver_type(),

Preconditioner_type(),

set_preconditioner_type(),

attach_preconditioner(),

1 {
2
if(this->_is_initialized)

3 libmesh_error_msg("Preconditioner must be attached before the s
olver is initialized");
4
5
_preconditioner_type = SHELL_PRECOND;
6
_preconditioner = preconditioner //function input;

7 }

virtual solver() 重载方法，一般给矩阵a，rhs列向量，求解列向量，pc矩阵可以给也可以有矩阵a计算。

static AutoPtr<LinearSolver<T> > build(), build a LinearSolver using linear solver package specified by solver_package

9 {

10 switch (solver_package)

11 {

12 #ifdef LIBMESH_HAVE_LASPACK

13
case LASPACK_SOLVERS:

14 {

15 AutoPtr<LinearSolver<T> > ap(new LaspackLinearSolver<T>(comm) );

16
return ap;

17 }

18 #endif

19

20 #ifdef LIBMESH_HAVE_PETSC

21
case PETSC_SOLVERS:

22 {

23 AutoPtr<LinearSolver<T> > ap(new PetscLinearSolver<T>(comm) );

24
return ap;

25 }

26 #endif

27

28 #ifdef LIBMESH_HAVE_TRILINOS

29
case TRILINOS_SOLVERS:

30 {

31 AutoPtr<LinearSolver<T> > ap(new AztecLinearSolver<T>(comm) );

32
return ap;

33 }

34 #endif

35

36 #ifdef LIBMESH_HAVE_EIGEN

37
case EIGEN_SOLVERS:

38 {

39 AutoPtr<LinearSolver<T> > ap(new EigenSparseLinearSolver<T>(comm) );

40
return ap;

41 }

42 #endif

43
default:

44 libmesh_error_msg("ERROR:Unrecognized solver package:" << solver_package);

45 }

46 AutoPtr<LinearSolver<T> > ap(NULL);

47 return ap;

48 }

PetscLinearSolver 重载 除构造／析构函数以外的 linear solver方法， 另外还有方法：

pc(),
ksp(), return PETSc ksp context pointer
get_residual_history()
get_initial_residual()
set_preconditioner_type(),