Hello,
I'm trying to reproduce the results of the band diagram calculation on
the example page using the c interface.
holey-wvg-bands example:
http://ab-initio.mit.edu/wiki/index.php/Meep_Tutorial/Band_diagram%2C_resonant_modes%2C_and_transmission_in_a_holey_waveguide#Band_diagram
I noticed that my c results are not quite the same as the libctl
results, and I was trying to figure out why.
h5topng showed no obvious differences in dielectric functions, but a
closer inspection using h5diff revealed some significant differences,
apparently near the boundary of the cylindrical hole.
I checked to see if the conditions were possibly incorrect (perhaps I
should use >= and not >) but this did not seem to have an effect.
--------
c snippet:
double resolution = 20; // pixels per distance
const int sy = 12;
const double cellwidth = 1.0;
const vec center = vec((cellwidth / 2),(sy/2));
const double width = 1.2;
const double rad = 0.36;
const double epsr = 13.0;
double eps(const vec &p) {
vec r = p - center;
if (abs(r.y()) <= (width / 2) && abs(r) > rad)
return epsr;
else return 1.0;
}
--------
Best Regards,
Matt
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
