def intersect(m1,m2, use_average): maptbx.intersection( map_data_1 = m1, map_data_2 = m2, thresholds = flex.double([0,0.1,0.2,0.3,0.4,0.5]), average = False) if(use_average): return (m1+m2)/2 else: return m1
def intersect(m1, m2, use_average): maptbx.intersection(map_data_1=m1, map_data_2=m2, thresholds=flex.double( [0, 0.1, 0.2, 0.3, 0.4, 0.5]), average=False) if (use_average): return (m1 + m2) / 2 else: return m1
def compute_map_and_combine(map_coeffs, crystal_gridding, map_data, thresholds=[0, 0.1, 0.2, 0.3, 0.4, 0.5]): m = get_map(map_coeffs=map_coeffs, crystal_gridding=crystal_gridding) if (map_data is None): map_data = m else: maptbx.intersection(map_data_1=m, map_data_2=map_data, thresholds=flex.double(thresholds), average=True) return map_data
def compute_map_and_combine( map_coeffs, crystal_gridding, map_data, thresholds = [0,0.1,0.2,0.3,0.4,0.5]): m = get_map(map_coeffs=map_coeffs, crystal_gridding=crystal_gridding) if(map_data is None): map_data = m else: maptbx.intersection( map_data_1 = m, map_data_2 = map_data, thresholds = flex.double(thresholds), average = True) return map_data
def exercise_intersection(): thresholds = flex.double([0,0.1,0.2,0.3,0.4,0.5, 0.6,0.7,0.8,0.8, 1.0]) def get_map(): av = [random.random() for i in xrange(10*20*30)] m = flex.double(av) m.resize(flex.grid((10,20,30))) return m m1 = get_map() m2 = get_map() for average in [True, False]: maptbx.intersection( map_data_1 = m1, map_data_2 = m2, thresholds = thresholds, average = average)