You can use this method, defined in polarization.cpp:
structure::add_polarizability(double sigma(const vec &), 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 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, matt wrote:
>>> Return a value from your eps function which is large and negative.
>>> -1e20 should work.
>>
>> Or better yet, -HUGE_VAL (which is defined in math.h, and should give -
>> Inf).
>>
>> Steven
>>
>> _______________________________________________
>> meep-discuss mailing list
>> [email protected]
>> http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
>>
>
> _______________________________________________
> meep-discuss mailing list
> [email protected]
> http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
>
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
