Index: mttroot/mtt/bin/trans/mtt_optimise.sh
==================================================================
--- mttroot/mtt/bin/trans/mtt_optimise.sh
+++ mttroot/mtt/bin/trans/mtt_optimise.sh
@@ -1,37 +1,50 @@
 #! /bin/sh
 
 set -e
+set -x
 
 filename="$1"
 
-# get expression to optimise
-
-lhs=`cat ${filename} | gawk '{ print $1 }'`
-rhs=`cat ${filename} | gawk '{ print $2 }'`
-
 # set work file names
 
 tmp="$0.tmp"
 log="$0.log"
 
-
-cat <<EOF | reduce > $log
+optimise_with_reduce ()
+{
+    cat <<EOF | reduce > $log
 off echo$
 load scope$
 on double$
 on noconvert$
 on rounded$
 off int$
 off nat$
 
 out "${tmp}"$
-optimize mtt_tmp0 :=: ${rhs} iname mtt_tmp$
-write "${lhs} := mtt_tmp0"$
+optimize mtt_tmp0 :=: `cat ${filename} | gawk '{ print $2 }'` iname mtt_tmp$
+write "`cat ${filename} | gawk '{ print $1 }'` := mtt_tmp0"$
 shut "${tmp}"$
 
 ;end;
 EOF
+}
+
+optimise_with_maxima ()
+{
+    cat ${filename}					|\
+	tr "\t\n" "  "					|\
+	sed -e 's/\ //g'				|\
+	sed -e 's/\(MTTYz([0-9.]*,[0-9.]*)\)/\1 /'     	|\
+	sed -e 's/[A-Za-z0-9_.:()]*(D2)\(.*\)/\1/'	|\
+	sed -e 's/\([A-Za-z0-9_.:()]*\)(C3)/\1/g'	|\
+	optimise_max2r.awk				|\
+	tr "A-Z" "a-z"					> ${tmp}
+}
+
 
+optimise_with_maxima
+#cat $filename
 mv ${tmp} ${filename}
 
 exit 0

ADDED   mttroot/mtt/bin/trans/optimise_max2r.awk
Index: mttroot/mtt/bin/trans/optimise_max2r.awk
==================================================================
--- /dev/null
+++ mttroot/mtt/bin/trans/optimise_max2r.awk
@@ -0,0 +1,49 @@
+#! /usr/bin/awk -f
+
+# Maxima optimized expression format:
+# (C1) DISPLAY2D:FALSE$
+# (C2) OPTIMIZE(x);
+# (D2) x
+# (C3) OPTIMIZE((1.0+a*b)/(2.0+a*b));
+# (D3) BLOCK([%1],%1:a*b,(%1+1.0)/(%1+2.0))
+
+function decipher_block(lhs,rhs)
+{
+  s=rhs;
+  s=gensub(/BLOCK/, "", "1", s);		# strip BLOCK keyword
+  s=gensub(/\(/, "", "1", s);			# strip leading (
+  s=gensub(/\)$/, "", "1", s);			# strip trailing )
+  s=gensub(/\[.*\],/, "", "g", s);		# strip declarations
+  s=gensub(/%/, "mtt_tmp", "g", s);		# rename tmp variables
+  s=gensub(/:/, " := ", "g", s);	        # fix assignment operator
+
+  # commas delimit statements in maxima, so
+  # preserve commas within function calls
+  # replace remaining commas with semi-colons
+  # then put commas back in functions and reinstate parentheses
+  while (match(s, (/\(([^()]*),([^()]*)\)/))) {
+    s=gensub(/\(([^()]*),([^()]*)\)/, "(\\1__MTT_COMMA__\\2)", "g", s);
+    s=gensub(/\(([^,()]*)\)/, "__MTT_LPAREN__\\1__MTT_RPAREN__", "g", s);
+  }
+  s=gensub(/,/, ";\n", "g", s);
+  s=gensub(/__MTT_COMMA__/, ",", "g", s);
+  s=gensub(/__MTT_LPAREN__/, "(", "g", s);
+  s=gensub(/__MTT_RPAREN__/, ")", "g", s);
+
+  # separate statements
+  rhs=gensub(/.*;([^;]*)\n/, "\\1", 1, s);
+  s=gensub(/;([^;]*)$/, ";", 1, s);
+
+  print s;					# print intermediate statements
+  printf "%s := %s;\n", lhs, rhs;		# print final assignment
+}  
+
+# pattern matching:
+
+(! /.* BLOCK.*/ ) { 				# unoptimised statement
+  print $1" := "$2";";
+}
+
+(  /.* BLOCK.*/ ) {				# optimised statement
+  decipher_block($1,$2)
+}