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