- cctbxbb - phenix-online.org

Re: [cctbxbb] Patterson map data structure
by Ten Eyck, Lynn 16 Oct '12

16 Oct '12

Any Fourier map without an F000 term will analytically have a mean value of 0. Lynn Connected by DROID on Verizon Wireless -----Original message----- From: Nathaniel Echols <nechols(a)lbl.gov> To: cctbx mailing list <cctbxbb(a)phenix-online.org> Sent: Tue, Oct 16, 2012 17:22:17 GMT+00:00 Subject: Re: [cctbxbb] Patterson map data structure On Tue, Oct 16, 2012 at 6:47 AM, Dmytro Guzenko <Dmytro.Guzenko(a)pharm.kuleuven.be> wrote: > rmap = pmap.real_map_unpadded() > n_real = pmap.n_real() > unit_cell = pmap.unit_cell() > grid2cart = maptbx.grid2cart(n_real,unit_cell.orthogonalization_matrix()) > > n_points_total = n_real[0]*n_real[1]*n_real[2] > cart_coord = flex.double(flex.grid(3,n_points_total)) > p_weight = flex.double(n_points_total) > > for plain_ind,grid_point in enumerate(flex.nested_loop(n_real)): > cart_coord[0,plain_ind],cart_coord[1,plain_ind],cart_coord[2,plain_ind] = grid2cart(grid_point) > p_weight[plain_ind] = rmap[grid_point] There are better arrays for storing the coordinates, e.g.: cart_coord = flex.vec3_double(flex.grid(n_real)) p_weight = flex.double(rmap.size()) for plain_ind,grid_point in enumerate(flex.nested_loop(n_real)): cart_coord[grid_point] = grid2cart(grid_point) p_weight[plain_ind] = rmap[grid_point] > By the way, I noticed that many weights are negative, is this possible in a patterson map or am I doing something wrong? This is normal - I believe it's a Fourier series termination error, but I still don't fully understand the reason. An Fobs map will have the same problem; the mean will always be approximately zero. Technically, if the map were perfect with no missing data, it would be uniformly positive, but that's never actually the case. There are ways to "fix" the map so it's all-positive, such as maximum entropy cleanup, but in practice I think we just ignore the negative values (at least for Fobs/2mFo-DFc maps - not sure how Patterson maps are usually handled). -Nat _______________________________________________ cctbxbb mailing list cctbxbb(a)phenix-online.org http://phenix-online.org/mailman/listinfo/cctbxbb

1 0

Patterson map data structure
by Dmytro Guzenko 16 Oct '12

16 Oct '12

Hello all, I am playing with the Patterson function as calculated by cctbx and I am confused as to how I can extract and manipulate the data. More specifically, miller.patterson_map(...).real_map() returns a one-dimensional array of values. How do I get the corresponding 3D coordinates? I suspect it has to do with the gridding and some kind of fixed order enumeration, but I haven't found an example (the peak_search().sites() seems to be relevant, but I would like the whole thing). On a related note - is there a way to debug the python code of cctbx and step into the c++ code? Debugging is my favorite way to figure things out, but with cctbx most interesting things are hidden. I am using Eclipse + PyDev now. Thank you! Kind regards, Dmytro.

3 4

iotbx.xds
by Graeme.Winter＠diamond.ac.uk 16 Oct '12

16 Oct '12

Hi Folks, Would anyone have any objection if I added some code to iotbx.xds to read things like experimental geometry etc in a straightforward way from the various files that the program produces? I see something to read reflection files in there already. Thanks, Graeme -- This e-mail and any attachments may contain confidential, copyright and or privileged material, and are for the use of the intended addressee only. If you are not the intended addressee or an authorised recipient of the addressee please notify us of receipt by returning the e-mail and do not use, copy, retain, distribute or disclose the information in or attached to the e-mail. Any opinions expressed within this e-mail are those of the individual and not necessarily of Diamond Light Source Ltd. Diamond Light Source Ltd. cannot guarantee that this e-mail or any attachments are free from viruses and we cannot accept liability for any damage which you may sustain as a result of software viruses which may be transmitted in or with the message. Diamond Light Source Limited (company no. 4375679). Registered in England and Wales with its registered office at Diamond House, Harwell Science and Innovation Campus, Didcot, Oxfordshire, OX11 0DE, United Kingdom

3 4

no more Python 2.4
by Nathaniel Echols 15 Oct '12

15 Oct '12

FYI for those of you who didn't get my original email: the minimum Python version required for CCTBX is now 2.5, so for anyone contributing to CCTBX, you can now use more advanced syntax such as "with" and "try/except/else/finally" (and decorators, which were already supported in 2.4, but we haven't used them so far). Full list of changes in 2.5 here: http://docs.python.org/whatsnew/2.5.html -Nat

1 0

Python documentation
by David Waterman 10 Oct '12

10 Oct '12

Hi, I've got a few questions about the automatically generated Python documentation. 1. Is this done by doxygen like the C++ documentation? 2. What are the rules that determine whether text in source files is included as documentation or not? Clearly docstrings are, but also some comments do (only those at module scope?) 3. Similarly, what determines which packages are documented or not? For example, from http://cctbx.sourceforge.net/current/python/rstbx.html, cftbx is documented but bpcx is not. Is there at some point a manual decision about what should be documented? Cheers -- David

4 7

Going from a pdb file to an x_ray_structure object
by Claudia Millán Nebot 09 Oct '12

09 Oct '12

Hi everyone, i'm trying to get an xray.structure object from a pdb file, as i'm interested on using both this and a Miller array coming from a reflection file to calculate correlation coefficients between Fobs and Fcalc. For this purpose, i was searching for something that will do the same as the iotbx.pdb.xray_structure.as_pdb_file but the other way around. Is there such a tool? I've look for a while but as i'm new to all that maybe i missed something. Thanks in advance and my excuses if is a too evident question, Claudia Millán (cmncri(a)ibmb.csic.es) Institut de Biologia Molecular de Barcelona (IBMB-CSIC) Barcelona, Spain

3 8

Re: [cctbxbb] RE : cctbx debian package: new commit
by Radostan Riedel 04 Oct '12

04 Oct '12

On Wed, 03. Oct 11:11, Baptiste Carvello wrote: > 2) there are a few C++ applications which depend on cctbx, the one I > know is FOX ( http://vincefn.net/Fox/ , objcryst-fox in Debian). For > those use cases, should we still build a libcctbx-dev package with the > headers and static libs (for those we don't need to care about ABI > stability)? I think this is a very important point here. Let's assume I want to extent the cctbx by adding a module, for let's say absorbtion correction and I want to implement a method called "multiscan" and I want to make a fast python boost extension to integrate with my python code. I need some symbols from any of the shared libs in cctbx. How could this be done in Debian or in any other distribution? About the check for version: In my patch I did this with SCons[1]. What this does is to compile and link a test program and check for success. I saw in your code that you are doing similar things and we could write some simple c++ program that for example includes boost/version and check some "BOOST_VERSION"? This way it checks globally for the availability of a shared library. [1] http://www.scons.org/doc/2.0.1/HTML/scons-user/x4146.html

2 1

Explicit import for iotbx.pdb.amino_acid_codes
by Claudia Millán Nebot 01 Oct '12

01 Oct '12

Hello everyone :) I've just started to look into the cctbx facilities, and as i'm a newbie in both object oriented programming and in general to the system, i'm tryng to play around and find by myself how things are working. Using the cctbx.python, i realised that it is necessary to import explicitly the module amino_acid_codes, while for others like atom_name_interpretation this is done straigthforward if you do an import from a lower level (like from iotbx.pdb import *). Is there some good reason because this has to be like that? Thanks in advance for the answer, and excuse me if is it is a question too simple! Claudia Millán (cmncri(a)ibmb.csic.es) Institut de Biologia Molecular de Barcelona (IBMB-CSIC) Barcelona, Spain

2 2

Build failing
by Graeme.Winter＠diamond.ac.uk 01 Oct '12

01 Oct '12

Hi Folks, Following an svn update and several attempts to build, I am still getting this compile error /Users/graeme/svn/cctbx/sources/boost/boost/math/special_functions/detail/bessel_jy.hpp:346:11: error: 'bessel_y_small_z_series' was not declared in this scope, and no declarations were found by argument-dependent lookup at the point of instantiation [-fpermissive] In file included from /Users/graeme/svn/cctbx/sources/boost/boost/math/special_functions/detail/bessel_jn.hpp:17:0, from /Users/graeme/svn/cctbx/sources/boost/boost/math/special_functions/bessel.hpp:18, from /Users/graeme/svn/cctbx/sources/cctbx_project/scitbx/math/bessel.h:9, from /Users/graeme/svn/cctbx/sources/cctbx_project/scitbx/math/boost_python/math_ext.cpp:3: /Users/graeme/svn/cctbx/sources/boost/boost/math/special_functions/detail/bessel_jy_series.hpp:139:10: note: 'template<class T, class Policy> T boost::math::detail::bessel_y_small_z_series(T, T, T*, const Policy&)' declared here, later in the translation unit In file included from /Users/graeme/svn/cctbx/sources/boost/boost/math/special_functions/bessel.hpp:17:0, from /Users/graeme/svn/cctbx/sources/cctbx_project/scitbx/math/bessel.h:9, from /Users/graeme/svn/cctbx/sources/cctbx_project/scitbx/math/boost_python/math_ext.cpp:3: /Users/graeme/svn/cctbx/sources/boost/boost/math/special_functions/detail/bessel_jy.hpp:356:11: error: 'bessel_j_small_z_series' was not declared in this scope, and no declarations were found by argument-dependent lookup at the point of instantiation [-fpermissive] In file included from /Users/graeme/svn/cctbx/sources/boost/boost/math/special_functions/detail/bessel_jn.hpp:17:0, from /Users/graeme/svn/cctbx/sources/boost/boost/math/special_functions/bessel.hpp:18, from /Users/graeme/svn/cctbx/sources/cctbx_project/scitbx/math/bessel.h:9, from /Users/graeme/svn/cctbx/sources/cctbx_project/scitbx/math/boost_python/math_ext.cpp:3: /Users/graeme/svn/cctbx/sources/boost/boost/math/special_functions/detail/bessel_jy_series.hpp:47:10: note: 'template<class T, class Policy> T boost::math::detail::bessel_j_small_z_series(T, T, const Policy&)' declared here, later in the translation unit In file included from /Users/graeme/svn/cctbx/sources/boost/boost/math/special_functions/bessel.hpp:17:0, from /Users/graeme/svn/cctbx/sources/cctbx_project/scitbx/math/bessel.h:9, from /Users/graeme/svn/cctbx/sources/cctbx_project/scitbx/math/boost_python/math_ext.cpp:3: /Users/graeme/svn/cctbx/sources/boost/boost/math/special_functions/detail/bessel_jy.hpp:363:11: error: 'bessel_yn_small_z' was not declared in this scope, and no declarations were found by argument-dependent lookup at the point of instantiation [-fpermissive] In file included from /Users/graeme/svn/cctbx/sources/boost/boost/math/special_functions/detail/bessel_jn.hpp:17:0, from /Users/graeme/svn/cctbx/sources/boost/boost/math/special_functions/bessel.hpp:18, from /Users/graeme/svn/cctbx/sources/cctbx_project/scitbx/math/bessel.h:9, from /Users/graeme/svn/cctbx/sources/cctbx_project/scitbx/math/boost_python/math_ext.cpp:3: /Users/graeme/svn/cctbx/sources/boost/boost/math/special_functions/detail/bessel_jy_series.hpp:213:3: note: 'template<class T, class Policy> T boost::math::detail::bessel_yn_small_z(int, T, T*, const Policy&)' declared here, later in the translation unit scons: *** [scitbx/math/boost_python/math_ext.o] Error 1 scons: building terminated because of errors. usr+sys time: 26.64 seconds, ticks: 51849408, micro-seconds/tick: 0.514 wall clock time: 8 minutes 5.55 seconds (485.55 seconds total) make: *** [default] Error 2 OS X Lion machine Graemes-MacBook-Pro:build graeme$ uname -a Darwin Graemes-MacBook-Pro.local 11.4.0 Darwin Kernel Version 11.4.0: Mon Apr 9 19:32:15 PDT 2012; root:xnu-1699.26.8~1/RELEASE_X86_64 x86_64 Graemes-MacBook-Pro:build graeme$ gcc --version gcc (GCC) 4.7.0 20120314 (prerelease) Copyright (C) 2012 Free Software Foundation, Inc. This is free software; see the source for copying conditions. There is NO warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. Have I missed something? Nightly builds seem to be OK. Commands ../../cctbx_plus/base/Python.framework/Versions/2.7/Resources/Python.app/Contents/MacOS/Python ../sources/cctbx_project/libtbx/configure.py --compiler=gcc labelit cctbx bpcx_regression rstbx . setpaths.sh make Thanks, Graeme -- This e-mail and any attachments may contain confidential, copyright and or privileged material, and are for the use of the intended addressee only. If you are not the intended addressee or an authorised recipient of the addressee please notify us of receipt by returning the e-mail and do not use, copy, retain, distribute or disclose the information in or attached to the e-mail. Any opinions expressed within this e-mail are those of the individual and not necessarily of Diamond Light Source Ltd. Diamond Light Source Ltd. cannot guarantee that this e-mail or any attachments are free from viruses and we cannot accept liability for any damage which you may sustain as a result of software viruses which may be transmitted in or with the message. Diamond Light Source Limited (company no. 4375679). Registered in England and Wales with its registered office at Diamond House, Harwell Science and Innovation Campus, Didcot, Oxfordshire, OX11 0DE, United Kingdom

3 8

dropping Python 2.4 support?
by Nathaniel Echols 27 Sep '12

27 Sep '12

Hi all-- Will it ruin anyone's work if we make Python 2.5 the minimum platform for using CCTBX? 2.4 is now almost eight years old and the lack of certain newer features (decorators, 'with') is starting to hold us back. thanks, Nat

2 2