[ADOL-C] sparse jacobian: inf's and nan's in non-zeros values array

[Antoine De Blois]

Hi everyone,

I created a ticket on the trac system, but I don't know which is the preferred way (trac or mailing list) please advise. Here is my problem:


I am trying to solve the jacobian of a residual function in my structured  CFD solver. In my application, I have 5 dependent variables and roughly

1700 independent variables. The number of non-zero values in the jacobian  matrix is roughly 200, hence the usage of the sparse_jac driver.



The jacobian driver works perfectly fine, but I obtain inf's and nan's in  the non-zeros values array. Note that the row and column index arrays are  fine when compared to the entire matrix output of the jacobian driver. I  tried changing the option array without success.



I was using version 2.3.0 and recently upgraded to 4.3.1 with the same  behavior. I am using colpack version 1.0.9.



Can you please help me with my problem?

Thank you for your time,

Antoine
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://list.coin-or.org/pipermail/adol-c/attachments/20130322/4dbf5807/attachment.html>