ADDED   mttroot/mtt/bin/trans/lang_header
Index: mttroot/mtt/bin/trans/lang_header
==================================================================
--- /dev/null
+++ mttroot/mtt/bin/trans/lang_header
@@ -0,0 +1,139 @@
+#!/bin/sh
+
+# Generates function header for a computer language
+# PJ Gawthrop May 1998
+
+###############################################################
+## Version control history
+###############################################################
+## $Id$
+## $Log$
+###############################################################
+
+
+# Set up the language specific stuff
+system=$1
+rep=$2
+language=$3
+args=$4
+output=$5
+extras=$6
+
+case $language in
+    m)
+	ext='m';
+	Lc='#';
+	Rc='#';
+	Lb='(';
+	Rb=')';
+	;;
+    c)
+	ext='c';
+	Lc='/*';
+	Rc='*/';
+	Lb='[';
+	Rb=']';
+	;;
+    *)
+	echo Language $language not supported - sorry; exit 1
+esac
+
+######### Octave (matlab) code
+if [ "$language" = m ]; then
+cat <<EOF
+function $output = $1_$rep($args);
+% $output = $1_ode($args);
+%System $system, representation $rep, language $language;
+%File $1_$rep.$ext;
+%Generated by MTT on `date`;
+%
+EOF
+fi
+
+######### c code
+if [ "$language" = c ]; then
+cat <<EOF
+void $1_$rep() {
+/*
+System $system, representation $rep, language $language;
+File $1_$rep.$ext;
+Generated by MTT on `date`;
+
+NB Arrays are be defined to be one larger than expected
+ - the 0 element is not used.
+*/
+
+EOF
+
+
+# Declarations
+$MATRIX -q  <<EOF 
+  %System structure
+  [nx,ny,nu,nz,nyz] = $1_def;
+
+  printf("$Lc Declare standard arrays $Rc\n");
+  printf("  extern double y[%i]; \t $Lc $1_ode output $Rc\n", ny+1);
+  printf("  extern double dx[%i]; \t $Lc $1_ode state derivative $Rc\n", nx+1);
+  printf("  extern double x[%i]; \t $Lc $1_ode state $Rc\n", nx+1);
+  printf("  extern double u[%i]; \t $Lc $1_ode input $Rc\n\n", nu+1);
+
+  printf("$Lc Dummy variable list $Rc\n");
+  printf("  double t0;\n\n");
+
+
+  printf("$Lc State variable list $Rc\n");
+  for i = 1:nx
+    printf("  double mttx%i\n",i);
+  end;
+printf("\n");
+
+printf("$Lc State variable input list $Rc\n");
+for i = 1:nu
+  printf("  double mttu%i\n",i);
+end;
+printf("\n");
+
+printf("$Lc Counter $Rc\n");
+printf("  int i;\n");
+printf("\n");
+
+EOF
+
+fi
+
+# Any extra bits ..
+if [ -n "$extras" ]; then
+cat <<EOF
+$Lc Extra bits here $Rc
+$extras
+EOF
+fi
+
+# Globals
+sympar2global_txt2lang $1 $language
+
+# Common part - ode
+# Use octave for this bit - needs the definition file
+$MATRIX -q  <<EOF 
+  %System structure
+  [nx,ny,nu,nz,nyz] = $1_def;
+
+
+printf("$Lc======  Read in the input ======$Rc\n");
+printf("mttu = $1_input(mttx,t)\n");
+printf("\n");
+
+
+printf("$Lc====== Set up the state variables ======$Rc\n");
+for i = 1:nx
+  printf("  mttx%i = mttx$Lb%i$Rb;\n",i,i);
+end;
+printf("\n");
+
+printf("$Lc====== Set up the input variables ======$Rc\n");
+for i = 1:nu
+  printf("  mttu%i = mttu$Lb%i$Rb;\n",i,i);
+end;
+printf("\n");
+
+EOF

ADDED   mttroot/mtt/bin/trans/ode_r2lang
Index: mttroot/mtt/bin/trans/ode_r2lang
==================================================================
--- /dev/null
+++ mttroot/mtt/bin/trans/ode_r2lang
@@ -0,0 +1,214 @@
+#! /bin/sh
+
+     ###################################### 
+     ##### Model Transformation Tools #####
+     ######################################
+
+# Bourne shell script: ode_r2lang
+# Reduce ODE to general language ODE
+# Based on obsolete ode_r2m  
+# P.J.Gawthrop 14 June 1991, 12 Jan 1994, April 1994, Jan 95, May 1998
+# Copyright (c) P.J.Gawthrop 1991, 1994, 1995, 1996, 1998
+
+###############################################################
+## Version control history
+###############################################################
+## $Id$
+## $Log$
+###############################################################
+
+# Set up the language specific stuff
+language=$2
+case $language in
+    m)
+	ext='m';
+	codegenerator='FORTRAN'
+	Lc='#';
+	Rc='#';
+	Lb='(';
+	Rb=')';
+	;;
+    c)
+	ext='c';
+	codegenerator='C'
+	Lc='/*';
+	Rc='*/';
+	Lb='[';
+	Rb=']';
+	;;
+    *)
+	echo Language $2 not supported - sorry; exit 1
+esac
+
+#Inform user
+echo Creating $1_ode.$ext
+echo Creating $1_odea.$ext
+echo Creating $1_odeo.$ext
+
+
+# Remove the old log file
+rm -f ode_r2lang.log
+
+#Remove the temporary files
+rm -f $1_ode.$ext.1
+rm -f $1_odea.$ext.1;
+rm -f $1_odeo.$ext.1;
+
+# Use reduce to accomplish the transformation
+$SYMBOLIC >ode_r2lang.log << EOF
+
+%Read the reduce definitions file
+in "$1_def.r";
+
+%Read the reduce  ODE   file
+in "$1_ode.r";
+
+
+% Set up the code generator
+% Load the general translator package
+LOAD GENTRAN;
+GENTRANLANG!* := '$codegenerator;
+ON GENTRANSEG;
+MAXEXPPRINTLEN!* := 40;
+
+% Matrix output function
+in"$MTTPATH/trans/lang_matrix.r";
+
+%Set up the number of argument variables to zero in case the user has forgotten
+MTTNVar := 0;
+
+%Read the parameter file
+in "$1_sympar.r";
+
+% The body of the ode function
+  GENTRANOUT "$1_ode.$ext.1";
+  mtt_matrix := MTTdX$
+  mtt_matrix_n := MTTNx$
+  mtt_matrix_m := 1$
+  mtt_matrix_name := MTTdX$
+  lang_matrix(); 
+  GENTRANSHUT "$1_ode.$ext.1";
+  
+% 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$
+  mtt_matrix_name := MTTy$
+  lang_matrix(); 
+  GENTRANSHUT "$1_odeo.$ext.1";
+
+EOF
+
+# Create the ode.$ext function
+lang_header $1 ode $ext 'mttx,t' mttdx > $1_ode.$ext
+
+# Algebraic bits
+$MATRIX -q  <<EOF  >> $1_ode.$ext
+  %System structure
+  [nx,ny,nu,nz,nyz] = $1_def;
+if nyz>0
+  printf("global mttxx mttt;\n");
+  printf("mttxx = mttx; mttt=t;\n");
+
+  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("mttui = fsolve('$1_odea',mttui);\n");
+  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");
+
+  printf("$Lc====== Put internal input variables at end of state derivative ======$Rc\n");
+  for i = 1:nyz
+    printf(" mttdx$Lb%i$Rb = mttui%i;\n",i+nx,i);
+  end;
+  printf("\n");
+end;
+
+EOF
+
+cat <<EOF >> $1_ode.$ext
+
+% The differential equations
+EOF
+
+sed 's/mtt_matrix/mttdx/' $1_ode.$ext.1 >> $1_ode.$ext
+
+
+
+# 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;
+if nyz>0
+  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
+
+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
+lang_header $1 odeo $ext 'mttx,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;
+
+EOF
+
+
+cat <<EOF >> $1_odeo.$ext
+
+% The output equations
+EOF
+sed 's/mtt_matrix/mtty/' $1_odeo.$ext.1 >> $1_odeo.$ext
+
+
+
+
+
+
+