142 |
epsx = eps |
epsx = eps |
143 |
epsg = eps |
epsg = eps |
144 |
endif |
endif |
145 |
|
c always force cold start for the 0th cycle |
146 |
|
if ( optimcycle .eq. 0 ) coldStart=.true. |
147 |
|
|
148 |
call optim_readdata ( nn, ctrlname, .true., ff, vv) |
c read header from costname rather than ctrlname, because the |
149 |
|
c cost function value in costname is what we need, ctrlname only |
150 |
|
c contains a valid cost function value for the 0th iteration |
151 |
|
call optim_readdata ( nn, costname, .true., ff, vv) |
152 |
|
|
153 |
if ( dfminFrac.lt.0.0 .or. dfminFrac.ge.1.0) then |
if ( dfminFrac.lt.0.0 .or. dfminFrac.ge.1.0) then |
154 |
print*, ' OPTIM_READPARMS: dfminFrac = ', dfminFrac, |
print*, ' OPTIM_READPARMS: dfminFrac = ', dfminFrac, |
157 |
endif |
endif |
158 |
if ( dfminFrac.ne.0.0 ) dfminFrac = 1.0 - dfminFrac |
if ( dfminFrac.ne.0.0 ) dfminFrac = 1.0 - dfminFrac |
159 |
if ( fmin.eq.UNSET_RL ) then |
if ( fmin.eq.UNSET_RL ) then |
160 |
if ( optimcycle .eq. 0 ) then |
if ( coldStart ) then |
161 |
c only in this case does ff contain the actual cost function value |
c only in this case does ff contain the actual cost function value |
162 |
fmin = dfminFrac*ff |
fmin = dfminFrac*ff |
163 |
else |
else |
165 |
fmin = ff |
fmin = ff |
166 |
endif |
endif |
167 |
endif |
endif |
168 |
if ( optimcycle.eq.0 ) then |
if ( coldStart ) then |
169 |
if ( fmin.le.0.0 ) then |
if ( fmin.le.0.0 ) then |
170 |
print '(A,E12.6,A)', ' OPTIM_READPARMS: fmin = ', fmin, |
print '(A,E12.6,A)', ' OPTIM_READPARMS: fmin = ', fmin, |
171 |
& ' should be > 0' |
& ' should be > 0' |
183 |
print*, ' OPTIM_READPARMS: Iteration number = ', optimcycle |
print*, ' OPTIM_READPARMS: Iteration number = ', optimcycle |
184 |
print*, ' number of control variables = ', nn |
print*, ' number of control variables = ', nn |
185 |
print*, ' cost function value in ', ctrlname, ' = ', ff |
print*, ' cost function value in ', ctrlname, ' = ', ff |
186 |
|
if ( coldStart ) then |
187 |
|
print*, ' expected cost function minimum = ', fmin |
188 |
|
print*, ' expected cost function decrease = ', ff-fmin |
189 |
|
endif |
190 |
print '(a,4a,i4.4)', |
print '(a,4a,i4.4)', |
191 |
& ' Data will be read from the following file: ', |
& ' Data will be read from the following file: ', |
192 |
& ctrlname,'_',yctrlid(1:10),'.opt', optimcycle |
& ctrlname,'_',yctrlid(1:10),'.opt', optimcycle |