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 : }
|