Changeset 4688 for trunk

Timestamp:

Dec 30, 2020 2:11:42 PM (3 years ago)

Author:

toby

Message:

correct setting of singular vars in old HessianLSQ

File:

• trunk/GSASIImath.py (modified) (1 diff)

trunk/GSASIImath.py

@@ -176 +176 @@
         Yvec,Amat = Hess(x0,*args)
         Adiag = np.sqrt(np.diag(Amat))
-        psing = np.where(np.abs(Adiag) < 1.e-14,True,False)
+        #psing = np.where(np.abs(Adiag) < 1.e-14,True,False)
+        psing = np.where(np.abs(Adiag) < 1.e-14)[0]
         if np.any(psing):                #hard singularity in matrix
             return [x0,None,{'num cyc':icycle,'fvec':M,'nfev':nfev,'lamMax':lamMax,'psing':psing,'SVD0':-1}]