matrixops.hpp

Go to the documentation of this file.

//==============================================================================
//==============================================================================
#ifndef MATRIXOPS_HPP_
#define MATRIXOPS_HPP_

#include <dune/common/fmatrix.hh> 
#include <dune/common/dynmatrix.hh>
#include <dune/istl/bcrsmatrix.hh>

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