Changeset 1092 for trunk/GSASIImath.py

Timestamp:

Oct 9, 2013 10:27:52 AM (10 years ago)

Author:

vondreele

Message:

remove GSASIIsolve - never will be used
put back matrix rescaling for v-cov matrix
add detail for penalty fxn output - still something wrong with this

File:

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

trunk/GSASIImath.py

@@ -148 +148 @@
     nfev += 1
     Yvec,Amat = Hess(x0,*args)
-    Amatlam = Amat          #*(One+Lam)/Anorm
+    Amatlam = Amat*(One+Lam)/Anorm              #scale Amat to Marquardt array
     try:
-        Bmat = nl.inv(Amatlam)          #*(One+Lam)/Anorm
+        Bmat = nl.inv(Amatlam)*(One+Lam)/Anorm      #rescale Bmat to Marquardt array
         return [x0,Bmat,{'num cyc':icycle,'fvec':M,'nfev':nfev,'lamMax':lamMax,'psing':[]}]
     except nl.LinAlgError: