analysis
All user defined applications using HAD must include a subroutine called "analysis"
and another called "analysis_pt".
One of these two routines is called after each evolution step. Quantities used to monitor the
solution may be calculed here, such as constraints, energies, etc. The "analysis_pt" routine
is for pointwise analysis. Default behavior is to call "analysis".
Notes
- Parameter analysis_[field_name] enables the user to output global min, max, l2norm, and spatial
integrals of any gridfuntion.
- If the parameter point_wise_analysis
is set, "analysis_pt" will be called instead of "analysis".
- If the parameter
initial_analysis is set, "analysis" or "analysis_pt" will be called before the first evolution step.
- Header: subroutine analysis(u2, u0, dxu, dyu, dzu, v, dxv, dyv, dzv, w, par)
type(gridfunction), dimension(NV) :: v, dxv, dyv, dzv
type(gridfunction), dimension(NU) :: u2, u0, dxu, dyu, dzu
type(gridfunction), dimension(NW) :: w
CCTK_REAL, dimension(:) :: par
- Header: subroutine analysis_pt(u_pt, u0_pt, dxu_pt, dyu_pt, dzu_pt,v_pt , dxv_pt, dyv_pt, dzv_pt, w_pt, par)
CCTK_REAL, dimension(NU), intent(inout) :: u_pt,u0_pt,dxu_pt,dyu_pt,dzu_pt
CCTK_REAL, dimension(NV), intent(inout) :: v_pt,dxv_pt,dyv_pt,dzv_pt
CCTK_REAL, dimension(NW), intent(in) :: w_pt
CCTK_REAL, dimension(NPAR),intent(inout) :: par
- u2 points to the updated ufuncs.
- u0 points to the previous timestep ufuncs.
- dxu points to the x-spatial derivatives of u2.
- dyu points to the y-spatial derivatives of u2.
- dzu points to the z-spatial derivatives of u2.
- v points to the current vfuncs.
- dxv points to the x-spatial derivatives of v.
- dyv points to the y-spatial derivatives of v.
- dzv points to the z-spatial derivatives of v.
- w points to the wfuncs
- "_pt" indicates a pointwise value rather than a gridfunction.
See Also
setup, point_wise_analysis, initial_analysis,
analysis_[field_name],
List of required routines
Index of all manual pages
Examples
had/src/hyperCurvedWave/analysis.f90