Hi all, I am using mpb to reproduce the results in Andreani-Gerace: "Photonic crystal slabs with a triangular lattice of triangular holes investigated using a guided-mode expansion method", (Physical Review B 73). I therefore wrote the program I attached to this email which uses slabs of the cone class to create triangles.
When outputting epsilon, after mpb's message "Initializing dielectric
function", I keep getting the error message "Speicherzugriffsfehler", which
in German is literally '"memory access failure" and probably means a
segmentation fault. Since it is German I assume that the message does not
come from your program directly but from the system (SuSE Linux 10, Kernel
2.6.13-15.13-smp ([EMAIL PROTECTED]), gcc version 4.0.2) which replies to some
illegal memory access from mpb. Varying the constants "innermaterialheight"
and "outermaterialheight" of my program which control the thickness of the
outer homogeneous and inner photonic lattice, I can get the program to work
but the most interesting combinations cause the crashs and the only
regularity as to what combinations cause crashes and which do not is, that
for low resolutions, small layer sizes and low meshsizes the error does not
occur. I tried to debug the program, which gives me:
Program received signal SIGSEGV, Segmentation fault.
[Switching to Thread 46912514210560 (LWP 26447)]
0x00002aaaaabf94f8 in scm_apply () from /usr/lib64/libguile.so.12
and occurs when the function gh_call1() is fed with my material function in
epsilon.c 103.
The version of guile on the system is 1.6.7-8 and was already installed.
Has anybody experienced similar problems or does anyone have a clue to
what the reason for this behaviour is? I have used mpb a lot since
installing it and did not get any errors, so I do not assume that a wrong
installation is the cause. Do you agree?
Best regards
Julia
tri-triangles_cones.ctl
Description: Binary data
_______________________________________________ mpb-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
