Changeset 1097 for trunk/GSASIIstrMath.py
- Timestamp:
- Oct 10, 2013 10:35:50 AM (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified trunk/GSASIIstrMath.py ¶
r1092 r1097 356 356 phi,psi = G2mth.getRestRama(XYZ,Amat) 357 357 restr,calc = G2mth.calcRamaEnergy(phi,psi,coeffDict[cofName]) 358 pVals.append( obs-calc)358 pVals.append(restr) 359 359 pWt.append(wt/esd**2) 360 pWsum[name] += wt*( (obs-calc)/esd)**2360 pWsum[name] += wt*(restr/esd)**2 361 361 elif name == 'ChemComp': 362 362 for i,[indx,factors,obs,esd] in enumerate(itemRest[rest]): … … 2003 2003 pSum = np.sum(pWt*pVals**2) 2004 2004 for name in pWsum: 2005 print ' Penalty function for % s = %.3f'%(name,pWsum[name])2005 print ' Penalty function for %8s = %.3f'%(name,pWsum[name]) 2006 2006 print 'Total penalty function: %.3f on %d terms'%(pSum,len(pVals)) 2007 2007 Nobs += len(pVals)
Note: See TracChangeset
for help on using the changeset viewer.