Changeset 3149

Timestamp:

Nov 3, 2017 9:39:05 AM (6 years ago)

Author:

vondreele

Message:

Fix Jacobian refine bug

File:

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

trunk/GSASIIstrMath.py

@@ -4059 +4059 @@
     Histograms,Phases,restraintDict,rigidbodyDict = HistoPhases
     dependentVars = G2mv.GetDependentVars()
-    dMdv = np.empty(0)
     histoList = list(Histograms.keys())
     histoList.sort()
+    First = True
     for histogram in histoList:
         if 'PWDR' in histogram[:4]:
@@ -4087 +4087 @@
         else:
             continue        #skip non-histogram entries
-        if len(dMdv):
+        if First:
+            dMdv = np.sqrt(wtFactor)*dMdvh
+            First = False
+        else:
             dMdv = np.concatenate((dMdv.T,np.sqrt(wtFactor)*dMdvh.T)).T
-        else:
-            dMdv = np.sqrt(wtFactor)*dMdvh
             
     GetFobsSq(Histograms,Phases,parmDict,calcControls)