Hi All! I would like to preset small utility, useful for exploring some PhC devices with meep. It converts 2D PhC device from pseudo graphic file into a list of coordinates of periodic structure. Only hexagonal and square lattices are supported and auto-detected. Source code, linux binary, examples (input, output and resulting epsilon output from meep) and available here: http://zalil.ru/30716530 (or ask me directly by email). Usage of utility is simple as:
$ ./patt2coord.out input-pseudo-graphic-file > output-coordinates-file Any remarks and ideas for future are welcome. ===== As the result i have a question: There is an example of using output-coordinates-file from meep cpp interface in file "test-fc2.cpp" . The final result of its execution is presented on the plot "hex-fc-total-r0-0.1-0.01a.jpg" . I was trying to excite plane wave with Gaussian source from the left of structures (in inserts) and compute transmission spectra from the right. The only changing parameter for 3 curves is hole radius (about zero, 0.016 mkm and 0.1 mkm) of hexagonal PhC . The total flux from the system was plotted (dtf_flux box was set on the edge of PML). If there is no PhC (r == 0) the Gauss is here. But if there is PhC, the result is stange. At some wavelength the output is bigger, then source can give. As far as integrated area is about the same (0.00712673 0.00712576 and 0.00712732) there is some redistribution of energy between frequencies. I hoped that I was setting linear system, but the output is not. Any ideas how to deal with it? WBR, Kosty _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
