added pml and new subroutines to calculate boundary and initial conditions for aux and q
created pml and no pml case in fdtd_pml
corrected PML components
changes and cases introduced for general solver. COrrections to python main routine to accomodate old example and typos
implemented none in boundary line conditions
typos corrected
Fix boundary conditions
Many changes to PML and main routines
a sanity check
etar, a.k.a set_aux has not vacuum mode, uses the da to build the local refractive index map
Add reader, write and probe
Add VTK writter