Changeset 4447


Ignore:
Timestamp:
May 27, 2020 11:56:18 AM (16 months ago)
Author:
vondreele
Message:

change RMCProfie chi2 to log10 chi2 to match fulrmc chi2 plot
fix angle restraint generator for fullrmc - remove duplicates

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIphsGUI.py

    r4442 r4447  
    58775877                X = np.array(X)
    58785878                Chi = np.array(Chi).T
    5879                 XY = [[X.T,Y.T] for Y in Chi]
     5879                XY = [[X.T,np.log10(Y.T)] for Y in Chi]
    58805880                G2plt.PlotXY(G2frame,XY,labelX='no. generated',
    5881                     labelY=r'$\mathsf{\chi^2}$',newPlot=True,Title='RMCP Chi^2 for '+pName,
     5881                    labelY=r'$log_{10}$ (reduced $\mathsf{\chi^2})$',newPlot=True,Title='RMCP Chi^2 for '+pName,
    58825882                    lines=True,names=Names[3:])
    58835883
  • trunk/GSASIIpwd.py

    r4440 r4447  
    25742574                XYZ = np.array(G2mth.GetAtomItemsById(atomData,AtLookup,IdList,cx,numItems=3))
    25752575                calAngle = G2mth.getRestAngle(XYZ,Amat)
    2576                 if angData[0] < calAngle < angData[1]:
    2577                     AngleList.append(str((MidName,nbr1[0],nbr3[0])+tuple(angData))+',\n')
     2576                if angData[0] <= calAngle <= angData[1]:
     2577                    angStr = str((MidName,nbr1[0],nbr3[0])+tuple(angData))+',\n'
     2578                    revangStr = str((MidName,nbr3[0],nbr1[0])+tuple(angData))+',\n'
     2579                    if angStr not in AngleList and revangStr not in AngleList:
     2580                        AngleList.append(angStr)
    25782581    return AngleList
    25792582
Note: See TracChangeset for help on using the changeset viewer.