Hi (Pavel),
I have a dataset (P31) which I am processing with the pairwise combinations
of reindexing operators (4 incl. no reindexing therefore h,k,l) and twin
laws (4 incl. no twinning as h,k,l). While some of these combinations are
fine during refinement (Vers 1.11.1), for some combinations I get the
following error, which ultimately kills the refinement. Thanks in advance
and apologies if I am missing something stupid.
Incidentally, a different dataset seems to tolerate all combinations.
*approx_equal eps: 1e-05*
*approx_equal multiplier: 10000000000.0*
*0.512707954357 approx_equal ERROR*
*0.514972857238 approx_equal ERROR*
Traceback (most recent call last):
File "/Applications/phenix-1.11.1-2575/build/../modules/phenix/
phenix/command_line/refine.py", line 11, in <module>
command_line.run(command_name="phenix.refine", args=sys.argv[1:])
File "/Applications/phenix-1.11.1-2575/modules/phenix/phenix/refinement/command_line.py",
line 154, in run
call_back_handler=call_back_handler)
File "/Applications/phenix-1.11.1-2575/modules/phenix/phenix/refinement/driver.py",
line 1529, in run
call_back_handler = call_back_handler).result
File "/Applications/phenix-1.11.1-2575/modules/phenix/phenix/refinement/strategies.py",
line 121, in __init__
log = log)
File "/Applications/phenix-1.11.1-2575/modules/phenix/phenix/refinement/fmodels.py",
line 159, in run
fmodel.update_solvent_and_scale(initialise=True)
File "/Applications/phenix-1.11.1-2575/modules/cctbx_project/
mmtbx/twinning/twin_f_model.py", line 797, in update_solvent_and_scale
params=params, log=log, apply_back_trace=apply_back_trace)
File "/Applications/phenix-1.11.1-2575/modules/cctbx_project/mmtbx/f_model/f_model.py",
line 1457, in update_all_scales
params = params, log = log, refine_hd_scattering=refine_hd_scattering)
File "/Applications/phenix-1.11.1-2575/modules/cctbx_project/
mmtbx/bulk_solvent/f_model_all_scales.py", line 38, in __init__
refine_hd_scattering = refine_hd_scattering, log = log)
File "/Applications/phenix-1.11.1-2575/modules/cctbx_project/
mmtbx/bulk_solvent/f_model_all_scales.py", line 61, in compute
refine_hd_scattering = refine_hd_scattering)
File "/Applications/phenix-1.11.1-2575/modules/cctbx_project/
mmtbx/bulk_solvent/f_model_all_scales.py", line 329, in
update_solvent_and_scale_twin
* assert approx_equal(self.r_all(), o_kbu_sol.r(), 1.e-5)*
*AssertionError*