/usr/include/dolfin/la/PETScLUSolver.h is in libdolfin-dev 2017.2.0.post0-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 | // Copyright (C) 2005-2017 Anders Logg and Garth N. Wells
//
// 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/>.
#ifndef __DOLFIN_PETSC_LU_SOLVER_H
#define __DOLFIN_PETSC_LU_SOLVER_H
#ifdef HAS_PETSC
#include <map>
#include <memory>
#include <string>
#include <petscmat.h>
#include <petscpc.h>
#include <dolfin/common/MPI.h>
#include "PETScKrylovSolver.h"
namespace dolfin
{
/// Forward declarations
class GenericLinearOperator;
class GenericVector;
class PETScLinearOperator;
class PETScMatrix;
class PETScVector;
/// 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 GenericLinearSolver
{
public:
/// Constructor
PETScLUSolver(MPI_Comm comm, std::string method="default");
/// Constructor
PETScLUSolver(std::string method="default");
/// Constructor
PETScLUSolver(MPI_Comm comm,
std::shared_ptr<const PETScMatrix> A,
std::string method="default");
/// Constructor
PETScLUSolver(std::shared_ptr<const PETScMatrix> A,
std::string method="default");
/// Destructor
~PETScLUSolver();
/// Set operator (matrix)
void set_operator(std::shared_ptr<const GenericLinearOperator> A);
/// Set operator (matrix)
void set_operator(const PETScMatrix& A);
/// Solve linear system Ax = b
std::size_t solve(GenericVector& x, const GenericVector& b);
/// Solve linear system Ax = b (A^t x = b if transpose is true)
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);
/// Sets the prefix used by PETSc when searching the options
/// database
void set_options_prefix(std::string options_prefix);
/// Returns the prefix used by PETSc when searching the options
/// database
std::string get_options_prefix() const;
/// Set options from the PETSc options database
void set_from_options() const;
/// Returns the MPI communicator
MPI_Comm mpi_comm() const;
/// Return informal string representation (pretty-print)
std::string str(bool verbose) const;
/// Return PETSc KSP pointer
KSP ksp() const;
/// Return a list of available solver methods
static std::map<std::string, std::string> methods();
/// Default parameter values
static Parameters default_parameters();
/// Return parameter type: "krylov_solver" or "lu_solver"
std::string parameter_type() const
{ return "lu_solver"; }
// FIXME: These should not be friend classes
friend class PETScSNESSolver;
friend class PETScTAOSolver;
private:
// FIXME: Remove
// Available LU solvers
static std::map<std::string, const MatSolverPackage> lumethods;
// Select LU solver type
static const MatSolverPackage select_solver(MPI_Comm comm,
std::string method);
PETScKrylovSolver _solver;
};
}
#endif
#endif
|