Overview
Comment:Fixed bug to do with multiple outputs
Downloads: Tarball | ZIP archive
Timelines: family | ancestors | descendants | both | origin/master | trunk
Files: files | file ages | folders
SHA3-256: a230d6637dc194895b57dc32539957878d7d7d323ee003f512e6e8cc5fb2ed98
User & Date: gawthrop@users.sourceforge.net on 2001-04-11 15:34:56.000
Other Links: branch diff | manifest | tags
Context
2001-04-11
15:41:45
New identification example - friction components only check-in: 6e8c859b73 user: gawthrop@users.sourceforge.net tags: origin/master, trunk
15:34:56
Fixed bug to do with multiple outputs check-in: a230d6637d user: gawthrop@users.sourceforge.net tags: origin/master, trunk
09:44:26
Fixed cc and c problems to do with pow(x,y) and integers
mtt/lib/reduce/fix_c.r is included in rdae2dae and cse2smx_lang for
-c, -cc and -oct options
check-in: 47742bde53 user: gawthrop@users.sourceforge.net tags: origin/master, trunk
Changes
18
19
20
21
22
23
24






25
26
27
28
29
30
31
32
33
34
35











36

37
38
39
40
41
42
43
18
19
20
21
22
23
24
25
26
27
28
29
30











31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50







+
+
+
+
+
+
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+

+








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

  if nargin<5			# No index given
    MTT_input_index = 0;
    [mtt_data] = ${sys}_ode2odes(x0,par,simpar);
    y = mtt_data(:,2:2:1+2*ny);
    ypar = [];
  else				# Compute sensitivities as well
  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
    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:2:1+2*ny);
      endif
      y_par = [y_par, mtt_data(:,3:2:2+2*ny)];
    endfor
  endif

endfunction
EOF

}

case ${lang} in
    m)
	make_m

MTT: Model Transformation Tools
GitHub | SourceHut | Sourceforge | Fossil RSS ]