#include "Interpolate.H"
#include "Interpolate_K.H"

namespace Interpolate
{
    using namespace amrex;

    std::unique_ptr<MultiFab>
    getInterpolatedScalar(
        const MultiFab& F_cp, const MultiFab& F_fp,
        const DistributionMapping& dm, const int r_ratio,
        const Real* /*dx*/, const int ngrow )
    {
        // Prepare the structure that will contain the returned fields
        std::unique_ptr<MultiFab> interpolated_F;
        interpolated_F = std::make_unique<MultiFab>(F_fp.boxArray(), dm, 1, ngrow);
        interpolated_F->setVal(0.);

        // Loop through the boxes and interpolate the values from the _cp data
#ifdef _OPENMP
#pragma omp parallel
#endif
        {
            FArrayBox ffab; // Temporary array ; contains interpolated fields
            for (MFIter mfi(*interpolated_F); mfi.isValid(); ++mfi)
            {
                Box finebx = mfi.fabbox();
                finebx.coarsen(r_ratio).refine(r_ratio); // so that finebx is coarsenable

                const FArrayBox& cfab = (F_cp)[mfi];
                ffab.resize(finebx);

                // - Fully nodal
                if ( F_fp.is_nodal() ){
                    IntVect refinement_vector{AMREX_D_DECL(r_ratio, r_ratio, r_ratio)};
                    node_bilinear_interp.interp(cfab, 0, ffab, 0, 1,
                                                finebx, refinement_vector, {}, {}, {}, 0, 0, RunOn::Cpu);
                } else {
                    amrex::Abort("Unknown field staggering.");
                }

                // Add temporary array to the returned structure
                const Box& bx = (*interpolated_F)[mfi].box();
                (*interpolated_F)[mfi].plus<RunOn::Host>(ffab, bx, bx, 0, 0, 1);
            }
        }
        return interpolated_F;
    }

    std::array<std::unique_ptr<MultiFab>, 3>
    getInterpolatedVector(
        const MultiFab* Fx_cp,
        const MultiFab* Fy_cp,
        const MultiFab* Fz_cp,
        const MultiFab* Fx_fp,
        const MultiFab* Fy_fp,
        const MultiFab* Fz_fp,
        const DistributionMapping& dm, const int r_ratio,
        const Real* /*dx*/, const int ngrow )
    {

        // Prepare the structure that will contain the returned fields
        std::array<std::unique_ptr<MultiFab>, 3> interpolated_F;
        interpolated_F[0] = std::make_unique<MultiFab>(Fx_fp->boxArray(), dm, 1, ngrow);
        interpolated_F[1] = std::make_unique<MultiFab>(Fy_fp->boxArray(), dm, 1, ngrow);
        interpolated_F[2] = std::make_unique<MultiFab>(Fz_fp->boxArray(), dm, 1, ngrow);

        IntVect fx_type = interpolated_F[0]->ixType().toIntVect();
        IntVect fy_type = interpolated_F[1]->ixType().toIntVect();
        IntVect fz_type = interpolated_F[2]->ixType().toIntVect();

#ifdef _OPENMP
#pragma omp parallel
#endif
        for (MFIter mfi(*interpolated_F[0], TilingIfNotGPU()); mfi.isValid(); ++mfi)
        {
            Box const& boxx = mfi.growntilebox(fx_type);
            Box const& boxy = mfi.growntilebox(fy_type);
            Box const& boxz = mfi.growntilebox(fz_type);

            Array4<Real      > const& fx = interpolated_F[0]->array(mfi);
            Array4<Real      > const& fy = interpolated_F[1]->array(mfi);
            Array4<Real      > const& fz = interpolated_F[2]->array(mfi);
            Array4<Real const> const& cx = Fx_cp->const_array(mfi);
            Array4<Real const> const& cy = Fy_cp->const_array(mfi);
            Array4<Real const> const& cz = Fz_cp->const_array(mfi);

            amrex::ParallelFor(boxx, boxy, boxz,
                [=] AMREX_GPU_DEVICE (int j, int k, int l) noexcept
                {
                    interp(j,k,l,fx,cx,r_ratio,fx_type);
                },
                [=] AMREX_GPU_DEVICE (int j, int k, int l) noexcept
                {
                    interp(j,k,l,fy,cy,r_ratio,fy_type);
                },
                [=] AMREX_GPU_DEVICE (int j, int k, int l) noexcept
                {
                    interp(j,k,l,fz,cz,r_ratio,fz_type);
                });
        }
        return interpolated_F;
    }
}
id=17d807f05fddbcd58656cd6888c7d5b85d52e441&amp;follow=1'>golang.org</a>/<a href='/ansg191/coredns/log/vendor/golang.org/x?id=17d807f05fddbcd58656cd6888c7d5b85d52e441&amp;follow=1'>x</a>/<a href='/ansg191/coredns/log/vendor/golang.org/x/net?id=17d807f05fddbcd58656cd6888c7d5b85d52e441&amp;follow=1'>net</a>/<a href='/ansg191/coredns/log/vendor/golang.org/x/net/AUTHORS?id=17d807f05fddbcd58656cd6888c7d5b85d52e441&amp;follow=1'>AUTHORS</a> (<a href='/ansg191/coredns/log/vendor/golang.org/x/net/AUTHORS?id=17d807f05fddbcd58656cd6888c7d5b85d52e441'>unfollow</a>)</div><div class='content'><table class='list nowrap'><tr class='nohover'><th class='left'>Age</th><th class='left'>Commit message (<a href='/ansg191/coredns/log/vendor/golang.org/x/net/AUTHORS?id=17d807f05fddbcd58656cd6888c7d5b85d52e441&amp;showmsg=1&amp;follow=1'>Expand</a>)</th><th class='left'>Author</th><th class='left'>Files</th><th class='left'>Lines</th></tr>
<tr><td><span title='2018-06-27 21:12:27 +0100'>2018-06-27</span></td><td><a href='/ansg191/coredns/commit/plugin/pkg/dnsutil/ttl_test.go?id=dae506b5638c7309399cb273d7f76bc20ee518dd&amp;follow=1'>Fix max-age in http server (#1890)</a></td><td><img src='//www.gravatar.com/avatar/25ba01eee5a1c0e007e810afe805766c?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Miek Gieben
</td><td>5</td><td><span class='deletions'>-68</span>/<span class='insertions'>+91</span></td></tr>
<tr><td><span title='2018-06-27 07:45:32 -0700'>2018-06-27</span></td><td><a href='/ansg191/coredns/commit/plugin/pkg/watch/watcher.go?id=99287d091c2db4028e54782fd4de43f63ca4b040&amp;follow=1'>Watch feature (#1527)</a></td><td><img src='//www.gravatar.com/avatar/92aec22a2f92e81152610b8b40f49a21?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> John Belamaric
</td><td>17</td><td><span class='deletions'>-57</span>/<span class='insertions'>+977</span></td></tr>
<tr><td><span title='2018-06-26 09:11:17 +0200'>2018-06-26</span></td><td><a href='/ansg191/coredns/commit/plugin/cache/handler.go?id=b7480d5d1216aa87d80d240a31de750079eba904&amp;follow=1'>plugin/cache: Fix: metric `cache miss total`  shall include ttl case (#1897)</a></td><td><img src='//www.gravatar.com/avatar/a826be33bdce4e546dac998498ae9c7f?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Jun Li
</td><td>1</td><td><span class='deletions'>-8</span>/<span class='insertions'>+9</span></td></tr>
<tr><td><span title='2018-06-24 19:21:58 +0100'>2018-06-24</span></td><td><a href='/ansg191/coredns/commit/Makefile.release?id=c25a2e0fac398c4454578abef8ea0870099fac97&amp;follow=1'>Replace gh-release with curl (#1896)</a></td><td><img src='//www.gravatar.com/avatar/10c4cad3ba13f0378dbd67185314d3d4?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Yong Tang
</td><td>1</td><td><span class='deletions'>-3</span>/<span class='insertions'>+10</span></td></tr>
<tr><td><span title='2018-06-22 10:54:40 +0100'>2018-06-22</span></td><td><a href='/ansg191/coredns/commit/plugin/pkg/dnstest/server.go?id=1ef9e31262b6dff796d35ee667cb9b93e723b13d&amp;follow=1'>retry creating dnstest.Server if tcp or udp port is in use (#1799)</a></td><td><img src='//www.gravatar.com/avatar/13cf34712525e2848d4468120b95674e?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Ruslan Drozhdzh
</td><td>1</td><td><span class='deletions'>-10</span>/<span class='insertions'>+20</span></td></tr>
<tr><td><span title='2018-06-21 14:38:29 -0400'>2018-06-21</span></td><td><a href='/ansg191/coredns/commit/test/template_upstream_test.go?id=716791aa094b54ea17241112a14eeecac10dcfa5&amp;follow=1'>plugin/template: fix CNAME upstream handling (#1886)</a></td><td><img src='//www.gravatar.com/avatar/c4e1dff64e2c0210e86f613a22916fc1?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Chris O'Haver
</td><td>4</td><td><span class='deletions'>-4</span>/<span class='insertions'>+97</span></td></tr>
<tr><td><span title='2018-06-21 04:34:02 -0700'>2018-06-21</span></td><td><a href='/ansg191/coredns/commit/plugin/kubernetes/handler_ignore_emptyservice_test.go?id=ad8021230c775f129556fe7f681f0b09a59d480e&amp;follow=1'>Makefile fixes: (#1892)</a></td><td><img src='//www.gravatar.com/avatar/25ba01eee5a1c0e007e810afe805766c?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Miek Gieben
</td><td>4</td><td><span class='deletions'>-17</span>/<span class='insertions'>+12</span></td></tr>
<tr><td><span title='2018-06-21 11:40:19 +0100'>2018-06-21</span></td><td><a href='/ansg191/coredns/commit/plugin/forward/proxy_test.go?id=422aec5f5fc9d1f4d4017357f86098ec34eb1676&amp;follow=1'>plugin/forward: Increase minimum read timeout to 200ms (#1889)</a></td><td><img src='//www.gravatar.com/avatar/90f0e512184d202445616446467a3da9?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Tobias Schmidt
</td><td>5</td><td><span class='deletions'>-10</span>/<span class='insertions'>+10</span></td></tr>
<tr><td><span title='2018-06-20 11:35:57 +0100'>2018-06-20</span></td><td><a href='/ansg191/coredns/commit/plugin/forward/metrics.go?id=e3534205c7c3a89b09c9a9e2aca1a2159eb9ecde&amp;follow=1'>Rename forward metrics socket_count_total to sockets_open (#1885)</a></td><td><img src='//www.gravatar.com/avatar/90f0e512184d202445616446467a3da9?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Tobias Schmidt
</td><td>1</td><td><span class='deletions'>-1</span>/<span class='insertions'>+1</span></td></tr>
<tr><td><span title='2018-06-20 11:35:07 +0100'>2018-06-20</span></td><td><a href='/ansg191/coredns/commit/plugin/proxy/upstream.go?id=41c724780f6a4e9e8461b56472660d9ad3e49540&amp;follow=1'>remove deprecated code (#1888)</a></td><td><img src='//www.gravatar.com/avatar/25ba01eee5a1c0e007e810afe805766c?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Miek Gieben
</td><td>3</td><td><span class='deletions'>-4</span>/<span class='insertions'>+0</span></td></tr>
<tr><td><span title='2018-06-19 19:50:08 +0100'>2018-06-19</span></td><td><a href='/ansg191/coredns/commit/plugin/cache/handler.go?id=9c2dc7a1564a4cfc5c197d99c828364fc02bfbae&amp;follow=1'>Ensure cache.ResponseWriter can be used asynchronously during prefetch (#1884)</a></td><td><img src='//www.gravatar.com/avatar/90f0e512184d202445616446467a3da9?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Tobias Schmidt
</td><td>2</td><td><span class='deletions'>-8</span>/<span class='insertions'>+40</span></td></tr>
<tr><td><span title='2018-06-19 17:57:35 +0100'>2018-06-19</span></td><td><a href='/ansg191/coredns/commit/man/coredns-rewrite.7?id=f78f30231df90da6184d5f811ecf9c06b0160c2b&amp;follow=1'>Release 1.1.4</a><span class='decoration'><a class='tag-deco' href='/ansg191/coredns/tag/?h=v1.1.4'>v1.1.4</a></span></td><td><img src='//www.gravatar.com/avatar/25ba01eee5a1c0e007e810afe805766c?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Miek Gieben
</td><td>3</td><td><span class='deletions'>-4</span>/<span class='insertions'>+7</span></td></tr>
<tr><td><span title='2018-06-18 11:16:56 +0100'>2018-06-18</span></td><td><a href='/ansg191/coredns/commit/plugin/backend_lookup.go?id=ab9efdcac4914b8bed486d32128ae00adba0dc90&amp;follow=1'>Follow up on #1867 : use one map in backend_lookup (#1869)</a></td><td><img src='//www.gravatar.com/avatar/25ba01eee5a1c0e007e810afe805766c?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Miek Gieben
</td><td>1</td><td><span class='deletions'>-32</span>/<span class='insertions'>+33</span></td></tr>
<tr><td><span title='2018-06-18 10:02:05 +0100'>2018-06-18</span></td><td><a href='/ansg191/coredns/commit/plugin/rewrite/README.md?id=6fbc1f8990601c322255f78476b281ba967dbd24&amp;follow=1'>Plugin/Rewrite - review documentation of "append" (#1877)</a></td><td><img src='//www.gravatar.com/avatar/85ef5ede2ac268872e757897d59a4021?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Francois Tur
</td><td>1</td><td><span class='deletions'>-2</span>/<span class='insertions'>+2</span></td></tr>
<tr><td><span title='2018-06-15 08:12:56 -0700'>2018-06-15</span></td><td><a href='/ansg191/coredns/commit/plugin/template/setup.go?id=177e32b62e1bd937317ec3fd7c7f3a8114c52d77&amp;follow=1'>plugin/forward: add REFUSED test (#1878)</a></td><td><img src='//www.gravatar.com/avatar/25ba01eee5a1c0e007e810afe805766c?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Miek Gieben
</td><td>10</td><td><span class='deletions'>-10</span>/<span class='insertions'>+62</span></td></tr>
<tr><td><span title='2018-06-15 07:37:22 +0100'>2018-06-15</span></td><td><a href='/ansg191/coredns/commit/plugin/forward/proxy.go?id=70c957d885a6c6b33a0828c680e451a9e3a3994a&amp;follow=1'>Plugin/Forward - autotune the dialTimeout for connection (#1852)</a></td><td><img src='//www.gravatar.com/avatar/85ef5ede2ac268872e757897d59a4021?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Francois Tur
</td><td>5</td><td><span class='deletions'>-52</span>/<span class='insertions'>+84</span></td></tr>
<tr><td><span title='2018-06-13 11:02:14 +0100'>2018-06-13</span></td><td><a href='/ansg191/coredns/commit/man/coredns.1?id=6f865a9de7d5ddc3daf52512016a5bbf72e3fb9a&amp;follow=1'>Run make -f Makefile.doc (#1875)</a></td><td><img src='//www.gravatar.com/avatar/25ba01eee5a1c0e007e810afe805766c?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Miek Gieben
</td><td>9</td><td><span class='deletions'>-71</span>/<span class='insertions'>+50</span></td></tr>
<tr><td><span title='2018-06-13 10:38:24 +0100'>2018-06-13</span></td><td><a href='/ansg191/coredns/commit/Makefile?id=d6288c25ba4d2ce06e9c7c2ac63fc6b01379a533&amp;follow=1'>Bump dns version to 1.0.8 (#1873)</a></td><td><img src='//www.gravatar.com/avatar/25ba01eee5a1c0e007e810afe805766c?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Miek Gieben
</td><td>1</td><td><span class='deletions'>-1</span>/<span class='insertions'>+1</span></td></tr>
<tr><td><span title='2018-06-12 14:54:37 +0100'>2018-06-12</span></td><td><a href='/ansg191/coredns/commit/test/file_upstream_test.go?id=26c41a0c177fbe89f1f81214634aed49790581f4&amp;follow=1'>plugin/file: fix local CNAME lookup (#1866)</a></td><td><img src='//www.gravatar.com/avatar/25ba01eee5a1c0e007e810afe805766c?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Miek Gieben
</td><td>5</td><td><span class='deletions'>-28</span>/<span class='insertions'>+83</span></td></tr>
<tr><td><span title='2018-06-11 19:23:25 -0700'>2018-06-11</span></td><td><a href='/ansg191/coredns/commit/plugin/pkg/dnsutil/dedup.go?id=6e466d509281953fdf0209a5b50611e89b4689ae&amp;follow=1'>Remove dnsutil.Dedup (#1867)</a></td><td><img src='//www.gravatar.com/avatar/25ba01eee5a1c0e007e810afe805766c?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Miek Gieben
</td><td>6</td><td><span class='deletions'>-61</span>/<span class='insertions'>+66</span></td></tr>
<tr><td><span title='2018-06-07 13:48:23 -0700'>2018-06-07</span></td><td><a href='/ansg191/coredns/commit/.github/PULL_REQUEST_TEMPLATE.md?id=58d69913639b7b210ed32afacb447ba53c7487a3&amp;follow=1'>Update the pr template (#1862)</a></td><td><img src='//www.gravatar.com/avatar/25ba01eee5a1c0e007e810afe805766c?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Miek Gieben
</td><td>1</td><td><span class='deletions'>-5</span>/<span class='insertions'>+1</span></td></tr>
<tr><td><span title='2018-06-07 16:21:17 +0100'>2018-06-07</span></td><td><a href='/ansg191/coredns/commit/plugin/pkg/replacer/replacer_test.go?id=18b11fc85156be8b027c61665c80e239324e2ebf&amp;follow=1'>Normalizing Response Duration in Log Plugin (#1860)</a></td><td><img src='//www.gravatar.com/avatar/d88558178af5dfa7a346c9a0f737a4f6?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Malcolm Akinje
</td><td>2</td><td><span class='deletions'>-1</span>/<span class='insertions'>+5</span></td></tr>
<tr><td><span title='2018-06-05 17:21:09 +0100'>2018-06-05</span></td><td><a href='/ansg191/coredns/commit/plugin/forward/proxy.go?id=751a08d6a2e3ace62e40ed8170015e7051d0020c&amp;follow=1'>plugin/forward: fix alignment for sync.Atomic (#1855)</a></td><td><img src='//www.gravatar.com/avatar/25ba01eee5a1c0e007e810afe805766c?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Miek Gieben
</td><td>1</td><td><span class='deletions'>-3</span>/<span class='insertions'>+3</span></td></tr>
<tr><td><span title='2018-06-05 17:20:38 +0100'>2018-06-05</span></td><td><a href='/ansg191/coredns/commit/Makefile.release?id=ec334014d73222981ef84667f0e82d04d1033537&amp;follow=1'>Update Makefile.release for docker image build (#1856)</a></td><td><img src='//www.gravatar.com/avatar/10c4cad3ba13f0378dbd67185314d3d4?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Yong Tang
</td><td>1</td><td><span class='deletions'>-15</span>/<span class='insertions'>+0</span></td></tr>
<tr><td><span title='2018-06-04 19:20:21 +0100'>2018-06-04</span></td><td><a href='/ansg191/coredns/commit/Dockerfile?id=a77d9834d867b44a7f10d1481916498332c586ff&amp;follow=1'>Docker: drop alpine (#1843)</a></td><td><img src='//www.gravatar.com/avatar/25ba01eee5a1c0e007e810afe805766c?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Miek Gieben
</td><td>1</td><td><span class='deletions'>-3</span>/<span class='insertions'>+6</span></td></tr>
<tr><td><span title='2018-06-02 20:28:56 +0100'>2018-06-02</span></td><td><a href='/ansg191/coredns/commit/test/plugin_dnssec_test.go?id=70a53e25c942d95de98e30d9f5f5636d0d315d25&amp;follow=1'>Fix presubmit (#1850)</a></td><td><img src='//www.gravatar.com/avatar/25ba01eee5a1c0e007e810afe805766c?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Miek Gieben
</td><td>2</td><td><span class='deletions'>-15</span>/<span class='insertions'>+15</span></td></tr>
<tr><td><span title='2018-06-02 11:48:39 -0700'>2018-06-02</span></td><td><a href='/ansg191/coredns/commit/plugin/test/helpers.go?id=22c0b30d5f7468c2477291419bee3e88dfba490c&amp;follow=1'>presubmit: Check errorf as well (#1845)</a></td><td><img src='//www.gravatar.com/avatar/25ba01eee5a1c0e007e810afe805766c?s=13&amp;d=retro' width='13' height='13' alt='Gravatar' /> Miek Gieben
</td><td>19</td><td><span class='deletions'>-52</span>/<span class='insertions'>+57</span></td></tr>