/usr/include/dolfin/la/EigenFactory.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 | // Copyright (C) 2015 Chris Richardson and Garth 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/>.
//
// First added: 2015-02-01
#ifndef __EIGEN_FACTORY_H
#define __EIGEN_FACTORY_H
#include <map>
#include <memory>
#include <string>
#include <dolfin/common/MPI.h>
#include <dolfin/log/log.h>
#include "EigenKrylovSolver.h"
#include "EigenLUSolver.h"
#include "EigenMatrix.h"
#include "EigenVector.h"
#include "TensorLayout.h"
#include "GenericLinearAlgebraFactory.h"
namespace dolfin
{
// Forward declaration
class GenericLinearSolver;
/// Eigen linear algebra factory
class EigenFactory : public GenericLinearAlgebraFactory
{
public:
/// Destructor
virtual ~EigenFactory() {}
/// Create empty matrix
std::shared_ptr<GenericMatrix> create_matrix(MPI_Comm comm) const
{ return std::make_shared<EigenMatrix>(); }
/// Create empty vector
std::shared_ptr<GenericVector> create_vector(MPI_Comm comm) const
{ return std::make_shared<EigenVector>(comm); }
/// Create empty tensor layout
std::shared_ptr<TensorLayout> create_layout(MPI_Comm comm,
std::size_t rank) const
{
TensorLayout::Sparsity sparsity = TensorLayout::Sparsity::DENSE;
if (rank > 1)
sparsity = TensorLayout::Sparsity::SPARSE;
return std::make_shared<TensorLayout>(comm, 0, sparsity);
}
/// Create empty linear operator
std::shared_ptr<GenericLinearOperator> create_linear_operator(MPI_Comm comm) const
{
dolfin_not_implemented();
std::shared_ptr<GenericLinearOperator> A;
return A;
}
/// Create LU solver
std::shared_ptr<GenericLinearSolver>
create_lu_solver(MPI_Comm comm, std::string method) const
{
return std::make_shared<EigenLUSolver>(method);
}
/// Create Krylov solver
std::shared_ptr<GenericLinearSolver>
create_krylov_solver(MPI_Comm comm,
std::string method,
std::string preconditioner) const
{
std::shared_ptr<GenericLinearSolver>
solver(new EigenKrylovSolver(method, preconditioner));
return solver;
}
/// Return a list of available LU solver methods
std::map<std::string, std::string> lu_solver_methods() const
{ return EigenLUSolver::methods(); }
/// Return a list of available Krylov solver methods
std::map<std::string, std::string> krylov_solver_methods() const
{ return EigenKrylovSolver::methods(); }
/// Return a list of available preconditioners
std::map<std::string, std::string> krylov_solver_preconditioners() const
{ return EigenKrylovSolver::preconditioners(); }
/// Return singleton instance
static EigenFactory& instance()
{ return factory; }
private:
// Private Constructor
EigenFactory() {}
// Singleton instance
static EigenFactory factory;
};
}
#endif
|