optima

来自「COOOL：CWP面向对象最优化库（CWP Object Oriented Op」· 代码 · 共 44 行

#! /bin/sh
####################################
#
#	This shell illustrates how the optimization of the stacking
#	power objective function can be implemented by the COOOL library.
#
#	In the file optima.cc all the initialization of the conjugate 
#	gradient procedure is performed. The "expect" mechanism is 
#	responsible for the communication between the processes optima 
#	and statics, in which the stacking power is computed
#
####################################

objfcn=estimates	# executable code that computes the stacking power
dimension=62		# number of surface consistent statics parameters
lower=-20		# lower boundary of model parameters
upper=20		# upper boundary of model parameters
iterCG=$2		# number of iteration performed in the conjugate 
			# gradient
tol=0.0001  		# tolerance of the gradient error
iterLine=$1		# number of iteration performed in the line search 
verbose=$3		# dialogue flag

###################
# set flags for little/big-edian machine
# uncomment either one of the following lines for your machine
#ENDIAN=1;  	# for NeXT, IBMRS6000, ...
ENDIAN=0;  	# for PC, ...
###################

endian=$ENDIAN		# =0 flag for little-edian machine

set -x

# prepare the right kind of data for little-edian machine 
#segyread > data tape=data.segy verbose=$verbose endian=$ENDIAN 
#segyread > xcor.su tape=xcor.segy verbose=$verbose endian=$ENDIAN 
#sustrip < xcor.su > xcor

optima -S $objfcn -N $dimension -M $iterLine -C $iterCG -T $tol\
	-L $lower -U $upper -V $verbose > junk

#/bin/rm xcor.su data xcor binary header