# Copyright (c) P.J.Gawthrop 1991, 1994, 1995, 1996, 1998

###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.5  1998/08/13 16:25:07  peterg
## Removed some documentation form odeo.m
##
## Revision 1.4  1998/08/13 15:10:47  peterg
## Change temp var name to mtt_t1 etc
##
## Revision 1.3  1998/08/13 12:43:19  peterg
## Cahnged tmp variable to mtt_1 etc
##
## Revision 1.2  1998/07/27 20:26:44  peterg

ext='m';
codegenerator='PASCAL'
Lc='#';
Rc='#';
Lb='(';
Rb=')';

outfile=$1"_"$rep.$ext
outfile=$1"_"$rep

rm -f $outfile.*

echo $outfile
rm -fv $outfile.[123]
#Set up gentran
cat <<EOF >mtt_setreduce.r
% Set up the code generator
% Load the general translator package
LOAD GENTRAN;
GENTRANLANG!* := '$codegenerator;
ON GENTRANSEG;

%Read the reduce definitions file
in "$1_def.r";

%Read the reduce  $REP   file
in "$1_$rep.r";

%Read the parameter file
in "$1_sympar.r";

EOF

#Inform user
echo Creating $outfile
echo Creating $outfile.m



# Remove the old log file
rm -f mtt_r2lang.log
rm -f mtt_r2m.log

#Remove the temporary files
rm -f $outfile.*

if [ "$rep" = "ode" ]; then
  #echo Creating $1_odea.$ext
  outfileo=$outfile"o"
  echo Creating $1_odeo.$ext
  rm -f $1_odeo.$ext.1;
  echo Creating $outfileo.m
   rm ode odeo
  #rm -f $1_odea.$ext.1;

  #Header
  lang_header $1 $rep $ext 'mttx,mttu,t' mttdx > $outfile.1
  lang_header $1 $rep $ext 'mttx,mttu,t' mttdx > $1_ode.p
  lang_header $1 odeo $ext 'mttx,mttu,t' mtty > $1_odeo.p

  # Use reduce to accomplish the transformation
  $SYMBOLIC >mtt_r2lang.log << EOF
    
$SYMBOLIC  > mtt_r2m.log <<EOF
  %Set up the code generator
  in"mtt_setreduce.r";
      
  % The body of the ode function
    GENTRANOUT "$outfile.3";
    GENTRANOUT "ode";
    mtt_matrix := MTTdX$
    mtt_matrix_n := MTTNx$
    mtt_matrix_m := 1$
    mtt_matrix_name := 'MTTdX$
    lang_matrix(); 
    GENTRANSHUT "$outfile.3";
    lang_matrix();
    GENTRANSHUT "ode";
    
%%      % The algebraic equations (if any)
%%    GENTRANOUT "$1_odea.$ext.1";
%%    mtt_matrix := MTTYz$
%%    mtt_matrix_n := MTTNYz$
%%    mtt_matrix_m := 1$
%%    mtt_matrix_name := MTTYz$
%%    lang_matrix(); 
%%    GENTRANSHUT "$1_odea.$ext.1;
      
 
      
      % Now do the y = g(x,t) function.
      % The body of the odeo function
    GENTRANOUT "$1_odeo.$ext.1";
    mtt_matrix := MTTy$
    mtt_matrix_n := MTTNy$
    mtt_matrix_m := 1$
     GENTRANOUT "odeo";
     mtt_matrix := MTTy$
     mtt_matrix_n := MTTNy$
     mtt_matrix_m := 1$
    mtt_matrix_name := MTTy$
    lang_matrix(); 
    GENTRANSHUT "$1_odeo.$ext.1";
      
     lang_matrix(); 
     GENTRANSHUT "odeo";

END;;


EOF

     
#  cat $1_odeo.$ext.1

##  # Algebraic bits
##  $MATRIX -q  <<EOF  >> $1_ode.$ext
##    %System structure
##    [nx,ny,nu,nz,nyz] = $1_def;
##  if nyz>0

##      printf(" mttdx$Lb%i$Rb = mttui%i;\n",i+nx,i);
##    end;
##    printf("\n");
##  end;
##  
##EOF
##  
  cat <<EOF >> $outfile.2
  
  $Lc The differential equations $Rc
EOF
  
  sed 's/mtt_matrix/mttdx/' $outfile.3 > $outfile.4
  rm -f $outfile.3
  
##  # Create the odea.$ext function
##  lang_header $1 odea $ext 'mttui' mttyz 'global mttxx mttt; mttx = mttxx; t=mttt;' > $1_odea.$ext
##  
##  $MATRIX -q  <<EOF  >> $1_odea.$ext
##    %System structure
##    [nx,ny,nu,nz,nyz] = $1_def;

##  cat <<EOF >> $1_odea.$ext
##  
##  % The algebraic equations
##EOF
##  sed 's/mtt_matrix/mttyz/' $1_odea.$ext.1 >> $1_odea.$ext
  
  
  # Create the odeo.$ext function
  # Create the odeo function
  lang_header $1 odeo $ext 'mttx,mttu,t' mtty > $1_odeo.$ext
  
$MATRIX -q  <<EOF  >> $1_odeo.$ext
    %System structure
    [nx,ny,nu,nz,nyz] = $1_def;
  if nyz>0
    printf("$Lc====== Extract the internal input variables ======$Rc\n");
    for i = 1:nyz
      printf("  mttui$Lb%i$Rb = mttx$Lb%i$Rb;\n",i,i+nx);
    end;
    printf("\n");
  
    printf("$Lc====== Set up the internal input variables ======$Rc\n");
    for i = 1:nyz
      printf("  mttui%i = mttui$Lb%i$Rb;\n",i,i);
    end;
    printf("\n");
  end;
##$MATRIX -q  <<EOF  >> $1_odeo.p
##    %System structure
##    [nx,ny,nu,nz,nyz] = $1_def;
##  if nyz>0
##    printf("$Lc====== Extract the internal input variables ======$Rc\n");
##    for i = 1:nyz
##      printf("  mttui$Lb%i$Rb = mttx$Lb%i$Rb;\n",i,i+nx);
##    end;
##    printf("\n");
##  
##    printf("$Lc====== Set up the internal input variables ======$Rc\n");
##    for i = 1:nyz
##      printf("  mttui%i = mttui$Lb%i$Rb;\n",i,i);
##    end;
##    printf("\n");
##  end;
  
EOF
##EOF

sed 's/mtt_matrix/mttdx/' <ode  >> $1_ode.p
sed 's/mtt_matrix/mtty/' $1_odeo.$ext.1 >> $1_odeo.$ext
sed 's/mtt_matrix/mtty/'  <odeo >> $1_odeo.p

# Convert from Pascal syntax to Octave syntax
  mv $1_odeo.m mtt_junk; mtt_p2m<mtt_junk > $1_odeo.m
# Convert from Pascal syntax to Octave syntax 
   mtt_p2m<$1_odeo.p > $1_odeo.m 
   mtt_p2m<$1_ode.p > $1_ode.m 
#  mv $1_odea.m mtt_junk; mtt_p2m<mtt_junk > $1_odea.m

fi
  
if [ "$rep" = "sm" ]||[ "$rep" = "dm" ]; then

if [ "$rep" = "sm" ]; then
   #Header
  lang_header $1 $rep $ext '' [mtta,mttb,mttc,mttd] > $outfile.A

    mtt_matrix := MTT$Symbol;
    mtt_matrix_n := $N;
    mtt_matrix_m := $M;
    lang_matrix();
    GENTRANSHUT "$outfile.1";
EOF

echo "% $Symbol matrix" >$outfile.$Symbol
cat "$outfile.2" >> $outfile.$Symbol
echo "%$Symbol matrix#" > $outfile.$Symbol
cat "$outfile.2"  >> $outfile.$Symbol
sed "s/mtt_matrix/mtt$Symbol/" < $outfile.1 |\
grep -v '=0\.0;' >> $outfile.$Symbol
grep -v '=0\.0;'  >> $outfile.$Symbol

done

rm -f $outfile.[12]
fi

fi 

# Convert from Pascal syntax to Octave syntax
cat $outfile.?  | mtt_p2m > $outfile

rm -f $outfile.?  
#rm -f $outfile.?