Ok so this is Sasha's README file - taking notes about how to use peter's filamentation code.

Waveion is the fdtd code for ionization off a surface - includes Peter's new and fancy ionization equations (i.e. non-Keldysh to allow for sub-cycle pulse changes).

The main filamentation code is contained in gsolver and the matlab code to load the results is in m. Bin contains the output of runs and run configurations (pulse power, output results, etc.) It is recommended to run the code from inside the bin directory using: mpirun -np <# of CPUS> ./gsolver_solid_qdht

Physical processes taken into account:

1. single mode Raman (with given frequency and decay time).
2. instantaneous Kerr; self-steepening
3. Plasma; Keldysh ionization, decay with constant time, avalanche ionization (goverend by collision time), plasma absorption and defocusing.

All of the code needs gsl, gslcblas, cuda compiler, openmpi, and fftw3 with openmpi flags to run.

Different gsolvers: Use ./gsolver_solid_qdht for axially symmetric beams (i.e. using r, z coordinates instead of x, y, z). QDHT stands for Quasi Discrete Henkel Transform - this is fast only for small problems, N_r < 512, 1024. Need a total grid of <16384 for reversibility.

Matlab functions: Use load_UDF function to load the results of the simulation. See seperate comments there for reading/loading info.