# Copyright (c) P.J.Gawthrop, 1989, 1990, 1991, 1993, 1994, 1995, 1996.

###############################################################
## Version control history
###############################################################
## $Header$
## $Log$
# Revision 1.54  1997/02/07  13:33:03  peterg
# Added ; to _input.c file - bug fix.
#
# Revision 1.53  1997/02/07  13:25:04  peterg
# Version 2.3
#
# Revision 1.52  1997/01/21  22:55:59  peterg
# Minor bug fixes
#
## Revision 1.51  1997/01/21 13:16:40  peterg
## Included transformations to language c for:
## ode
## params
## numpar
##
## Revision 1.50  1996/12/21 19:46:52  peterg

## Revision 1.1  1996/08/04 17:29:42  peter
## Initial revision
##
###############################################################

# MTT recursion level is zero unless explicitly set
level=0

#Computation mode is octave by default
computation=octave

#Look for a command line argument
while [ -n "`echo $1 | grep '-'`" ]; do
  case $1 in
	-q )
		quiet=quiet ;;
	-c )
		computation=c ;;
	-l )
		level=$2; shift ;;
  esac
  shift
done

#Print header if not in quiet (-q) mode.
if [ "$quiet" != "quiet" ]; then
  echo
  echo 'MTT (Model Transformation Tools) version 2.3'
  echo 'This is free software with ABSOLUTELY NO WARRANTY.'
  echo 'Type `mtt warranty'\' 'for details.'
  echo
fi

if [ -z "$1" ]; then
  echo 'Usage: mtt help'

        )> $1_simp.r

#SUMMARY numpar	numerical parameter declaration (m) 
$1_numpar.m: $1_sympar.r $1_params.m
	sympar2numpar_r2m $1

#SUMMARY numpar	numerical parameter declaration (c) 
$1_numpar.c: $1_sympar.r
	sympar2numpar_r2c $1

#SUMMARY params	numerical parameter setting (m)
 $1_params.m: 
	echo Creating $1_params.m
	( \
	echo "%% Parameter file for system $1 ($1_params.m)"; \

	echo Creating $1_params.c
	( \
	echo '/*'; \
	echo "%% Parameter file for system $1 ($1_params.c)"; \
	echo "%% This file provides the system params for simulation:";\
	echo '*/'; \
	echo '$1_params()'; \
        echo '{'; \
	echo '};'; \
        )> $1_params.c

#SUMMARY odes ODE simulation header file (h)
 $1_odes.h: 
	echo Creating $1_odes.h
	( \
	echo '/*'; \
	echo "%% Parameter file for system $1 ($1_odes.h)"; \
	echo "%% This file provides the params for simulation:";\
	echo '*/'; \
	echo '#define DT 0.1 /* Time step (for printing) */'; \
	echo '#define LAST 10.0 /* Last time */'; \
	echo '#define STEPFACTOR 1000 /* Integration steps per time step */'; \
        )> $1_odes.h

#SUMMARY params	symbolic parameter setting (r)
 $1_params.r: 
	echo Creating $1_params.r
	( \
	echo "%% Parameter file for system $1 ($1_params.r)"; \
	echo "%% This file provides symbolic parameters for simplification";\

	cat $MTTPATH/trans/m/rcs_header.txt; \
	echo; \
	echo "% Set up system definitions"; \
	echo "[nx,ny,nu,nz,nyz] = $1_def;"; \
	echo "% Change the following lines as you wish"; \
	echo "u = ones(nu,1);"; \
        )> $1_input.m

#SUMMARY input	input signal definition (c)
$1_input.c:
	echo Creating $1_input.c
	( \
	echo "$1_input(t)"; \
	echo '  float t;'; \
	echo '/*'; \
	echo "%% Input file for system $1 ($1_input.c)"; \
	echo "%% This file provides the system input for simulation:";\
	echo '*/'; \
	cat $MTTPATH/trans/rcs_header.c; \
	echo; \
	echo "/* Change the following lines as you wish */"; \
        echo '{'; \
	echo '};'; \
        )> $1_input.c

#SUMMARY sspar	steady-state definition (r)
$1_sspar.r:
	echo Creating $1_sspar.r
	(\
	echo "%% Reduce steady-state parameter file for $1 ($1_sspar.r)"; \
	cat $MTTPATH/trans/m/rcs_header.txt; \

#SUMMARY ode	ordinary differential equations (view)
#SUMMARY ode	ordinary differential equations (ps)
#Ordinary differential equations
$1_ode.r: $1_cse.r $1_def.r $1_sympar.r
	cse2ode_r $1; tidy $1_ode.r
$1_ode.m: $1_def.r $1_ode.r $1_sympar.r $1_numpar.m
	ode_r2m $1; matlab_tidy $1_ode.m; matlab_tidy $1_odeo.m
$1_ode.c: $1_def.r $1_ode.r $1_sympar.r
	ode_r2c $1
$1_ode.tex: $1_ode.r  $1_sympar.r $1_simp.r
	ode_r2tex $1; latex_tidy $1_ode.tex

#SUMMARY obs	observer equations for CGPC (r)
#SUMMARY obs	observer equations for CGPC (m)
#SUMMARY obs	observer equations for CGPC (tex)

#SUMMARY sro	step response - output (ps)
#SUMMARY sro	step response - output (view)

#SUMMARY odes	ode solution - state (m)
#SUMMARY odes	ode solution - state (dat)
#SUMMARY odes	ode solution - state (ps)
#SUMMARY odes	ode solution - state (view)
#SUMMARY odes	ode solution - state & output(c)

#SUMMARY odeso	ode solution - output (m)
#SUMMARY odeso	ode solution - output (dat)
#SUMMARY odeso	ode solution - output (ps)
#SUMMARY odeso	ode solution - output (view)

#SUMMARY daes	dae solution - state (m)

	sm2ir_m $1 '$ARGS'
$1_iro.m: $1_ir.m
	touch $1_iro.m
$1_sr.m: $1_sm.m $1_numpar.m $1_params.m $1_args.m $1_def.m
	sm2sr_m $1 '$ARGS'
$1_sro.m: $1_sr.m
	touch $1_sro.m
$1_odes.c: $1_ode.c $1_def.r $1_sympar.r
	ode2odes_r2c $1
ifeq ($computation,octave)
$1_odes.m: $1_ode.m $1_numpar.m $1_params.m $1_args.m $1_def.m $1_input.m
	ode2odes_m $1 '$ARGS'
endif
ifeq ($computation,c)
$1_odes.m: $1_odes.c $1_ode.c $1_numpar.c $1_odes.h $1_params.c $1_input.c
	echo Creating $1_odes.m
	echo 'Compiling ...'
	gcc $1_odes.c $1_ode.c $1_numpar.c $1_params.c $1_input.c
	echo 'Running  ...'; a.out; echo 'Done.'
endif
$1_odeso.m: $1_odes.m
	touch $1_odeso.m
$1_daes.m: $1_dae.m $1_numpar.m $1_params.m $1_args.m $1_def.m $1_input.m
	dae2daes_m $1 '$ARGS'
$1_daeso.m: $1_daes.m
	touch $1_daeso.m