def makeRandGrains(numGrains, omegaMM, detectorGeom):
spotsDummy = Spots(planeData, None, detectorGeom, omegaMM)
if numGrains == 1:
rMats = num.atleast_3d(num.eye(3)).T
else:
quats = num.atleast_2d(ors.Quat.getRandQuat(n=numGrains)).T
rMats = ors.quatToMat(quats)
gList = []
for rMat in rMats:
g = grain.Grain(spotsDummy, rMat=rMat)
gList.append(g)
return gList
def makeRandGrains(numGrains, omegaMM, detectorGeom):
spotsDummy = Spots(planeData, None, detectorGeom, omegaMM)
if numGrains == 1:
rMats = num.atleast_3d(num.eye(3)).T
else:
quats = num.atleast_2d(ors.Quat.getRandQuat(n=numGrains)).T
rMats = ors.quatToMat(quats)
gList = []
for rMat in rMats:
g = grain.Grain(spotsDummy, rMat=rMat)
gList.append(g)
return gList
def makeRandGrains(numGrains, omegaMM, detectorGeom):
spotsDummy = Spots(planeData, None, detectorGeom, omegaMM)
if numGrains == 1:
import orientations as ors
r = ors.RotInv(3.08805198 - num.pi, 0., 1., 0.).toMatrix()
# r = num.eye(3)
rMats = num.atleast_3d(r).T
else:
quats = num.atleast_2d(ors.Quat.getRandQuat(n=numGrains)).T
rMats = ors.quatToMat(quats)
gList = []
for rMat in rMats:
g = grain.Grain(spotsDummy, rMat=rMat)
gList.append(g)
return gList
def makeRandGrains(numGrains, omegaMM, detectorGeom):
spotsDummy = Spots(planeData, None, detectorGeom, omegaMM)
if numGrains == 1:
#import orientations as ors
r = ors.RotInv(3.08805198-num.pi, 0., 1., 0.).toMatrix()
# r = num.eye(3)
rMats = num.atleast_3d(r).T
else:
quats = num.atleast_2d(ors.Quat.getRandQuat(n=numGrains)).T
rMats = ors.quatToMat(quats)
gList = []
for rMat in rMats:
g = grain.Grain(spotsDummy, rMat=rMat)
gList.append(g)
return gList
