I am solving a large system of ODEs with the DifferentialEquation.jl package. Since I am using a number of for loops to set up the RHS to the du, so I would like to see the RHS matrix before solving. Is there a way I can do that?
Lets say I have the following
for i in 1:n
for j in 1:n
prob = ODEProblem(flux!, u0, tspan)
sol = solve(prob)
So basically I am trying to put a convolution like function on the right hand side, and trying to see if the correct ODE system is generated or not.
Thanks for all the help, the issue was when I was writing the flux as a sum on a new line for each component, only the 1st component was being added somehow. So when I changed the coefficient for the 1st component to zero to see how the others affected the solution, the RHS essentially became zero, as the rest of the components were not being added. That’s why I wanted to look at the RHS.
It took me a while to realize that whatever I sum to the flux after putting a newline is not being added to the RHS.