123 if (fflag(ii,jj ,k).isRegular() && fflag(ii+1,jj ,k).isRegular() &&
124 fflag(ii,jj+1,k).isRegular() && fflag(ii+1,jj+1,k).isRegular())
126 cflag(i,j,k).setRegular();
127 cvol(i,j,k) = 1.0_rt;
128 ccent(i,j,k,0) = 0.0_rt;
129 ccent(i,j,k,1) = 0.0_rt;
131 cbc(i,j,k,0) = -1.0_rt;
132 cbc(i,j,k,1) = -1.0_rt;
133 cbn(i,j,k,0) = 0.0_rt;
134 cbn(i,j,k,1) = 0.0_rt;
136 else if (fflag(ii,jj ,k).isCovered() && fflag(ii+1,jj ,k).isCovered() &&
137 fflag(ii,jj+1,k).isCovered() && fflag(ii+1,jj+1,k).isCovered())
139 cflag(i,j,k).setCovered();
140 cvol(i,j,k) = 0.0_rt;
141 ccent(i,j,k,0) = 0.0_rt;
142 ccent(i,j,k,1) = 0.0_rt;
144 cbc(i,j,k,0) = -1.0_rt;
145 cbc(i,j,k,1) = -1.0_rt;
146 cbn(i,j,k,0) = 0.0_rt;
147 cbn(i,j,k,1) = 0.0_rt;
151 cflag(i,j,k).setSingleValued();
153 cvol(i,j,k) = 0.25_rt*(fvol(ii,jj ,k) + fvol(ii+1,jj ,k) +
154 fvol(ii,jj+1,k) + fvol(ii+1,jj+1,k));
155 Real cvolinv = 1.0_rt/cvol(i,j,k);
157 ccent(i,j,k,0) = 0.25_rt * cvolinv *
158 (fvol(ii ,jj ,k)*(0.5_rt*fcent(ii ,jj ,k,0)-0.25_rt) +
159 fvol(ii+1,jj ,k)*(0.5_rt*fcent(ii+1,jj ,k,0)+0.25_rt) +
160 fvol(ii ,jj+1,k)*(0.5_rt*fcent(ii ,jj+1,k,0)-0.25_rt) +
161 fvol(ii+1,jj+1,k)*(0.5_rt*fcent(ii+1,jj+1,k,0)+0.25_rt));
162 ccent(i,j,k,1) = 0.25_rt * cvolinv *
163 (fvol(ii ,jj ,k)*(0.5_rt*fcent(ii ,jj ,k,1)-0.25_rt) +
164 fvol(ii+1,jj ,k)*(0.5_rt*fcent(ii+1,jj ,k,1)-0.25_rt) +
165 fvol(ii ,jj+1,k)*(0.5_rt*fcent(ii ,jj+1,k,1)+0.25_rt) +
166 fvol(ii+1,jj+1,k)*(0.5_rt*fcent(ii+1,jj+1,k,1)+0.25_rt));
168 cba(i,j,k) = 0.5_rt*(fba(ii,jj ,k) + fba(ii+1,jj ,k) +
169 fba(ii,jj+1,k) + fba(ii+1,jj+1,k));
170 Real cbainv = 1.0_rt/cba(i,j,k);
172 cbc(i,j,k,0) = 0.5_rt * cbainv *
173 (fba(ii ,jj ,k)*(0.5_rt*fbc(ii ,jj ,k,0)-0.25_rt) +
174 fba(ii+1,jj ,k)*(0.5_rt*fbc(ii+1,jj ,k,0)+0.25_rt) +
175 fba(ii ,jj+1,k)*(0.5_rt*fbc(ii ,jj+1,k,0)-0.25_rt) +
176 fba(ii+1,jj+1,k)*(0.5_rt*fbc(ii+1,jj+1,k,0)+0.25_rt));
177 cbc(i,j,k,1) = 0.5_rt * cbainv *
178 (fba(ii ,jj ,k)*(0.5_rt*fbc(ii ,jj ,k,1)-0.25_rt) +
179 fba(ii+1,jj ,k)*(0.5_rt*fbc(ii+1,jj ,k,1)-0.25_rt) +
180 fba(ii ,jj+1,k)*(0.5_rt*fbc(ii ,jj+1,k,1)+0.25_rt) +
181 fba(ii+1,jj+1,k)*(0.5_rt*fbc(ii+1,jj+1,k,1)+0.25_rt));
183 Real nx = fbn(ii ,jj ,k,0)*fba(ii ,jj ,k)
184 + fbn(ii+1,jj ,k,0)*fba(ii+1,jj ,k)
185 + fbn(ii ,jj+1,k,0)*fba(ii ,jj+1,k)
186 + fbn(ii+1,jj+1,k,0)*fba(ii+1,jj+1,k);
187 Real ny = fbn(ii ,jj ,k,1)*fba(ii ,jj ,k)
188 + fbn(ii+1,jj ,k,1)*fba(ii+1,jj ,k)
189 + fbn(ii ,jj+1,k,1)*fba(ii ,jj+1,k)
190 + fbn(ii+1,jj+1,k,1)*fba(ii+1,jj+1,k);
191 Real nfac = 1.0_rt/std::sqrt(nx*nx+ny*ny+1.e-30_rt);
192 cbn(i,j,k,0) = nx*nfac;
193 cbn(i,j,k,1) = ny*nfac;
194 ierr = (nx == 0.0_rt && ny == 0.0_rt)
197 || ( fapx(ii,jj ,k)==1.0_rt && fapx(ii+2,jj ,k)==1.0_rt
198 && fapx(ii,jj+1,k)==1.0_rt && fapx(ii+2,jj+1,k)==1.0_rt
199 && fapy(ii,jj ,k)==1.0_rt && fapy(ii+1,jj ,k)==1.0_rt
200 && fapy(ii,jj+2,k)==1.0_rt && fapy(ii+1,jj+2,k)==1.0_rt);
205 cvol(i,j,k) = 1.0_rt;
206 ccent(i,j,k,0) = 0.0_rt;
207 ccent(i,j,k,1) = 0.0_rt;
209 cbc(i,j,k,0) = -1.0_rt;
210 cbc(i,j,k,1) = -1.0_rt;
211 cbn(i,j,k,0) = 0.0_rt;
212 cbn(i,j,k,1) = 0.0_rt;
217 capx(i,j,k) = 0.5_rt*(fapx(ii,jj,k) + fapx(ii,jj+1,k));
218 if (capx(i,j,k) != 0.0_rt) {
219 cfcx(i,j,k) = 0.5_rt / capx(i,j,k) *
220 (fapx(ii,jj ,k)*(0.5_rt*ffcx(ii,jj ,k)-0.25_rt) +
221 fapx(ii,jj+1,k)*(0.5_rt*ffcx(ii,jj+1,k)+0.25_rt));
222 if (fecy(ii,jj,k) == Real(1.0) && fecy(ii,jj+1,k) == Real(1.0)) {
223 cecy(i,j,k) = Real(1.0);
225 cecy(i,j,k) = cfcx(i,j,k);
229 cfcx(i,j,k) = 0.0_rt;
230 cecy(i,j,k) = -1.0_rt;
235 capx(i,j,k) = 1.0_rt;
236 cfcx(i,j,k) = 0.0_rt;
237 cecy(i,j,k) = 1.0_rt;
242 capy(i,j,k) = 0.5_rt*(fapy(ii,jj,k) + fapy(ii+1,jj,k));
243 if (capy(i,j,k) != 0.0_rt) {
244 cfcy(i,j,k) = 0.5_rt / capy(i,j,k) *
245 (fapy(ii ,jj,k)*(0.5_rt*ffcy(ii ,jj,k)-0.25_rt) +
246 fapy(ii+1,jj,k)*(0.5_rt*ffcy(ii+1,jj,k)+0.25_rt));
247 if (fecx(ii,jj,k) == Real(1.0) && fecx(ii+1,jj,k) == Real(1.0)) {
248 cecx(i,j,k) = Real(1.0);
250 cecx(i,j,k) = cfcy(i,j,k);
253 cfcy(i,j,k) = 0.0_rt;
254 cecx(i,j,k) = -1.0_rt;
259 capy(i,j,k) = 1.0_rt;
260 cfcy(i,j,k) = 0.0_rt;
261 cecx(i,j,k) = 1.0_rt;
274 auto flg = cflag(i,j,k);
276 if (apx(i ,j ,k) == 0.0_rt) { flg.setDisconnected(-1, 0, 0); }
277 if (apx(i+1,j ,k) == 0.0_rt) { flg.setDisconnected( 1, 0, 0); }
278 if (apy(i ,j ,k) == 0.0_rt) { flg.setDisconnected( 0,-1, 0); }
279 if (apy(i ,j+1,k) == 0.0_rt) { flg.setDisconnected( 0, 1, 0); }
281 if ((apx(i,j ,k) == 0.0_rt || apy(i-1,j,k) == 0.0_rt) &&
282 (apx(i,j-1,k) == 0.0_rt || apy(i ,j,k) == 0.0_rt))
284 flg.setDisconnected(-1,-1,0);
287 if ((apx(i+1,j ,k) == 0.0_rt || apy(i+1,j,k) == 0.0_rt) &&
288 (apx(i+1,j-1,k) == 0.0_rt || apy(i ,j,k) == 0.0_rt))
290 flg.setDisconnected(1,-1,0);
293 if ((apx(i,j ,k) == 0.0_rt || apy(i-1,j+1,k) == 0.0_rt) &&
294 (apx(i,j+1,k) == 0.0_rt || apy(i ,j+1,k) == 0.0_rt))
296 flg.setDisconnected(-1,1,0);
299 if ((apx(i+1,j ,k) == 0.0_rt || apy(i+1,j+1,k) == 0.0_rt) &&
300 (apx(i+1,j+1,k) == 0.0_rt || apy(i ,j+1,k) == 0.0_rt))
302 flg.setDisconnected(1,1,0);
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE int coarsen_from_fine(int i, int j, Box const &bx, int ngrow, Array4< Real > const &cvol, Array4< Real > const &ccent, Array4< Real > const &cba, Array4< Real > const &cbc, Array4< Real > const &cbn, Array4< Real > const &capx, Array4< Real > const &capy, Array4< Real > const &cfcx, Array4< Real > const &cfcy, Array4< Real > const &cecx, Array4< Real > const &cecy, Array4< EBCellFlag > const &cflag, Array4< Real const > const &fvol, Array4< Real const > const &fcent, Array4< Real const > const &fba, Array4< Real const > const &fbc, Array4< Real const > const &fbn, Array4< Real const > const &fapx, Array4< Real const > const &fapy, Array4< Real const > const &ffcx, Array4< Real const > const &ffcy, Array4< Real const > const &fecx, Array4< Real const > const &fecy, Array4< EBCellFlag const > const &fflag)
Definition AMReX_EB2_2D_C.H:96