NanoComp / libctl

Guile-based library implementing flexible control files for scientific simulations
GNU General Public License v2.0
18 stars 23 forks source link

Skip empty dimensions in geom.c #43

Closed ChristopherHogan closed 4 years ago

ChristopherHogan commented 4 years ago

Replaces #42.

oskooi commented 4 years ago

This PR is causing python/tests/refl_angular.py of the Meep make check suite to fail:

Using MPI version 3.1, 2 processes
-----------
Initializing structure...
Splitting into 2 chunks evenly
time for choose_chunkdivision = 0.00153303 s
Working in 1D dimensions.
Computational cell is 0 x 0 x 12 with resolution 100
time for set_epsilon = 0.00033456 s
-----------
Meep: using complex fields.
field decay(t = 50.005): 0.2508223760350792 / 0.2508223760350792 = 1.0
field decay(t = 100.01): 2.909213103141057e-15 / 0.2508223760350792 = 1.1598698446003814e-14
run 0 finished at t = 100.01 (20002 timesteps)
-----------
Initializing structure...
Splitting into 2 chunks evenly
BUG in divide_geom_box_tree
time for choose_chunkdivision = 0.000618045 s
Working in 1D dimensions.
Computational cell is 0 x 0 x 12 with resolution 100
     block, center = (0,0,3)
          size (1e+20,1e+20,6)
          axes (1,0,0), (0,1,0), (0,0,1)
          dielectric constant epsilon diagonal = (12.25,12.25,12.25)
BUG in divide_geom_box_tree
-------------------------------------------------------
Primary job  terminated normally, but 1 process returned
a non-zero exit code.. Per user-direction, the job has been aborted.
-------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:

  Process name: [[8489,1],0]
  Exit code:    1
--------------------------------------------------------------------------
FAIL tests/refl_angular.py (exit status: 1)

All the other tests are fine.