1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
|
include <math/nlfit.h>
include "nlfitdefr.h"
# NLSOLVE -- Procedure to solve nonlinear system
procedure nlsolver (nl, ier)
pointer nl # pointer to the nlfit structure
int ier # error code
int nfree
begin
# Make temporary arrays.
call amovr (ALPHA(NL_ALPHA(nl)), COVAR(NL_COVAR(nl)),
NL_NFPARAMS(nl) ** 2)
call amovr (BETA(NL_BETA(nl)), DPARAM(NL_DPARAM(nl)), NL_NFPARAMS(nl))
# Add the lambda damping factor.
call nl_dampr (COVAR(NL_COVAR(nl)), COVAR(NL_COVAR(nl)),
(1.0 + NL_LAMBDA(nl)), NL_NFPARAMS(nl), NL_NFPARAMS(nl))
ier = OK
nfree = NL_NPTS(nl) - NL_NFPARAMS(nl)
if (nfree < 0) {
ier = NO_DEG_FREEDOM
return
}
# Compute the factorization of the matrix.
call nl_chfacr (COVAR(NL_COVAR(nl)), NL_NFPARAMS(nl), NL_NFPARAMS(nl),
CHOFAC(NL_CHOFAC(nl)), ier)
# Solve the equations for the parameter increments.
call nl_chslvr (CHOFAC(NL_CHOFAC(nl)), NL_NFPARAMS(nl), NL_NFPARAMS(nl),
DPARAM(NL_DPARAM(nl)), DPARAM(NL_DPARAM(nl)))
end
|
