examples.fvm.light_prop Namespace Reference

Variables
int	torder = 3
int	sorder = 2
int	omega_c_norm = 3
int	omega_p_norm = 5
int	c_0_norm = 4
warpy	mesh
warpy	field
int	E0 = omega_c_norm*c_0_norm/omega_p_norm**2
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	ic_applicators
list	ss1
list	ss_sets = [ss1]
warpy	writer
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_cell = B_f.read_data()[1].flatten()
	shading
	cmap

Variable Documentation

B0

int examples.fvm.light_prop.B0 = 1

B_f

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

By_cell

examples.fvm.light_prop.By_cell = B_f.read_data()[1].flatten()

By_ic

warpy examples.fvm.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])

Bz_ic

warpy examples.fvm.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])

c_0_norm

int examples.fvm.light_prop.c_0_norm = 4

cmap

examples.fvm.light_prop.cmap

configs

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

connectivity

examples.fvm.light_prop.connectivity

coords

examples.fvm.light_prop.coords

detect_nonscalar

examples.fvm.light_prop.detect_nonscalar

dset

warpy examples.fvm.light_prop.dset

Initial value:

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

E0

int examples.fvm.light_prop.E0 = omega_c_norm*c_0_norm/omega_p_norm**2

else

examples.fvm.light_prop.else :

Ey_ic

warpy examples.fvm.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.fvm.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.fvm.light_prop.field

Initial value:

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

gen_xdmf

examples.fvm.light_prop.gen_xdmf

ic_applicators

list examples.fvm.light_prop.ic_applicators

Initial value:

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

mesh

warpy examples.fvm.light_prop.mesh

Initial value:

=  warpy.mesh.block(Bounds=[-1, 1]*2,
                        NumCells=[50]*2,
                        NodeSets=['Left', 'Right','Bottom','Top'],
                        NumLayers=sorder,
                        PeriodicBoundaries=['Left','Right','Bottom','Top'])

null_func

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

omega_c_norm

int examples.fvm.light_prop.omega_c_norm = 3

omega_p_norm

int examples.fvm.light_prop.omega_p_norm = 5

pmesh

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

sd

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

shading

examples.fvm.light_prop.shading

sim

warpy examples.fvm.light_prop.sim

Initial value:

=  warpy.fvm_sim(name='light_prop',
                    torder=torder,
                    sorder=sorder,
                    meshes=[mesh],
                    initial_conditions=ic_applicators,
                    boundary_conditions=None,
                    writers=[writer],
                    subsolver_sets=ss_sets,
                    time=[0,1],
                    write_steps=100,
                    verbosity='info')

sorder

int examples.fvm.light_prop.sorder = 2

ss1

list examples.fvm.light_prop.ss1

Initial value:

=  [warpy.hyperapps.pcmaxwell(name='pcmaxwell_flux',
                                 field_var=field,
                                 Chi_E = 0,
                                 Chi_B = 0,
                                 omega_c_norm = omega_c_norm,
                                 omega_p_norm = omega_p_norm,
                                 c_0_norm = c_0_norm,
                                 fluxes_only=True)]

ss_sets

list examples.fvm.light_prop.ss_sets = [ss1]

torder

int examples.fvm.light_prop.torder = 3

True

examples.fvm.light_prop.True

try

examples.fvm.light_prop.try :

writer

warpy examples.fvm.light_prop.writer

Initial value:

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