Sun 08 Jan 2017 03:19:11 PM UTC, comment #5:
"Are you trying to make Octave work with both IDA 2.8.1 and 2.9.0, or do you want to require 2.9.0?"
The intention was to require the latest version of IDA, the one that ships with sundials 2.7.
The data format and API for sparse matrices have changed between 2.6 and 2.7 and supporting both seemed like a lot of work.
"Which library defines SparseSetMatToZero? "
It is in the sparse matrix module of sundials 2.7
"Why exactly do you need to check for SparseSetMatToZero? What feature will be enabled/disabled if it is present/missing?"
I used that as a way to detect the API version of sundials_sparse
the same function had a different name and syntax in versions prior to 2.7
"Note that it is now possible to build without IDAKLU, but sparse Jacobian handling is simply disabled."
What do you mean by "disabled": will Octave error out if the Jacobian is sparse or just silently convert it to full?
As I mentioned in a mailing list thread, I would not like to silently ignore missing IDA_KLU and let Octave work with sparse jacobians as if they were full.
I have shown a very simple example where the difference in performance is a factor of 150, more meaningful examples have even more striking differences. I would like users in general and downstream packagers in particular to be well aware of this issue, so that they are a bit motivated to improve the sundials package for debian at least ...
|