@@ -3502,7 +3502,6 @@ def GetReflPosDerv(refl,im,wave,A,pfx,hfx,phfx,calcControls,parmDict):
35023502def GetHStrainShift (refl ,im ,SGData ,phfx ,hfx ,calcControls ,parmDict ):
35033503 '''Computes the shifts in peak position due to the Hydrostatic strain
35043504 (HStrain, Dij terms).
3505- This routine is not used anywhere
35063505 '''
35073506 laue = SGData ['SGLaue' ]
35083507 uniq = SGData ['SGUniq' ]
@@ -4034,7 +4033,7 @@ def SavePartial(phase,y):
40344033 Uniq = np .inner (refl [:3 ],SGMT )
40354034 refl [5 + im ] = GetReflPos (refl ,im ,0.0 ,A ,pfx ,hfx ,phfx ,calcControls ,parmDict ) #corrected reflection position - #TODO - what about tabluated offset?
40364035 Lorenz = sind (abs (parmDict [hfx + '2-theta' ])/ 2 )* refl [4 + im ]** 4 #TOF Lorentz correction
4037- # refl[5+im] += GetHStrainShift(refl,im,SGData,phfx,hfx,calcControls,parmDict) #apply hydrostatic strain shift
4036+ refl [5 + im ] += GetHStrainShift (refl ,im ,SGData ,phfx ,hfx ,calcControls ,parmDict ) #apply hydrostatic strain shift
40384037 refl [6 + im :8 + im ] = GetReflSigGamTOF (refl ,im ,G ,GB ,phfx ,calcControls ,parmDict ) #peak sig & gam
40394038 refl [12 + im :14 + im ] = GetReflAlpBet (refl ,im ,hfx ,parmDict ) #TODO - skip if alp, bet tabulated?
40404039 refl [11 + im ],refl [15 + im ],refl [16 + im ],refl [17 + im ] = GetIntensityCorr (refl ,im ,Uniq ,G ,g ,pfx ,phfx ,hfx ,SGData ,calcControls ,parmDict )
0 commit comments