File mttroot/mtt/bin/trans/make_ssim artifact 3a0f893d5a part of check-in 3b8cbf87d1

#!/bin/sh

## Makes the sensitivity simulation
## Copyright (C) 2001 by Peter J. Gawthrop

sys=$1
lang=$2

outfile=${sys}_ssim.${lang}

echo Creating $outfile
make_m() {
mtt_header ${sys} ssim m >  ${outfile}
cat >> ${outfile} <<EOF

  ## Pass input as a global
  global MTT_input MTT_input_index MTT_input_last

  [nx2,ny2] = ${sys}_def;
  ny = ny2/2;
  y_par = [];
  MTT_input = u;
  [MTT_input_last,m] = size(u);

  for i=index
    MTT_input_index = 0;
    p = par;           # Reset parameters
    p(i) = 1;          # Set sensitivity index to 1
    [mtt_data] = ${sys}_ode2odes(x0,p,simpar);
    if (i==index(1))
      y = mtt_data(:,2:1+ny);
    endif
    y_par = [y_par, mtt_data(:,2+ny:1+2*ny)];
  endfor
endfunction

EOF

}

case ${lang} in
    m)
	make_m
	;;
    *)
	echo Language ${lang} not supported - sorry
esac