From 91cf7ebbcc84c3f4f846461f4706842e825be2b9 Mon Sep 17 00:00:00 2001 From: Marco Cammarata Date: Fri, 13 Jan 2017 14:50:43 +0100 Subject: [PATCH] forgot to delete old file... --- xray/dspace.py | 70 -------------------------------------------------- 1 file changed, 70 deletions(-) delete mode 100644 xray/dspace.py diff --git a/xray/dspace.py b/xray/dspace.py deleted file mode 100644 index 40591fd..0000000 --- a/xray/dspace.py +++ /dev/null @@ -1,70 +0,0 @@ -from __future__ import division,print_function -import numpy as np -from numpy import sin,cos - - - -class Triclinic(object): - def __init__(self,a=1,b=1,c=1,alpha=90,beta=90,gamma=90): - self.a = a - self.b = b - self.c = c - alpha = alpha*np.pi/180 - beta = beta*np.pi/180 - gamma = gamma*np.pi/180 - self.alpha = alpha - self.beta = beta - self.gamma = gamma - - self._s11 = b**2 * c**2 * sin(alpha)**2 - self._s22 = a**2 * c**2 * sin(beta)**2 - self._s33 = a**2 * b**2 * sin(gamma)**2 - self._s12 = a*b*c**2*(cos(alpha) * cos(beta) - cos(gamma)) - self._s23 = a**2*b*c*(cos(beta) * cos(gamma) - cos(alpha)) - self._s13 = a*b**2*c*(cos(gamma) * cos(alpha) - cos(beta)) - self.V = (a*b*c)*np.sqrt(1-cos(alpha)**2 - cos(beta)**2 - cos(gamma)**2 + 2*cos(alpha)*cos(beta)*cos(gamma)) - - def __call__(self,h,k,l): return self.q(h,k,l) - - def d(self,h,k,l): - temp = self._s11*h**2 + \ - self._s22*k**2 + \ - self._s33*l**2 + \ - 2*self._s12*h*k+ \ - 2*self._s23*k*l+ \ - 2*self._s13*h*l - d = self.V/np.sqrt(temp) - return d - - def q(self,h,k,l): - return 2*np.pi/self.d(h,k,l) - -class Orthorombic(Triclinic): - def __init__(self,a=1,b=1,c=1): - Triclinic.__init__(self,a=a,b=b,c=c,alpha=90,beta=90,gamma=90) - -class Monoclinic(object): - def __init__(self,a=1,b=1,c=1,beta=90.): - self.a = a - self.b = b - self.c = c - beta = beta/np.pi*180 - self.beta = beta - - self.V = (a*b*c) - - def __call__(self,h,k,l): return self.Q(h,k,l) - - def Q(self,h,k,l): - temp = h**2/self.a**2 + (k*sin(self.beta))**2/self.b**2+l**2/self.c**2+2*h*l*cos(self.beta)/self.a/self.c - d = 1/np.sqrt(temp) - print(d) - return 2*np.pi/d - - - -ti3o5_lambda = Triclinic(a = 9.83776, b = 3.78674, c = 9.97069, beta = 91.2567) -ti3o5_beta = Triclinic(a = 9.7382 , b = 3.8005 , c = 9.4333 , beta = 91.496) -#ti3o5_beta = Monoclinic(a = 9.7382 , b = 3.8005 , c = 9.4333 , beta = 91.496) -ti3o5_alpha = Triclinic(a = 9.8372, b = 3.7921, c = 9.9717) -#ti3o5_alpha1 = Orthorombic(a = 9.8372, b = 3.7921, c = 9.9717)