In the case where reflections are a significant part of the dynamics (thin samples, nonlinear reflection from a surface) we will sometimes want to solve Maxwell's equations rather than the unidirectional propagation equation. This solution could be integrated within a sequence as a function for calling the Maxwell's equation solver and returning a new space vs. time grid from a given plane within the Maxwell grid.
First version will assume planar sample in vacuum for simplicity. Basic structure:
input field (taken from current value in solution) injected from lefthand side of maxwell grid (3D or 2D space only). Material pixels have running calculation of polarization/current which will only work for Lorentzian/Gaussian materials.
Inputs of function:
material
z extent of grid
start z of material on grid
end z of material on grid
z location of output plane
integer divisor of time step
integer divisor of z step
In the case where reflections are a significant part of the dynamics (thin samples, nonlinear reflection from a surface) we will sometimes want to solve Maxwell's equations rather than the unidirectional propagation equation. This solution could be integrated within a sequence as a function for calling the Maxwell's equation solver and returning a new space vs. time grid from a given plane within the Maxwell grid.
First version will assume planar sample in vacuum for simplicity. Basic structure: input field (taken from current value in solution) injected from lefthand side of maxwell grid (3D or 2D space only). Material pixels have running calculation of polarization/current which will only work for Lorentzian/Gaussian materials.
Inputs of function: material z extent of grid start z of material on grid end z of material on grid z location of output plane integer divisor of time step integer divisor of z step