Changeset 3603

Timestamp:

Sep 15, 2018 7:20:55 PM (5 years ago)

Author:

vondreele

Message:

small mods to SVD stuff

File:

• trunk/GSASIImath.py (modified) (2 diffs)

trunk/GSASIImath.py

@@ -91 +91 @@
 
     """
-    u, s, vt = nl.svd(a, 0)
+    u, s, vt = nl.svd(a)
     cutoff = rcond*np.maximum.reduce(s)
     s = np.where(s>cutoff,1./s,0.)
@@ -183 +183 @@
                 Ainv,Nzeros = pinv(Amatlam,xtol)    #do Moore-Penrose inversion (via SVD)
             except nl.LinAlgError:
-                print ('ouch #1 bad SVD inversion; change parameterization')
                 psing = list(np.where(np.diag(nl.qr(Amatlam)[1]) < 1.e-14)[0])
+                print ('ouch #1 bad SVD inversion; change parameterization for ',psing)
                 return [x0,None,{'num cyc':icycle,'fvec':M,'nfev':nfev,'lamMax':lamMax,'psing':psing,'SVD0':-1}]
             Xvec = np.inner(Ainv,Yvec)      #solve