def do_mosaic(self, alg):
print("-" * 79, file=self.log)
print("Refine k_masks", file=self.log)
result = mosaic.refinery(
fmodel=self.fmodel_0,
fv=self.mm.FV,
anomaly=True, # Refine all contributions at once!
alg=alg,
log=self.log)
print("", file=self.log)
result.fmodel.show(show_header=False, show_approx=False, log=self.log)
print(result.fmodel.r_factors(prefix=" "), file=self.log)
return result
def do_mosaic(self, alg):
print("-" * 79, file=self.log)
print("Refine k_masks", file=self.log)
assert approx_equal(self.fmodel_2013_opt.f_obs().data(),
self.fmodel_0.f_obs().data())
assert approx_equal(self.fmodel_2013_opt.f_calc().data(),
self.fmodel_0.f_calc().data())
if (self.fmodel_0.r_two_bins_lowest() <
self.fmodel_2013_opt.r_two_bins_lowest()):
fmodel = self.fmodel_0
else:
fmodel = self.fmodel_2013_opt
result = mosaic.refinery(fmodel=fmodel,
fv=self.mm.FV,
alg=alg,
log=self.log)
print("", file=self.log)
result.fmodel.show(show_header=False, show_approx=False, log=self.log)
print(result.fmodel.r_factors(prefix=" "), file=self.log)
return result