[Meep-discuss] Meep installation notes for macOS Sierra, July 2017

2017-07-22 Thread Filip Dominec
[On behalf of John Weiner, I forward his experience with installation of MEEP on Mac OS. John's kind contribution follows. -- Filip] Dear Meep Users and Developers, I have been able to install meep-1.3 on two different iMacs (one in Washington DC( Sierra 10.12.5) and the other in Paris (Sierra 10

Re: [Meep-discuss] TE polarized photonic band structure of metallic photonic crystals

2017-01-19 Thread Filip Dominec
Hi, would you please describe to all what is the reported issue with the TE mode? Filip 2017-01-19 11:00 GMT+01:00, Aliasghar Sedghi : > Dear meep users, > I have calculated the band structure of PCs composed of metallic rods > (Drude model) embedded in an air background using meep. There is no >

Re: [Meep-discuss] Compiling Meep 1.3 on Debian 8.5 (uniprocessor, not MPI); First three tutorials

2017-01-08 Thread Filip Dominec
. I should be able to run the interactive > commands at the start of the tutorial: > > http://ab-initio.mit.edu/wiki/index.php/Meep_Tutorial > > - - - - > > (set! geometry-lattice (make lattice (size 16 8 no-size))) > > (set! geometry (list > (ma

Re: [Meep-discuss] trying to find charge density of metamaterials

2016-10-26 Thread Filip Dominec
Hi, assuming you are able to compute a time-domain simulation of the metamaterial unit cell, there should be no big problem to 1) resolve the resonance frequencies of each mode, then to 2) switch to a narrow-band source with the actual frequency of the mode you are interested in, 3) record the elec

Re: [Meep-discuss] meep-discuss Digest, Vol 127, Issue 1

2016-10-07 Thread Filip Dominec
Excuse my little advertisement, but here I have written script that compiles the freshest version of MEEP on your Ubuntu system: https://github.com/FilipDominec/python-meep-install/blob/master/python-meep-install.sh It was tested on many systems and it should work out-of-the-box also on your versi

Re: [Meep-discuss] Question about permittivity and permeability

2016-09-19 Thread Filip Dominec
Hi, since MEEP uses relative permittivity and permeability which have no physical dimension, there is no need to convert units. Filip 2016-09-16 22:43 GMT+02:00, Heligman, Daniel M. : > Hi, > > > I have a simple question. Is the a conversion between SI and meep units for > permittivity and perm

Re: [Meep-discuss] MEEP slow after Ubuntu 16.04 upgrade

2016-09-01 Thread Filip Dominec
time. > Other computational tools are running normally, for example matlab, I have > not seen any performance drop after the Ubuntu upgrade. > Best regards, > Jisha > > > On Thu, Sep 1, 2016 at 2:13 PM, Filip Dominec > wrote: > >> Hi, >> 1) Do you use multiprocessi

Re: [Meep-discuss] MEEP slow after Ubuntu 16.04 upgrade

2016-09-01 Thread Filip Dominec
Hi, 1) Do you use multiprocessing? If so, which implementation, MPICH or OpenMPI? 2) Do you use HDF5 export? If so, which kind of HDF libraries, MPI-enabled or "serial"? The latter is known to be very slow when used in multiprocessing setup, though they should not alter subpixel averaging. 3) Did

Re: [Meep-discuss] LDOS calculation in Python-Meep

2016-08-31 Thread Filip Dominec
omething wrong?. Is there > any workaround to get the values of the ldos calculation in Python-Meep or > Scheme or C++ interface for MEEP are needed to get them?. Anyone tried this > before?. > > > Thank you in advanced for any help! > > Eduardo >

Re: [Meep-discuss] Plane Wave Source

2016-08-26 Thread Filip Dominec
Dear Craig, if I do not miss something, you may easily define a (perpendicular or even oblique) plane wave source as shown in https://www.mail-archive.com/meep-discuss@ab-initio.mit.edu/msg03088.html F. 2016-08-26 7:51 GMT+02:00, Craig Milde : > Hi meep discussion list, > > > I am attempting to s

Re: [Meep-discuss] MEEP output interpretation

2016-08-11 Thread Filip Dominec
Dear Craig, from your e-mail it is not clear what steps you actually tried to extract the field data from the HDF file. Can you run, e.g., this? h5tovtk input.h5 -o output.vtk paraview output.vtk You should be able to add some monitor plane in paraview and browse the field records as you wish. Fi

Re: [Meep-discuss] 3D Plane Wave source

2016-07-28 Thread Filip Dominec
Dear Saket, just define the source amplitude as a harmonic function, with the periodicity corresponding by the desired transverse component of the wave vector. Filip 2016-07-27 21:03 GMT+02:00, saket wankhade : > Dear meep users, > I need a tilted plane wave source with both E

Re: [Meep-discuss] Getting grid points/time points of meep simulation

2016-06-29 Thread Filip Dominec
Maybe it helps somebody if I mention that the python-meep interface solves the previous questions, since 1) each individual time step is accessible within a cycle (see https://github.com/FilipDominec/python-meep-utils/blob/master/scatter.py#L88), and 2) for every point in the grid, MEEP calls a

Re: [Meep-discuss] Compuation of electric field

2016-06-05 Thread Filip Dominec
Hi, I would rather guess the images Saket gets are really full black, which would mean that there are "NaN" values in the output of the computation. This may result from some numerical instability. I suggest to try first to halve the simulation time until the field gets finite. Then one can det

Re: [Meep-discuss] numerical instability in vacuum - no materials, no pml

2016-03-04 Thread Filip Dominec
Hi, whenever you set the phase in either the Y- or Z- symmetry to -1, the simulation eventually becomes unstable. I have tested out all four combinations for you, and I can tell that only (set! symmetries (list (make mirror-sym (direction Y) (phase 1)) (make mirror-sym (direction Z) (phase 1)))

Re: [Meep-discuss] Mesh size

2016-02-10 Thread Filip Dominec
Hi, it is directly in http://ab-initio.mit.edu/wiki/index.php/Meep_Reference, search for "resolution". The higher, the finer mesh. Filip On 11.2.2016 06:07, saket wankhade wrote: Dear meep users, I know meep uses uniform meshing but want to change size of it.Is it possible to change m

Re: [Meep-discuss] doubt regarding subpixel averaging behaviour

2015-09-22 Thread Filip Dominec
works and for 8 processors it > doesn't appear, so I'm guessing it is not working. I couldn't check if > it's affecting my results yet. > > Pablo. > > El 21/09/15 a las 19:36, Filip Dominec escribió: >> Dear Pablo, how exactly does the bug manifest? Can it

Re: [Meep-discuss] doubt regarding subpixel averaging behaviour

2015-09-21 Thread Filip Dominec
Dear Pablo, how exactly does the bug manifest? Can it be observed on the exported permittivity, e.g. in Paraview? Or does the simulation lose its sub-pixel precision with regards to some structure parameter? Filip 2015-09-21 19:28 GMT+02:00, Pablo Díaz Núñez : > Dear meep users, > > I have encount

Re: [Meep-discuss] diffraction limit

2015-09-15 Thread Filip Dominec
Dear Amin, I am afraid your question should be better searched/posed e.g. in http://physics.stackexchange.com/ rather than to the MEEP community. Your question has multiple answers depending on the point of view. 1) The most straightforward is that Maxwell equations have some solutions and extreme

Re: [Meep-discuss] Can source amplitude be changed?

2015-09-13 Thread Filip Dominec
Dear Steve, 1) before we involve in changing the source amplitude: Analysing the high-quality cavity ringdown by time-domain method running the whole time is *extremely* wasteful. In fact, the "harminv", or "filter diagonalization method" was designed to solve exactly this problem: from a short sa

Re: [Meep-discuss] Source code: Nonlinear polarization

2015-07-22 Thread Filip Dominec
Dear Vladimir, search for the variable 'chi2' in the code, it is used in multiple files. The actual computation of its physical effect should be in step_generic.cpp Filip 2015-07-22 23:05 GMT+02:00, Владимир Борисович Новиков : > Dear all, > Please show code place in https://github.com/stevengj/me

Re: [Meep-discuss] Reconfigurable medium meep computation

2015-06-18 Thread Filip Dominec
Dear Nicolas, I guess the optical and mechanical processes happen at completely different time scales (picoseconds to microseconds at least). Perhaps the simplest way is to re-iterate multiple separate simulations, with each simulation's structure being defined by the result of the previous simulat

Re: [Meep-discuss] bragg structure with kerr effect ....

2015-06-15 Thread Filip Dominec
Hi, Alshikh, from what I understand, you shall introduce another callback function like your eps_bragg(const vec &pt), which will return the value of the Kerr coefficient. It should be called in a similar manner as e.g. the susceptibilities are. I always found it useful to test out the Kerr nonli

Re: [Meep-discuss] Compiling Meep 1.3 on Debian 8 (testing) - HDF5 library troubles

2015-06-05 Thread Filip Dominec
Dear all, in 2015, there have been numerous reports of MEEP not compiling with HDF5 on Debian Jessie and newer versions of Ubuntu [1,2]. This happened due to a known change in the system [3]. After a longer trial-and-error, I found the required change in the LDFLAGS variable so that I can compile

Re: [Meep-discuss] missing _H5P_CLS_DATASET_CREATE_g

2015-06-02 Thread Filip Dominec
Dear all, I am, too experimenting with this issue at the moment. You can clone my repository with the installation scripts for the latest python-meep from: https://github.com/FilipDominec/python-meep-install. The advantage of using git is in clarity: Every modification of the installation procedur

Re: [Meep-discuss] Precompiled package for meep 1.3

2015-06-01 Thread Filip Dominec
Dear Vladimir and all interested, to my knowledge there is no precompiled MEEP of the version 1.3 for Debian-based systems. Version 1.2.1 has been readily available for over a year thanks to Thorsten Alteholz. I prepared a compilation script, and planned to package the very latest MEEP and Python-

Re: [Meep-discuss] FW: Si material parameters for MEEP C++

2015-05-11 Thread Filip Dominec
Dear Shin-Woong, I did not compile your code, but I have few thoughts: 1) First, the material should be tested on a half-space, not a slab. The latter combines "too much physics", i.e. the Lorentz-Drude model with a pronounced Fabry-Pérot interference (which is moreover modulated by the dispersion

Re: [Meep-discuss] flux spectrum of continuous source

2015-05-05 Thread Filip Dominec
ds or something, why do the transients take so long to die out.? >> >> The attached repo contains a _test.ctl file with my code and graphs in the >> _test-out/ folder. >> >> https://github.com/eoinmurray/meep-test.git >> <https://github.com/eoinmurray/meep-test.g

[Meep-discuss] Fwd: Plane wave in meep (C++)

2015-05-02 Thread Filip Dominec
would like to have it propagate in only one particular > direction (think of an experiment where you sign a laser light on a metallic > nanoparticle). > > > > Thanks for your time and help. > > Amrit > > > >> On May 1, 2015, at 3:02 PM, Filip Dominec >> wro

Re: [Meep-discuss] Plane wave in meep (C++)

2015-05-01 Thread Filip Dominec
Dear Amrit, I deduce you are interested in _oblique_ incidence of the plane wave. Perhaps the optimum solution would be to freely rotate the source object in space; however to my knowledge this is not possible. It is however relatively easy to make a source emit with a transverse component of the

Re: [Meep-discuss] Confusions

2015-04-27 Thread Filip Dominec
Dear Hamza, from what I know you need to run a loop of the simulations; this is probably accomplished the best through a loop in a bash script; the scheme interface accepts command-line parameters as can be found in numerous examples. As an example in python-meep, I just uploaded https://github.com

Re: [Meep-discuss] CAD tool for Meep

2015-04-21 Thread Filip Dominec
>From my personal point of view, MEEP is mostly used for scientific research rather than engineering. Most simulations therefore deal with complicated physical interpretation of rather simple structures. The highlights of MEEP are the possibility to defive various physical phenomena and its scripti

Re: [Meep-discuss] Python-Meep Investigation

2015-04-02 Thread Filip Dominec
re, with its full log output on the current cleanly installed system, is at: https://github.com/FilipDominec/python-meep-install Filip 2015-03-25 4:05 GMT-04:00, Filip Dominec : > Hi, > it is as you write, and therefore I would like to make the blue-boxed > installation failsafe on as

Re: [Meep-discuss] Resonant Modes in Microdisks

2015-04-02 Thread Filip Dominec
Dear Pierre, I observed similar mode splitting (from intended circular symmetric structures) when my simulations had too coarse a resolution. How the results change if you halve or double the resolution? You could also launch the frequency-domain solver to identify the detected two modes. Filip

Re: [Meep-discuss] Resonant Modes in Microdisks

2015-04-02 Thread Filip Dominec
Dear Pierre, I observed similar mode splitting (from intended circular symmetric structures) when my simulations had too coarse a resolution. How the results change if you halve or double the resolution? You could also launch the frequency-domain solver to identify the detected two modes. Filip

Re: [Meep-discuss] Low frequency modulation

2015-03-31 Thread Filip Dominec
ure > a smooth startup.I am about to test that. > > I did try to use the frequency domain solver, but it ran for hours and > failed to converge. I don't know what that signifies. > > Cheers > Gib > > From: Filip Dominec

Re: [Meep-discuss] Low frequency modulation

2015-03-31 Thread Filip Dominec
Hi, Gib, it appears there are two resonant modes relatively close in frequency (separated by 1/4 of their center frequency, as you wrote). You may wish to set the 'width' to some reasonable value so that the spectrum of the source is narrower. http://ab-initio.mit.edu/wiki/index.php/Meep_Reference

Re: [Meep-discuss] Numerical dispersion?

2015-03-30 Thread Filip Dominec
Hi, Glenn, in my experience, the speed of light and generally the results in MEEP are quite accurate (compared to theory, other methods or experiment) until the wavelength becomes comparable to several voxels - in which case however also the structure ceases to be described correctly. Thus I believ

Re: [Meep-discuss] What is "strides" in Meep sources?

2015-03-26 Thread Filip Dominec
case NO_DIRECTION: break; // no stride here, either > } > } > > 2015-03-26 20:07 GMT+03:00 Filip Dominec : >> As far as I understood the source, the 'stride' represents how much >> bytes in the computer memory shall be skipped to get to the voxel in >&g

Re: [Meep-discuss] What is "strides" in Meep sources?

2015-03-26 Thread Filip Dominec
As far as I understood the source, the 'stride' represents how much bytes in the computer memory shall be skipped to get to the voxel in the neighboring row. More readable, but perhaps less efficient, implementation would be to address the voxels as a 2-D or 3-D array. Is this correct? Filip 2015-

Re: [Meep-discuss] Calculation of response time of a device in meep

2015-03-26 Thread Filip Dominec
uide. It is like a Y shaped PC > waveguide. Now I want to calculate , how much time the device will take > responds to the input signals ? > > > Nirmala > > On Thu, Mar 26, 2015 at 12:35 PM, Filip Dominec > wrote: > >> Hi, Nimala, >> perhaps it is possib

Re: [Meep-discuss] Calculation of response time of a device in meep

2015-03-26 Thread Filip Dominec
Hi, Nimala, perhaps it is possible. Could you specify your task in more specific words, please? The custom source can be easily found on the net, e.g. http://comments.gmane.org/gmane.comp.science.electromagnetism.meep.general/3936 Filip 2015-03-26 3:57 GMT+01:00, nirmala dsouza : > I want to calcu

Re: [Meep-discuss] Meep can't find HDF5 libraries during ./configure

2015-03-25 Thread Filip Dominec
b/x86_64-linux-gnu/libhdf5_cpp.a > /usr/lib/x86_64-linux-gnu/libhdf5_openmpihl_fortran.so.8 > /usr/lib/x86_64-linux-gnu/libhdf5_serial_fortran.a > /usr/lib/jni/libjhdf5.so > > I also tried to add the libraries separately with: > export > LD_RUN_PATH="/usr/local/lib/:/usr/li

Re: [Meep-discuss] Frequency-domain solver

2015-03-19 Thread Filip Dominec
arying the frequency. It would be great > if you could point me to an example .ctl file that executes multiple > simulations, creating file names appropriately. > > Thanks > Gib > ________ > From: Filip Dominec [filip.domi...@gmail.com] > Sen

Re: [Meep-discuss] Frequency-domain solver

2015-03-18 Thread Filip Dominec
evengj/meep > > Gib > ________ > From: Filip Dominec [filip.domi...@gmail.com] > Sent: Thursday, 19 March 2015 12:05 p.m. > To: Gib Bogle > Cc: meep-discuss > Subject: Re: [Meep-discuss] Frequency-domain solver > > Hi, Gib, perhaps I was a bit confused. The fact is

Re: [Meep-discuss] Frequency-domain solver

2015-03-18 Thread Filip Dominec
1.2.1.tar.gz from > http://ab-initio.mit.edu/wiki/index.php/Meep_download in January this year. > This claims to be the latest release - is this not the case? > > Cheers > Gib > ________ > From: Filip Dominec [filip.domi...@gmail.com] > Sen

Re: [Meep-discuss] Frequency-domain solver

2015-03-18 Thread Filip Dominec
Hi, Gib, it looks exactly like the bug that has been resolved in git commit https://github.com/stevengj/meep/commit/0640e9f6a2ccfa2a40375a89b5692ebb1fb87efb Could you check a recent version of MEEP either compiling it from source (blue box in http://f.dominec.eu/meep/index.html#install), or insta

Re: [Meep-discuss] h5ls: unable to open file. python -meep

2015-03-13 Thread Filip Dominec
Hi, Julián, your script works for me. I get the same result from 1- and 2-process simulation. I have h5utils in version 1.12.1, and freshly-compiled meep 1.2.1 & python-meep according to the blue box at http://f.dominec.eu/meep/index.html#install Could you somehow examine the sample.h5 file? Its

Re: [Meep-discuss] Seeking advice for measuring transmission spectra

2015-02-18 Thread Filip Dominec
Dear John, I can not confirm not deny your approach is correct. I use another way to get the incident/transmitted/reflected fields from a structure in free space. Most importantly, it provides also the phase information, and also allows to get all the reflection and transmission coefficients in si

Re: [Meep-discuss] help for calculation of transmission and reflection of photnic crystals slabs

2015-02-17 Thread Filip Dominec
Hi, Saeid, your task has been perhaps done by somebody else, please take a look at https://timsaucer.wordpress.com/2013/04/30/building-l3-photonic-crystal-cavity-in-meep-1-finding-mode/, http://comments.gmane.org/gmane.comp.science.electromagnetism.meep.general/10, http://osdir.com/ml/science.elect

Re: [Meep-discuss] Error while using perfect-metal

2015-02-11 Thread Filip Dominec
Dear Jim, this seems to be a result of instability, resulting in saturated (i.e. infinite or NaN) values as you suggest. To verify this, you may shorten the simulation until you get some image that is not plain black. As a matter of fact, the instability does not occur instantly, but it takes it se

Re: [Meep-discuss] Gaissian Distribution of Phase

2015-02-10 Thread Filip Dominec
Hi, currently I do not have the example you mention, but surely you can make a focused beam from an oblique beam, by introducing a *quadratic* term into the imaginary part of the exponential that defines the source amplitude. If you wish to build a focused *gaussian* beam, you will also wish to mod

Re: [Meep-discuss] Puzzling results for reflection spectrum

2015-02-10 Thread Filip Dominec
Hi, if you properly divide the reflected amplitude by the incident amplitude, the reflection of a linear system should indeed be the same, independent of the source duration and spectrum. If it is not, one possible cause may be that you clip the time record while the source was still running, or th

Re: [Meep-discuss] how to define a invariable complex refractive index?

2015-02-10 Thread Filip Dominec
Dear Liu, from the principle of FDTD computation, there is no other way how to define complex refractive index than to include dispersion. Either you may introduce a Drude term (conductivity), or some high-frequency Lorentzian. More details are here: http://ab-initio.mit.edu/wiki/index.php/Material

Re: [Meep-discuss] meep and python-meep on OS-X

2015-02-10 Thread Filip Dominec
Thank you, Glenn, great job. I will be happy to add a link to your post! Filip ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss

Re: [Meep-discuss] Slow calculation with meep parallelisation

2014-12-02 Thread Filip Dominec
, I tried multiprocessing on single node too. And the problem persists. > Is it a common problem with meep multiprocessing? > > Thanks and Regards, > Ratnesh > > -Original Message- > From: Filip Dominec [mailto:filip.domi...@gmail.com] > Sent: Tuesday, December 02, 2014 3:41 P

Re: [Meep-discuss] Scaling - dielectric-epsilon, metal-copper, lattice/resolution and Symmetry

2014-11-25 Thread Filip Dominec
If my current run is successful I want to calculate forces > around the outside of the modeled resonate cavity. I am hopeful that > Meep can detect forces caused by evanescent wave momentum and if so, I > will be a long way ahead of game. I appreciate your continued support. > > S

Re: [Meep-discuss] Scaling - dielectric-epsilon, metal-copper, lattice/resolution and Symmetry

2014-11-24 Thread Filip Dominec
Dear Steve, 1) Use the value of 2.3. The static permittivity is invariant to scaling. 2) If you wanted to model realistic copper, with scaling of 1:1e6, you would have to feed MEEP with values normalized to the scaling (1e6) and also by the speed of light (2.998e8). For instance, the physical Lore

Re: [Meep-discuss] Dispersive anisotropic media

2014-11-19 Thread Filip Dominec
ence, Pilani > > On Wed, Nov 19, 2014 at 9:04 PM, Filip Dominec > wrote: > >> Dear Rajath, >> it pretty much looks like missing- or incompatible version of libhdf. >> Have you tried the compilation procedure from this site? >> http://fzu.cz/~dominecf/meep/index.h

Re: [Meep-discuss] Units in Meep

2014-11-10 Thread Filip Dominec
I am convinced that it would be more convenient if MEEP allowed one to use the true SI units. (Thanks to the terahertz time-domain spectroscopy, I can directly compare data from simulation and from our laboratory.) The current approach of normalized units is however just a formal complication and i

Re: [Meep-discuss] monochromatic source

2014-08-25 Thread Filip Dominec
Hi, Farhad, generally both the continuous and the narrow-band-gaussian sources are good choices. Note however that for mathematical reasons a source short in time will always have a broad spectrum. Spectral selectivity has to be paid with computation time. Fortunately, the Fourier transform of a g

Re: [Meep-discuss] Part of PML

2014-08-17 Thread Filip Dominec
Hi, Sophia, as far as I know, it is possible to overlap any material with PML. The PML condition however applies for the whole boundary, so it can not be defined in its part only. Note I had some numerical stability issues when the PML was defined on one axis only while the other boundaries were p

Re: [Meep-discuss] Complex transmission and reflection coefficients

2014-07-14 Thread Filip Dominec
Hi, Joshua, I am not sure if it will be helpful to you, but two years ago I wrote a similar simulation of complex scattering coefficients (reflection+transmission) from a SRR or any other metamaterial unit cell. The codes are in python-meep and are publicly available at https://github.com/FilipDom

Re: [Meep-discuss] Multiple wavelengths and light direction

2014-06-19 Thread Filip Dominec
Hi, 1) spectrum: although the gaussian source with frequency around (1/1.55) would give an approximate result, you probably want to use the band_src_time which is described here: http://fzu.cz/~dominecf/misc/meep/index.html#bandsource I am going to contribute the code to github today. 2) If you wi

Re: [Meep-discuss] real metal material

2014-06-19 Thread Filip Dominec
uge positive imaginary part, it is a inductance-free conductor with tiny resistivity, again --> PEC. Complex infinity has no sign. Hope this helps, Filip Dominec IOP ASCR, Prague 2014-06-19 14:07 GMT+02:00, Weiner John : > Dear meep users and experts: > > I am really stuck, and I need your

Re: [Meep-discuss] (no subject)

2014-06-12 Thread Filip Dominec
NSHA#". The capitals look like shouting, and the hashmarks look like something as banging your fist on the table before and after saying it... Regards, Filip Dominec 2014-06-12 18:01 GMT+02:00, #SHAMPY MANSHA# : > Dear Steven and other Meep users, > > I am doing a simulation in whic

Re: [Meep-discuss] Beam splitters in meep?

2014-05-26 Thread Filip Dominec
Hi, as a simplest approach, I would suggest to reduce the 'sigma' parameter of the metal, first. It should allow the transmission to be adjusted even in the case of a low resolution. Filip 2014-05-26 7:06 GMT+02:00, jayv : > I need some beam splitters with varying T/R ratios, and I've had some > s

Re: [Meep-discuss] Permittivity of Silicon

2014-05-19 Thread Filip Dominec
Dear Thomas and all, I wrote most of my thoughts on the referred page; while sometimes the discrete Lorentzian oscillators suffice to precisely describe the dielectric function in the frequency range of interest, the optical spectra of semiconductors are usually of the harder cases. You can manual

Re: [Meep-discuss] Read refractive index data from file

2014-05-13 Thread Filip Dominec
If I correctly understand your question, you wish to define material properties as an arbitrary function of frequency. This is generally not possible in FDTD, see e.g. http://ab-initio.mit.edu/wiki/index.php/Dielectric_materials_in_Meep Filip 2014-05-13 19:07 GMT+02:00, Charalampos Mavidis : > Hel

Re: [Meep-discuss] PML overlapping dispersive material

2014-05-07 Thread Filip Dominec
Hi, Sophia, from my personal experience, the dispersive material can overlap with PML in the simulations. The question of stability under such conditions may be quite hard, though. You probably have to find out by trying which parameters are stable - and which not. Regards, Filip 2014-05-07 12:07

Re: [Meep-discuss] Drude material in MEEP

2014-04-29 Thread Filip Dominec
Dear Allen, surely the best guide for defining the Lorentzian model for materials in MEEP is http://ab-initio.mit.edu/wiki/index.php/Meep_Tutorial/Material_dispersion where also a working example is presented. I must apologize for the statement that you have to use a low-frequency Lorentzian oscil

Re: [Meep-discuss] Drude material in MEEP

2014-04-28 Thread Filip Dominec
Dear Allen, If I correctly understand your question, the Drude model can be approximated by the polarisation coming from a Lorentz oscillator with as extremely low resonance frequency. General approach is to select a realistic plasma frequency (for instance, omega_p := 2*pi*3.64e15 rad/s for alumin

Re: [Meep-discuss] Meaning of Complex fields in hd5 files

2014-04-28 Thread Filip Dominec
Hi, jayv, the real and imaginary fields should be the same except for that the imaginary part if shifted by 1/4 of its period. Thus computing the modulus of the complex fields you get a smooth, non-oscillating function, describing the magnitude of the field amplitude. I think you can disregard the

Re: [Meep-discuss] Zero Amplitude source

2014-04-24 Thread Filip Dominec
Hi, you can probably set the source cutoff to 0 (http://ab-initio.mit.edu/wiki/index.php/Meep_Reference#source). (But I must admit that I do not understand the purpose to do this.) F. 2014-04-23 22:13 GMT+02:00, James Green : > Dear MEEP users, > > Does anyone have any idea about how to make a ga

Re: [Meep-discuss] why the energy spectrum radiated by a point source doesn't match with the spectrum of source itself?

2014-04-22 Thread Filip Dominec
This appears exactly as the result of the said effect. The good point is that the gaussian source shape is reasonably preserved upon differentiation, so you may pre-compensate it by "fcen=0.40". Or, if you have tighter requirements on the source, you may use the newly implemented feature of band-so

Re: [Meep-discuss] why the energy spectrum radiated by a point source doesn't match with the spectrum of source itself?

2014-04-22 Thread Filip Dominec
May it be attributed to that the radiated field is a derivative of the current, so that high-frequency components are enhanced compared to the lower frequency ones? Or is it a more complicated deviation? Regards, Filip 2014-04-22 16:01 GMT+02:00, Ehsan Saei : > Dear Steven and MEEP users, > > I si

Re: [Meep-discuss] Doubt about concept of flux spectrum in MEEP

2014-04-15 Thread Filip Dominec
the Fourier-transformed fields in the same cross-section. > This is certainly do-able to implement (at least for dispersionless > materials). > > On Apr 14, 2014, at 5:21 PM, Filip Dominec wrote: > >> This is true - I mostly simulate the behaviour of an infinite periodic >> arr

Re: [Meep-discuss] Doubt about concept of flux spectrum in MEEP

2014-04-14 Thread Filip Dominec
Another option would be to define an amplitude-recording-plane: for instance, knowing the Ex-field and Hy-field, one can easily separate the amplitudes of both forward- and backward-propagating waves along the z-axis. The usual scattering problems with properly absorbing boundaries then require onl

Re: [Meep-discuss] meep-mpi doesn't print anything in terminal

2014-04-10 Thread Filip Dominec
Dear Ehsan, on Linux, you can use the 'tee' command: mpirun -np 8 meep-mpi test.ctl| tee test.out Did it help? Filip 2014-04-10 8:39 GMT+02:00, Ehsan Saei : > Dear all, > > I try to run a ctl file using meep-mpi on a ubuntu machine with the > following command line: > > mpirun -np

Re: [Meep-discuss] Meep frequency snapshots and near to far field transform

2014-03-28 Thread Filip Dominec
Dear Arthur, thank you for sharing your code! I cloned your source, downloaded the meep 1.2.1 sources and copied your files over the original directory. Then I used my old compilation procedure on Ubuntu [1] to compile meep and rebuilt the python-meep interface. The new functions of snapshot() an

Re: [Meep-discuss] meep mpi harminv problem

2014-03-13 Thread Filip Dominec
Hi, it appears to me that harminv may not use the MPI synchronisation at its end. If this is the case, it could result in a race condition, when one instance of harminv ends the multithreaded program, while the second one is still running. I suggest to add some field-manipulating function at the

Re: [Meep-discuss] silica dispersion model

2014-02-26 Thread Filip Dominec
Dear Pablo, I plot the Lorentz model for silica at http://fzu.cz/~dominecf/misc/eps/#SiO2 The page also contains the numerical data for the Lorentzian oscillators, along with additional information about the sources and application of the data. Hope this helps, Fillip 2014-02-26 13:14 GMT+01:00,

Re: [Meep-discuss] wave guide coupling

2014-02-13 Thread Filip Dominec
simple task. If there is coupling of freely propagating wave to any waveguide, the task becomes a bit more complicated and a sketch of your geometry would be needed. Best regards, Filip Dominec 2014-02-13 16:44 GMT+01:00, ISAAC MATERE : > Hi meep users > am trying to come up with a script for

Re: [Meep-discuss] Gaussian pulses in Meep

2014-02-07 Thread Filip Dominec
Dear Simon, I think your conjectures are all correct. The pseudo-gaussian spectral profile (I guess the spectrum is more precisely ~ f*exp[-(f-fcen)**2 / fwidth**2] or something like that) requires a gaussian time evolution of the source. The cutoff defines sharp edges in time when the source start

Re: [Meep-discuss] (no subject)

2014-01-06 Thread Filip Dominec
t how can I include the effects of > temperature on the emission of the EM wave from a surface with a finite > temperature? > > Thank you again. sorry if my message is long. > > Farhad > > > > On Monday, January 6, 2014 5:50 PM, Filip Dominec > wrote: > > D

Re: [Meep-discuss] (no subject)

2014-01-06 Thread Filip Dominec
implemented in MEEP successfully so even this can be done. The article touches a quite interesting topic, by the way. Best wishes, Filip Dominec 2014/1/6, FARHAD KAZEMI : > Dear all, > > I would like to model the radiation transfer between two plates > separated by nano-gap in which we h

Re: [Meep-discuss] Designing inclined geometry

2013-12-04 Thread Filip Dominec
Dear Prem, as far as I know, this task would require either 1) to build the desired structure by overlapping a big (potentially infinite) number of overlapping blocks, or 2) load the structure pixel-by-pixel either from a data file, or from a mathematical expression. I spent several days implementi

Re: [Meep-discuss] Near-field radiation between two parallel plates in MEEP

2013-11-27 Thread Filip Dominec
Dear Farhad, I must admit that maybe I do not fully understand your question, but anything is predicted by Maxwell equations should be possible in MEEP. In many cases my computations required MEEP to take into account the non-radiative electric or magnetic fields (such as when building a split-r

Re: [Meep-discuss] refractive index of silicon

2013-09-21 Thread Filip Dominec
consumes CPU time! Some remarks on the materials in MEEP and FDTD stability is also here: http://fzu.cz/~dominecf/misc/meep/#mat_mod Regards, Filip Dominec On Fri, 20 Sep 2013 15:55:50 -0700 Rui La wrote: > Dear meep users, > > I want to study the short wavelength light(300-700nm) pr

Re: [Meep-discuss] How to set frequency and fwidth

2013-09-20 Thread Filip Dominec
Hi, Rui, I also always need to compute it from the basics. 25 nm corresponds to 0.25 units, therefore one unit corresponds to 1e-7 m, 400 nm is roughly 700 THz in reality, therefore you would rescale the frequency to 700e12 * 1e-7 = 7e7, but MEEP defines speed of light as 1 instead of 3e8 m/s, so

Re: [Meep-discuss] Fw: question about meep

2013-07-09 Thread Filip Dominec
Dear Lin, for better illustration of what actually happens, could you generate the images again with the source shifted downwards, e. g. to 1/4 of the picture height? Thanks, Filip 2013/7/9, nw...@mail.sim.ac.cn : > > > > > > > > Dear meep users > > I'm a student of CAS(Chinese Academy of Scie

Re: [Meep-discuss] add-field, save-field, load-field, etc.

2013-06-21 Thread Filip Dominec
Hi, John, other people have asked a similar question and I am afraid it is not implemented yet. It is one of the features I proposed to implement in the mail few weeks ago: http://permalink.gmane.org/gmane.comp.science.electromagnetism.meep.general/4959 However, I did not get any reply so far, the

Re: [Meep-discuss] New feature: band-pass source

2013-06-19 Thread Filip Dominec
Just for comparison, these are some spectra recorded using the old Gaussian source: http://fzu.cz/~dominecf/misc/meep/img/ampli_debug_gaussian.png The structure couples to the low frequencies and induces ugly artifacts by spectral leakage. And here are those obtained using the band source with wid

[Meep-discuss] New feature: band-pass source

2013-06-19 Thread Filip Dominec
Dear MEEP users, let me share a little code snippet with you. Sometimes it happens that one wants to run a FDTD simulation with a broadband source, which however avoids some frequencies. A single gaussian source is not useful in such case, as its spectrum has always broad wings. Composing the sourc

[Meep-discuss] FDTD packages convergence

2013-06-12 Thread Filip Dominec
Dear all, this is an additional thought following my previous e-mail. Many interesting functions have been already implemented in different free FDTD codes. Just few examples of "special capabilities": * MEEP implements interesting physics and frequency-domain solver, * EMTL uses sub-pixel smoothin

[Meep-discuss] Amplitude monitor plane for MEEP 1.3

2013-06-12 Thread Filip Dominec
Dear MEEP developers and users, in our team, we employ MEEP to solve various problems and it is very useful. However, now and then we encounter some important functions to be missing. Most of them may be somehow written in Python code, but they would be the most efficient if hard-coded in the C++ s

Re: [Meep-discuss] organic material in meep

2013-06-04 Thread Filip Dominec
Dear Gholamhosain, the permittivity has to be established from experimental (or quantum-mechanical) sources, see for instance http://xplqa30.ieee.org/ielx5/5503869/6055237/6004802/html/img/6004802-fig-2-large.gif This graph may be easily converted to numerical data using Engauge or similar program

Re: [Meep-discuss] Get (n, k) for a frequency from a dispersive material.

2013-05-24 Thread Filip Dominec
And I forgot to add that n = real(sqrt(eps)) k = imag(sqrt(eps)) F. D. ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss

Re: [Meep-discuss] Get (n, k) for a frequency from a dispersive material.

2013-05-24 Thread Filip Dominec
Hi, Oscar, this is rather straightforward. The dielectric function used by MEEP is given by the well-known Lorentz formula. (I have verified it numerically, but this is a bit longer to explain.) In Python, I use the following function for the material. As an example I will give the SiO2 definition

Re: [Meep-discuss] Porting MEEP to the Web

2013-05-24 Thread Filip Dominec
with Physical simulations using NaCl, > but I'm no expert in this technology. > We can discuss issues with Google developers which are very supportive, if > you got interested in the project. > I can't port MEEP to this technology by myself, but with some help, I think &g

  1   2   >