LM_nonlinear_fit Subroutine

private subroutine LM_nonlinear_fit(is, g, n_params, nJT, params, param_mask, iperp, npar, ipparbar_lower, quality)


  • proc~~lm_nonlinear_fit~~UsesGraph proc~lm_nonlinear_fit LM_nonlinear_fit module~alps_io alps_io proc~lm_nonlinear_fit->module~alps_io module~alps_var alps_var proc~lm_nonlinear_fit->module~alps_var mpi mpi proc~lm_nonlinear_fit->mpi module~alps_var->mpi

This subroutine processes the nonlinear Levenberg-Marquart algorithm and returns the one-dimensional array params at a given iperp. The variable quality is the sum of the squares of all residuals.


Type IntentOptional Attributes Name
integer, intent(in) :: is

Index of species for which LM_nonlinear_fit is executed.

double precision, intent(in) :: g(0:npar)

Array of function to be fitted.

integer, intent(in) :: n_params

Total number of fit parameters for a given species.

integer, intent(in) :: nJT

First dimension of matrix JT (see determine_JT).

double precision, intent(inout) :: params(n_params)

Array of fit parameters.

logical, intent(in) :: param_mask(n_params)

Bit mask for required fit parameters.

integer, intent(in) :: iperp

Index of perpendicular momentum.

integer, intent(in) :: npar

Number of steps in parallel momentum

integer, intent(in) :: ipparbar_lower

Lower index of parallel momentum (relativistic).

double precision, intent(out) :: quality

Quality of the individual fit result.


proc~~lm_nonlinear_fit~~CallsGraph proc~lm_nonlinear_fit LM_nonlinear_fit dgetrf dgetrf proc~lm_nonlinear_fit->dgetrf dgetri dgetri proc~lm_nonlinear_fit->dgetri proc~alps_error alps_error proc~lm_nonlinear_fit->proc~alps_error proc~determine_jt determine_JT proc~lm_nonlinear_fit->proc~determine_jt proc~fit_function fit_function proc~lm_nonlinear_fit->proc~fit_function mpi_abort mpi_abort proc~alps_error->mpi_abort proc~output_time output_time proc~alps_error->proc~output_time

Called by

proc~~lm_nonlinear_fit~~CalledByGraph proc~lm_nonlinear_fit LM_nonlinear_fit proc~determine_param_fit determine_param_fit proc~determine_param_fit->proc~lm_nonlinear_fit program~alps alps program~alps->proc~determine_param_fit