diff options
| author | Erik Schnetter <schnetter@cct.lsu.edu> | 2012-02-24 14:43:13 -0500 |
|---|---|---|
| committer | Barry Wardell <barry.wardell@gmail.com> | 2012-09-11 18:23:03 +0100 |
| commit | a50b94db8d36ffd6c9aca449299e8148c2b8fd33 (patch) | |
| tree | c45c57f43e14e22807bc954dc3bd89270d22f643 /Carpet/CarpetLib/src/data.cc | |
| parent | 181db01dd6579975d5b6d9a56a5ad4792012f1dd (diff) | |
CarpetLib: Add preliminary support for DGFE operators
Diffstat (limited to 'Carpet/CarpetLib/src/data.cc')
| -rw-r--r-- | Carpet/CarpetLib/src/data.cc | 26 |
1 files changed, 26 insertions, 0 deletions
diff --git a/Carpet/CarpetLib/src/data.cc b/Carpet/CarpetLib/src/data.cc index 24534b376..1ee6c94e7 100644 --- a/Carpet/CarpetLib/src/data.cc +++ b/Carpet/CarpetLib/src/data.cc @@ -711,6 +711,8 @@ transfer_prolongate (data const * const src, ibbox const & srcbox, int const order_space) { + DECLARE_CCTK_PARAMETERS; + static Timer total ("prolongate"); total.start (); @@ -773,6 +775,17 @@ transfer_prolongate (data const * const src, break; } case cell_centered: { + if (use_dgfe) { + call_operator<T>(prolongate_3d_dgfe_rf2<T,5>, + static_cast<T const *>(src->storage()), + src->shape(), + static_cast<T *>(this->storage()), + this->shape(), + src->extent(), + this->extent(), + srcbox, dstbox, NULL); + break; + } static void (* the_operators[]) (T const * restrict const src, ivect3 const & restrict srcext, @@ -1084,6 +1097,8 @@ transfer_restrict (data const * const src, islab const * restrict const slabinfo, int const /*order_space*/) { + DECLARE_CCTK_PARAMETERS; + static Timer total ("restrict"); total.start (); @@ -1119,6 +1134,17 @@ transfer_restrict (data const * const src, ibbox const& dstbox = this->extent(); if (all(is_centered == ivect(1,1,1))) { + if (use_dgfe) { + call_operator<T>(restrict_3d_dgfe_rf2<T,5>, + static_cast<T const *>(src->storage()), + src->shape(), + static_cast<T *>(this->storage()), + this->shape(), + srcbox, + dstbox, + srcregbox, dstregbox, NULL); + break; + } call_operator<T> (& restrict_3d_cc_rf2, static_cast <T const *> (src->storage()), src->shape(), |