Skip to content

Commit 2027ac3

Browse files
JamesMcClungJamesMcClung
committed
psc_bnd_fields_impl: restore if debugs
This reverts commit f5765b8.
1 parent f5765b8 commit 2027ac3

1 file changed

Lines changed: 98 additions & 0 deletions

File tree

‎src/libpsc/psc_bnd_fields/psc_bnd_fields_impl.hxx‎

Lines changed: 98 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -143,6 +143,19 @@ struct BndFields_ : BndFieldsBase
143143
Int3 ib = mflds.ib(), im = mflds.im();
144144

145145
if (d == 1) {
146+
#ifdef DEBUG
147+
for (int iz = -2; iz < ldims[2] + 2; iz++) {
148+
for (int ix = MAX(-2, F.ib_[0]);
149+
ix < MIN(ldims[0] + 2, F.ib_[0] + F.im_[0]); ix++) {
150+
fields_t_set_nan(&F(EX, ix, -1, iz));
151+
fields_t_set_nan(&F(EX, ix, -2, iz));
152+
fields_t_set_nan(&F(EY, ix, -1, iz));
153+
fields_t_set_nan(&F(EY, ix, -2, iz));
154+
fields_t_set_nan(&F(EZ, ix, -1, iz));
155+
fields_t_set_nan(&F(EZ, ix, -2, iz));
156+
}
157+
}
158+
#endif
146159
for (int iz = -2; iz < ldims[2] + 2; iz++) {
147160
// FIXME, needs to be for other dir, too, and it's ugly
148161
for (int ix = MAX(-2, ib[0]); ix < MIN(ldims[0] + 2, ib[0] + im[0]);
@@ -157,6 +170,18 @@ struct BndFields_ : BndFieldsBase
157170
}
158171
}
159172
} else if (d == 2) {
173+
#ifdef DEBUG
174+
for (int iy = -2; iy < ldims[1] + 2; iy++) {
175+
for (int ix = -2; ix < ldims[0] + 2; ix++) {
176+
fields_t_set_nan(&F(EX, ix, iy, -1));
177+
fields_t_set_nan(&F(EX, ix, iy, -2));
178+
fields_t_set_nan(&F(EY, ix, iy, -1));
179+
fields_t_set_nan(&F(EY, ix, iy, -2));
180+
fields_t_set_nan(&F(EZ, ix, iy, -1));
181+
fields_t_set_nan(&F(EZ, ix, iy, -2));
182+
}
183+
}
184+
#endif
160185
for (int iy = -2; iy < ldims[1] + 2; iy++) {
161186
for (int ix = MAX(-2, ib[0]); ix < MIN(ldims[0] + 2, ib[0] + im[0]);
162187
ix++) {
@@ -182,6 +207,19 @@ struct BndFields_ : BndFieldsBase
182207

183208
if (d == 1) {
184209
int my _mrc_unused = ldims[1];
210+
#ifdef DEBUG
211+
for (int iz = -2; iz < ldims[2] + 2; iz++) {
212+
for (int ix = MAX(-2, F.ib_[0]);
213+
ix < MIN(ldims[0] + 2, F.ib_[0] + F.im_[0]); ix++) {
214+
fields_t_set_nan(&F(EX, ix, my, iz));
215+
fields_t_set_nan(&F(EX, ix, my + 1, iz));
216+
fields_t_set_nan(&F(EY, ix, my, iz));
217+
fields_t_set_nan(&F(EY, ix, my + 1, iz));
218+
fields_t_set_nan(&F(EZ, ix, my, iz));
219+
fields_t_set_nan(&F(EZ, ix, my + 1, iz));
220+
}
221+
}
222+
#endif
185223
for (int iz = -2; iz < ldims[2] + 2; iz++) {
186224
for (int ix = MAX(-2, ib[0]); ix < MIN(ldims[0] + 2, ib[0] + im[0]);
187225
ix++) {
@@ -196,6 +234,18 @@ struct BndFields_ : BndFieldsBase
196234
}
197235
} else if (d == 2) {
198236
int mz = ldims[2];
237+
#ifdef DEBUG
238+
for (int iy = -2; iy < ldims[1] + 2; iy++) {
239+
for (int ix = -2; ix < ldims[0] + 2; ix++) {
240+
fields_t_set_nan(&F(EX, ix, iy, mz));
241+
fields_t_set_nan(&F(EX, ix, iy, mz + 1));
242+
fields_t_set_nan(&F(EY, ix, iy, mz));
243+
fields_t_set_nan(&F(EY, ix, iy, mz + 1));
244+
fields_t_set_nan(&F(EZ, ix, iy, mz));
245+
fields_t_set_nan(&F(EZ, ix, iy, mz + 1));
246+
}
247+
}
248+
#endif
199249
for (int iy = -2; iy < ldims[1] + 2; iy++) {
200250
for (int ix = MAX(-2, ib[0]); ix < MIN(ldims[0] + 2, ib[0] + im[0]);
201251
ix++) {
@@ -220,6 +270,19 @@ struct BndFields_ : BndFieldsBase
220270
Int3 ib = mflds.ib(), im = mflds.im();
221271

222272
if (d == 1) {
273+
#ifdef DEBUG
274+
for (int iz = -2; iz < ldims[2] + 2; iz++) {
275+
for (int ix = MAX(-2, F.ib_[0]);
276+
ix < MIN(ldims[0] + 2, F.ib_[0] + F.im_[0]); ix++) {
277+
fields_t_set_nan(&F(HX, ix, -1, iz));
278+
fields_t_set_nan(&F(HX, ix, -2, iz));
279+
fields_t_set_nan(&F(HY, ix, -1, iz));
280+
fields_t_set_nan(&F(HY, ix, -2, iz));
281+
fields_t_set_nan(&F(HZ, ix, -1, iz));
282+
fields_t_set_nan(&F(HZ, ix, -2, iz));
283+
}
284+
}
285+
#endif
223286
for (int iz = -1; iz < ldims[2] + 2; iz++) {
224287
for (int ix = MAX(-2, ib[0]); ix < MIN(ldims[0] + 2, ib[0] + im[0]);
225288
ix++) {
@@ -231,6 +294,18 @@ struct BndFields_ : BndFieldsBase
231294
}
232295
}
233296
} else if (d == 2) {
297+
#ifdef DEBUG
298+
for (int iy = -2; iy < ldims[1] + 2; iy++) {
299+
for (int ix = -2; ix < ldims[0] + 2; ix++) {
300+
fields_t_set_nan(&F(HX, ix, iy, -1));
301+
fields_t_set_nan(&F(HX, ix, iy, -2));
302+
fields_t_set_nan(&F(HY, ix, iy, -1));
303+
fields_t_set_nan(&F(HY, ix, iy, -2));
304+
fields_t_set_nan(&F(HZ, ix, iy, -1));
305+
fields_t_set_nan(&F(HZ, ix, iy, -2));
306+
}
307+
}
308+
#endif
234309
for (int iy = -2; iy < ldims[1] + 2; iy++) {
235310
for (int ix = MAX(-2, ib[0]); ix < MIN(ldims[0] + 2, ib[0] + im[0]);
236311
ix++) {
@@ -255,6 +330,18 @@ struct BndFields_ : BndFieldsBase
255330

256331
if (d == 1) {
257332
int my _mrc_unused = ldims[1];
333+
#ifdef DEBUG
334+
for (int iz = -2; iz < ldims[2] + 2; iz++) {
335+
for (int ix = MAX(-2, F.ib_[0]);
336+
ix < MIN(ldims[0] + 2, F.ib_[0] + F.im_[0]); ix++) {
337+
fields_t_set_nan(&F(HX, ix, my, iz));
338+
fields_t_set_nan(&F(HX, ix, my + 1, iz));
339+
fields_t_set_nan(&F(HY, ix, my + 1, iz));
340+
fields_t_set_nan(&F(HZ, ix, my, iz));
341+
fields_t_set_nan(&F(HZ, ix, my + 1, iz));
342+
}
343+
}
344+
#endif
258345
for (int iz = -2; iz < ldims[2] + 2; iz++) {
259346
for (int ix = MAX(-2, ib[0]); ix < MIN(ldims[0] + 2, ib[0] + im[0]);
260347
ix++) {
@@ -267,6 +354,17 @@ struct BndFields_ : BndFieldsBase
267354
}
268355
} else if (d == 2) {
269356
int mz = ldims[2];
357+
#ifdef DEBUG
358+
for (int iy = -2; iy < ldims[1] + 2; iy++) {
359+
for (int ix = -2; ix < ldims[0] + 2; ix++) {
360+
fields_t_set_nan(&F(HX, ix, iy, mz));
361+
fields_t_set_nan(&F(HX, ix, iy, mz + 1));
362+
fields_t_set_nan(&F(HY, ix, iy, mz));
363+
fields_t_set_nan(&F(HY, ix, iy, mz + 1));
364+
fields_t_set_nan(&F(HZ, ix, iy, mz + 1));
365+
}
366+
}
367+
#endif
270368
for (int iy = -2; iy < ldims[1] + 2; iy++) {
271369
for (int ix = MAX(-2, ib[0]); ix < MIN(ldims[0] + 2, ib[0] + im[0]);
272370
ix++) {

0 commit comments

Comments
 (0)