/usr/include/dolfin/la/ublas.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 | // Copyright (C) 2006-2011 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/>.
//
// First added: 2006-
// Last changed: 2011-12-04
#ifndef __DOLFIN_UBLAS_H
#define __DOLFIN_UBLAS_H
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/vector.hpp>
//#include <boost/numeric/ublas/vector_of_vector.hpp>
#include <boost/numeric/ublas/matrix_sparse.hpp>
// These two files must be included due to a bug in Boost version < 1.33.
#include <boost/numeric/ublas/vector_proxy.hpp>
#include <boost/numeric/ublas/triangular.hpp>
#include <boost/numeric/ublas/operation.hpp>
#include <boost/numeric/ublas/lu.hpp>
#include <boost/numeric/ublas/io.hpp>
namespace dolfin
{
/// Various typedefs for uBLAS data types
namespace ublas = boost::numeric::ublas;
// uBLAS vector
typedef ublas::vector<double> ublas_vector;
typedef ublas::vector_range<ublas_vector> ublas_vector_range;
// uBLAS dense matrix
typedef ublas::matrix<double> ublas_dense_matrix;
typedef ublas::matrix_range<ublas_dense_matrix> ublas_matrix_range;
// uBLAS dense matrix (column major format)
typedef ublas::matrix<double, ublas::column_major> ublas_matrix_cmajor;
typedef ublas::matrix_range<ublas_matrix_cmajor> ublas_matrix_cmajor_range;
// uBLAS sparse matrix
typedef ublas::compressed_matrix<double, ublas::row_major> ublas_sparse_matrix;
// uBLAS sparse matrix (column major format)
typedef ublas::compressed_matrix<double, ublas::column_major> ublas_sparse_matrix_cmajor;
// uBLAS sparse matrix for temporary assembly
typedef ublas::generalized_vector_of_vector<double, ublas::row_major,
ublas::vector<ublas::compressed_vector<double> > > ublas_assembly_matrix;
// uBLAS sparse matrix for temporary assembly (column major format)
typedef ublas::generalized_vector_of_vector<double, ublas::column_major,
ublas::vector<ublas::compressed_vector<double> > > ublas_assembly_matrix_cmajor;
// uBLAS upper triangular matrix (column major format)
typedef ublas::triangular_matrix<double, ublas::upper, ublas::column_major> ublas_matrix_cmajor_tri;
typedef ublas::matrix_range<ublas_matrix_cmajor_tri> ublas_matrix_cmajor_tri_range;
typedef ublas::matrix_column<ublas_matrix_cmajor_tri> ublas_matrix_cmajor_tri_column;
}
#endif
|