examples.dg.maxwell.light_prop Namespace Reference

Variables
str	torder = 'SSPRK3'
str	sorder = 'thirdOrder'
int	skin_depth_norm = 1/3
int	omega_p_norm = 5
warpy	mesh
warpy	field
tuple	energy_factor = (omega_p_norm*skin_depth_norm)
tuple	E0 = energy_factor
int	B0 = 1
	else :
warpy	Ey_ic
warpy	Bz_ic
warpy	Ez_ic
warpy	By_ic
warpy	null_func = warpy.functions.set_to(name='null', value=0)
list	applicators
list	apps
warpy	writer
warpy	spatial_solver = warpy.spatial_solvers.dg(name="dg", spatial_order=sorder, applications=apps)
warpy	temporal_solver
int	dt = 1e-4
warpy	dt_controller = warpy.dt_calc.fixed_dt(init_dt=dt)
warpy	sim
	gen_xdmf
	True
	detect_nonscalar
	try :
warpy	configs = warpy.load_configs('light_prop_config.py')
warpy	dset
list	sd = [sd for sd in dset.subdomains if sd.name=='domain'][0]
list	B_f = [f for f in dset.fields[0] if f.name=='B'][0]
	coords
	connectivity
Triangulation	pmesh = Triangulation(coords[:,0],coords[:,1],connectivity)
list	By_node = B_f.read_data()[1].flatten()
	shading
	cmap

Variable Documentation

applicators

list examples.dg.maxwell.light_prop.applicators

Initial value:

=  [
    warpy.applicator(spatial_order=sorder,
                     fun=Ey_ic,
                     var=field,
                     components=['Ey'],
                     spatial_scheme='Nodal'),
    warpy.applicator(spatial_order=sorder,
                     fun=Bz_ic,
                     var=field,
                     components=['Bz'],
                     spatial_scheme='Nodal'),
    warpy.applicator(spatial_order=sorder,
                     fun=Ez_ic,
                     var=field,
                     components=['Ez'],
                     spatial_scheme='Nodal'),
    warpy.applicator(spatial_order=sorder,
                     fun=By_ic,
                     var=field,
                     components=['By'],
                     spatial_scheme='Nodal'),
    warpy.applicator(spatial_order=sorder,
                     fun=null_func,
                     var=field,
                     components=['Ex','Bx'],
                     spatial_scheme='Nodal')
]

apps

list examples.dg.maxwell.light_prop.apps

Initial value:

=  [warpy.apps.maxwell.flux(name='maxwell', field=field,
                                skin_depth_norm=skin_depth_norm,
                                omega_p_norm=omega_p_norm)]

B0

int examples.dg.maxwell.light_prop.B0 = 1

B_f

list examples.dg.maxwell.light_prop.B_f = [f for f in dset.fields[0] if f.name=='B'][0]

By_ic

warpy examples.dg.maxwell.light_prop.By_ic

Initial value:

=  warpy.functions.fourier(name='By_ic',base=0,
                                direction=[1,0,0],
                                center=[0,0,0],
                                sine_amplitudes=[B0],
                                sine_frequencies=[pi],
                                sine_phases=[0])

By_node

examples.dg.maxwell.light_prop.By_node = B_f.read_data()[1].flatten()

Bz_ic

warpy examples.dg.maxwell.light_prop.Bz_ic

Initial value:

=  warpy.functions.fourier(name='Bz_ic',base=0,
                                direction=[1,0,0],
                                center=[0,0,0],
                                sine_amplitudes=[B0],
                                sine_frequencies=[pi],
                                sine_phases=[0])

cmap

examples.dg.maxwell.light_prop.cmap

configs

warpy examples.dg.maxwell.light_prop.configs = warpy.load_configs('light_prop_config.py')

connectivity

examples.dg.maxwell.light_prop.connectivity

coords

examples.dg.maxwell.light_prop.coords

detect_nonscalar

examples.dg.maxwell.light_prop.detect_nonscalar

dset

warpy examples.dg.maxwell.light_prop.dset

Initial value:

= warpy.dataset(mesh_file_prefixes=['meshes/plotMesh'],
                           var_prefixes=['data/light_prop'],
                           dset_path=configs['test_rundir']('maxwell',sim.name))

dt

int examples.dg.maxwell.light_prop.dt = 1e-4

dt_controller

warpy examples.dg.maxwell.light_prop.dt_controller = warpy.dt_calc.fixed_dt(init_dt=dt)

E0

int examples.dg.maxwell.light_prop.E0 = energy_factor

else

examples.dg.maxwell.light_prop.else :

energy_factor

tuple examples.dg.maxwell.light_prop.energy_factor = (omega_p_norm*skin_depth_norm)

Ey_ic

warpy examples.dg.maxwell.light_prop.Ey_ic

Initial value:

=  warpy.functions.fourier(name='Ey_ic',
                                base=0,
                                direction=[1,0,0],
                                center=[0,0,0],
                                sine_amplitudes=[E0],
                                sine_frequencies=[pi],
                                sine_phases=[0])

Ez_ic

warpy examples.dg.maxwell.light_prop.Ez_ic

Initial value:

=  warpy.functions.fourier(name='Ez_ic',
                                base=0,
                                direction=[1,0,0],
                                center=[0,0,0],
                                sine_amplitudes=[E0],
                                sine_frequencies=[pi],
                                sine_phases=[0])

field

warpy examples.dg.maxwell.light_prop.field

Initial value:

=  warpy.variable(name='field',
                       components=['Ex', 'Ey', 'Ez', 'Bx', 'By', 'Bz'],
                       basis_array_set=sorder)

gen_xdmf

examples.dg.maxwell.light_prop.gen_xdmf

mesh

warpy examples.dg.maxwell.light_prop.mesh

Initial value:

=  warpy.mesh.block(Bounds=[-1, 1],
                        NumCells=[200],
                        NodeSets=['Left', 'Right'],
                        NumLayers=1,
                        PeriodicBoundaries=['Left','Right'],
                        basis_array_set=sorder)

null_func

warpy examples.dg.maxwell.light_prop.null_func = warpy.functions.set_to(name='null', value=0)

omega_p_norm

int examples.dg.maxwell.light_prop.omega_p_norm = 5

pmesh

examples.dg.maxwell.light_prop.pmesh = Triangulation(coords[:,0],coords[:,1],connectivity)

sd

list examples.dg.maxwell.light_prop.sd = [sd for sd in dset.subdomains if sd.name=='domain'][0]

shading

examples.dg.maxwell.light_prop.shading

sim

warpy examples.dg.maxwell.light_prop.sim

Initial value:

=  warpy.dg_sim(name='light_prop',meshes=[mesh],
                   initial_conditions=applicators,
                   temporal_solvers=[temporal_solver],
                   writers=[writer],
                   dt_controller = dt_controller,
                   flexible_writeout=True,
                   time=[0,2/(skin_depth_norm*omega_p_norm)],
                   write_steps=100,
                   verbosity='info',
                   global_verbosity='info')

skin_depth_norm

int examples.dg.maxwell.light_prop.skin_depth_norm = 1/3

sorder

str examples.dg.maxwell.light_prop.sorder = 'thirdOrder'

spatial_solver

warpy examples.dg.maxwell.light_prop.spatial_solver = warpy.spatial_solvers.dg(name="dg", spatial_order=sorder, applications=apps)

temporal_solver

warpy examples.dg.maxwell.light_prop.temporal_solver

Initial value:

=  warpy.host_actions.erk(name="rk", scheme=torder,
                                         spatial_solvers=[spatial_solver])

torder

str examples.dg.maxwell.light_prop.torder = 'SSPRK3'

True

examples.dg.maxwell.light_prop.True

try

examples.dg.maxwell.light_prop.try :

writer

warpy examples.dg.maxwell.light_prop.writer

Initial value:

=  warpy.host_actions.writer(name='writer',
                                   ReadVars=[field])