8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
# Copyright (C) 2000 by Peter J. Gawthrop
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.37 2001/07/13 04:54:04 geraint
## Branch merge: numerical-algebraic-solution back to main.
##
## Revision 1.36 2001/07/12 04:00:51 gawthrop
## Now zeros y correctly - ie Ny NOT Nx elements
##
## Revision 1.35 2001/06/13 10:39:51 gawthrop
|
>
>
>
|
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
# Copyright (C) 2000 by Peter J. Gawthrop
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.38 2001/07/23 23:43:15 gawthrop
## header only version does not need to compute sizes from _def.r
##
## Revision 1.37 2001/07/13 04:54:04 geraint
## Branch merge: numerical-algebraic-solution back to main.
##
## Revision 1.36 2001/07/12 04:00:51 gawthrop
## Now zeros y correctly - ie Ny NOT Nx elements
##
## Revision 1.35 2001/06/13 10:39:51 gawthrop
|
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
|
# Argument specific stuff
get_arg_specific_stuff ()
{
arg_name=${1:-""}
case ${arg_name} in
mtta | mtte)
arg_type="Matrix"
arg_size="(MTTNX,MTTNX)"
;;
mttb)
arg_type="Matrix"
arg_size="(MTTNX,MTTNU)"
;;
mttc)
arg_type="Matrix"
arg_size="(MTTNY,MTTNX)"
;;
mttd)
arg_type="Matrix"
arg_size="(MTTNY,MTTNU)"
;;
mttax | mttdx | mttedx | mttx | mttopen)
arg_type="ColumnVector"
arg_size="(MTTNX)"
;;
mttpar)
arg_type="ColumnVector"
arg_size="(MTTNPAR)"
;;
mttsimpar_map)
arg_type="Octave_map"
arg_size="";
;;
mttu)
arg_type="ColumnVector"
arg_size="(MTTNU)"
;;
mttyz)
arg_type="ColumnVector"
arg_size="(MTTNYZ)"
;;
mtty)
arg_type="ColumnVector"
arg_size="(MTTNY)"
;;
mttt)
arg_type="const double"
arg_size=""
;;
nil)
arg_type="void *"
arg_size="(0x0)"
;;
*)
echo "Argument ${arg} not supported - sorry"; exit 1
esac
}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
|
# Argument specific stuff
get_arg_specific_stuff ()
{
arg_name=${1:-""}
case ${arg_name} in
mtta | mtte)
arg_type="Matrix"
arg_init="(MTTNX,MTTNX,0.0)"
;;
mttb)
arg_type="Matrix"
arg_init="(MTTNX,MTTNU,0.0)"
;;
mttc)
arg_type="Matrix"
arg_init="(MTTNY,MTTNX,0.0)"
;;
mttd)
arg_type="Matrix"
arg_init="(MTTNY,MTTNU,0.0)"
;;
mttax | mttdx | mttedx | mttx | mttopen)
arg_type="ColumnVector"
arg_init="(MTTNX,0.0)"
;;
mttpar)
arg_type="ColumnVector"
arg_init="(MTTNPAR,0.0)"
;;
mttsimpar_map)
arg_type="Octave_map"
arg_init="";
;;
mttu)
arg_type="ColumnVector"
arg_init="(MTTNU,0.0)"
;;
mttyz)
arg_type="ColumnVector"
arg_init="(MTTNYZ,0.0)"
;;
mtty)
arg_type="ColumnVector"
arg_init="(MTTNY,0.0)"
;;
mttt)
arg_type="const double"
arg_init=""
;;
nil)
arg_type="void *"
arg_init="(0x0)"
;;
*)
echo "Argument ${arg} not supported - sorry"; exit 1
esac
}
|
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
|
s=${1:-""} # comma separated output list
c=`get_field ${s} 0` # count of outputs
i=0
while [ ${i} -lt ${c} ]; do
i=`expr ${i} + 1`
w=`get_field ${s} ${i}` # argument name
get_arg_specific_stuff ${w}
printf " ${arg_type}\t${w}\t${arg_size};\n"
done
}
return_DLD_outputs ()
{
s=${1:-""} # comma separated output list
|
|
|
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
|
s=${1:-""} # comma separated output list
c=`get_field ${s} 0` # count of outputs
i=0
while [ ${i} -lt ${c} ]; do
i=`expr ${i} + 1`
w=`get_field ${s} ${i}` # argument name
get_arg_specific_stuff ${w}
printf " static ${arg_type}\t${w}\t${arg_init};\n"
done
}
return_DLD_outputs ()
{
s=${1:-""} # comma separated output list
|
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
|
write_DLD_header
map_DLD_inputs ${args}
undeclared=`get_extra_fields ${args:-"nil"} ${output:-"nil"}`
declare_DLD_outputs ${undeclared}
array2constant
case ${arg_type} in
Matrix)
printf " ${w}\t\t\t= zeros ${arg_size};\n"
;;
*)
;;
esac
cat <<EOF
|
|
|
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
|
write_DLD_header
map_DLD_inputs ${args}
undeclared=`get_extra_fields ${args:-"nil"} ${output:-"nil"}`
declare_DLD_outputs ${undeclared}
array2constant
case ${arg_type} in
Matrix)
printf " ${w}\t\t\t= zeros ${arg_init};\n"
;;
*)
;;
esac
cat <<EOF
|