idaes-pse icon indicating copy to clipboard operation
idaes-pse copied to clipboard

Report Numerical Issues does not issue a warning for a singular Jacobian

Open dallan-keylogic opened this issue 1 year ago • 1 comments
trafficstars

When troubleshooting issues revealed by #1395, I encountered this:

====================================================================================
Model Statistics

    Jacobian Condition Number: 9.195E+14

------------------------------------------------------------------------------------
1 WARNINGS

    WARNING: 54 pairs of variables are parallel (to tolerance 1.0E-08)

------------------------------------------------------------------------------------
6 Cautions

    Caution: 72 Variables with value close to zero (tol=1.0E-08)
    Caution: 84 Variables with extreme value (<1.0E-04 or >1.0E+04)
    Caution: 8 Variables with None value
    Caution: 74 Variables with extreme Jacobian values (<1.0E-04 or >1.0E+04)
    Caution: 54 Constraints with extreme Jacobian values (<1.0E-04 or >1.0E+04)
    Caution: 76 extreme Jacobian Entries (<1.0E-04 or >1.0E+04)

------------------------------------------------------------------------------------
Suggested next steps:

    display_near_parallel_variables()

====================================================================================

This matrix is more-or-less singular, which is why it's shocking that we're not issuing a warning for it.

dallan-keylogic avatar Jun 15 '24 21:06 dallan-keylogic

This requires #1429 and will also require many models to improve their scaling. A working session with all the owners of the impacted models is being organized.

ksbeattie avatar Oct 03 '24 18:10 ksbeattie