@@ -143,19 +143,6 @@ 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
159146 for (int iz = -2 ; iz < ldims[2 ] + 2 ; iz++) {
160147 // FIXME, needs to be for other dir, too, and it's ugly
161148 for (int ix = MAX (-2 , ib[0 ]); ix < MIN (ldims[0 ] + 2 , ib[0 ] + im[0 ]);
@@ -170,18 +157,6 @@ struct BndFields_ : BndFieldsBase
170157 }
171158 }
172159 } 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
185160 for (int iy = -2 ; iy < ldims[1 ] + 2 ; iy++) {
186161 for (int ix = MAX (-2 , ib[0 ]); ix < MIN (ldims[0 ] + 2 , ib[0 ] + im[0 ]);
187162 ix++) {
@@ -207,19 +182,6 @@ struct BndFields_ : BndFieldsBase
207182
208183 if (d == 1 ) {
209184 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
223185 for (int iz = -2 ; iz < ldims[2 ] + 2 ; iz++) {
224186 for (int ix = MAX (-2 , ib[0 ]); ix < MIN (ldims[0 ] + 2 , ib[0 ] + im[0 ]);
225187 ix++) {
@@ -234,18 +196,6 @@ struct BndFields_ : BndFieldsBase
234196 }
235197 } else if (d == 2 ) {
236198 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
249199 for (int iy = -2 ; iy < ldims[1 ] + 2 ; iy++) {
250200 for (int ix = MAX (-2 , ib[0 ]); ix < MIN (ldims[0 ] + 2 , ib[0 ] + im[0 ]);
251201 ix++) {
@@ -270,19 +220,6 @@ struct BndFields_ : BndFieldsBase
270220 Int3 ib = mflds.ib (), im = mflds.im ();
271221
272222 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
286223 for (int iz = -1 ; iz < ldims[2 ] + 2 ; iz++) {
287224 for (int ix = MAX (-2 , ib[0 ]); ix < MIN (ldims[0 ] + 2 , ib[0 ] + im[0 ]);
288225 ix++) {
@@ -294,18 +231,6 @@ struct BndFields_ : BndFieldsBase
294231 }
295232 }
296233 } 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
309234 for (int iy = -2 ; iy < ldims[1 ] + 2 ; iy++) {
310235 for (int ix = MAX (-2 , ib[0 ]); ix < MIN (ldims[0 ] + 2 , ib[0 ] + im[0 ]);
311236 ix++) {
@@ -330,18 +255,6 @@ struct BndFields_ : BndFieldsBase
330255
331256 if (d == 1 ) {
332257 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
345258 for (int iz = -2 ; iz < ldims[2 ] + 2 ; iz++) {
346259 for (int ix = MAX (-2 , ib[0 ]); ix < MIN (ldims[0 ] + 2 , ib[0 ] + im[0 ]);
347260 ix++) {
@@ -354,17 +267,6 @@ struct BndFields_ : BndFieldsBase
354267 }
355268 } else if (d == 2 ) {
356269 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
368270 for (int iy = -2 ; iy < ldims[1 ] + 2 ; iy++) {
369271 for (int ix = MAX (-2 , ib[0 ]); ix < MIN (ldims[0 ] + 2 , ib[0 ] + im[0 ]);
370272 ix++) {
0 commit comments