125 |
mm = nn |
mm = nn |
126 |
mupd = nupdate |
mupd = nupdate |
127 |
|
|
128 |
REAL_BYTE = 8 |
REAL_BYTE = 4 |
129 |
isize = REAL_BYTE |
isize = REAL_BYTE |
130 |
|
|
131 |
c----------------------------------------- |
c----------------------------------------- |
255 |
c--- start if cold start --- |
c--- start if cold start --- |
256 |
if (lphprint) then |
if (lphprint) then |
257 |
print '(a)', 'pathei-lsopt: cold start' |
print '(a)', 'pathei-lsopt: cold start' |
|
print * |
|
258 |
end if |
end if |
259 |
|
|
260 |
|
print *, 'pathei-lsopt vor simul', nn |
261 |
|
print *, 'pathei-lsopt xx(1), gg(1) ', xx(1), gg(1) |
262 |
|
|
263 |
call simul( indic, nn, xx, ff, gg ) |
call simul( indic, nn, xx, ff, gg ) |
264 |
cph( |
|
265 |
print *, 'pathei: nach simul: nn, ff = ', nn, ff |
print *, 'pathei: nach simul: nn, ff = ', nn, ff |
266 |
print *, 'pathei: nach simul: xx, gg = ', xx(1), gg(1) |
print *, 'pathei: nach simul: xx(1), gg(1) = ', xx(1), gg(1) |
|
cph) |
|
267 |
|
|
268 |
do i = 1, nn |
do i = 1, nn |
269 |
xdiff(i) = 1. |
xdiff(i) = 1. |
279 |
cph( |
cph( |
280 |
print *, 'pathei: vor lswri ' |
print *, 'pathei: vor lswri ' |
281 |
cph) |
cph) |
282 |
call lswri( iiter, nn, xx, gg, lphprint ) |
call lswri( isize, iiter, nn, xx, gg, lphprint ) |
283 |
|
|
284 |
cph( |
cph( |
285 |
print *, 'pathei: vor gnorm0 ' |
print *, 'pathei: vor gnorm0 ' |
450 |
call dostore( nn, xx, .true., 1 ) |
call dostore( nn, xx, .true., 1 ) |
451 |
call dostore( nn, gg, .true., 2 ) |
call dostore( nn, gg, .true., 2 ) |
452 |
cph( |
cph( |
453 |
call lswri( iiter, nn, xx, gg, lphprint ) |
call lswri( isize, iiter, nn, xx, gg, lphprint ) |
454 |
cph) |
cph) |
455 |
|
|
456 |
gnorm = DNRM2( nn, gg, 1 ) |
gnorm = DNRM2( nn, gg, 1 ) |