AMReX_MultiFabUtilI.H

Go to the documentation of this file.

 
// Used by EB_LSCore only.  Might be removed in the future.  So don't use.
 
namespace amrex::MFUtil {
 
 
        template<typename T>
        struct SymmetricGhost {
            static void copy (T & target, const T & source, int nc, int ng) {
                target.ParallelCopy(source, 0, 0, nc, ng, ng);
            }
        };
 
 
        template<typename T> // T can be either MultiFab or iMultiFab
        struct AsymmetricGhost {
            static void copy (T & target, const T & source, int nc, int ng) {
                target.ParallelCopy(source, 0, 0, nc, 0, ng);
            }
        };
 
 
        template<typename T, template<typename> class Cpy>
        T duplicate(const BoxArray & ba, const DistributionMapping & dm, const T & mf_in) {
 
            int ng = mf_in.nGrow();
            int nc = mf_in.nComp();
 
            T mf_out(ba, dm, nc, ng);
 
            Cpy<T>::copy(mf_out, mf_in, nc, ng);
 
            return mf_out;
        }
 
 
        template<typename T, template<typename> class Cpy>
        T duplicate(const BoxArray & ba, const DistributionMapping & dm, const T & mf_in,
                    const Periodicity & periodicity) {
 
            T mf_out = duplicate<T, Cpy>(ba, dm, mf_in);
            mf_out.FillBoundary(periodicity);
 
            return mf_out;
        }
 
 
        template<typename T>
        T regrid(const BoxArray & ba, const DistributionMapping & dm,
                 const T & mf_in, bool regrid_ghost = false ) {
 
            if (regrid_ghost) {
                T mf_out = duplicate<T, SymmetricGhost>(ba, dm, mf_in);
                return mf_out;
            } else {
                T mf_out = duplicate<T, AsymmetricGhost>(ba, dm, mf_in);
                return mf_out;
            }
        }
 
 
        template<typename T>
        void regrid( T & mf_out, const BoxArray & ba, const DistributionMapping & dm,
                     const T & mf_in, bool regrid_ghost = false ) {
 
            int ng = mf_in.nGrow();
            int nc = mf_in.nComp();
 
            mf_out.define(ba, dm, nc, ng);
            if (regrid_ghost) { SymmetricGhost<T>::copy(mf_out, mf_in, nc, ng); }
            else              { AsymmetricGhost<T>::copy(mf_out, mf_in, nc, ng); }
        }
 
 
        template<typename T>
        void regrid( T & mf_out, const BoxArray & ba, const DistributionMapping & dm,
                     const FabFactory<FArrayBox> & eb_factory,
                     const T & mf_in, bool regrid_ghost = false ) {
 
            int ng = mf_in.nGrow();
            int nc = mf_in.nComp();
 
            mf_out.define(ba, dm, nc, ng, MFInfo(), eb_factory);
            if (regrid_ghost) {  SymmetricGhost<T>::copy(mf_out, mf_in, nc, ng); }
            else              { AsymmetricGhost<T>::copy(mf_out, mf_in, nc, ng); }
        }
 
}