[Nek5000-users] Using CVODE with a vector of absolute tolerances

[nek5000-users at lists.mcs.anl.gov]

Dear all,

            I am trying to solve a catalytic combustion problem with nek5000 and I am using CVODE to integrate the temperature and species equations. I was doing some tests with the tolerances and I wanted to try different absolute tolerances for the equations.

In the nek folder I changed some parameters in CVODE:

- IATOL = 2;
- real cv_atol(15)

Then in cv_init I assign a value to each component of the vector cv_atol.

When fcvmalloc is called, I get the error message:

[CVODE ERROR] CVodeMalloc
abstol has a negative component(s) (illegal).

If I debug the code, before and after the call to fcvmalloc the vector cv_atol is fine.

Does anybody experienced something similar or has a suggestion about that?
In which folder is the source code of CVODE?

Thanks in advance for your answer.

Regards,
Andreas.