petsc_linear_solver.h File Reference

Go to the source code of this file.

Classes
class	libMesh::PetscMatrixBase< T >
	This class provides a nice interface to the PETSc C-based data structures for parallel, sparse matrices. More...
class	libMesh::PetscLinearSolver< T >
	This class provides an interface to PETSc iterative solvers that is compatible with the libMesh LinearSolver<> More...

Namespaces
	libMesh
	The libMesh namespace provides an interface to certain functionality in the library.

Functions
PetscErrorCode	libmesh_petsc_preconditioner_setup (PC)
	This function is called by PETSc to initialize the preconditioner. More...
PetscErrorCode	libmesh_petsc_preconditioner_apply (PC, Vec x, Vec y)
	This function is called by PETSc to actually apply the preconditioner. More...
PetscErrorCode	__libmesh_petsc_preconditioner_setup (PC)
	This function is called by PETSc to initialize the preconditioner. More...
PetscErrorCode	__libmesh_petsc_preconditioner_apply (PC, Vec x, Vec y)
	This function is called by PETSc to actually apply the preconditioner. More...

Function Documentation

__libmesh_petsc_preconditioner_apply()

PetscErrorCode __libmesh_petsc_preconditioner_apply	(	PC	,
		Vec	x,
		Vec	y
	)

This function is called by PETSc to actually apply the preconditioner.

ctx will hold the Preconditioner.

Definition at line 93 of file petsc_linear_solver.C.

References libMesh::libmesh_petsc_preconditioner_apply(), and libMesh::PetscFunctionReturn().

  {
    PetscFunctionBegin;
    libmesh_deprecated();
    PetscFunctionReturn(libmesh_petsc_preconditioner_apply(pc, x, y));
  }

__libmesh_petsc_preconditioner_setup()

PetscErrorCode __libmesh_petsc_preconditioner_setup

(

PC

)

This function is called by PETSc to initialize the preconditioner.

ctx will hold the Preconditioner.

Definition at line 86 of file petsc_linear_solver.C.

References libMesh::libmesh_petsc_preconditioner_setup(), and libMesh::PetscFunctionReturn().

  {
    PetscFunctionBegin;
    libmesh_deprecated();
    PetscFunctionReturn(libmesh_petsc_preconditioner_setup(pc));
  }

libmesh_petsc_preconditioner_apply()

PetscErrorCode libmesh_petsc_preconditioner_apply	(	PC	,
		Vec	x,
		Vec	y
	)

This function is called by PETSc to actually apply the preconditioner.

ctx will hold the Preconditioner.

Definition at line 69 of file petsc_linear_solver.C.

References libMesh::Preconditioner< T >::apply(), libMesh::ParallelObject::comm(), libMesh::ctx, libMesh::LibmeshPetscCallQ(), and libMesh::PetscFunctionReturn().

Referenced by libMesh::__libmesh_petsc_preconditioner_apply(), libMesh::PetscNonlinearSolver< Number >::init(), and libMesh::PetscLinearSolver< Number >::init().

  {
    PetscFunctionBegin;

    void * ctx;
    LibmeshPetscCallQ(PCShellGetContext(pc,&ctx));
    Preconditioner<Number> * preconditioner = static_cast<Preconditioner<Number> *>(ctx);

    PetscVector<Number> x_vec(x, preconditioner->comm());
    PetscVector<Number> y_vec(y, preconditioner->comm());

    preconditioner->apply(x_vec,y_vec);

    PetscFunctionReturn(LIBMESH_PETSC_SUCCESS);
  }

libmesh_petsc_preconditioner_setup()

PetscErrorCode libmesh_petsc_preconditioner_setup

(

PC

)

This function is called by PETSc to initialize the preconditioner.

ctx will hold the Preconditioner.

Definition at line 53 of file petsc_linear_solver.C.

References libMesh::ctx, libMesh::Preconditioner< T >::initialized(), libMesh::LibmeshPetscCallQ(), libMesh::PetscFunctionReturn(), and libMesh::Preconditioner< T >::setup().

Referenced by libMesh::__libmesh_petsc_preconditioner_setup(), libMesh::PetscNonlinearSolver< Number >::init(), and libMesh::PetscLinearSolver< Number >::init().

  {
    PetscFunctionBegin;

    void * ctx;
    LibmeshPetscCallQ(PCShellGetContext(pc,&ctx));
    Preconditioner<Number> * preconditioner = static_cast<Preconditioner<Number> *>(ctx);

    libmesh_error_msg_if(!preconditioner->initialized(),
                         "Preconditioner not initialized!  Make sure you call init() before solve!");

    preconditioner->setup();

    PetscFunctionReturn(LIBMESH_PETSC_SUCCESS);
  }