Hi,
I guess it's not the correct place to report a bug in harminv here but I didn't
find any harminv dedicated mailing list so here I post.
After compiling meep with harminv on debian testing/powerpc, I do a make check
which fail like this:
make[2]: entrant dans le repertoire /home/loic/meep/meep-0.10/tests
./convergence_cyl_waveguide.dac
Testing convergence of a waveguide mode frequency...
harminv: failure on line 853 of harminv.c: argument out of range in
harminv_get_amplitude
make[2]: *** [convergence_cyl_waveguide.done] Erreur 1
make[2]: quittant le repertoire /home/loic/meep/meep-0.10/tests
make[1]: *** [check-am] Erreur 2
make[1]: quittant le repertoire /home/loic/meep/meep-0.10/tests
make: *** [check-recursive] Erreur 1
If I run this test with gdb, I got:
Starting program:
/home/loic/meep/meep-0.10/tests/convergence_cyl_waveguide.dac [Thread
debugging using libthread_db enabled]
[New Thread 805421152 (LWP 20454)]
Testing convergence of a waveguide mode frequency...
Program received signal SIGSEGV, Segmentation fault.
[Switching to Thread 805421152 (LWP 20454)]
0x1007713c in harminv_get_amplitude (d=0x7fb94938, k=269113576)
at harminv.c:856
856 return d->amps[k];
Current language: auto; currently c
(gdb) bt
#0 0x1007713c in harminv_get_amplitude (d=0x7fb94938, k=269113576)
at harminv.c:856
#1 0x10005814 in meep::do_harminv (data=0x1009bc48, n=2499,
dt=0.050000000000000003, fmin=0.10000000000000001, fmax=0.5,
maxbands=10, amps=0x1009a3d8, freq_re=0x1009a480,
freq_im=0x1009a4d8, errors=0x100a5890,
spectral_density=1.1000000000000001, Q_thresh=50,
rel_err_thresh=1e+20, err_thresh=0.01, rel_amp_thresh=-1, amp_thresh=-1)
at bands.cpp:469
#2 0x10003f40 in test_convergence_without_averaging ()
at convergence_cyl_waveguide.cpp:61
#3 0x10004660 in main (argc=1, argv=0x7fb95814)
at convergence_cyl_waveguide.cpp:155
This problem also append in the example material-dispersion.ctl
[EMAIL PROTECTED] ~/bin$ gdb meep
GNU gdb 6.4.90-debian
Copyright (C) 2006 Free Software Foundation, Inc.
GDB is free software, covered by the GNU General Public License, and you are
welcome to change it and/or distribute copies of it under certain conditions.
Type "show copying" to see the conditions.
There is absolutely no warranty for GDB. Type "show warranty" for details.
This GDB was configured as "powerpc-linux-gnu"...Using host
libthread_db library "/lib/tls/libthread_db.so.1".
(gdb) set args ~/meep/meep-0.10/examples/material-dispersion.ctl
(gdb) run
Starting program: /home/loic/bin/meep
~/meep/meep-0.10/examples/material-dispersion.ctl
[Thread debugging using libthread_db enabled]
[New Thread 805426976 (LWP 20477)]
-----------
Initializing structure...
Working in 2D dimensions.
time for set_epsilon = 0.00037384 s
polarizability: omega=0.5, gamma=0.1, deps=2e-05, esat=0
polarizability: omega=1.1, gamma=1e-05, deps=0.5, esat=0
-----------
creating output file "./material-dispersion-eps-000000.00.h5"...
Field time usage:
outputting fields: 0.151938 s
everything else: 0.00169206 s
harminv0:, frequency, imag. freq., Q, |amp|, amplitude, error
harminv: failure on line 853 of harminv.c: argument out of range in
harminv_get_amplitude
Program exited with code 01.
I don't get any problem on i386.
Hope this help
Cheers.
--
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
