Changeset 3068 for trunk/GSASIIstrMain.py

Timestamp:

Sep 12, 2017 10:54:45 AM (4 years ago)

Author:

vondreele

Message:

fix parameter display to show nan (crashed before)
trap with error if any parameter is nan
skip postprocessing for failed refinements (IfOK=False)
fix Jacobian to use GetProfileDervMP

File:

• trunk/GSASIIstrMain.py (modified) (4 diffs)

trunk/GSASIIstrMain.py

@@ -54 +54 @@
         begin = time.time()
         values =  np.array(G2stMth.Dict2Values(parmDict, varyList))
+        if np.any(np.isnan(values)):
+            raise G2obj.G2Exception('ERROR - nan found in LS parameters - use Calculate/View LS parms to locate')
         # test code to compute GOF and save for external repeat
         #args = ([Histograms,Phases,restraintDict,rigidbodyDict],parmDict,varyList,calcControls,pawleyLookup,dlg)
@@ -135 +137 @@
                         del(varyList[ipvt-1])
                         break
-    G2stMth.GetFobsSq(Histograms,Phases,parmDict,calcControls)
+    if IfOK:
+        G2stMth.GetFobsSq(Histograms,Phases,parmDict,calcControls)
     return IfOK,Rvals,result,covMatrix,sig
 
@@ -212 +215 @@
         IfOK,Rvals,result,covMatrix,sig = RefineCore(Controls,Histograms,Phases,restraintDict,
             rigidbodyDict,parmDict,varyList,calcControls,pawleyLookup,ifPrint,printFile,dlg)
-        sigDict = dict(zip(varyList,sig))
-        newCellDict = G2stMth.GetNewCellParms(parmDict,varyList)
-        newAtomDict = G2stMth.ApplyXYZshifts(parmDict,varyList)
-        covData = {'variables':result[0],'varyList':varyList,'sig':sig,'Rvals':Rvals,
-                   'varyListStart':varyListStart,
-                   'covMatrix':covMatrix,'title':GPXfile,'newAtomDict':newAtomDict,
-                   'newCellDict':newCellDict,'freshCOV':True}
-        # add the uncertainties into the esd dictionary (sigDict)
-        sigDict.update(G2mv.ComputeDepESD(covMatrix,varyList,parmDict))
-        G2mv.PrintIndependentVars(parmDict,varyList,sigDict,pFile=printFile)
-        G2stMth.ApplyRBModels(parmDict,Phases,rigidbodyDict,True)
-        G2stIO.SetRigidBodyModels(parmDict,sigDict,rigidbodyDict,printFile)
-        G2stIO.SetPhaseData(parmDict,sigDict,Phases,rbIds,covData,restraintDict,printFile)
-        G2stIO.SetHistogramPhaseData(parmDict,sigDict,Phases,Histograms,calcControls['FFtables'],pFile=printFile)
-        G2stIO.SetHistogramData(parmDict,sigDict,Histograms,calcControls['FFtables'],pFile=printFile)
-        G2stIO.SetUsedHistogramsAndPhases(GPXfile,Histograms,Phases,rigidbodyDict,covData,makeBack)
-        printFile.close()
-        print ' Refinement results are in file: '+ospath.splitext(GPXfile)[0]+'.lst'
-        print ' ***** Refinement successful *****'
+        if IfOK:
+            sigDict = dict(zip(varyList,sig))
+            newCellDict = G2stMth.GetNewCellParms(parmDict,varyList)
+            newAtomDict = G2stMth.ApplyXYZshifts(parmDict,varyList)
+            covData = {'variables':result[0],'varyList':varyList,'sig':sig,'Rvals':Rvals,
+                       'varyListStart':varyListStart,
+                       'covMatrix':covMatrix,'title':GPXfile,'newAtomDict':newAtomDict,
+                       'newCellDict':newCellDict,'freshCOV':True}
+            # add the uncertainties into the esd dictionary (sigDict)
+            sigDict.update(G2mv.ComputeDepESD(covMatrix,varyList,parmDict))
+            G2mv.PrintIndependentVars(parmDict,varyList,sigDict,pFile=printFile)
+            G2stMth.ApplyRBModels(parmDict,Phases,rigidbodyDict,True)
+            G2stIO.SetRigidBodyModels(parmDict,sigDict,rigidbodyDict,printFile)
+            G2stIO.SetPhaseData(parmDict,sigDict,Phases,rbIds,covData,restraintDict,printFile)
+            G2stIO.SetHistogramPhaseData(parmDict,sigDict,Phases,Histograms,calcControls['FFtables'],pFile=printFile)
+            G2stIO.SetHistogramData(parmDict,sigDict,Histograms,calcControls['FFtables'],pFile=printFile)
+            G2stIO.SetUsedHistogramsAndPhases(GPXfile,Histograms,Phases,rigidbodyDict,covData,makeBack)
+            printFile.close()
+            print ' Refinement results are in file: '+ospath.splitext(GPXfile)[0]+'.lst'
+            print ' ***** Refinement successful *****'
+        else:
+            print '****ERROR - Refinement failed'
+            raise G2obj.G2Exception('****ERROR - Refinement failed')
     except G2obj.G2Exception,Msg:
         printFile.close()
@@ -236 +243 @@
 
 #for testing purposes!!!
-    if DEBUG:
+    if DEBUG and IfOK:
 #needs: values,HistoPhases,parmDict,varylist,calcControls,pawleyLookup
         import cPickle