GlobalMatrixVectorTypes.h

Go to the documentation of this file.

 
#pragma once
 
#include "SparsityPattern.h"
 
//
// Global vector/matrix types and linear solver.
//
#if defined(USE_LIS)
 
#include "MathLib/LinAlg/Eigen/EigenMatrix.h"
#include "MathLib/LinAlg/Eigen/EigenVector.h"
#include "MathLib/LinAlg/EigenLis/EigenLisLinearSolver.h"
 
using GlobalVector = MathLib::EigenVector;
using GlobalMatrix = MathLib::EigenMatrix;
 
using GlobalLinearSolver = MathLib::EigenLisLinearSolver;
 
#elif defined(USE_PETSC)
#include "MathLib/LinAlg/PETSc/PETScLinearSolver.h"
#include "MathLib/LinAlg/PETSc/PETScMatrix.h"
#include "MathLib/LinAlg/PETSc/PETScVector.h"
 
using GlobalVector = MathLib::PETScVector;
using GlobalMatrix = MathLib::PETScMatrix;
 
using GlobalLinearSolver = MathLib::PETScLinearSolver;
 
#else
#include "MathLib/LinAlg/Eigen/EigenLinearSolver.h"
#include "MathLib/LinAlg/Eigen/EigenMatrix.h"
#include "MathLib/LinAlg/Eigen/EigenVector.h"
 
using GlobalVector = MathLib::EigenVector;
using GlobalMatrix = MathLib::EigenMatrix;
 
using GlobalLinearSolver = MathLib::EigenLinearSolver;
 
#endif
 
static_assert(std::is_integral_v<GlobalMatrix::IndexType>,
              "The index type for global matrices is not an integral type.");
static_assert(std::is_integral_v<GlobalVector::IndexType>,
              "The index type for global vectors is not an integral type.");
static_assert(std::is_same_v<GlobalMatrix::IndexType, GlobalVector::IndexType>,
              "The global matrix and vector index types do not match.");
// Both types are integral types and equal, define a single GlobalIndexType.
using GlobalIndexType = GlobalMatrix::IndexType;
 
using GlobalSparsityPattern = MathLib::SparsityPattern<GlobalIndexType>;