I am trying to code an structure with a 'jump' in epsilon value during the execution of the code. I know the built-in function phase_in_materia But, I don't want a linear interpolation. I need a 'jump' in epsilon. For example, for time < t0 -> epsilon = 2.7 and for time > t0 -> epsilon = 3.5.
Below, an example of a straightforward code to verify the structure of the desired code. Please, forgive me for such a simple question. -- import meep as mp cell_size = mp.Vector3(6,6,0) geometry1 = [mp.Cylinder(center=mp.Vector3(),radius=1.0,material=mp.Medium(index=3.5))] geometry2 = [mp.Cylinder(center=mp.Vector3(1,1),radius=1.0,material=mp.Medium(index=3.5))] sim = mp.Simulation(cell_size=cell_size, geometry=geometry1, resolution=20) sim.init_sim() sim.run(mp.at_beginning(mp.output_epsilon),mp.at_time(4.0,sim.set_materials(geometry=geometry2)),mp.at_time(8.0,mp.output_epsilon),until=10) -- Cícero Julião da Silva Junior ____________________________________ Instituto Federal de Alagoas - IFAL Campus Murici http://www.ifal.edu.br/
_______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss