/usr/include/dolfin/la/uBLASFactory.h is in libdolfin-dev 1.4.0+dfsg-4.
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 | // Copyright (C) 2007 Ola Skavhaug
//
// 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 2008
// Modified by Anders Logg 2011-2012
//
// First added: 2007-12-06
// Last changed: 2012-08-20
#ifndef __UBLAS_FACTORY_H
#define __UBLAS_FACTORY_H
#include <string>
#include <memory>
#include "uBLASKrylovSolver.h"
#include "uBLASMatrix.h"
#include "uBLASVector.h"
#include "TensorLayout.h"
#include "UmfpackLUSolver.h"
#include "GenericLinearAlgebraFactory.h"
namespace dolfin
{
// Forward declaration
class GenericLinearSolver;
template<typename Mat = ublas_sparse_matrix>
class uBLASFactory : public GenericLinearAlgebraFactory
{
public:
/// Destructor
virtual ~uBLASFactory() {}
/// Create empty matrix
std::shared_ptr<GenericMatrix> create_matrix() const
{
std::shared_ptr<GenericMatrix> A(new uBLASMatrix<Mat>);
return A;
}
/// Create empty vector
std::shared_ptr<GenericVector> create_vector() const
{
std::shared_ptr<GenericVector> x(new uBLASVector);
return x;
}
/// Create empty tensor layout
std::shared_ptr<TensorLayout> create_layout(std::size_t rank) const
{
bool sparsity = false;
if (rank > 1)
sparsity = true;
std::shared_ptr<TensorLayout> pattern(new TensorLayout(0, sparsity));
return pattern;
}
/// Create empty linear operator
std::shared_ptr<GenericLinearOperator> create_linear_operator() const
{
std::shared_ptr<GenericLinearOperator> A(new uBLASLinearOperator);
return A;
}
/// Create LU solver
std::shared_ptr<GenericLUSolver> create_lu_solver(std::string method) const
{
std::shared_ptr<GenericLUSolver> solver(new UmfpackLUSolver);
return solver;
}
/// Create Krylov solver
std::shared_ptr<GenericLinearSolver> create_krylov_solver(std::string method,
std::string preconditioner) const
{
std::shared_ptr<GenericLinearSolver>
solver(new uBLASKrylovSolver(method, preconditioner));
return solver;
}
/// Return a list of available LU solver methods
std::vector<std::pair<std::string, std::string> >
lu_solver_methods() const
{
std::vector<std::pair<std::string, std::string> > methods;
methods.push_back(std::make_pair("default",
"default LU solver"));
methods.push_back(std::make_pair("umfpack",
"UMFPACK (Unsymmetric MultiFrontal sparse LU factorization)"));
return methods;
}
/// Return a list of available Krylov solver methods
std::vector<std::pair<std::string, std::string> >
krylov_solver_methods() const
{
return uBLASKrylovSolver::methods();
}
/// Return a list of available preconditioners
std::vector<std::pair<std::string, std::string> >
krylov_solver_preconditioners() const
{
return uBLASKrylovSolver::preconditioners();
}
/// Return singleton instance
static uBLASFactory<Mat>& instance()
{ return factory; }
private:
// Private Constructor
uBLASFactory() {}
// Singleton instance
static uBLASFactory<Mat> factory;
};
}
// Initialise static data
template<typename Mat> dolfin::uBLASFactory<Mat> dolfin::uBLASFactory<Mat>::factory;
#endif
|