[bug #64303] Statistics fit --fitestimate crash when covariance matrix has NaN

[Mohammad Akhlaghi]

URL:
  <https://savannah.gnu.org/bugs/?64303>

                 Summary: Statistics fit --fitestimate crash when covariance
matrix has NaN
                   Group: GNU Astronomy Utilities
               Submitter: makhlaghi
               Submitted: Thu 15 Jun 2023 11:24:27 AM BST
                Category: Statistics
                Severity: 3 - Normal
              Item Group: Crash
                  Status: In Progress
                 Privacy: Public
             Assigned to: makhlaghi
             Open/Closed: Open
         Discussion Lock: Any


    _______________________________________________________

Follow-up Comments:


-------------------------------------------------------
Date: Thu 15 Jun 2023 11:24:27 AM BST By: Mohammad Akhlaghi <makhlaghi>
When the covariance matrix of the fit contains a NaN value, Statistics will
abort with a CFITSIO error. This can be reproduced with the simple table and
command below:


$ cat table.txt
1  3
2  4


$ aststatistics table.txt -c1,2 --fit=polynomial --fitmaxpower=1
--fitestimate=self -oout.fits

Statistics (GNU Astronomy Utilities) 0.20
-------
Fitting results (remove extra info with '--quiet' or '-q)
  Input file:    table.txt with 2 non-blank rows.
  X column: 1
  Y column: 2

Fit function: Y = c0 + (c1 * X^1) + (c2 * X^2) + ... (cN * X^N)
  N:  1
  c0:  +2.0000000000
  c1:  +1.0000000000

Covariance matrix:
  -nan                 -nan                 
  -nan                 -nan                 

Reduced chi^2 of fit:
  -nan

Requested estimation:
  Written to: out.fits
aststatistics: gal_fits_key_write_in_ptr_nan_check: (WARNING) value of
'FRDCHISQ' is NaN and FITS doesn't recognize a NaN key value

FITSIO status = 402: bad float to string conversion
Error reading data buffer from file:
table.txt
ffopen could not interpret primary array header of file: 
table.txt
Error reading data buffer from file:
table.txt
ffopen could not interpret primary array header of file: 
table.txt
Error reading data buffer from file:
table.txt
ffopen could not interpret primary array header of file: 
table.txt
Error in ffd2e: double value is a NaN or INDEF
aststatistics: Error in CFITSIO, see above.


This bug was reported by Sepideh Eskandarlou.







    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?64303>

_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/