| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=pr), | intent(in) | :: | J(:,:) |
Jacobian at x |
||
| real(kind=pr), | intent(in) | :: | F(:) |
Residual at x |
||
| real(kind=pr), | intent(inout) | :: | lambda |
LM damping (absolute) |
||
| real(kind=pr), | intent(in) | :: | jscale |
||J||_F^2/n |
||
| type(newton_settings), | intent(in) | :: | s |
Solver settings |
||
| real(kind=pr), | intent(out) | :: | dX(:) |
Computed Newton/LM step |
||
| real(kind=pr), | intent(out) | :: | cond_est |
cond(A) estimate |
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| real(kind=pr), | private, | allocatable | :: | A(:,:) |
Working copy of J (or J^T J + lam I) |
||
| real(kind=pr), | private | :: | anorm |
1-norm of A (input to dgecon) |
|||
| integer, | private | :: | info |
LAPACK return code (0 = success) |
|||
| integer, | private, | allocatable | :: | ipiv(:) |
Pivot indices from dgetrf, size n |
||
| integer, | private, | allocatable | :: | iwork(:) |
Integer workspace for dgecon, size n |
||
| integer, | private | :: | n |
Problem size |
|||
| real(kind=pr), | private | :: | rcond |
Reciprocal condition number from dgecon |
|||
| real(kind=pr), | private, | allocatable | :: | rhs(:) |
Right-hand side (-F or -J^T F) |
||
| logical, | private | :: | use_lm |
.true. when the LM branch is used |
|||
| real(kind=pr), | private, | allocatable | :: | work(:) |
LAPACK workspace, size 4n |