Matlab framework to compute projected preconditioned gradient descent method and projected multigrid descent method. Experiments optimizing the residues of the Eikonal equation for fire arrival time in a coupled atmosphere-fire data (from WRF-SFIRE simulations). File startup.m includes into the path all the source routines.
Routines to compute sine FFT for the minus Laplacian operator with Dirichlet boundary conditions.
Routines to compute fire dynamics as rate of spread, fuels, and interpolation of wind components and fuel moisture.
Routines to interpolate and condense operator H in the linear system of constraints Hu=g.
Routines to read data from Netcdf files (WRF-SFIRE simulation output files).
Routines in the optimization process: computation of the objective function, binary search, linesearch, projected preconditioned gradient descent method and projected multigrid descent method.
Routines to plot results.
Routines to setup cases: ideal, ideal from WRF-SFIRE simulation and real cases.
Other routines as: abstract KML files information to matlab, record the graphics in a gif file, etc.
To test the functions:
-
Start Matlab inside Multigrid folder, so inside Multigrid folder type in terminal
$ matlab -nodesktop -nodisplay
-
If you started Matlab in another folder, run inside Matlab
>> startup
-
Run any of the tests in the tests folder. For example:
>> cd tests
>> pdirichlet_constr_test