Changeset 2815


Ignore:
Timestamp:
Apr 28, 2017 10:38:40 AM (5 years ago)
Author:
vondreele
Message:

sign changes & fix error in texture & Pawley restraints

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrMath.py

    r2808 r2815  
    378378                        ODFln = G2lat.Flnh(False,SHCoef,phi,beta,SGData)
    379379                        R,P,Z = G2mth.getRestPolefig(ODFln,SamSym[textureData['Model']],grid)
    380                         Z1 = -ma.masked_greater(Z,0.0)
     380                        Z1 = ma.masked_greater(Z,0.0)           #is this + or -?
    381381                        IndZ1 = np.array(ma.nonzero(Z1))
    382382                        for ind in IndZ1.T:
     
    389389                            for ind in np.ndindex(grid,grid):
    390390                                pNames.append('%d:%s:%d:%.2f:%.2f'%(pId,name+'-unit',i,R[ind[0],ind[1]],P[ind[0],ind[1]]))
    391                                 pVals.append(Z1[ind[0]][ind[1]])
     391                                pVals.append(Z2[ind[0]][ind[1]])
    392392                                pWt.append(wt/esd2**2)
    393393                                pWsum[name] += wt*(Z2/esd2)**2
     
    420420                        ODFln = G2lat.Flnh(False,SH3Coef,phi,beta,SGData)
    421421                        X = np.linspace(0,90.0,26)
    422                         Y = -ma.masked_greater(G2lat.polfcal(ODFln,'0',X,0.0),0.0)
     422                        Y = ma.masked_greater(G2lat.polfcal(ODFln,'0',X,0.0),0.0)       #+ or -?
    423423                        IndY = ma.nonzero(Y)
    424424                        for ind in IndY[0]:
     
    433433            if negWt[pId]:
    434434                pNames.append(item)
    435                 pVals.append(-parmDict[item])
     435                pVals.append(parmDict[item])
    436436                pWt.append(negWt[pId])
    437                 pWsum['PWLref'] += negWt[pId]*(-parmDict[item])**2
     437                pWsum['PWLref'] += negWt[pId]*(parmDict[item])**2
    438438    pVals = np.array(pVals)
    439439    pWt = np.array(pWt)         #should this be np.sqrt?
Note: See TracChangeset for help on using the changeset viewer.