Using Sundials cvode solver

Hi All,

I'm currently using ode15s to solve a set of stiff differential equations.

[t, species]  = ode15s(@(t,s) factory(t,s), tspan , species0, options);

I would like to try solving the same set of odes using the CVODE solver in Sundials.

I checked some examples which show how to specify the solver type and use Sundials.

configsetObj = getconfigset(modelObj);
configsetObj.SolverType = 'sundials';
configsetObj.MaximumWallClock = 30;
configsetObj.MaximumNumberOfLogs = 100;
configsetObj.StopTime = 10;
[time, states] = sbiosimulate(modelObj);
plot(time,states)

I would like to know how modelObj should be defined.

[Deepa Maheshvare]

Hi All,

Any suggestions?

[Deepa Maheshvare]

Any suggestions?

[Walter Roberson]

No suggestions from me. I had no idea that facility existed until you asked your question.

[Deepa Maheshvare]

Thank you for the reponse.

I hope I could get some suggestions from the others.

[Star Strider]

I had to look that up. I found SUNDIALS Solvers as the only relevant documentation. I get the impression that these are stand-alone routines, however I have not explored the LLNL documentation for the software.

Answer 1

[Aditya]

Hi Deepa,

I understand that you want to solve set of stiff differential equations using the CVODE solver in Sundials and need help defining the “modelObj” for it.

To convert your ODE model to SimBiology, you can follow a direct approach by converting each differential equation into a rate rule. You can find an example of that here:

• https://www.mathworks.com/help/simbio/ug/rule-object.html#mw_fef86e0c-1589-4ff5-b866-0c8dfd230895

In addition, I recommend checking out the following resources to understand how to define a model object according to your specific needs.

• https://www.mathworks.com/help/simbio/ref/modelobject.html#bqu195b-7
• https://www.mathworks.com/help/simbio/gs/construct-a-simple-model.html

Hope this helps!

Answer 2

[Mike Croucher]

As of MATLAB R2024a, SUNDIALS is shipped with and supported directly in MATLAB. See my blog post for details

Faster Ordinary Differential Equations (ODEs) solvers and Sensitivity Analysis of Parameters: Introducing SUNDIALS support in MATLAB » The MATLAB Blog - MATLAB & Simulink (mathworks.com)