def rotate_around_rotax_to_get_vecs_in_plane(self,v0,v1): v0 = cross(np.array([0,0,1]),v0) v1 = cross(np.array([0,0,1]),v1) #print rotmat3Dfrom2vectors(v0,v1) self.crystalRotMat = rotmat3Dfrom2vectors(v0,v1)*self.crystalRotMat
def rotate_vec_parralel_to_rotax(self,vec): self.crystalRotMat = rotmat3Dfrom2vectors(vec,np.array([0,0,1]))