examples.dg.maxwell.light_reflect Namespace Reference

Variables
str	torder = 'SSPRK3'
str	sorder = 'thirdOrder'
int	skin_depth_norm = 1
int	omega_p_norm = 1
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	bcs = []
list	apps
warpy	writer
warpy	variable_adjusters
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_reflect.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', 'Ey', 'Ez', 'Bx', 'By', 'Bz'],
                     spatial_scheme='Nodal')
]

apps

list examples.dg.maxwell.light_reflect.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_reflect.B0 = 1

B_f

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

bcs

list examples.dg.maxwell.light_reflect.bcs = []

By_ic

warpy examples.dg.maxwell.light_reflect.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_reflect.By_node = B_f.read_data()[1].flatten()

Bz_ic

warpy examples.dg.maxwell.light_reflect.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_reflect.cmap

configs

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

connectivity

examples.dg.maxwell.light_reflect.connectivity

coords

examples.dg.maxwell.light_reflect.coords

detect_nonscalar

examples.dg.maxwell.light_reflect.detect_nonscalar

dset

warpy examples.dg.maxwell.light_reflect.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_reflect.dt = 1e-4

dt_controller

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

E0

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

else

examples.dg.maxwell.light_reflect.else :

energy_factor

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

Ey_ic

warpy examples.dg.maxwell.light_reflect.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_reflect.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_reflect.field

Initial value:

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

gen_xdmf

examples.dg.maxwell.light_reflect.gen_xdmf

mesh

warpy examples.dg.maxwell.light_reflect.mesh

Initial value:

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

null_func

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

omega_p_norm

int examples.dg.maxwell.light_reflect.omega_p_norm = 1

pmesh

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

sd

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

shading

examples.dg.maxwell.light_reflect.shading

sim

warpy examples.dg.maxwell.light_reflect.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,6/(skin_depth_norm*omega_p_norm)],
                   write_steps=100,
                   verbosity='info',
                   global_verbosity='info')

skin_depth_norm

int examples.dg.maxwell.light_reflect.skin_depth_norm = 1

sorder

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

spatial_solver

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

temporal_solver

warpy examples.dg.maxwell.light_reflect.temporal_solver

Initial value:

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

torder

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

True

examples.dg.maxwell.light_reflect.True

try

examples.dg.maxwell.light_reflect.try :

variable_adjusters

warpy examples.dg.maxwell.light_reflect.variable_adjusters

Initial value:

=  warpy.variable_adjusters.boundary_conditions(
        name = 'bc_va', 
        priority = 0,
        spatial_order = sorder,
        applications = bcs,
        is_virtual_bc = False,
        copy_all=True,
    )

writer

warpy examples.dg.maxwell.light_reflect.writer

Initial value:

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