Dear Bin, Steven, meep users I did not solve the problem in C++. I do not understand what is going on there. I guess there is a problem of implementation somehow, while it works perfectly with the scheme interface!!! May be the developpers could have a look on this recurrent problem, at some of their available time. Best,
REnaud Le mercredi 27 octobre 2010 à 23:32 -0400, Bin Huang a écrit : > 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 -- Dr. Renaud Vallée Centre de Recherche Paul Pascal (CRPP, CNRS, UPR8641) 115, avenue du docteur Albert Schweitzer 33600 Pessac France _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss

