#! /bin/sh
######################################
##### Model Transformation Tools #####
######################################
# Bourne shell script: lde2ode.r
# Reduce differential-algebraic equations to Lagranges equations
# Copyright (C) 2000 by Peter J. Gawthrop
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.2.4.3 2002/09/13 10:10:42 geraint
## Commented out optimisation.
##
## Revision 1.2.4.2 2002/09/10 23:24:19 geraint
## Rationalised local and global optimisations.
## Fixes presentation of locally optimised code (ode view).
## Much more elegant :-)
##
## Revision 1.2.4.1 2002/09/03 23:44:43 geraint
## adding global optimisation (-optg).
##
## Revision 1.2 2001/07/27 23:29:10 geraint
## Optimises only when requested (-opt).
##
## Revision 1.1 2000/12/28 12:25:13 peterg
## Initial revision
##
###############################################################
optimise=''; optimise_msg=''
while [ -n "`echo $1 | grep '^-'`" ]; do
case $1 in
-optimise_global )
# optimise='-optimise_global'
# optimise_msg=' with global optimisation' ;;
optimise=''
optimise_msg='' ;;
-optimise_local )
# optimise='-optimise_local'
# optimise_msg=' with local optimisation' ;;
optimise=''
optimise_msg='' ;;
*)
echo "$1 is an invalid argument - ignoring" ;;
esac
shift
done
# Create the reduce output code
def2write_r $1 lde
def2write_r $1 ldeo
#Inform user
echo Creating $1_lde.r $optimise_msg
echo Creating $1_ldeo.r $optimise_msg
# Remove the old log file
rm -f cse2lde_r.log
# Use reduce to accomplish the transformation
$SYMBOLIC << EOF >cse2lde_r.log
%Read the formatting function
in "$MTTPATH/trans/reduce_matrix.r";
%Read the definitions file
in "$1_def.r";
%Read the substitution file
in "$1_subs.r";
%Read the dae file
in "$1_dae.r";
% Zap bits we dont need
clear mttui, mttu, mttdu, mttx, mttdz;
%% The lagrangian eqns are those set to zero.
MTTL := mttyz;
% Create the second derivative of z
MATRIX mttddz(MTTNz,1);
FOR i := 1:MTTNz DO
BEGIN
mttuii := mkid(mttui,i);
mttdduii := mkid(mttddui,i);
mttddz(i,1) := sub(mttuii=mttdduii,mttz(i,1));
END;
mttddz;
% Create Lagranges Equations
FOR j := 1:MTTNyz DO
BEGIN
%% Create derivative of the equations by substituting:
%du for u
mttuj := mkid(mttu,j);
mttduj := mkid(mttdu,j);
MTTL := sub(mttuj=mttduj,MTTL);
% dv for v
mttuij := mkid(mttui,j);
mttduij := mkid(mttdui,j);
MTTL := sub(mttuij=mttduij,MTTL);
% dx for x
mttxj := mkid(mttx,j);
MTTL := sub(mttxj=mttdx(j,1),MTTL);
% ddz for dz
mttdzj := mkid(mttdz,j);
MTTL := sub(mttdzj=mttddz(j,1),MTTL);
END;
L = MTTL;
OFF Echo;
OFF Nat;
%Write out the ordinary differential equations.
OUT "$1_lde.r";
write "%File: $1_lde.r";
in ("$1_lde_write.r");
write "in ""$1_ldeo.r"";";
write "END;";
SHUT "$1_lde.r";
%Write out the output equations
OUT "$1_ldeo.r";
write "%File: $1_ldeo.r";
in ("$1_ldeo_write.r");
write "END;";
SHUT "$1_ldeo.r";
quit;
EOF
# TODO: edit mtt_optimise to accept lde and ldeo
# if [ ${optimise:-""} = "-optimise_global" ]; then
# mtt_optimise global $1 lde
# mtt_optimise global $1 ldeo
# elif [ ${optimise:-""} = "-optimise_local" ]; then
# mtt_optimise local $1 lde
# mtt_optimise local $1 ldeo
# fi
# Now invoke the standard error handling.
mtt_error_r cse2lde_r.log