matrixops.hpp

Go to the documentation of this file.

//==============================================================================
//==============================================================================
#ifndef MATRIXOPS_HPP_
#define MATRIXOPS_HPP_
 
#include <opm/common/utility/platform_dependent/disable_warnings.h>
 
#include <dune/common/fmatrix.hh> 
#include <dune/common/dynmatrix.hh>
#include <dune/istl/bcrsmatrix.hh>
 
#include <opm/common/utility/platform_dependent/reenable_warnings.h>
 
namespace Opm {
namespace Elasticity {
 
typedef Dune::BCRSMatrix<Dune::FieldMatrix<double,1,1> > Matrix;
 
typedef std::vector< std::set<int> > AdjacencyPattern;
 
typedef Dune::BlockVector<Dune::FieldVector<double,1> > Vector;
 
class MatrixOps {
  public:
    static void fromAdjacency(Matrix& A, const AdjacencyPattern& adj,
                              int rows, int cols);
 
    static Matrix fromDense(const Dune::DynamicMatrix<double>& T);
 
    static void print(const Matrix& A);
 
    static Matrix Axpy(const Matrix& A, const Matrix& B, double alpha);
 
    static Matrix augment(const Matrix& A, const Matrix& B,
                          size_t r0, size_t c0, bool symmetric);
 
    static Matrix extractDiagonal(const Matrix& A);
 
    static Matrix diagonal(size_t N);
 
    static Matrix extractBlock(const Matrix& A,
                               size_t r0, size_t N, size_t c0, size_t M);
 
    static void saveAsc(const Matrix& A, const std::string& file);
};
 
}
}
 
#endif