def test1(self):
v = wripaca.new_doubleArray(9)
for i, x in enumerate([1, 0, 0, 0, 1, 0, 0, 0, 1]):
wripaca.doubleArray_setitem(v, i, x)
assert wripaca.determinant3(v) == 1.0
import sys
sys.path.append("build/lib.win-amd64-2.7")
import wripaca
print wripaca.__file__
import numpy as np
from ImageD11 import transform, gv_general
print "det(eye(3))",wripaca.determinant3( np.eye(3).ravel() )
pre = np.eye(3).ravel()
post = np.eye(3).ravel()
axis = np.array( (0.,0.,-1.))
tth = np.array( [60, 10, 15, 20,25, 1.9555],np.float)
wvln = 0.5
eta = np.array( [90, 10,120,-20,340, -73 ],np.float)
omega = np.array([60,90,180, 60,97, 131],np.float)
def testaxisang():
v = np.array([11.,12.,13.],np.float)
cv = v.copy()
axis = np.array([np.sqrt(1./3),-np.sqrt(1./3),np.sqrt(1./3)],np.float)
for i in range(len(omega)):
a = gv_general.rotation_axis(axis, angle=omega[i])
rv = a.rotate_vectors( v )
wripaca.rotate_vector_axis_angle(