125 |
mm = nn |
mm = nn |
126 |
mupd = nupdate |
mupd = nupdate |
127 |
|
|
128 |
REAL_BYTE = 8 |
REAL_BYTE = 8 |
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 ' |
377 |
print iform, iiter, tact, ifunc, ff, gnorm0 |
print iform, iiter, tact, ifunc, ff, gnorm0 |
378 |
$ , DNRM2( nn, xx, 1 ), 0. |
$ , DNRM2( nn, xx, 1 ), 0. |
379 |
|
|
380 |
write(94,'(i5,2x,1pe11.4,4(2x,1pe8.1))') |
Cml write(94,'(i5,2x,1pe11.4,4(2x,1pe8.1))') |
381 |
& iiter, ff, gnorm0, tact, |
Cml & iiter, ff, gnorm0, tact, |
382 |
& DNRM2( nn, xx, 1 ), 0. |
Cml & DNRM2( nn, xx, 1 ), 0. |
383 |
|
|
384 |
if ( itmax .EQ. 0 ) then |
if ( itmax .EQ. 0 ) then |
385 |
ifail = 10 |
ifail = 10 |
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 ) |
461 |
print iform, iiter, tact, ifunc, ff, gnorm |
print iform, iiter, tact, ifunc, ff, gnorm |
462 |
$ , DNRM2( nn, xx, 1 ), tact*DNRM2( nn, dd, 1 ) |
$ , DNRM2( nn, xx, 1 ), tact*DNRM2( nn, dd, 1 ) |
463 |
|
|
464 |
write(94,'(i5,2x,1pe11.4,4(2x,1pe8.1))') |
Cml write(94,'(i5,2x,1pe11.4,4(2x,1pe8.1))') |
465 |
& iiter, ff, gnorm, tact, |
Cml & iiter, ff, gnorm, tact, |
466 |
& DNRM2( nn, xx, 1 ), tact*DNRM2( nn, dd, 1 ) |
Cml & DNRM2( nn, xx, 1 ), tact*DNRM2( nn, dd, 1 ) |
467 |
|
|
468 |
c----------------------------------------- |
c----------------------------------------- |
469 |
c check output mode of ifail |
c check output mode of ifail |