EqnsStruct.hpp

//
// Created by Abhinav Singh on 17.04.20.
//
//Example file to be used only when DCPSE_Solver is defined before.
//Create your own structures for your equations after including the solver.
#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 constexpr bool boundary[]={NON_PERIODIC, NON_PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations2d2 {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 2;

    static constexpr bool boundary[]={NON_PERIODIC, NON_PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};


struct equations2d1p {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 1;

    static constexpr bool boundary[]={PERIODIC, PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations2d2p {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 2;

    static constexpr bool boundary[]={PERIODIC, PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};


struct equations2d3p {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 3;

    static constexpr bool boundary[]={PERIODIC, PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations2d3 {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 3;

    static constexpr bool boundary[]={NON_PERIODIC, NON_PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations2d4 {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 4;

    static constexpr bool boundary[]={NON_PERIODIC, NON_PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};


struct equations3d3 {
    static const unsigned int dims = 3;
    static const unsigned int nvar = 3;

    static constexpr bool boundary[]={NON_PERIODIC, NON_PERIODIC,NON_PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations3d1 {
    static const unsigned int dims = 3;
    static const unsigned int nvar = 1;

    static constexpr bool boundary[]={NON_PERIODIC, NON_PERIODIC,NON_PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations3d3Pz {
    static const unsigned int dims = 3;
    static const unsigned int nvar = 3;

    static constexpr bool boundary[]={NON_PERIODIC, NON_PERIODIC,PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations3d3Pyz {
    static const unsigned int dims = 3;
    static const unsigned int nvar = 3;

    static constexpr bool boundary[]={NON_PERIODIC, PERIODIC,PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations3d3EPxz {
    static const unsigned int dims = 3;
    static const unsigned int nvar = 3;

    static constexpr bool boundary[]={PERIODIC, NON_PERIODIC,PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations3d3EPz {
    static const unsigned int dims = 3;
    static const unsigned int nvar = 3;

    static constexpr bool boundary[]={PERIODIC, PERIODIC,PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations3d3Pxz {
    static const unsigned int dims = 3;
    static const unsigned int nvar = 3;

    static constexpr bool boundary[]={PERIODIC, NON_PERIODIC,PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations3d1Pz {
    static const unsigned int dims = 3;
    static const unsigned int nvar = 1;

    static constexpr bool boundary[]={NON_PERIODIC, NON_PERIODIC,PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};


struct equations2d1E {

    static const unsigned int dims=2;
    static const unsigned int nvar=1;

    static constexpr bool boundary[]={NON_PERIODIC, NON_PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations2d2E {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 2;

    static constexpr bool boundary[]={NON_PERIODIC, NON_PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations2d3E {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 3;

    static constexpr bool boundary[]={NON_PERIODIC, NON_PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations2d4E {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 4;

    static constexpr bool boundary[]={NON_PERIODIC, NON_PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};


struct equations2d1pE {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 1;

    static constexpr bool boundary[]={PERIODIC, PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations2d2pE {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 2;

    static constexpr bool boundary[]={PERIODIC, PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations2d3pE {
    static const unsigned int dims = 2;
    static const unsigned int nvar = 3;

    static constexpr bool boundary[]={PERIODIC, PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations3d3E {
    static const unsigned int dims = 3;
    static const unsigned int nvar = 3;

    static constexpr bool boundary[]={NON_PERIODIC, NON_PERIODIC,NON_PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};

struct equations3d1E {
    static const unsigned int dims = 3;
    static const unsigned int nvar = 1;

    static constexpr bool boundary[]={NON_PERIODIC, NON_PERIODIC,NON_PERIODIC};

    typedef double stype;

    typedef vector_dist<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;
};


#endif //OPENFPM_PDATA_EQNSSTRUCT_HPP