Changeset 1676 for trunk/GSASIImapvars.py
- Timestamp:
- Feb 27, 2015 9:46:37 PM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/GSASIImapvars.py
r1502 r1676 570 570 if notvaried: notvaried += ', ' 571 571 notvaried += mv 572 if mv not in indepVarList: indepVarList.append(mv)573 572 if parmDict is not None and mv not in parmDict: 574 msg += "\nCannot equivalence to variable "+str(mv)+". Not defined in this refinement" 575 continue 573 print "Dropping equivalence for variable "+str(mv)+". Not defined in this refinement" 574 #msg += "\nCannot equivalence to variable "+str(mv)+". Not defined in this refinement" 575 #continue 576 else: 577 if mv not in indepVarList: indepVarList.append(mv) 576 578 for v,m in zip(varlist,invmultarr): 577 579 if parmDict is not None and v not in parmDict: … … 873 875 for varlist,mapvars,invmultarr in zip(dependentParmList,indParmList,invarrayList): 874 876 #if invmultarr is None: continue # probably not needed 875 valuelist = [parmDict[var] for var in mapvars] 877 try: 878 valuelist = [parmDict[var] for var in mapvars] 879 except KeyError: 880 continue 876 881 # get the v-covar matrix for independent parameters 877 882 vcov = np.zeros((len(mapvars),len(mapvars))) … … 1066 1071 for varlist,mapvars,invmultarr in zip(dependentParmList,indParmList,invarrayList): 1067 1072 #if invmultarr is None: continue 1068 valuelist = [parmDict[var] for var in mapvars] 1073 try: 1074 valuelist = [parmDict[var] for var in mapvars] 1075 except KeyError: 1076 continue 1069 1077 parmDict.update(zip(varlist, 1070 1078 np.dot(invmultarr,np.array(valuelist)))
Note: See TracChangeset
for help on using the changeset viewer.