155 |
print *, ' OPTIM_SUB xx(1) = ', xx(1) |
print *, ' OPTIM_SUB xx(1) = ', xx(1) |
156 |
print *, ' OPTIM_SUB adxx(1) = ', adxx(1) |
print *, ' OPTIM_SUB adxx(1) = ', adxx(1) |
157 |
|
|
|
c compute expected decrease of cost function from objf and fmin; |
|
|
c this value is only used for a cold start of m1qn3_offline, for a |
|
|
c warm start df1 is overwritten with data from a restart file |
|
|
df1=objf-fmin |
|
|
if ( df1 .le. 0. ) then |
|
|
print *, ' OPTIM_SUB: df1 = objf-fmin = ', df1, ' should be > 0.' |
|
|
stop 'ABNORMAL in S/R OPTIM_SUB' |
|
|
endif |
|
|
|
|
|
c global variable coldStart is set in s/r optim_readparms |
|
|
c the default is false, always set it to true for the 0th cycle |
|
|
if ( optimcycle .eq. 0 ) coldStart=.true. |
|
158 |
if ( coldStart ) then |
if ( coldStart ) then |
159 |
c-- cold start |
c-- cold start |
160 |
print *, ' OPTIM_SUB: cold start, optimcycle =', optimcycle |
print *, ' OPTIM_SUB: cold start, optimcycle =', optimcycle |
163 |
c that needs to be initialized here to make sure that we have a |
c that needs to be initialized here to make sure that we have a |
164 |
c clean start |
c clean start |
165 |
reentry = 0 |
reentry = 0 |
166 |
|
c compute expected decrease of cost function from objf and fmin; |
167 |
|
c this value is only used for a cold start of m1qn3_offline, for a |
168 |
|
c warm start df1 is overwritten with data from a restart file |
169 |
|
df1=objf-fmin |
170 |
|
if ( df1 .le. 0. ) then |
171 |
|
print *, ' OPTIM_SUB: df1 = objf-fmin = ', df1, |
172 |
|
& ' but should be > 0.' |
173 |
|
stop 'ABNORMAL in S/R OPTIM_SUB' |
174 |
|
endif |
175 |
|
|
176 |
c open output file for m1qn3 |
c open output file for m1qn3 |
177 |
open(io,file=fname_m1qn3,status='unknown') |
open(io,file=fname_m1qn3,status='unknown') |
178 |
else |
else |