@@ -8,8 +8,8 @@ eqs = [D(x) ~ a*x - x*y,
8
8
"""
9
9
real[] diffeqf(real t,real[] internal_var___u,real[] internal_var___p,real[] x_r,int[] x_i) {
10
10
real internal_var___du[2];
11
- internal_var___du[1] = internal_var___p[1] * internal_var___u[1] - internal_var___u[1] * internal_var___u[2];
12
- internal_var___du[2] = -3 * internal_var___u[2] + internal_var___u[1] * internal_var___u[2];
11
+ internal_var___du[1] = internal_var___p[1] * internal_var___u[1] + -1 * internal_var___u[1] * internal_var___u[2];
12
+ internal_var___du[2] = internal_var___u[1] * internal_var___u[2] + -3 * internal_var___u[2];
13
13
return internal_var___du;
14
14
}
15
15
"""
@@ -19,14 +19,14 @@ eqs = [D(x) ~ a*x - x*y,
19
19
rhsnames= [:internal_var___u ,:internal_var___p ,:t ]) ==
20
20
"""
21
21
void diffeqf(double* internal_var___du, double* internal_var___u, double* internal_var___p, double t) {
22
- internal_var___du[0] = internal_var___p[0] * internal_var___u[0] - internal_var___u[0] * internal_var___u[1];
23
- internal_var___du[1] = -3 * internal_var___u[1] + internal_var___u[0] * internal_var___u[1];
22
+ internal_var___du[0] = internal_var___p[0] * internal_var___u[0] + -1 * internal_var___u[0] * internal_var___u[1];
23
+ internal_var___du[1] = internal_var___u[0] * internal_var___u[1] + -3 * internal_var___u[1];
24
24
}
25
25
"""
26
26
27
27
@test ModelingToolkit. build_function (eqs,[x,y],[a],t,target = ModelingToolkit. MATLABTarget ()) ==
28
28
"""
29
- diffeqf = @(t,internal_var___u) [internal_var___p(1) * internal_var___u(1) - internal_var___u(1) * internal_var___u(2); -3 * internal_var___u(2) + internal_var___u(1) * internal_var___u(2)];"""
29
+ diffeqf = @(t,internal_var___u) [internal_var___p(1) * internal_var___u(1) + -1 * internal_var___u(1) * internal_var___u(2); internal_var___u(1) * internal_var___u(2) + -3 * internal_var___u(2)];"""
30
30
31
31
sys = ODESystem (eqs,t,[x,y],[a])
32
32
0 commit comments