search

cmillion / gPhoton

The GALEX photon database and tools project.
Other
21 stars 11 forks source link

e38150 crashes gPipeline #255

Closed cmillion closed 7 years ago

cmillion commented 7 years ago

Data available from: http://galex.stsci.edu/gPhoton/RAW6/e38100_38199/38150/e38150-nd-raw6.fits.gz http://galex.stsci.edu/gPhoton/RAW6/e38100_38199/38150/e38150-scst.fits.gz

Reproduce w/:

import gPhoton
import gPhoton.gPipeline
gPhoton.PhotonPipe.photonpipe('e38150-nd-raw6.fits','e38150-scst.fits','NUV','test-nd-x',verbose=2)

Crashes w/:

linalg.py in lstsq(a, b, rcond) 1917 work = zeros((lwork,), t) 1918 results = lapack_routine(m, n, n_rhs, a, m, bstar, ldb, s, rcond, -> 1919 0, work, lwork, iwork, 0) 1920 if results['info'] > 0: 1921 raise LinAlgError('SVD did not converge in Linear Least Squares')

ValueError: On entry to DGELSD parameter number 6 had an illegal value

cmillion commented 7 years ago

The issues seems to be that when hunting for stim positions to calculate stim distances for the distortion correction, two of the four stims return zero hits. So these propagate as nan values which then make the linear trend fitting explode.