-
Notifications
You must be signed in to change notification settings - Fork 60
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
ludc #1700
Comments
@myravian I've checked the code, the INTERCHANGES keyword is silently ignored in GDL. |
mat is singular, it's determinant is exactly zero. Probably an error message should be printed, as with LA_LUDC in IDL. |
@fawltylanguage IDL does not complain about this mat. But yes, another issue, LUDC should return a warning when encountering a singular matrix. |
The matrix values were arbitrary and there just to illustrate the issue of the keyword that I wanted to raise. It's easy enough to check the determinant beforehand. |
I just corrected the bug in the trunk.
GDL should probably signal the singular matrix as IDL does, even if the unity matrix (wrongly?) passes for IDL. After LU decomposition, it would suffice to check the diagonal for any zero, agreed? |
For my records : |
Hi there,
While trying GDL 1.0.4 on a script, I found this issue:
GDL> mat = [[1,1,1],[1,1,1],[1,1,1]]
GDL> ludc, mat, index, interchanges=parity
GDL> help, parity
PARITY UNDEFINED =
expected behavior:
PARITY FLOAT = 1.00000
Does it have to do with the way the array is being declared?
Cheers,
Vian
The text was updated successfully, but these errors were encountered: