Matrix.h

Go to the documentation of this file.

// This file is part of the Peano project. For conditions of distribution and
// use, please see the copyright notice at www.peano-framework.org
#pragma once
 
 
namespace tarch {
  namespace la {
    template<int Rows, int Cols, typename Scalar>
    class Matrix;
 
    template<int Rows, int Cols, typename Scalar>
    Matrix<Rows,Cols,Scalar> buildMatrix(const Scalar* values);
  }
}
 
 
#include "tarch/la/MatrixOperations.h"
#include "tarch/la/MatrixVectorOperations.h"
#include "tarch/la/MatrixMatrixOperations.h"
#include "tarch/la/MatrixScalarOperations.h"
 
#include "tarch/compiler/CompilerSpecificSettings.h"
 
 
#include <string>
#include <initializer_list>
 
 
 
template<int Rows, int Cols, typename Scalar>
class tarch::la::Matrix {
  private:
    Scalar _values[Rows*Cols];
 
    friend Matrix<Rows,Cols,Scalar> buildMatrix<Rows,Cols,Scalar>(const Scalar* values);
 
  public:
    Matrix();
    Matrix( const Scalar&  value);
    Matrix( std::initializer_list<Scalar> values );
 
    int rows() const;
 
    int cols() const;
 
    int size() const;
 
    inline Scalar & operator() (
      int rowIndex,
      int colIndex
    )
      #ifdef UseManualInlining
      __attribute__((always_inline))
      #endif
      {
        assertion5( rowIndex >= 0, Rows, Cols, rowIndex, colIndex, toString() );
        assertion5( colIndex >= 0, Rows, Cols, rowIndex, colIndex, toString() );
        assertion5( rowIndex < Rows, Rows, Cols, rowIndex, colIndex, toString() );
        assertion5( colIndex < Cols, Rows, Cols, rowIndex, colIndex, toString() );
        return _values[rowIndex * Cols + colIndex];
      }
 
    inline const Scalar & operator() (
      int rowIndex,
      int colIndex
    ) const
      #ifdef UseManualInlining
      __attribute__((always_inline))
      #endif
      {
        assertion5( rowIndex >= 0, Rows, Cols, rowIndex, colIndex, toString() );
        assertion5( colIndex >= 0, Rows, Cols, rowIndex, colIndex, toString() );
        assertion5( rowIndex < Rows, Rows, Cols, rowIndex, colIndex, toString() );
        assertion5( colIndex < Cols, Rows, Cols, rowIndex, colIndex, toString() );
        return _values[rowIndex * Cols + colIndex];
      }
 
    std::string toString() const;
 
    std::string toPrettyString(int numberOfDigits=4) const;
 
    template <typename NewScalarType>
    tarch::la::Matrix<Rows,Cols,NewScalarType> convertScalar() const;
 
    Scalar* data() {
      return _values;
    }
 
    const Scalar *  data() const {
      return _values;
    }
};
 
 
 
#include "tarch/la/Matrix.cpph"