| File |
Line number |
Procedure |
Code |
|
./lsopt/dgscale.F |
24 |
DGSCALE |
double precision gold(nn), xdiff(nn), diag(nn) |
|
./lsopt/dgscale.F |
2 |
DGSCALE |
subroutine DGSCALE( nn, gold, xdiff, diag, rmin ) |
|
./lsopt/dgscale.F |
33 |
DGSCALE |
call DOSTORE( nn, diag, .false., 3 ) |
|
./lsopt/dgscale.F |
37 |
DGSCALE |
r1 = r1 + gold(i)*gold(i)*diag(i) |
|
./lsopt/dgscale.F |
41 |
DGSCALE |
call SSCAL( nn, r1, diag, 1 ) |
|
./lsopt/dgscale.F |
52 |
DGSCALE |
if (diag(i).LE.0) then |
|
./lsopt/dgscale.F |
54 |
DGSCALE |
diag(i) = rmin |
|
./lsopt/dgscale.F |
57 |
DGSCALE |
den = den + xdiff(i)*xdiff(i) / diag(i) |
|
./lsopt/dgscale.F |
61 |
DGSCALE |
diag(i) = 1./ |
|
./lsopt/dgscale.F |
62 |
DGSCALE |
$ (1./diag(i)+gold(i)**2-(xdiff(i)/diag(i))**2/den) |
|
./lsopt/dgscale.F |
63 |
DGSCALE |
if (diag(i).le.0.) then |
|
./lsopt/dgscale.F |
64 |
DGSCALE |
diag(i) = rmin |
|
./lsopt/dgscale.F |
71 |
DGSCALE |
call DOSTORE( nn, diag, .true., 3 ) |