RationalMatrix.h

Go to the documentation of this file.

/* This file is part of gbLAB.
 *
 * gbLAB is distributed without any warranty under the MIT License.
 */
 
 
#ifndef gbLAB_RationalMatrix_h_
#define gbLAB_RationalMatrix_h_
 
#include <Eigen/Dense>
namespace oILAB {
 
template <int dim> class RationalMatrix {
 
  typedef Eigen::Matrix<double, dim, 1> VectorDimD;
  typedef Eigen::Matrix<double, dim, dim> MatrixDimD;
  typedef long long int IntScalarType;
  typedef Eigen::Matrix<IntScalarType, dim, 1> VectorDimI;
  typedef Eigen::Matrix<IntScalarType, dim, dim> MatrixDimI;
 
  // static constexpr int64_t maxDen=10000000;
  static constexpr long long int maxDen = 1000000;
  static std::pair<MatrixDimI, IntScalarType> compute(const MatrixDimD &R);
  static std::pair<MatrixDimI, IntScalarType> reduce(const MatrixDimI &Rn,
                                                     const MatrixDimI &Rd);
 
  const std::pair<MatrixDimI, IntScalarType> returnPair;
 
public:
  const MatrixDimI &integerMatrix;
  const IntScalarType &mu;
 
  RationalMatrix(const MatrixDimD &R);
  RationalMatrix(const MatrixDimI &Rn, const IntScalarType &Rd);
  RationalMatrix(const MatrixDimI &Rn, const MatrixDimI &Rd);
  MatrixDimD asMatrix() const;
 
    };
 
    } // namespace oILAB
#endif