b.as_scitbx_matrix()
Out[50]: matrix.rec(elems=(0.0, ..., 0.0), n=(2,2))
'''
from scitbx import matrix
from dials.scratch.luiso_s import tst_prod , write_2d
from scitbx.array_family import flex
a = flex.double(flex.grid(3, 3))
#b = flex.double(flex.grid(3, 3))
import random
for xpos in range(3):
for ypos in range(3):
a[ypos, xpos] = random.random()
# b[ypos, xpos] = 1.0
print "a ="
write_2d(a)
#print "b ="
#write_2d(b)
a_mat = a.as_scitbx_matrix()
print a_mat.inverse()
x_mat = a_mat.inverse()
print "inverse(a)= "
x = x_mat.as_flex_double_matrix()
write_2d(x)
from __future__ import division
from dials.scratch.luiso_s import tst_prod , write_2d
from scitbx.array_family import flex
a = flex.double(flex.grid(3, 3))
b = flex.double(flex.grid(3, 3))
import random
contr = 0
for xpos in range(3):
for ypos in range(3):
a[ypos, xpos] = random.random()
b[ypos, xpos] = 1.0
contr += 1.0
print "a ="
write_2d(a)
print "b ="
write_2d(b)
a_mat = a.as_scitbx_matrix()
b_mat = b.as_scitbx_matrix()
x_mat = a_mat * b_mat
x = x_mat.as_flex_double_matrix()
print "x = a * b ="
x_np = x.as_numpy_array()
from __future__ import division from scitbx.array_family import flex from dials.scratch.luiso_s import model_2d, write_2d, test_compare_2d from dials.algorithms.integration import add_2d nrow = 10 ncol = 10 sumation = flex.double(flex.grid(21, 21)) descr = flex.double(flex.grid(1, 3)) descr[0, 0] = 0.5 descr[0, 1] = 5 descr[0, 2] = 5 ref2d = model_2d(5, 5, 2, 1, 0.2, 55, 0.5) sumation = add_2d(descr, ref2d, sumation) # sumation = test_compare_2d(descr, ref2d, sumation) write_2d(sumation)
from dials.scratch.luiso_s import tst_prod, write_2d
from scitbx.array_family import flex
a = flex.double(flex.grid(3, 3))
b = flex.double(flex.grid(3, 3))
import random
contr = 0
for xpos in range(3):
for ypos in range(3):
a[ypos, xpos] = random.random()
b[ypos, xpos] = 1.0
contr += 1.0
print("a =")
write_2d(a)
print("b =")
write_2d(b)
a_mat = a.as_scitbx_matrix()
b_mat = b.as_scitbx_matrix()
x_mat = a_mat * b_mat
x = x_mat.as_flex_double_matrix()
print("x = a * b =")
x_np = x.as_numpy_array()
print(x_np)
print("= as_flex")
from __future__ import division from scitbx.array_family import flex from dials.scratch.luiso_s import model_2d, write_2d, test_compare_2d from dials.algorithms.integration import add_2d nrow = 10 ncol = 10 sumation = flex.double(flex.grid(21, 21)) descr = flex.double(flex.grid(1, 3)) descr[0, 0] = .5 descr[0, 1] = 5 descr[0, 2] = 5 ref2d = model_2d(5, 5, 2, 1, 0.2, 55, 0.5) sumation = add_2d(descr, ref2d, sumation) #sumation = test_compare_2d(descr, ref2d, sumation) write_2d(sumation)
from __future__ import print_function
from scitbx.array_family import flex
from dials.scratch.luiso_s import model_2d, write_2d
from dials.algorithms.integration import add_2d
nrow = 10
ncol = 10
sumation = flex.double(flex.grid(9, 9), 1)
descr = flex.double(flex.grid(1, 3))
descr[0, 0] = 2.5
descr[0, 1] = 2.5
descr[0, 2] = 1
# double centr_col = descriptor(0,0);
# double centr_row = descriptor(0,1);
# double scale = descriptor(0,2);
ref2d = model_2d(5, 5, 2, 1, 0.2, 55, 0.5)
write_2d(ref2d)
write_2d(sumation)
print("<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<")
sumation01 = add_2d(descr, ref2d, sumation)
write_2d(ref2d)
write_2d(sumation01)
write_2d(sumation)
from scitbx.array_family import flex from dials.scratch.luiso_s import model_2d, write_2d from dials.algorithms.integration import add_2d nrow = 10 ncol = 10 sumation = flex.double(flex.grid(9, 9),1) descr = flex.double(flex.grid(1, 3)) descr[0, 0] = 2.5 descr[0, 1] = 2.5 descr[0, 2] = 1 #double centr_col = descriptor(0,0); #double centr_row = descriptor(0,1); #double scale = descriptor(0,2); ref2d = model_2d(5, 5, 2, 1, 0.2, 55, 0.5) write_2d(ref2d) write_2d(sumation) print "<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<" sumation01 = add_2d(descr, ref2d, sumation) write_2d(ref2d) write_2d(sumation01) write_2d(sumation)
from __future__ import division
from dials.scratch.luiso_s import tst_ref_prod , write_2d
from scitbx.array_family import flex
a = flex.double(flex.grid(3, 3))
b = flex.double(flex.grid(3, 3))
import random
for xpos in range(3):
for ypos in range(3):
a[ypos, xpos] = random.random()
b[ypos, xpos] = 1.0
print "a ="
write_2d(a)
print "b ="
write_2d(b)
c = tst_ref_prod(a, b)
print "c ="
write_2d(c)
print "______________________________________"
print "a ="
write_2d(a)
print "b ="
write_2d(b)
from __future__ import division
from __future__ import print_function
from dials.scratch.luiso_s import tst_ref_prod, write_2d
from scitbx.array_family import flex
a = flex.double(flex.grid(3, 3))
b = flex.double(flex.grid(3, 3))
import random
for xpos in range(3):
for ypos in range(3):
a[ypos, xpos] = random.random()
b[ypos, xpos] = 1.0
print("a =")
write_2d(a)
print("b =")
write_2d(b)
c = tst_ref_prod(a, b)
print("c =")
write_2d(c)
print("______________________________________")
print("a =")
write_2d(a)
