Hi all,
I figured out the lapack compile problems finally. Now everything
compiles fine, and runs well... except I never get any output from
harminv.
For example, if I run the ring example from the tutorial on the super
computer either with meep, meep-mpi (with mpirun) or meep (with mpirun
-nonmpi), I get this output:
-----------
Initializing structure...
Working in 2D dimensions.
cylinder, center = (0,0,0)
radius 2, height 1e+20, axis (0, 0, 1)
dielectric constant epsilon = 11.56
cylinder, center = (0,0,0)
radius 1, height 1e+20, axis (0, 0, 1)
dielectric constant epsilon = 1
time for set_epsilon = 0.058553 s
-----------
creating output file "./ring-eps-000000.00.h5"...
harminv0:, frequency, imag. freq., Q, |amp|, amplitude, error
run 0 finished at t = 300.0 (6000 timesteps)
Elapsed run time = 2.35695 s
As you know, I should get this output (result obtained on laptop, using
gcc throughout):
-----------
Initializing structure...
Working in 2D dimensions.
cylinder, center = (0,0,0)
radius 2, height 1e+20, axis (0, 0, 1)
dielectric constant epsilon = 11.56
cylinder, center = (0,0,0)
radius 1, height 1e+20, axis (0, 0, 1)
dielectric constant epsilon = 1
time for set_epsilon = 0.295364 s
-----------
creating output file "./ring-eps-000000.00.h5"...
Meep progress: 290.65/300.0 = 96.9% done in 4.0s, 0.1s to go
on time step 5823 (time=291.15), 0.000686978 s/step
harminv0:, frequency, imag. freq., Q, |amp|, amplitude, error
harminv0:, 0.118450442650229, -6.54478181615834e-4, 90.4922776476582,
0.00323035195862018, -0.00237726034205063-0.00218719158801294i,
3.268691961116e-4
harminv0:, 0.147154094816768, -3.90827000142487e-4, 188.259888343333,
0.0315632941394842, 0.0197993916076897-0.0245810013811672i,
2.61477140061881e-6
harminv0:, 0.17524196913769, -6.54604672407747e-5, 1338.53283152655,
0.00730271007900924, -2.49573057274589e-4-0.00729844420319466i,
1.90523431775226e-6
run 0 finished at t = 300.0 (6000 timesteps)
Elapsed run time = 5.10775 s
Is there a way to debug harminv? Thanks.
Nic
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
