Re: [Meep-discuss] Fast access to the field array for output to VTK, slice images, animations...

2008-08-07 Thread Jefferson Thomas
that the node could then operate on it easily and save it to a desired format by itself, without even needing to use the parallel I/O. What do you think of that ? Thomas Jefferson 2008/8/7 Steven G. Johnson [EMAIL PROTECTED]: On Aug 6, 2008, at 9:52 PM, Jefferson Thomas wrote: The problem

Re: [Meep-discuss] Using perfect metal in C++ interface

2008-08-06 Thread Jefferson Thomas
for the first dispersive material, 2 for the second one, 3 for the third one and so on ? After all the dispersive material can not overlap, can they ? Thomas Jefferson 2008/8/6 Steven G. Johnson [EMAIL PROTECTED]: On Aug 5, 2008, at 4:30 AM, Jefferson Thomas wrote: Thank you for your answer matt. I

[Meep-discuss] Fast access to the field array for output to VTK, slice images, animations...

2008-08-06 Thread Jefferson Thomas
Hi, In the program I am developping now I have to be able to output a given type of field (Ex, Hz, ... ) into VTK files. It is a format useful for 3d-data storage and very efficient for their 3d-visualisation. As for now I am doing it with first outputtting the field through the

Re: [Meep-discuss] Using perfect metal in C++ interface

2008-08-05 Thread Jefferson Thomas
), double omega, double gamma, double delta_epsilon, double energy_saturation); There are some examples available in tests/known_results.cpp and tests/symmetry.cpp. On Tue, 5 Aug 2008, Jefferson Thomas wrote: Thank you all. It works. Could I also ask how to implement in C++ interface

[Meep-discuss] Using perfect metal in C++ interface

2008-08-04 Thread Jefferson Thomas
Hi, I would like to know how to use the perfect metal in the C++ interface. I can not find it in the documentation. Do I have to return some specific value from the eps function or it is done in a different way ? Thomas Jefferson ___ meep-discuss

Re: [Meep-discuss] Using perfect metal in C++ interface

2008-08-04 Thread Jefferson Thomas
Thank you all. It works. Could I also ask how to implement in C++ interface the dispersive materials like those described in http://ab-initio.mit.edu/wiki/index.php/Meep_Tutorial/Material_dispersion ? Thomas Jefferson 2008/8/5 Steven G. Johnson [EMAIL PROTECTED]: On Aug 4, 2008, at 4:42 PM,

Re: [Meep-discuss] Field divergence in meep after many steps for a photonic crystal

2008-07-31 Thread Jefferson Thomas
, 2008, at 12:15 AM, Jefferson Thomas wrote: I have also some other kind of divergences in other cases. They seem not to be connected to symmetries but arise from PML layers. I am starting to think if I am not doing any similar mistake there with those last pixels beyond the calculation volume

Re: [Meep-discuss] Field divergence in meep after many steps for a photonic crystal

2008-07-30 Thread Jefferson Thomas
the fact that neither your structure nor your sources obey the symmetry you specified. On Jul 22, 2008, at 9:40 PM, Jefferson Thomas wrote: meep::symmetry S=meep::identity(); S = S+ meep::mirror(meep::X, v)*(-1); /// Introducing symmetries on X and Y axis S = S+ meep::mirror(meep::Y, v

Re: [Meep-discuss] Field divergence in meep after many steps for a photonic crystal

2008-07-30 Thread Jefferson Thomas
for your help Steven. Thomas Jefferson 2008/7/31 Steven G. Johnson [EMAIL PROTECTED]: On Jul 30, 2008, at 9:20 PM, Jefferson Thomas wrote: What you are saying is really very helpful. In fact I had no idea I have to supply the symmetric source. I will try it right now and see how it goes. See

Re: [Meep-discuss] The fastest way to calculate the total energy

2008-07-28 Thread Jefferson Thomas
2008/7/29 Steven G. Johnson [EMAIL PROTECTED]: On Jul 28, 2008, at 12:19 AM, Jefferson Thomas wrote: double ex = fields-field_energy_in_box(meep::Ex, volume.surroundings()); double ey = fields-field_energy_in_box(meep::Ey, volume.surroundings()); double ez = fields-field_energy_in_box(meep

[Meep-discuss] Influencing fields in meep in a specific point

2008-07-28 Thread Jefferson Thomas
Hi, Is there a way to influence the field in meep in a specific location? I would like to set to zero all the field values in some specific 3d box inside the calculation volume. There is function fields::zero_fields() but it sets to zero all the calculation volume. Is there a way to do it for a

Re: [Meep-discuss] The fastest way to calculate the total energy

2008-07-27 Thread Jefferson Thomas
, volume.surroundings()); double e=ex+ey+ez+hx+hy+hz; Is it optimal ? It is still not as fast as the algorithm in Fortran I mentioned earlier. Thank you, Thomas Jefferson 2008/7/28 Zheng Li [EMAIL PROTECTED]: Jefferson Thomas wrote: Hi, I want to calculate the total energy in my simulation

Re: [Meep-discuss] Unable to output fields to HDF in MPI meep.

2008-07-14 Thread Jefferson Thomas
, v.surroundings()); - does not work It seems to be just a problem of FIELD output in HDF while using MPI. Thomas Jefferson 2008/7/13 Steven G. Johnson [EMAIL PROTECTED]: On Jul 6, 2008, at 9:52 PM, Jefferson Thomas wrote: BUT, while running it in MPI mode, with 2 processors for example it stucks

Re: [Meep-discuss] Synchronizing in mpi meep

2008-07-14 Thread Jefferson Thomas
Thank you Thomas Jefferson 2008/7/13 Steven G. Johnson [EMAIL PROTECTED]: On Jul 3, 2008, at 2:27 AM, Jefferson Thomas wrote: Hi, I am using meep in MPI. I can not find a logical way to synchronize my myself between all the processes. In mpi should be a command like mpi_barrier which

Re: [Meep-discuss] Speed comparison of meep running in Scheme and C++. C++ is several times slower ?

2008-07-14 Thread Jefferson Thomas
on this forum (Unable to output fields to HDF in MPI meep.) seems to be solved. When using real fields I can output HDF fields in MPI mode without problems. Thank you. Thomas Jefferson 2008/7/13 Steven G. Johnson [EMAIL PROTECTED]: On Jul 7, 2008, at 4:37 AM, Jefferson Thomas wrote: I am

[Meep-discuss] How to use metallic boundary condition in C++

2008-07-09 Thread Jefferson Thomas
Hi, How can I use the metallic boundary condition in the C++ interface ? With structure s(v, eps, pml(1.0)); I create a PML layer but in Scheme interface it is also possible to create a metallic boundary. How to do the same in C++ ? Thank you for any help, Thomas Jefferson

[Meep-discuss] Speed comparison of meep running in Scheme and C++. C++ is several times slower ?

2008-07-07 Thread Jefferson Thomas
Hi, I am developing an application using C++ meep. I have just done some speed tests in the same environment using meep in Scheme mode and as a C++ program. Here are the files i used: CTL file test.ctl: (set! geometry-lattice (make lattice (size 5 10 2))) (set! sources (list (make

[Meep-discuss] Unable to output fields to HDF in MPI meep.

2008-07-06 Thread Jefferson Thomas
Hi, I am having problems outputting fields like Ez, Hx etc. to HDF while running C++ meep in MPI. Here is a simple test program: // Test program starts #include meep.hpp #include iostream using namespace meep; double eps(const vec p) { return p.x()+p.y(); } int

Re: [Meep-discuss] Converting meep time units to step numbers and vice versa

2008-07-03 Thread Jefferson Thomas
Thank you, it works Thomas Jefferson 2008/6/28 Zheng Li [EMAIL PROTECTED]: Jefferson Thomas wrote: Hi, I want to convert time units which Im getting from meep::fields::time() to the number of steps and vice versa. To make it more clear: I have to convert a stepping loop of type: for (int

[Meep-discuss] Meep C++ documentation

2008-07-03 Thread Jefferson Thomas
Hi, Is there any detailed documentation on meep c++ interface ? http://ab-initio.mit.edu/wiki/index.php/Meep_C-plus-plus_Reference is definitely not enough. It would be nice if there was something like doxygen documentation for every class and all the functions but I can not find it and the