Hi Dr. Renaud Vallée,
I am having the same problem as you encountered. May I know how did you
solve it eventually? My code is very similar to yours.
I added this though.
s.add_polarizability(sigma,type(Ez),.5,.1);
fields f(&s,0,true);
The 3rd argument of fields constructor is store_pol_energy, which is by
default 0.
fields(structure *, double m=0, bool store_pol_energy=0, double beta=0,
bool zero_fields_near_cylorigin=true);
I thought this would work but the result is still the same as the code
without active region.
Please kindly advise if you know how to add polarizability now.
Best Regards,
Bin Huang (Bryan)
MIT Graduate Student '10
Computation for Design and Optimization
(+65)98947649
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
