Line data Source code
1 : #include <AMReX_MultiFabUtil.H> 2 : #include <AMReX_REAL.H> 3 : #include <AMReX_MLCGSolver.H> 4 : #include "Set/Set.H" 5 : 6 : #include <AMReX_ArrayLim.H> 7 : 8 : #include "Util/Util.H" 9 : #include "Set/Set.H" 10 : #include "Implicit.H" 11 : 12 : namespace Operator 13 : { 14 : namespace Implicit 15 : { 16 0 : Implicit::Implicit (const Vector<Geometry>& a_geom, 17 : const Vector<BoxArray>& a_grids, 18 : const Vector<DistributionMapping>& a_dmap, 19 : BC::BC<Set::Scalar>& a_bc, 20 0 : const LPInfo& a_info) 21 : { 22 0 : define(a_geom, a_grids, a_dmap, a_bc, a_info); 23 0 : } 24 : 25 : void 26 0 : Implicit::Fapply (int /*amrlev*/, ///<[in] AMR Level 27 : int /*mglev*/, ///<[in] 28 : MultiFab& /*f*/,///<[out] The force vector 29 : const MultiFab& /*u*/ ///<[in] The displacements vector 30 : ) const 31 : { 32 0 : Util::Message(INFO); 33 0 : } 34 : 35 : 36 : void 37 0 : Implicit::Fsmooth (int /*amrlev*/, ///<[in] AMR level 38 : int /*mglev*/, ///<[in] 39 : MultiFab& /*u*/, ///<[inout] Solution (displacement field) 40 : const MultiFab& /*rhs*/, ///<[in] Body force vectors (rhs=right hand side) 41 : int /*redblack*/ ///<[in] Smooth even vs. odd modes 42 : ) const 43 : { 44 0 : Util::Message(INFO); 45 0 : } 46 : 47 0 : void Implicit::FFlux (int /*amrlev*/, const MFIter& /*mfi*/, 48 : const Array<FArrayBox*,AMREX_SPACEDIM>& sigmafab, 49 : const FArrayBox& /*sol*/, Location /*loc*/, const int /*face_only*/) const 50 : { 51 0 : amrex::BaseFab<amrex::Real> AMREX_D_DECL( &fxfab = *sigmafab[0], 52 : &fyfab = *sigmafab[1], 53 : &fzfab = *sigmafab[2] ) ; 54 0 : AMREX_D_TERM(fxfab.setVal(0.0);, 55 : fyfab.setVal(0.0);, 56 : fzfab.setVal(0.0);); 57 : 58 0 : } 59 : 60 : } 61 : }