Ignore:
Timestamp:
Oct 10, 2013 10:35:50 AM (12 years ago)
Author:
vondreele
Message:

mods to restraint GUI & LS output for restraints
correct torsion & Ramachandran penalties

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TabularUnified trunk/GSASIIstrMath.py

    r1092 r1097  
    356356                            phi,psi = G2mth.getRestRama(XYZ,Amat)
    357357                            restr,calc = G2mth.calcRamaEnergy(phi,psi,coeffDict[cofName])                               
    358                         pVals.append(obs-calc)
     358                        pVals.append(restr)
    359359                        pWt.append(wt/esd**2)
    360                         pWsum[name] += wt*((obs-calc)/esd)**2
     360                        pWsum[name] += wt*(restr/esd)**2
    361361                elif name == 'ChemComp':
    362362                    for i,[indx,factors,obs,esd] in enumerate(itemRest[rest]):
     
    20032003        pSum = np.sum(pWt*pVals**2)
    20042004        for name in pWsum:
    2005             print '  Penalty function for %s = %.3f'%(name,pWsum[name])
     2005            print '  Penalty function for %8s = %.3f'%(name,pWsum[name])
    20062006        print 'Total penalty function: %.3f on %d terms'%(pSum,len(pVals))
    20072007        Nobs += len(pVals)
Note: See TracChangeset for help on using the changeset viewer.