Open mickyaero opened 6 years ago
The dimension should be generated from VoxelPools
object. Each VoxelPool has its won system of ODE ( BoostSys
); there you need to figure out how to extract the dimension of dydt
, dy
vector.
See the file ksolve/VoxelPool.cpp
https://github.com/BhallaLab/moose-core/blob/master/ksolve/VoxelPools.cpp .
In short, the problem is how to get dimension of ODE from VoxelPools
?
To start with
VoxelPools
process function. How to compute the size of matrix ? Number of species you can get by reading the size of vector of species.I am writing this from memory. Few names may be different.
The seg fault is most probably due to the wrong dimension
f[i] is being accessed but variable i is of type size_t and the for loop is iterating over dimension and the size of the dimension is huge (6729915426273446717) which is beyond the limit of size_t and hence seg fault. Therefore the dimension needs to be corrected and that will resolve the issue hopefully.