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