EqnsStructFD.hpp

//
// Created by Abhinav Singh on 17.04.20.
//
 
#ifndef OPENFPM_PDATA_EQNSSTRUCT_HPP
#define OPENFPM_PDATA_EQNSSTRUCT_HPP
 
#include "Solvers/umfpack_solver.hpp"
#include "Solvers/petsc_solver.hpp"
 
struct equations2d1 {
 
    static const unsigned int dims=2;
    static const unsigned int nvar=1;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double,double,double>> b_part;
 
    typedef SparseMatrix<double, int, PETSC_BASE> SparseMatrix_type;
 
    typedef Vector<double, PETSC_BASE> Vector_type;
 
    typedef petsc_solver<double> solver_type;
};
 
const bool equations2d1::boundary[] = {NON_PERIODIC,NON_PERIODIC};
 
struct equations2d2 {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 2;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, PETSC_BASE> SparseMatrix_type;
 
    typedef Vector<double, PETSC_BASE> Vector_type;
 
    typedef petsc_solver<double> solver_type;
};
const bool equations2d2::boundary[]={NON_PERIODIC,NON_PERIODIC};
 
struct equations2d1p {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 1;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, PETSC_BASE> SparseMatrix_type;
 
    typedef Vector<double, PETSC_BASE> Vector_type;
 
    typedef petsc_solver<double> solver_type;
};
const bool equations2d1p::boundary[]={PERIODIC,PERIODIC};
 
struct equations2d2p {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 2;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, PETSC_BASE> SparseMatrix_type;
 
    typedef Vector<double, PETSC_BASE> Vector_type;
 
    typedef petsc_solver<double> solver_type;
};
const bool equations2d2p::boundary[]={PERIODIC,PERIODIC};
 
struct equations2d3p {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 3;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, PETSC_BASE> SparseMatrix_type;
 
    typedef Vector<double, PETSC_BASE> Vector_type;
 
    typedef petsc_solver<double> solver_type;
};
const bool equations2d3p::boundary[]={PERIODIC,PERIODIC};
 
 
struct equations2d3 {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 3;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, PETSC_BASE> SparseMatrix_type;
 
    typedef Vector<double, PETSC_BASE> Vector_type;
 
    typedef petsc_solver<double> solver_type;
};
const bool equations2d3::boundary[]={NON_PERIODIC,NON_PERIODIC};
 
struct equations2d4 {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 4;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, PETSC_BASE> SparseMatrix_type;
 
    typedef Vector<double, PETSC_BASE> Vector_type;
 
    typedef petsc_solver<double> solver_type;
};
const bool equations2d4::boundary[]={NON_PERIODIC,NON_PERIODIC};
 
 
struct equations3d3 {
    static const unsigned int dims = 3;
    static const unsigned int nvar = 3;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, PETSC_BASE> SparseMatrix_type;
 
    typedef Vector<double, PETSC_BASE> Vector_type;
 
    typedef petsc_solver<double> solver_type;
};
const bool equations3d3::boundary[]={NON_PERIODIC,NON_PERIODIC,NON_PERIODIC};
 
struct equations3d1 {
    static const unsigned int dims = 3;
    static const unsigned int nvar = 1;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, PETSC_BASE> SparseMatrix_type;
 
    typedef Vector<double, PETSC_BASE> Vector_type;
 
    typedef petsc_solver<double> solver_type;
};
const bool equations3d1::boundary[]={NON_PERIODIC,NON_PERIODIC,NON_PERIODIC};
 
 
struct equations2d1E {
 
    static const unsigned int dims=2;
    static const unsigned int nvar=1;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, EIGEN_BASE> SparseMatrix_type;
 
    typedef Vector<double> Vector_type;
 
    typedef umfpack_solver<double> solver_type;
};
const bool equations2d1E::boundary[]={NON_PERIODIC,NON_PERIODIC};
 
struct equations2d2E {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 2;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, EIGEN_BASE> SparseMatrix_type;
 
    typedef Vector<double> Vector_type;
 
    typedef umfpack_solver<double> solver_type;
};
const bool equations2d2E::boundary[]={NON_PERIODIC,NON_PERIODIC};
 
struct equations2d3E {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 3;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, EIGEN_BASE> SparseMatrix_type;
 
    typedef Vector<double> Vector_type;
 
    typedef umfpack_solver<double> solver_type;
};
const bool equations2d3E::boundary[]={NON_PERIODIC,NON_PERIODIC};
 
struct equations2d4E {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 4;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, EIGEN_BASE> SparseMatrix_type;
 
    typedef Vector<double> Vector_type;
 
    typedef umfpack_solver<double> solver_type;
};
const bool equations2d4E::boundary[]={NON_PERIODIC,NON_PERIODIC};
 
 
 
struct equations2d1pE {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 1;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, EIGEN_BASE> SparseMatrix_type;
 
    typedef Vector<double> Vector_type;
 
    typedef umfpack_solver<double> solver_type;
};
const bool equations2d1pE::boundary[]={PERIODIC,PERIODIC};
 
 
struct equations2d2pE {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 2;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, EIGEN_BASE> SparseMatrix_type;
 
    typedef Vector<double> Vector_type;
 
    typedef umfpack_solver<double> solver_type;
};
 
const bool equations2d2pE::boundary[]={PERIODIC,PERIODIC};
 
 
struct equations2d3pE {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 3;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, EIGEN_BASE> SparseMatrix_type;
 
    typedef Vector<double> Vector_type;
 
    typedef umfpack_solver<double> solver_type;
};
const bool equations2d3pE::boundary[]={PERIODIC,PERIODIC,PERIODIC};
 
 
struct equations3d3E {
    static const unsigned int dims = 3;
    static const unsigned int nvar = 3;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, EIGEN_BASE> SparseMatrix_type;
 
    typedef Vector<double> Vector_type;
 
    typedef umfpack_solver<double> solver_type;
};
 
const bool equations3d3E::boundary[]={NON_PERIODIC, NON_PERIODIC,NON_PERIODIC};
 
struct equations3d1E {
    static const unsigned int dims = 3;
    static const unsigned int nvar = 1;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef grid_dist_id<dims, double, aggregate<double>> b_part;
 
    typedef SparseMatrix<double, int, EIGEN_BASE> SparseMatrix_type;
 
    typedef Vector<double> Vector_type;
 
    typedef umfpack_solver<double> solver_type;
};
const bool equations3d1E::boundary[]={NON_PERIODIC, NON_PERIODIC,NON_PERIODIC};
 
 
struct equations2d1_stag {
 
    static const unsigned int dims=2;
    static const unsigned int nvar=1;
 
    static const bool boundary[];
 
    typedef double stype;
 
    typedef staggered_grid_dist<dims, double, aggregate<double,double,double>> b_part;
 
    typedef SparseMatrix<double, int, PETSC_BASE> SparseMatrix_type;
 
    typedef Vector<double, PETSC_BASE> Vector_type;
 
    typedef petsc_solver<double> solver_type;
};
 
const bool equations2d1_stag::boundary[] = {NON_PERIODIC,NON_PERIODIC};
 
 
#endif //OPENFPM_PDATA_EQNSSTRUCT_HPP