Changeset 4903 for trunk/GSASIImath.py

Timestamp:

May 17, 2021 11:23:03 AM (3 years ago)

Author:

vondreele

Message:

fix 2 issues in HessianLSQ for SVN LinAlgErrors?
fix a typo in G2plot for SASD plot display
change GetTthAzmDsp? --> GeTthAzmDsp2 in G2plot & G2Image
comment old GetTthAzmDsp? code - fails in recent python/numpy versions; replaced by GetTthAzmDsp2
change GetTthAzmG --> GetTthAzm2 in G2image

File:

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

trunk/GSASIImath.py

@@ -263 +263 @@
             Amatlam = Amat*(1.+np.eye(Amat.shape[0])*lam)
             try:
+                Nzeros = 1
                 Ainv,Nzeros = pinv(Amatlam,xtol)    # Moore-Penrose SVD inversion
             except nl.LinAlgError:
@@ -412 +413 @@
         psing = list(np.where(np.abs(np.diag(nl.qr(Amat)[1])) < 1.e-14)[0])
         if not len(psing): # make sure at least the worst term is flagged
+            d = np.abs(np.diag(nl.qr(Amat)[1]))        
             psing = [np.argmin(d)]
         Amat, indices, Yvec = dropTerms(psing, Amat, indices, Yvec)