30 int iloc=-1, jloc=-1, kloc=-1;
31 for(
int k=0;k<nz;k++){
36 else if (zvec[k] == z) {
42 iloc =
static_cast<int>(std::floor((x - xvec[0]) / dxvec));
43 jloc =
static_cast<int>(std::floor((y - yvec[0]) / dyvec));
45 if(iloc > nx-1 or iloc < 0 or
46 jloc > ny-1 or iloc < 0 or
47 kloc > nz-1 or kloc < 0){
67 tmp_var = (1-w_x)*(1-w_y)*(1-w_z)*varvec[ind0] + w_x*(1-w_y)*(1-w_z)*varvec[ind1] +
68 w_x*w_y*(1-w_z)*varvec[ind2] + (1-w_x)*w_y*(1-w_z)*varvec[ind3] +
69 (1-w_x)*(1-w_y)*w_z*varvec[ind4] + w_x*(1-w_y)*w_z*varvec[ind5] +
70 w_x*w_y*w_z*varvec[ind6] + (1-w_x)*w_y*w_z*varvec[ind7];
AMREX_FORCE_INLINE AMREX_GPU_HOST_DEVICE int get_single_index(int i, int j, int k, int nx, int ny)
Definition: ERF_Interpolation_Bilinear.H:13
amrex::Real Real
Definition: ERF_ShocInterface.H:16