Analysis routines¶
In addition to the main pyro program, there are many analysis tools that we describe here. Note: some problems write a report at the end of the simulation specifying the analysis routines that can be used with their data.
compare.py
: this takes two simulation output files as input and compares zone-by-zone for exact agreement. This is used as part of the regression testing.usage:
./compare.py file1 file2
plot.py
: this takes an output file as input and plots the data using the solver’s dovis method. It deduces the solver from the attributes stored in the HDF5 file.usage:
./plot.py file
analysis/
convergence.py
: this compares two files with different resolutions (one a factor of 2 finer than the other). It coarsens the finer data and then computes the norm of the difference. This is used to test the convergence of solvers.dam_compare.py
: this takes an output file from the shallow water dam break problem and plots a slice through the domain together with the analytic solution (calculated in the script).usage:
./dam_compare.py file
gauss_diffusion_compare.py
: this is for the diffusion solver’s Gaussian diffusion problem. It takes a sequence of output files as arguments, computes the angle-average, and the plots the resulting points over the analytic solution for comparison with the exact result.usage:
./gauss_diffusion_compare.py file*
incomp_converge_error.py
: this is for the incompressible solver’s converge problem. This takes a single output file as input and compares the velocity field to the analytic solution, reporting the L2 norm of the error.usage:
./incomp_converge_error.py file
plotvar.py
: this takes a single output file and a variable name and plots the data for that variable.usage:
./plotvar.py file variable
sedov_compare.py
: this takes an output file from the compressible Sedov problem, computes the angle-average profile of the solution and plots it together with the analytic data (read in fromcylindrical-sedov.out
).usage:
./sedov_compare.py file
smooth_error.py
: this takes an output file from the advection solver’s smooth problem and compares to the analytic solution, outputing the L2 norm of the error.usage:
./smooth_error.py file
sod_compare.py
: this takes an output file from the compressible Sod problem and plots a slice through the domain over the analytic solution (read in fromsod-exact.out
).usage:
./sod_compare.py file