AMReX_SDCstruct.H

Go to the documentation of this file.

#ifndef SDCSTRUCT_H_
#define SDCSTRUCT_H_
#include <AMReX_Config.H>
 
#include <AMReX_Geometry.H>
#include <AMReX_MultiFab.H>
#include <AMReX_BCRec.H>
 
using namespace amrex;
 
extern "C"  {
 
  void SDC_quadrature(int* qtype_in, int* nnodes, int* nnodes0,
                      amrex_real* nodes, int* nflags,amrex_real* qmats);
}
class SDCstruct
 {
 
 public:
   Real* qnodes;      
   int*  Nflags;      
 
   int Nnodes;         
   int qtype=1;        
   int Nsweeps=8;      
   int Npieces;        
 
 
   // The quadrature matrices
   Real* Qall;
   Vector<Vector<Real>> Qgauss;     
   Vector<Vector<Real>> Qexp;       
   Vector<Vector<Real>> Qimp;       
   Vector<Vector<Real>> QLU;        
 
 
   // SDC storage
   Vector<MultiFab> sol;           
   Vector<Vector<MultiFab> > f;    
   Vector<MultiFab> res;           
   Vector<MultiFab> Ithree;        
 
 
   SDCstruct(int Nnodes_in,int Npieces_in, MultiFab& sol_in);
 
   //  Sweeper routines
   void SDC_rhs_integrals(Real dt);
 
   void SDC_rhs_k_plus_one(MultiFab& rhs, Real dt,int m);
 
   void SDC_rhs_misdc(MultiFab& rhs, Real dt,int m);
 
 };
 
#endif