/usr/include/dolfin/la/uBLASFactory.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 | // 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
//
// First added: 2007-12-06
// Last changed: 2011-10-19
#ifndef __UBLAS_FACTORY_H
#define __UBLAS_FACTORY_H
#include <string>
#include "uBLASKrylovSolver.h"
#include "uBLASMatrix.h"
#include "uBLASVector.h"
#include "SparsityPattern.h"
#include "UmfpackLUSolver.h"
#include "LinearAlgebraFactory.h"
namespace dolfin
{
// Forward declaration
class GenericLinearSolver;
template<typename Mat = ublas_sparse_matrix>
class uBLASFactory : public LinearAlgebraFactory
{
public:
/// Destructor
virtual ~uBLASFactory() {}
/// Create empty matrix
uBLASMatrix<Mat>* create_matrix() const
{ return new uBLASMatrix<Mat>(); }
/// Create empty vector
uBLASVector* create_vector() const
{ return new uBLASVector(); }
/// Create empty vector (local)
uBLASVector* create_local_vector() const
{ return new uBLASVector(); }
/// Create empty sparsity pattern
SparsityPattern* create_pattern() const
{ return new SparsityPattern(); }
/// Create LU solver
UmfpackLUSolver* create_lu_solver(std::string method) const
{ return new UmfpackLUSolver(); }
/// Create Krylov solver
GenericLinearSolver* create_krylov_solver(std::string method,
std::string preconditioner) const
{ return new uBLASKrylovSolver(method, preconditioner); }
/// 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
|