Le 24/04/2015 14:34, Stéphane Mottelet a écrit :
Hello,
this is not trivial indexing, in fact some terms are linear
combination of v's components
M1_v=[v(17)
v(104)
v(149)
-(v(18)+v(63)+v(103))
-(v(18)+v(63)+v(103))
v(17)
...
v(104)
v(149)
]
How do you take this into account in your proposed method ? These
combinations are sums of influxes in a metabolic network, and the code
is automatically generated.
In this way: from
1> M1_v=[v(17)
2> v(104)
3> v(149)
4> -(v(18)+v(63)+v(103))
5> -(v(18)+v(63)+v(103))
6> v(17)
...
838> v(104)
839> v(149)
]
set:
i = [1 2 3 6 ... 838 839]
iv = [17 104 149 17 ... 104 149]
M1_v(i) = v(iv)
im = [4 5 ] // or setdiff(1:839, i).
ivm = [18 18 ... ]
M1_v(im) = -v(ivm)
// + loop over the rank of elements in the linear combination
// + adapt according to the automatic generation of your code
Samuel
_______________________________________________
users mailing list
[email protected]
http://lists.scilab.org/mailman/listinfo/users
