11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
+
+
+
+
+
|
# Copyright (c) P.J.Gawthrop 1991, 1992, 1994.
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.25 2002/09/18 19:46:48 geraint
## Refined test to reduce shell loop overhead for unnecessary E matrix.
## Removed unnecessary attempt to optimise cse.
## Fixes sm representation.
##
## Revision 1.24 2002/09/16 08:08:00 geraint
## Merged changes from global-optimisation branch.
##
## Revision 1.23.2.3 2002/09/12 18:50:50 geraint
## Uncommented cse optimisations - they seem to work ok.
##
## Revision 1.23.2.2 2002/09/10 23:24:19 geraint
|
300
301
302
303
304
305
306
307
308
309
310
311
312
313
|
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
|
+
+
|
IF MTTNx>0 THEN
BEGIN
MTTEdX := MTTdX; %Ie MTTEdX is MTTdX with the dz terms deleted ie EdX.
MTTdX := MTTdXs; %Restore the symbolic dX
%% Add on input derivative terms
MTTEdX := MTTEdX + MTTExu*MTTdu;
%% Add on output derivative terms
MTTEdx := MTTEdX + MTTExx*(MTTE^(-1))*MTTEdX;
END;
%%%%%MTTY := MTTY + MTTEyx*MTTEdX;
%%% This causes the matrix mismatch
%%% MTTdXs and MTTdu need setting in _def.r file
MTTY := MTTY + MTTEyu*MTTdu;
|