Dear Phenixers, Is is possible for phenix.refine to output a ccp4 map? If so, how? Thanks, Paul.
Nope.
2008/10/24 Paul Emsley
Dear Phenixers,
Is is possible for phenix.refine to output a ccp4 map? If so, how?
Thanks,
Paul.
_______________________________________________ phenixbb mailing list [email protected] http://www.phenix-online.org/mailman/listinfo/phenixbb
-- ----------------------------------------------------------------- P.H. Zwart Beamline Scientist Berkeley Center for Structural Biology Lawrence Berkeley National Laboratories 1 Cyclotron Road, Berkeley, CA-94703, USA Cell: 510 289 9246 BCSB: http://bcsb.als.lbl.gov PHENIX: http://www.phenix-online.org CCTBX: http://cctbx.sf.net -----------------------------------------------------------------
typically I use coot and the mtz directly.
Have you heard of that program? It is rather nice.
P
2008/10/24 Peter Zwart
Nope.
2008/10/24 Paul Emsley
: Dear Phenixers,
Is is possible for phenix.refine to output a ccp4 map? If so, how?
Thanks,
Paul.
_______________________________________________ phenixbb mailing list [email protected] http://www.phenix-online.org/mailman/listinfo/phenixbb
-- ----------------------------------------------------------------- P.H. Zwart Beamline Scientist Berkeley Center for Structural Biology Lawrence Berkeley National Laboratories 1 Cyclotron Road, Berkeley, CA-94703, USA Cell: 510 289 9246 BCSB: http://bcsb.als.lbl.gov PHENIX: http://www.phenix-online.org CCTBX: http://cctbx.sf.net -----------------------------------------------------------------
-- ----------------------------------------------------------------- P.H. Zwart Beamline Scientist Berkeley Center for Structural Biology Lawrence Berkeley National Laboratories 1 Cyclotron Road, Berkeley, CA-94703, USA Cell: 510 289 9246 BCSB: http://bcsb.als.lbl.gov PHENIX: http://www.phenix-online.org CCTBX: http://cctbx.sf.net -----------------------------------------------------------------
Hi Paul, I do not think it can : *electron_density_maps* *Electron density maps calculation parameters* map_format= *xplor So you have to get the ccp4 map using FFT. Cheers! liu Paul Emsley wrote:
Dear Phenixers,
Is is possible for phenix.refine to output a ccp4 map? If so, how?
Thanks,
Paul.
_______________________________________________ phenixbb mailing list [email protected] http://www.phenix-online.org/mailman/listinfo/phenixbb
Hi Paul, phenix.refine outputs maps in two formats: - CNS/X-plor formatted maps (by default, mFo-DFc and 2mFo-DFc, but you can request any coefficients), and - MTZ file with map coefficients that is good to use in COOT (by default, the map coefficients are for mFo-DFc and 2mFo-DFc maps). Cheers, Pavel. On 10/24/2008 12:23 PM, Paul Emsley wrote:
Dear Phenixers,
Is is possible for phenix.refine to output a ccp4 map? If so, how?
Thanks,
Paul.
_______________________________________________ phenixbb mailing list [email protected] http://www.phenix-online.org/mailman/listinfo/phenixbb
Hi Pavel/Peter,
The reason why I ask is that 64-bit Coot fails to read MTZ files [1] and there
are problems with 'CNS-style' maps from phenix.refine in Coot (they don't cover
the asymmetric unit). A CCP4 map from phenix.refine would be a convenient
work-around.
Paul.
[1] For me, at least.
In message <[email protected]> PHENIX user mailing list
Hi Paul,
phenix.refine outputs maps in two formats:
- CNS/X-plor formatted maps (by default, mFo-DFc and 2mFo-DFc, but you can request any coefficients), and - MTZ file with map coefficients that is good to use in COOT (by default, the map coefficients are for mFo-DFc and 2mFo-DFc maps).
Cheers, Pavel.
On 10/24/2008 12:23 PM, Paul Emsley wrote:
Dear Phenixers,
Is is possible for phenix.refine to output a ccp4 map? If so, how?
Thanks,
Paul.
_______________________________________________ phenixbb mailing list [email protected] http://www.phenix-online.org/mailman/listinfo/phenixbb
_______________________________________________ phenixbb mailing list [email protected] http://www.phenix-online.org/mailman/listinfo/phenixbb
The reason why I ask is that 64-bit Coot fails to read MTZ files [1]
That's weird. We are using the cmtz library essentially directly from the current ccp4 cvs. Reading and writing mtz files is routinely and heavily exercised on all platforms we support, and is known to work flawlessly everywhere, including a bunch of 64-bit platforms. It is not too long ago that I ran the tests with valgrind, which reported 0 errors. Are you compiling with -fno-strict-aliasing? That could in theory be important. If that doesn't help, what are the symptoms of the mtz read failure? Ralf
cc: Kevin
In message <[email protected]> "Ralf W. Grosse-Kunstleve"
The reason why I ask is that 64-bit Coot fails to read MTZ files [1]
That's weird. We are using the cmtz library essentially directly from the current ccp4 cvs. Reading and writing mtz files is routinely and heavily exercised on all platforms we support, and is known to work flawlessly everywhere, including a bunch of 64-bit platforms. It is not too long ago that I ran the tests with valgrind, which reported 0 errors.
Are you compiling with -fno-strict-aliasing? That could in theory be important.
If that doesn't help, what are the symptoms of the mtz read failure?
Hi Ralf, thanks for getting back to me. This info is mostly second-hand of course, I don't have a usable 64 bit machine to hand. These are the symptoms so far (as far as I can gather). It fails (hangs) in fftw (IIRC) and libm (different cases) 0x0000003486b13885 in __ieee754_pow () from /lib64/tls/libm.so.6 (gdb) where #0 0x0000003486b13885 in __ieee754_pow () from /lib64/tls/libm.so.6 #1 0x0000003486b144b8 in __ieee754_pow () from /lib64/tls/libm.so.6 #2 0x0000003486b27903 in pow () from /lib64/tls/libm.so.6 #3 0x0000002a9641b18f in clipper::Grid_sampling::init (this=0x7fbfffc190, spacegroup=Variable "spacegroup" is not available. ) at coords.cpp:399 Yes, indeed I do compile clipper with -fno-strict-aliasing, for both the GTK1 and GTK2 versions. (not anything else though). I did hear that 64 bit GTK1 version works with MTZ files (or can be made to work) but GTK2 version does not (above is gtk2 version). I did try to track this down with Kevin a while ago, the clipper test program for fft works fine with the test data in 64 bits. But not in Coot. Then my disk died and I never got back to booting 64 bit OS. Cheers, Paul.
Hm... sorry, this doesn't ring any bells. I've used fftw on 64-bit architectures (fftw3tbx module in cctbx) without a problem. We compile everything with -fno-strict-aliasing. In the past I tried to cut corners by using -fno-strict-aliasing only for some translation units. When I switched to compiling everything with this option I didn't notice an impact on the runtime. If it hangs in the pow() function, maybe there is a floating-point problem? Unfortunately, gcc's default "non-stop" handling of floating-point exceptions is useless IMO. We install proper traps if possible. When we started doing this, we found a whole bunch of numerical problems that went unnoticed before. (The trap stuff is in boost_adaptbx/meta_ext.cpp.) Ralf
participants (5)
-
junfeng liu
-
Paul Emsley
-
Pavel Afonine
-
Peter Zwart
-
Ralf W. Grosse-Kunstleve