Block-Structured AMR Software Framework
AMReX_MFInterp_3D_C.H
Go to the documentation of this file.
1 #ifndef AMREX_MF_INTERP_3D_C_H_
2 #define AMREX_MF_INTERP_3D_C_H_
3 
4 namespace amrex {
5 
6 AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
7 void mf_cell_cons_lin_interp_limit_minmax_llslope (int i, int j, int k, Array4<Real> const& slope,
8  Array4<Real const> const& u, int scomp, int ncomp,
9  Box const& domain, IntVect const& ratio, BCRec const* bc) noexcept
10 {
11  Real sfx = Real(1.0);
12  Real sfy = Real(1.0);
13  Real sfz = Real(1.0);
14 
15  Real sx = Real(0.0);
16  Real sy = Real(0.0);
17  Real sz = Real(0.0);
18 
19  Real dcx = Real(0.0);
20  Real dcy = Real(0.0);
21  Real dcz = Real(0.0);
22 
23  for (int ns = 0; ns < ncomp; ++ns) {
24  int nu = ns + scomp;
25 
26  // x-direction
27  if (ratio[0] > 1) {
28  dcx = mf_compute_slopes_x(i, j, k, u, nu, domain, bc[ns]);
29  Real df = Real(2.0) * (u(i+1,j,k,nu) - u(i ,j,k,nu));
30  Real db = Real(2.0) * (u(i ,j,k,nu) - u(i-1,j,k,nu));
31  sx = (df*db >= Real(0.0)) ? amrex::min(std::abs(df),std::abs(db)) : Real(0.);
32  sx = std::copysign(Real(1.),dcx)*amrex::min(sx,std::abs(dcx));
33  slope(i,j,k,ns ) = dcx; // unlimited slope
34  } else {
35  slope(i,j,k,ns ) = Real(0.0);
36  }
37 
38  // y-direction
39  if (ratio[1] > 1) {
40  dcy = mf_compute_slopes_y(i, j, k, u, nu, domain, bc[ns]);
41  Real df = Real(2.0) * (u(i,j+1,k,nu) - u(i,j ,k,nu));
42  Real db = Real(2.0) * (u(i,j ,k,nu) - u(i,j-1,k,nu));
43  sy = (df*db >= Real(0.0)) ? amrex::min(std::abs(df),std::abs(db)) : Real(0.);
44  sy = std::copysign(Real(1.),dcy)*amrex::min(sy,std::abs(dcy));
45  slope(i,j,k,ns+ ncomp) = dcy; // unlimited slope
46  } else {
47  slope(i,j,k,ns+ ncomp) = Real(0.0);
48  }
49 
50  // z-direction
51  if (ratio[2] > 1) {
52  dcz = mf_compute_slopes_z(i, j, k, u, nu, domain, bc[ns]);
53  Real df = Real(2.0) * (u(i,j,k+1,nu) - u(i,j,k ,nu));
54  Real db = Real(2.0) * (u(i,j,k ,nu) - u(i,j,k-1,nu));
55  sz = (df*db >= Real(0.0)) ? amrex::min(std::abs(df),std::abs(db)) : Real(0.);
56  sz = std::copysign(Real(1.),dcz)*amrex::min(sz,std::abs(dcz));
57  slope(i,j,k,ns+2*ncomp) = dcz; // unlimited slope
58  } else {
59  slope(i,j,k,ns+2*ncomp) = Real(0.0);
60  }
61 
62  // adjust limited slopes to prevent new min/max for this component
63  Real alpha = 1.0;
64  if (sx != Real(0.0) || sy != Real(0.0) || sz != Real(0.0)) {
65  Real dumax = std::abs(sx) * Real(ratio[0]-1)/Real(2*ratio[0])
66  + std::abs(sy) * Real(ratio[1]-1)/Real(2*ratio[1])
67  + std::abs(sz) * Real(ratio[2]-1)/Real(2*ratio[2]);
68  Real umax = u(i,j,k,nu);
69  Real umin = u(i,j,k,nu);
70  int ilim = ratio[0] > 1 ? 1 : 0;
71  int jlim = ratio[1] > 1 ? 1 : 0;
72  int klim = ratio[2] > 1 ? 1 : 0;
73  for (int koff = -klim; koff <= klim; ++koff) {
74  for (int joff = -jlim; joff <= jlim; ++joff) {
75  for (int ioff = -ilim; ioff <= ilim; ++ioff) {
76  umin = amrex::min(umin, u(i+ioff,j+joff,k+koff,nu));
77  umax = amrex::max(umax, u(i+ioff,j+joff,k+koff,nu));
78  }}}
79  if (dumax * alpha > (umax - u(i,j,k,nu))) {
80  alpha = (umax - u(i,j,k,nu)) / dumax;
81  }
82  if (dumax * alpha > (u(i,j,k,nu) - umin)) {
83  alpha = (u(i,j,k,nu) - umin) / dumax;
84  }
85  }
86  sx *= alpha;
87  sy *= alpha;
88  sz *= alpha;
89 
90  // for each direction, compute minimum of the ratio of limited to unlimited slopes
91  if (dcx != Real(0.0)) {
92  sfx = amrex::min(sfx, sx / dcx);
93  }
94  if (dcy != Real(0.0)) {
95  sfy = amrex::min(sfy, sy / dcy);
96  }
97  if (dcz != Real(0.0)) {
98  sfz = amrex::min(sfz, sz / dcz);
99  }
100  }
101 
102  // multiply unlimited slopes by the minimum of the ratio of limited to unlimited slopes
103  for (int ns = 0; ns < ncomp; ++ns) {
104  slope(i,j,k,ns ) *= sfx;
105  slope(i,j,k,ns+ ncomp) *= sfy;
106  slope(i,j,k,ns+2*ncomp) *= sfz;
107  }
108 }
109 
110 AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
111 void mf_cell_cons_lin_interp_llslope (int i, int j, int k, Array4<Real> const& slope,
112  Array4<Real const> const& u, int scomp, int ncomp,
113  Box const& domain, IntVect const& ratio, BCRec const* bc) noexcept
114 {
115  Real sfx = Real(1.0);
116  Real sfy = Real(1.0);
117  Real sfz = Real(1.0);
118 
119  for (int ns = 0; ns < ncomp; ++ns) {
120  int nu = ns + scomp;
121 
122  // x-direction
123  if (ratio[0] > 1) {
124  Real dc = mf_compute_slopes_x(i, j, k, u, nu, domain, bc[ns]);
125  Real df = Real(2.0) * (u(i+1,j,k,nu) - u(i ,j,k,nu));
126  Real db = Real(2.0) * (u(i ,j,k,nu) - u(i-1,j,k,nu));
127  Real sx = (df*db >= Real(0.0)) ? amrex::min(std::abs(df),std::abs(db)) : Real(0.);
128  sx = std::copysign(Real(1.),dc)*amrex::min(sx,std::abs(dc));
129  if (dc != Real(0.0)) {
130  sfx = amrex::min(sfx, sx / dc);
131  }
132  slope(i,j,k,ns ) = dc;
133  } else {
134  slope(i,j,k,ns ) = Real(0.0);
135  }
136 
137  // y-direction
138  if (ratio[1] > 1) {
139  Real dc = mf_compute_slopes_y(i, j, k, u, nu, domain, bc[ns]);
140  Real df = Real(2.0) * (u(i,j+1,k,nu) - u(i,j ,k,nu));
141  Real db = Real(2.0) * (u(i,j ,k,nu) - u(i,j-1,k,nu));
142  Real sy = (df*db >= Real(0.0)) ? amrex::min(std::abs(df),std::abs(db)) : Real(0.);
143  sy = std::copysign(Real(1.),dc)*amrex::min(sy,std::abs(dc));
144  if (dc != Real(0.0)) {
145  sfy = amrex::min(sfy, sy / dc);
146  }
147  slope(i,j,k,ns+ ncomp) = dc;
148  } else {
149  slope(i,j,k,ns+ ncomp) = Real(0.0);
150  }
151 
152 
153  // z-direction
154  if (ratio[2] > 1) {
155  Real dc = mf_compute_slopes_z(i, j, k, u, nu, domain, bc[ns]);
156  Real df = Real(2.0) * (u(i,j,k+1,nu) - u(i,j,k ,nu));
157  Real db = Real(2.0) * (u(i,j,k ,nu) - u(i,j,k-1,nu));
158  Real sz = (df*db >= Real(0.0)) ? amrex::min(std::abs(df),std::abs(db)) : Real(0.);
159  sz = std::copysign(Real(1.),dc)*amrex::min(sz,std::abs(dc));
160  if (dc != Real(0.0)) {
161  sfz = amrex::min(sfz, sz / dc);
162  }
163  slope(i,j,k,ns+2*ncomp) = dc;
164  } else {
165  slope(i,j,k,ns+2*ncomp) = Real(0.0);
166  }
167 
168  }
169 
170  for (int ns = 0; ns < ncomp; ++ns) {
171  slope(i,j,k,ns ) *= sfx;
172  slope(i,j,k,ns+ ncomp) *= sfy;
173  slope(i,j,k,ns+2*ncomp) *= sfz;
174  }
175 }
176 
177 AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
178 void mf_cell_cons_lin_interp_mcslope (int i, int j, int k, int ns, Array4<Real> const& slope,
179  Array4<Real const> const& u, int scomp, int ncomp,
180  Box const& domain, IntVect const& ratio,
181  BCRec const* bc) noexcept
182 {
183  int nu = ns + scomp;
184 
185  Real sx = Real(0.);
186  Real sy = Real(0.);
187  Real sz = Real(0.);
188 
189  // x-direction
190  if (ratio[0] > 1) {
191  Real dc = mf_compute_slopes_x(i, j, k, u, nu, domain, bc[ns]);
192  Real df = Real(2.0) * (u(i+1,j,k,nu) - u(i ,j,k,nu));
193  Real db = Real(2.0) * (u(i ,j,k,nu) - u(i-1,j,k,nu));
194  sx = (df*db >= Real(0.0)) ? amrex::min(std::abs(df),std::abs(db)) : Real(0.);
195  sx = std::copysign(Real(1.),dc)*amrex::min(sx,std::abs(dc));
196  }
197 
198  // y-direction
199  if (ratio[1] > 1) {
200  Real dc = mf_compute_slopes_y(i, j, k, u, nu, domain, bc[ns]);
201  Real df = Real(2.0) * (u(i,j+1,k,nu) - u(i,j ,k,nu));
202  Real db = Real(2.0) * (u(i,j ,k,nu) - u(i,j-1,k,nu));
203  sy = (df*db >= Real(0.0)) ? amrex::min(std::abs(df),std::abs(db)) : Real(0.);
204  sy = std::copysign(Real(1.),dc)*amrex::min(sy,std::abs(dc));
205  }
206 
207  // z-direction
208  if (ratio[2] > 1) {
209  Real dc = mf_compute_slopes_z(i, j, k, u, nu, domain, bc[ns]);
210  Real df = Real(2.0) * (u(i,j,k+1,nu) - u(i,j,k ,nu));
211  Real db = Real(2.0) * (u(i,j,k ,nu) - u(i,j,k-1,nu));
212  sz = (df*db >= Real(0.0)) ? amrex::min(std::abs(df),std::abs(db)) : Real(0.);
213  sz = std::copysign(Real(1.),dc)*amrex::min(sz,std::abs(dc));
214  }
215 
216  Real alpha = 1.0;
217  if (sx != Real(0.0) || sy != Real(0.0) || sz != Real(0.0)) {
218  Real dumax = std::abs(sx) * Real(ratio[0]-1)/Real(2*ratio[0])
219  + std::abs(sy) * Real(ratio[1]-1)/Real(2*ratio[1])
220  + std::abs(sz) * Real(ratio[2]-1)/Real(2*ratio[2]);
221  Real umax = u(i,j,k,nu);
222  Real umin = u(i,j,k,nu);
223  int ilim = ratio[0] > 1 ? 1 : 0;
224  int jlim = ratio[1] > 1 ? 1 : 0;
225  int klim = ratio[2] > 1 ? 1 : 0;
226  for (int koff = -klim; koff <= klim; ++koff) {
227  for (int joff = -jlim; joff <= jlim; ++joff) {
228  for (int ioff = -ilim; ioff <= ilim; ++ioff) {
229  umin = amrex::min(umin, u(i+ioff,j+joff,k+koff,nu));
230  umax = amrex::max(umax, u(i+ioff,j+joff,k+koff,nu));
231  }}}
232  if (dumax * alpha > (umax - u(i,j,k,nu))) {
233  alpha = (umax - u(i,j,k,nu)) / dumax;
234  }
235  if (dumax * alpha > (u(i,j,k,nu) - umin)) {
236  alpha = (u(i,j,k,nu) - umin) / dumax;
237  }
238  }
239 
240  slope(i,j,k,ns ) = sx * alpha;
241  slope(i,j,k,ns+ ncomp) = sy * alpha;
242  slope(i,j,k,ns+2*ncomp) = sz * alpha;
243 }
244 
245 AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
246 void mf_cell_cons_lin_interp (int i, int j, int k, int ns, Array4<Real> const& fine, int fcomp,
247  Array4<Real const> const& slope, Array4<Real const> const& crse,
248  int ccomp, int ncomp, IntVect const& ratio) noexcept
249 {
250  const int ic = amrex::coarsen(i, ratio[0]);
251  const int jc = amrex::coarsen(j, ratio[1]);
252  const int kc = amrex::coarsen(k, ratio[2]);
253  const Real xoff = (static_cast<Real>(i - ic*ratio[0]) + Real(0.5)) / Real(ratio[0]) - Real(0.5);
254  const Real yoff = (static_cast<Real>(j - jc*ratio[1]) + Real(0.5)) / Real(ratio[1]) - Real(0.5);
255  const Real zoff = (static_cast<Real>(k - kc*ratio[2]) + Real(0.5)) / Real(ratio[2]) - Real(0.5);
256  fine(i,j,k,fcomp+ns) = crse(ic,jc,kc,ccomp+ns)
257  + xoff * slope(ic,jc,kc,ns)
258  + yoff * slope(ic,jc,kc,ns+ncomp)
259  + zoff * slope(ic,jc,kc,ns+ncomp*2);
260 }
261 
262 template<typename T>
263 AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
264 void mf_cell_bilin_interp (int i, int j, int k, int n, Array4<T> const& fine, int fcomp,
265  Array4<T const> const& crse, int ccomp, IntVect const& ratio) noexcept
266 {
267  int ic = amrex::coarsen(i,ratio[0]);
268  int jc = amrex::coarsen(j,ratio[1]);
269  int kc = amrex::coarsen(k,ratio[2]);
270  int ioff = i - ic*ratio[0];
271  int joff = j - jc*ratio[1];
272  int koff = k - kc*ratio[2];
273  int sx, sy, sz;
274  Real wx, wy, wz;
275  if (ioff*2 < ratio[0]) {
276  sx = -1;
277  wx = Real(ratio[0]+1+2*ioff) / Real(2*ratio[0]);
278  } else {
279  sx = 1;
280  wx = Real(3*ratio[0]-1-2*ioff) / Real(2*ratio[0]);
281  }
282  if (joff*2 < ratio[1]) {
283  sy = -1;
284  wy = Real(ratio[1]+1+2*joff) / Real(2*ratio[1]);
285  } else {
286  sy = 1;
287  wy = Real(3*ratio[1]-1-2*joff) / Real(2*ratio[1]);
288  }
289  if (koff*2 < ratio[2]) {
290  sz = -1;
291  wz = Real(ratio[2]+1+2*koff) / Real(2*ratio[2]);
292  } else {
293  sz = 1;
294  wz = Real(3*ratio[2]-1-2*koff) / Real(2*ratio[2]);
295  }
296  fine(i,j,k,n+fcomp) =
297  crse(ic ,jc ,kc ,n+ccomp)* wx * wy * wz +
298  crse(ic+sx,jc ,kc ,n+ccomp)*(Real(1.0)-wx)* wy * wz +
299  crse(ic ,jc+sy,kc ,n+ccomp)* wx *(Real(1.0)-wy)* wz +
300  crse(ic+sx,jc+sy,kc ,n+ccomp)*(Real(1.0)-wx)*(Real(1.0)-wy)* wz +
301  crse(ic ,jc ,kc+sz,n+ccomp)* wx * wy *(Real(1.0)-wz) +
302  crse(ic+sx,jc ,kc+sz,n+ccomp)*(Real(1.0)-wx)* wy *(Real(1.0)-wz) +
303  crse(ic ,jc+sy,kc+sz,n+ccomp)* wx *(Real(1.0)-wy)*(Real(1.0)-wz) +
304  crse(ic+sx,jc+sy,kc+sz,n+ccomp)*(Real(1.0)-wx)*(Real(1.0)-wy)*(Real(1.0)-wz);
305 }
306 
307 AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
308 void mf_nodebilin_interp (int i, int j, int k, int n, Array4<Real> const& fine, int fcomp,
309  Array4<Real const> const& crse, int ccomp, IntVect const& ratio) noexcept
310 {
311  int ic = amrex::coarsen(i,ratio[0]);
312  int jc = amrex::coarsen(j,ratio[1]);
313  int kc = amrex::coarsen(k,ratio[2]);
314  int ioff = i - ic*ratio[0];
315  int joff = j - jc*ratio[1];
316  int koff = k - kc*ratio[2];
317  Real rxinv = Real(1.0) / Real(ratio[0]);
318  Real ryinv = Real(1.0) / Real(ratio[1]);
319  Real rzinv = Real(1.0) / Real(ratio[2]);
320  if (ioff != 0 && joff != 0 && koff != 0) {
321  // Fine node at center of cell
322  fine(i,j,k,n+fcomp) = rxinv * ryinv * rzinv *
323  (crse(ic ,jc ,kc ,n+ccomp) * static_cast<Real>((ratio[0]-ioff) * (ratio[1]-joff) * (ratio[2]-koff)) +
324  crse(ic+1,jc ,kc ,n+ccomp) * static_cast<Real>(( ioff) * (ratio[1]-joff) * (ratio[2]-koff)) +
325  crse(ic ,jc+1,kc ,n+ccomp) * static_cast<Real>((ratio[0]-ioff) * ( joff) * (ratio[2]-koff)) +
326  crse(ic+1,jc+1,kc ,n+ccomp) * static_cast<Real>(( ioff) * ( joff) * (ratio[2]-koff)) +
327  crse(ic ,jc ,kc+1,n+ccomp) * static_cast<Real>((ratio[0]-ioff) * (ratio[1]-joff) * ( koff)) +
328  crse(ic+1,jc ,kc+1,n+ccomp) * static_cast<Real>(( ioff) * (ratio[1]-joff) * ( koff)) +
329  crse(ic ,jc+1,kc+1,n+ccomp) * static_cast<Real>((ratio[0]-ioff) * ( joff) * ( koff)) +
330  crse(ic+1,jc+1,kc+1,n+ccomp) * static_cast<Real>(( ioff) * ( joff) * ( koff)));
331  } else if (joff != 0 && koff != 0) {
332  // Node on a Y-Z face
333  fine(i,j,k,n+fcomp) = ryinv * rzinv *
334  (crse(ic,jc ,kc ,n+ccomp) * static_cast<Real>((ratio[1]-joff) * (ratio[2]-koff)) +
335  crse(ic,jc+1,kc ,n+ccomp) * static_cast<Real>(( joff) * (ratio[2]-koff)) +
336  crse(ic,jc ,kc+1,n+ccomp) * static_cast<Real>((ratio[1]-joff) * ( koff)) +
337  crse(ic,jc+1,kc+1,n+ccomp) * static_cast<Real>(( joff) * ( koff)));
338  } else if (ioff != 0 && koff != 0) {
339  // Node on a Z-X face
340  fine(i,j,k,n+fcomp) = rxinv * rzinv *
341  (crse(ic ,jc,kc ,n+ccomp) * static_cast<Real>((ratio[0]-ioff) * (ratio[2]-koff)) +
342  crse(ic+1,jc,kc ,n+ccomp) * static_cast<Real>(( ioff) * (ratio[2]-koff)) +
343  crse(ic ,jc,kc+1,n+ccomp) * static_cast<Real>((ratio[0]-ioff) * ( koff)) +
344  crse(ic+1,jc,kc+1,n+ccomp) * static_cast<Real>(( ioff) * ( koff)));
345  } else if (ioff != 0 && joff != 0) {
346  // Node on a X-Y face
347  fine(i,j,k,n+fcomp) = rxinv * ryinv *
348  (crse(ic ,jc ,kc,n+ccomp) * static_cast<Real>((ratio[0]-ioff) * (ratio[1]-joff)) +
349  crse(ic+1,jc ,kc,n+ccomp) * static_cast<Real>(( ioff) * (ratio[1]-joff)) +
350  crse(ic ,jc+1,kc,n+ccomp) * static_cast<Real>((ratio[0]-ioff) * ( joff)) +
351  crse(ic+1,jc+1,kc,n+ccomp) * static_cast<Real>(( ioff) * ( joff)));
352  } else if (ioff != 0) {
353  // Node on X line
354  fine(i,j,k,n+fcomp) = rxinv*(static_cast<Real>(ratio[0]-ioff)*crse(ic ,jc,kc,n+ccomp) +
355  static_cast<Real>( ioff)*crse(ic+1,jc,kc,n+ccomp));
356  } else if (joff != 0) {
357  // Node on Y line
358  fine(i,j,k,n+fcomp) = ryinv*(static_cast<Real>(ratio[1]-joff)*crse(ic,jc ,kc,n+ccomp) +
359  static_cast<Real>( joff)*crse(ic,jc+1,kc,n+ccomp));
360  } else if (koff != 0) {
361  // Node on Z line
362  fine(i,j,k,n+fcomp) = rzinv*(static_cast<Real>(ratio[2]-koff)*crse(ic,jc,kc ,n+ccomp) +
363  static_cast<Real>( koff)*crse(ic,jc,kc+1,n+ccomp));
364  } else {
365  // Node coincident with coarse node
366  fine(i,j,k,n+fcomp) = crse(ic,jc,kc,n+ccomp);
367  }
368 }
369 
370 AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
371 void mf_cell_quadratic_calcslope (int i, int j, int k, int n,
372  Array4<Real const> const& crse, int ccomp,
373  Array4<Real> const& slope,
374  Box const& domain,
375  BCRec const* bc) noexcept
376 {
377  int nu = ccomp + n;
378  Real sx = mf_compute_slopes_x(i, j, k, crse, nu, domain, bc[n]);
379  Real sy = mf_compute_slopes_y(i, j, k, crse, nu, domain, bc[n]);
380  Real sz = mf_compute_slopes_z(i, j, k, crse, nu, domain, bc[n]);
381  Real sxx = mf_cell_quadratic_compute_slopes_xx (i, j, k, crse, nu, domain, bc[n]);
382  Real syy = mf_cell_quadratic_compute_slopes_yy (i, j, k, crse, nu, domain, bc[n]);
383  Real szz = mf_cell_quadratic_compute_slopes_zz (i, j, k, crse, nu, domain, bc[n]);
384  Real sxy = mf_cell_quadratic_compute_slopes_xy (i, j, k, crse, nu, domain, bc[n]);
385  Real sxz = mf_cell_quadratic_compute_slopes_xz (i, j, k, crse, nu, domain, bc[n]);
386  Real syz = mf_cell_quadratic_compute_slopes_yz (i, j, k, crse, nu, domain, bc[n]);
387 
388  slope(i,j,k,9*n ) = sx; // x
389  slope(i,j,k,9*n+1) = sy; // y
390  slope(i,j,k,9*n+2) = sz; // z
391  slope(i,j,k,9*n+3) = sxx; // x^2
392  slope(i,j,k,9*n+4) = syy; // y^2
393  slope(i,j,k,9*n+5) = szz; // z^2
394  slope(i,j,k,9*n+6) = sxy; // xy
395  slope(i,j,k,9*n+7) = sxz; // xz
396  slope(i,j,k,9*n+8) = syz; // yz
397 }
398 
399 AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
400 void mf_cell_quadratic_interp (int i, int j, int k, int n,
401  Array4<Real> const& fine, int fcomp,
402  Array4<Real const> const& crse, int ccomp,
403  Array4<Real const> const& slope,
404  IntVect const& ratio) noexcept
405 {
406  int ic = amrex::coarsen(i, ratio[0]);
407  int jc = amrex::coarsen(j, ratio[1]);
408  int kc = amrex::coarsen(k, ratio[2]);
409  int irx = i - ic*ratio[0]; // = abs(i % ratio[0])
410  int jry = j - jc*ratio[1]; // = abs(j % ratio[1])
411  int krz = k - kc*ratio[2]; // = abs(k % ratio[2])
412 
413  // Compute offsets.
414  Real xoff = ( Real(irx) + 0.5_rt ) / Real(ratio[0]) - 0.5_rt;
415  Real yoff = ( Real(jry) + 0.5_rt ) / Real(ratio[1]) - 0.5_rt;
416  Real zoff = ( Real(krz) + 0.5_rt ) / Real(ratio[2]) - 0.5_rt;
417 
418  fine(i,j,k,fcomp+n) = crse(ic,jc,kc,ccomp+n)
419  + xoff * slope(ic,jc,kc,9*n ) // x
420  + yoff * slope(ic,jc,kc,9*n+1) // y
421  + zoff * slope(ic,jc,kc,9*n+2) // z
422  + 0.5_rt * xoff * xoff * slope(ic,jc,kc,9*n+3) // x^2
423  + 0.5_rt * yoff * yoff * slope(ic,jc,kc,9*n+4) // y^2
424  + 0.5_rt * zoff * zoff * slope(ic,jc,kc,9*n+5) // z^2
425  + xoff * yoff * slope(ic,jc,kc,9*n+6) // xy
426  + xoff * zoff * slope(ic,jc,kc,9*n+7) // xz
427  + yoff * zoff * slope(ic,jc,kc,9*n+8); // yz
428 }
429 
430 } // namespace amrex
431 
432 #endif
AMREX_FORCE_INLINE
#define AMREX_FORCE_INLINE
Definition: AMReX_Extension.H:119
AMREX_GPU_HOST_DEVICE
#define AMREX_GPU_HOST_DEVICE
Definition: AMReX_GpuQualifiers.H:20
fine
Array4< Real > fine
Definition: AMReX_InterpFaceRegister.cpp:90
slope
Array4< Real > slope
Definition: AMReX_InterpFaceRegister.cpp:91
crse
Array4< Real const > crse
Definition: AMReX_InterpFaceRegister.cpp:92
abs
#define abs(x)
Definition: complex-type.h:85
amrex
Definition: AMReX_Amr.cpp:49
amrex::mf_compute_slopes_y
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE Real mf_compute_slopes_y(int i, int j, int k, Array4< Real const > const &u, int nu, Box const &domain, BCRec const &bc)
Definition: AMReX_MFInterp_C.H:37
amrex::mf_cell_bilin_interp
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mf_cell_bilin_interp(int i, int, int, int n, Array4< T > const &fine, int fcomp, Array4< T const > const &crse, int ccomp, IntVect const &ratio) noexcept
Definition: AMReX_MFInterp_1D_C.H:186
amrex::mf_cell_quadratic_calcslope
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mf_cell_quadratic_calcslope(int i, int j, int, int n, Array4< Real const > const &crse, int ccomp, Array4< Real > const &slope, Box const &domain, BCRec const *bc) noexcept
Definition: AMReX_MFInterp_2D_C.H:363
amrex::Box
BoxND< AMREX_SPACEDIM > Box
Definition: AMReX_BaseFwd.H:27
amrex::mf_cell_cons_lin_interp_llslope
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mf_cell_cons_lin_interp_llslope(int i, int, int, Array4< Real > const &slope, Array4< Real const > const &u, int scomp, int ncomp, Box const &domain, IntVect const &, BCRec const *bc) noexcept
Definition: AMReX_MFInterp_1D_C.H:39
amrex::max
AMREX_GPU_HOST_DEVICE constexpr AMREX_FORCE_INLINE const T & max(const T &a, const T &b) noexcept
Definition: AMReX_Algorithm.H:35
amrex::min
AMREX_GPU_HOST_DEVICE constexpr AMREX_FORCE_INLINE const T & min(const T &a, const T &b) noexcept
Definition: AMReX_Algorithm.H:21
amrex::mf_cell_cons_lin_interp_limit_minmax_llslope
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mf_cell_cons_lin_interp_limit_minmax_llslope(int i, int, int, Array4< Real > const &slope, Array4< Real const > const &u, int scomp, int ncomp, Box const &domain, IntVect const &ratio, BCRec const *bc) noexcept
Definition: AMReX_MFInterp_1D_C.H:7
amrex::mf_cell_quadratic_compute_slopes_yy
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE Real mf_cell_quadratic_compute_slopes_yy(int i, int j, int k, Array4< Real const > const &u, int nu, Box const &domain, BCRec const &bc)
Definition: AMReX_MFInterp_C.H:110
amrex::mf_cell_quadratic_compute_slopes_xy
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE Real mf_cell_quadratic_compute_slopes_xy(int i, int j, int k, Array4< Real const > const &u, int nu, Box const &domain, BCRec const &bc)
Definition: AMReX_MFInterp_C.H:152
amrex::mf_cell_cons_lin_interp_mcslope
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mf_cell_cons_lin_interp_mcslope(int i, int, int, int ns, Array4< Real > const &slope, Array4< Real const > const &u, int scomp, int, Box const &domain, IntVect const &ratio, BCRec const *bc) noexcept
Definition: AMReX_MFInterp_1D_C.H:66
amrex::mf_cell_quadratic_compute_slopes_yz
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE Real mf_cell_quadratic_compute_slopes_yz(int i, int j, int k, Array4< Real const > const &u, int nu, Box const &domain, BCRec const &bc)
Definition: AMReX_MFInterp_C.H:220
amrex::mf_cell_quadratic_compute_slopes_zz
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE Real mf_cell_quadratic_compute_slopes_zz(int i, int j, int k, Array4< Real const > const &u, int nu, Box const &domain, BCRec const &bc)
Definition: AMReX_MFInterp_C.H:131
amrex::mf_cell_quadratic_compute_slopes_xx
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE Real mf_cell_quadratic_compute_slopes_xx(int i, int j, int k, Array4< Real const > const &u, int nu, Box const &domain, BCRec const &bc)
Definition: AMReX_MFInterp_C.H:89
amrex::IntVect
IntVectND< AMREX_SPACEDIM > IntVect
Definition: AMReX_BaseFwd.H:30
amrex::mf_cell_cons_lin_interp
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mf_cell_cons_lin_interp(int i, int, int, int ns, Array4< Real > const &fine, int fcomp, Array4< Real const > const &slope, Array4< Real const > const &crse, int ccomp, int, IntVect const &ratio) noexcept
Definition: AMReX_MFInterp_1D_C.H:102
amrex::mf_cell_quadratic_compute_slopes_xz
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE Real mf_cell_quadratic_compute_slopes_xz(int i, int j, int k, Array4< Real const > const &u, int nu, Box const &domain, BCRec const &bc)
Definition: AMReX_MFInterp_C.H:186
amrex::mf_nodebilin_interp
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mf_nodebilin_interp(int i, int, int, int n, Array4< Real > const &fine, int fcomp, Array4< Real const > const &crse, int ccomp, IntVect const &ratio) noexcept
Definition: AMReX_MFInterp_1D_C.H:206
amrex::mf_compute_slopes_x
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE Real mf_compute_slopes_x(int i, int j, int k, Array4< Real const > const &u, int nu, Box const &domain, BCRec const &bc)
Definition: AMReX_MFInterp_C.H:11
amrex::coarsen
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE BoxND< dim > coarsen(const BoxND< dim > &b, int ref_ratio) noexcept
Coarsen BoxND by given (positive) refinement ratio. NOTE: if type(dir) = CELL centered: lo <- lo/rati...
Definition: AMReX_Box.H:1304
amrex::mf_cell_quadratic_interp
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void mf_cell_quadratic_interp(int i, int j, int, int n, Array4< Real > const &fine, int fcomp, Array4< Real const > const &crse, int ccomp, Array4< Real const > const &slope, IntVect const &ratio) noexcept
Definition: AMReX_MFInterp_2D_C.H:384
amrex::mf_compute_slopes_z
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE Real mf_compute_slopes_z(int i, int j, int k, Array4< Real const > const &u, int nu, Box const &domain, BCRec const &bc)
Definition: AMReX_MFInterp_C.H:63