- GetHess(CenterParams, parsteps, Hess, CenterLike)
Get the curvature matrix (Hess) for the grid defined by CenterParams
and parsteps.
Return the answer in Hess.
Input only: CenterParams, parsteps
Ouput only: Hess
-
ParamSet
:: CenterParams
-
real(mcp) :: parsteps(num_params_used)
-
real(mcp) :: Hess(num_params_used,num_params_used)
-
real(mcp) :: CenterLike
- GetStepsForDChisq1(CenterParams, maxstepsizetries, parsteps, bestdll, LCenter)
Stays centered on CenterParams
But for each paramter, iterates for a maximum of maxstepsizetries to try to find
parsteps such that Delta chisq is of order 1 (when stepping only in that parameter)
Inputs: CenterParams, LCenter, maxstepsizetries
Outputs: parsteps
-
ParamSet
:: CenterParams
-
integer intent(in) :: maxstepsizetries
-
real(mcp) :: parsteps(num_params_used)
-
real(mcp) :: bestdll
-
real(mcp) :: LCenter
-
real(mcp)
EstCovmat(sucess, bestdll_in, Params)
-
integer intent(out) :: sucess
-
real(mcp) intent(in) :: bestdll_in
the dloglike value used to determine grid spacing eg. try 4
-
ParamSet
:: Params
-
real(mcp)
PlaceGrid(CenterParams, parsteps, bestdll)
Takes the input CenterParams and parsteps values and improves on them
While doing so, fills the likelhood grid
CenterParams and parsteps values
-
ParamSet
:: CenterParams
-
real(mcp) :: parsteps(num_params_used)
-
real(mcp) :: bestdll