11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
# Copyright (c) P.J.Gawthrop 1991, 1992, 1994.
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.26 2002/10/28 23:31:21 gawthrop
## Added additional term to MTTEdx to account for zdot terms on
## Right-Hand Side state equations
##
## 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.
|
>
>
>
|
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
|
# Copyright (c) P.J.Gawthrop 1991, 1992, 1994.
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.27 2005/03/21 11:50:39 gawthrop
## Don't write an empty cseo file
##
## Revision 1.26 2002/10/28 23:31:21 gawthrop
## Added additional term to MTTEdx to account for zdot terms on
## Right-Hand Side state equations
##
## 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.
|
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
|
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;
|
|
|
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
|
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 + MTTEyx*(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;
|