/usr/include/dolfin/la/PETScLUSolver.h is in libdolfin1.3-dev 1.3.0+dfsg-2.
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 | // Copyright (C) 2005-2011 Anders Logg
//
// 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 Garth N. Wells, 2009-2010.
//
// First added: 2005
// Last changed: 2011-10-19
#ifndef __DOLFIN_PETSC_LU_SOLVER_H
#define __DOLFIN_PETSC_LU_SOLVER_H
#ifdef HAS_PETSC
#include <map>
#include <boost/shared_ptr.hpp>
#include <petscksp.h>
#include <petscpc.h>
#include "GenericLUSolver.h"
#include "PETScObject.h"
namespace dolfin
{
/// Forward declarations
class GenericLinearOperator;
class GenericVector;
class PETScLinearOperator;
class PETScMatrix;
class PETScVector;
class PETScSNESSolver;
/// This class implements the direct solution (LU factorization) for
/// linear systems of the form Ax = b. It is a wrapper for the LU
/// solver of PETSc.
class PETScLUSolver : public GenericLUSolver, public PETScObject
{
public:
/// Constructor
PETScLUSolver(std::string method="default");
/// Constructor
PETScLUSolver(boost::shared_ptr<const PETScMatrix> A,
std::string method="default");
/// Destructor
~PETScLUSolver();
/// Set operator (matrix)
void set_operator(const boost::shared_ptr<const GenericLinearOperator> A);
/// Set operator (matrix)
void set_operator(const boost::shared_ptr<const PETScMatrix> A);
/// Get operator (matrix)
const GenericLinearOperator& get_operator() const;
/// Solve linear system Ax = b
std::size_t solve(GenericVector& x, const GenericVector& b);
/// Solve linear system Ax = b
std::size_t solve(GenericVector& x, const GenericVector& b,
bool transpose);
/// Solve linear system Ax = b
std::size_t solve(const GenericLinearOperator& A, GenericVector& x,
const GenericVector& b);
/// Solve linear system Ax = b
std::size_t solve(const PETScMatrix& A, PETScVector& x,
const PETScVector& b);
/// Solve linear system A^Tx = b
std::size_t solve_transpose(GenericVector& x, const GenericVector& b);
/// Solve linear system A^Tx = b
std::size_t solve_transpose(const GenericLinearOperator& A,
GenericVector& x, const GenericVector& b);
/// Solve linear system A^Tx = b
std::size_t solve_transpose(const PETScMatrix& 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();
/// Default parameter values
static Parameters default_parameters();
friend class PETScSNESSolver;
private:
const MatSolverPackage _solver_package;
// Available LU solvers
static const std::map<std::string, const MatSolverPackage> _methods;
// Whether those solvers support Cholesky
static const std::map<const MatSolverPackage, const bool> _methods_cholesky;
// Available LU solvers descriptions
static const std::vector<std::pair<std::string, std::string> >
_methods_descr;
// Select LU solver type
const MatSolverPackage select_solver(std::string& method) const;
// Does an LU solver support Cholesky?
bool solver_has_cholesky(const MatSolverPackage package) const;
// Initialise solver
void init_solver(std::string& method);
// Configure PETSc options
void configure_ksp(const MatSolverPackage solver_package);
// Set PETSc operators
void set_petsc_operators();
// Print pre-solve report
void pre_report(const PETScMatrix& A) const;
/// PETSc solver pointer
boost::shared_ptr<KSP> _ksp;
// Operator (the matrix)
boost::shared_ptr<const PETScMatrix> _A;
};
}
#endif
#endif
|