| ︙ | | | ︙ | |
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.57 2004/06/20 10:18:00 gawthrop
## Explicitly size vectors - avoid default row vector bug
##
## Revision 1.56 2003/04/17 20:07:33 geraint
## Added rule to create _sese.m
##
## Revision 1.55 2002/10/30 01:15:20 gawthrop
|
>
>
>
>
|
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
# Copyright (C) 2000 by Peter J. Gawthrop
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.58 2004/06/20 13:34:11 gawthrop
## Create vectors using zeros(N,1) - avoids new octave
## default to row vector - more reps now included
##
## Revision 1.57 2004/06/20 10:18:00 gawthrop
## Explicitly size vectors - avoid default row vector bug
##
## Revision 1.56 2003/04/17 20:07:33 geraint
## Added rule to create _sese.m
##
## Revision 1.55 2002/10/30 01:15:20 gawthrop
|
| ︙ | | | ︙ | |
348
349
350
351
352
353
354
355
356
357
358
359
360
361
|
zeromatrices='open';
;;
numpar)
states=no;
inputs=no;
parameters=no;
output='mttpar'
;;
ode)
states=yes;
inputs=yes;
parameters=yes;
output='mttdx'
args=$eqnargs
|
>
|
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
|
zeromatrices='open';
;;
numpar)
states=no;
inputs=no;
parameters=no;
output='mttpar'
zeromatrices='par'
;;
ode)
states=yes;
inputs=yes;
parameters=yes;
output='mttdx'
args=$eqnargs
|
| ︙ | | | ︙ | |
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
|
case ${rep} in
dm)
let Nx=Nx+2*Nz+Nyz # Matrices not usual size
;;
*)
esac
##Npar=`wc -l $system\_sympar.txt | gawk '{print $1}'`
}
zero_matrices()
{
## Set matrices to zero
echo
## echo '## Set matrices to zero'
|
|
|
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
|
case ${rep} in
dm)
let Nx=Nx+2*Nz+Nyz # Matrices not usual size
;;
*)
esac
Npar=`wc -l $system\_sympar.txt | gawk '{print $1}'`
}
zero_matrices()
{
## Set matrices to zero
echo
## echo '## Set matrices to zero'
|
| ︙ | | | ︙ | |
605
606
607
608
609
610
611
612
613
614
615
616
617
618
|
;;
open)
N=$Nx; M=1
;;
ax)
N=$Nx; M=1
;;
*)
esac
if [ "${language}" = "oct" ]; then
if [ "$M" = "1" ]; then
echo " mtt$name = zeros($N);"
|
>
>
>
|
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
|
;;
open)
N=$Nx; M=1
;;
ax)
N=$Nx; M=1
;;
par)
N=$Npar; M=1
;;
*)
esac
if [ "${language}" = "oct" ]; then
if [ "$M" = "1" ]; then
echo " mtt$name = zeros($N);"
|
| ︙ | | | ︙ | |