/usr/include/dolfin/la/PETScKrylovSolver.h is in libdolfin1.0-dev 1.0.0-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 | // Copyright (C) 2004-2005 Johan Jansson
//
// This file is part of DOLFIN.
//
// DOLFIN is free software: you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// DOLFIN is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with DOLFIN. If not, see <http://www.gnu.org/licenses/>.
//
// Modified by Anders Logg, 2005-2009.
// Modified by Johan Hoffman, 2005.
// Modified by Andy R. Terrel, 2005.
// Modified by Garth N. Wells, 2005-2010.
//
// First added: 2005-12-02
// Last changed: 2011-10-19
#ifndef __DOLFIN_PETSC_KRYLOV_SOLVER_H
#define __DOLFIN_PETSC_KRYLOV_SOLVER_H
#ifdef HAS_PETSC
#include <map>
#include <petscksp.h>
#include <boost/scoped_ptr.hpp>
#include <boost/shared_ptr.hpp>
#include <dolfin/common/types.h>
#include "GenericLinearSolver.h"
#include "PETScObject.h"
namespace dolfin
{
/// Forward declarations
class GenericMatrix;
class GenericVector;
class PETScBaseMatrix;
class PETScMatrix;
class PETScVector;
class PETScPreconditioner;
class PETScUserPreconditioner;
/// This class implements Krylov methods for linear systems
/// of the form Ax = b. It is a wrapper for the Krylov solvers
/// of PETSc.
class PETScKrylovSolver : public GenericLinearSolver, public PETScObject
{
public:
/// Create Krylov solver for a particular method and names preconditioner
PETScKrylovSolver(std::string method = "default",
std::string preconditioner = "default");
/// Create Krylov solver for a particular method and PETScPreconditioner
PETScKrylovSolver(std::string method, PETScPreconditioner& preconditioner);
/// Create Krylov solver for a particular method and PETScPreconditioner
PETScKrylovSolver(std::string method, PETScUserPreconditioner& preconditioner);
/// Create solver from given PETSc KSP pointer
explicit PETScKrylovSolver(boost::shared_ptr<KSP> ksp);
/// Destructor
~PETScKrylovSolver();
/// Set operator (matrix)
void set_operator(const boost::shared_ptr<const GenericMatrix> A);
/// Set operator (matrix)
void set_operator(const boost::shared_ptr<const PETScBaseMatrix> A);
/// Set operator (matrix) and preconditioner matrix
void set_operators(const boost::shared_ptr<const GenericMatrix> A,
const boost::shared_ptr<const GenericMatrix> P);
/// Set operator (matrix) and preconditioner matrix
void set_operators(const boost::shared_ptr<const PETScBaseMatrix> A,
const boost::shared_ptr<const PETScBaseMatrix> P);
/// Get operator (matrix)
const PETScBaseMatrix& get_operator() const;
/// Solve linear system Ax = b and return number of iterations
uint solve(GenericVector& x, const GenericVector& b);
/// Solve linear system Ax = b and return number of iterations
uint solve(PETScVector& x, const PETScVector& b);
/// Solve linear system Ax = b and return number of iterations
uint solve(const GenericMatrix& A, GenericVector& x, const GenericVector& b);
/// Solve linear system Ax = b and return number of iterations
uint solve(const PETScBaseMatrix& A, PETScVector& x, const PETScVector& b);
/// Return informal string representation (pretty-print)
std::string str(bool verbose) const;
/// Return PETSc KSP pointer
boost::shared_ptr<KSP> ksp() const;
/// Return a list of available solver methods
static std::vector<std::pair<std::string, std::string> > methods();
/// Return a list of available preconditioners
static std::vector<std::pair<std::string, std::string> > preconditioners();
/// Default parameter values
static Parameters default_parameters();
private:
/// Initialize KSP solver
void init(const std::string& method);
// Set PETSc operators
void set_petsc_operators();
// Set options
void set_petsc_options();
/// Report the number of iterations
void write_report(int num_iterations, KSPConvergedReason reason);
void check_dimensions(const PETScBaseMatrix& A, const GenericVector& x,
const GenericVector& b) const;
// Available solvers
static const std::map<std::string, const KSPType> _methods;
/// DOLFIN-defined PETScUserPreconditioner
PETScUserPreconditioner* pc_dolfin;
/// PETSc solver pointer
boost::shared_ptr<KSP> _ksp;
/// Preconditioner
boost::shared_ptr<PETScPreconditioner> preconditioner;
/// Operator (the matrix)
boost::shared_ptr<const PETScBaseMatrix> A;
/// Matrix used to construct the preconditoner
boost::shared_ptr<const PETScBaseMatrix> P;
bool preconditioner_set;
};
}
#endif
#endif
|