- cctbxbb - phenix-online.org

spec file for cctbx rpm (on F10)
by Cameron Mura 28 Jul '09

28 Jul '09

Hi all, In case it's of interest, I've attached a spec file that works for building a cctbx RPM on Fedora 10 (nothing really F10-specific, so may work on other distros/releases too). Please note that it works by playing the dangerous (and extremely ugly) game of moving the RPM build into the system-wide tree, so it builds "in situ" in /usr/local/cctbx instead of in the RPM default locale (/my/rpm/BUILDROOT/usr/local/cctbx). On the same note, the patch which is referenced in the spec file ('cctbx_bundle_StrucBio.patch') is just a small snippet that edits the head of 'cctbx_install_script.csh' to explicitly control the build path (i.e., override $cwd by explicitly setting something like 'set install_root="/usr/local/cctbx/2009_02_15_2320"', or wherever you'd like cctbx to live...). Probably a much cleaner procedure would be to use the libtbx/configure.py approach in a post-processing step; let me know if anyone pkg's an RPM that way. :-) Enjoy, Cam -- ------------------------------------------------------- Cameron Mura Assistant Professor University of Virginia Department of Chemistry McCormick Road, Box 400319 Charlottesville, VA 22904-4319 Email: cmura(a)virginia.edu Web: http://muralab.org Tel: 434.924.7824 Fax: 434.924.3710 ------------------------------------------------------- #------------------------------------------------------------------------------- # Cameron Mura <cmura(a)virginia.edu>; 27jul09; #------------------------------------------------------------------------------- %define _buildroot /usr/local/cctbx/2009_02_15_2320 %define debug_package %{nil} #%define _builddir /usr/local/cctbx/2009_02_15_2320 Name: cctbx Version: 2009_02_15_2320 Release: 1 Summary: Computational Crystallography Toolbox License: Open Source (BSD-like) Group: Applications/Engineering URL: http://cctbx.sourceforge.net Source0: %{name}_bundle_%{version}.tar.gz Source1: ccp4_mon_lib.tar.gz Packager: Cameron Mura <cmura(a)virginia.edu> Patch0: cctbx_bundle_StrucBio.patch #BuildRoot: /var/tmp/%{name}-buildroot BuildRoot: /usr/local/cctbx/2009_02_15_2320 %description The Computational Crystallography Toolbox (cctbx) is being developed as the open source component of the PHENIX system. The goal of the PHENIX project is to advance automation of macromolecular structure determination. PHENIX depends on the cctbx, but not vice versa. This hierarchical approach enforces a clean design as a reusable library. The cctbx is therefore also useful for small-molecule crystallography and even general scientific applications. The cctbx is designed with an open and flexible architecture to promote extendability and easy incorporation into other software environments. The package is organized as a set of ISO C++ classes with Python bindings. This organization combines the computational efficiency of a strongly typed compiled language with the convenience and flexibility of a dynamically typed scripting language in a strikingly uniform and very maintainable way. Note that CCTBX is developed under an open source (BSD-like) license; for addition information, see http://cctbx.svn.sourceforge.net/viewvc/cctbx/trunk/cctbx/LICENSE_2_0.txt. %prep # Require weird gymnastics to relocate cctbx outside the RPM build tree... #export MYCCTBXHOME="/usr/local/cctbx/2009_02_15_2320" #if [ ! -d $MYCCTBXHOME ]; then # mkdir -p $MYCCTBXHOME #fi %setup -a 0 -D -c /usr/local/cctbx/2009_02_15_2320 -n . cp %{_topdir}/SOURCES/ccp4_mon_lib.tar.gz /usr/local/cctbx/2009_02_15_2320/. #%setup -a 1 -D -T -c /usr/local/cctbx/2009_02_15_2320 -n . %patch0 -p1 %build /usr/local/cctbx/2009_02_15_2320/cctbx_install_script.csh /usr/bin/python 8 . "/usr/local/cctbx/2009_02_15_2320/cctbx_build/setpaths.sh" libtbx.unpack_in_sources ccp4_mon_lib.tar.gz %install mkdir -p $RPM_BUILD_ROOT/usr/local/cctbx/2009_02_15_2320 mv /usr/local/cctbx/2009_02_15_2320/* $RPM_BUILD_ROOT/usr/local/cctbx/2009_02_15_2320/. mkdir -p $RPM_BUILD_ROOT/etc/profile.d echo '. "/usr/local/cctbx/2009_02_15_2320/cctbx_build/setpaths.sh"' > $RPM_BUILD_ROOT/etc/profile.d/cctbx.sh echo 'source "/usr/local/cctbx/2009_02_15_2320/cctbx_build/setpaths.csh"' > $RPM_BUILD_ROOT/etc/profile.d/cctbx.csh #chmod 755 $RPM_BUILD_ROOT/etc/profile.d/cctbx.sh #chmod 755 $RPM_BUILD_ROOT/etc/profile.d/cctbx.csh #[ $RPM_BUILD_ROOT != "/" ] && rm -rf $RPM_BUILD_ROOT #install -m 755 src/rasmol_32BIT $RPM_BUILD_ROOT/usr/local/rasmol #install -m 755 src/rasmol_XFORMS_32BIT $RPM_BUILD_ROOT/usr/local/rasmol #install -m 644 doc/rasmol.hlp $RPM_BUILD_ROOT/usr/local/rasmol #install -m 644 INSTALL* $RPM_BUILD_ROOT/usr/local/rasmol #install -m 644 README* $RPM_BUILD_ROOT/usr/local/rasmol #install -m 644 TODO* $RPM_BUILD_ROOT/usr/local/rasmol #cp -R working/cctbx_sources $RPM_BUILD_ROOT/usr/local/cctbx/%{version} #cp -R working/cctbx_build $RPM_BUILD_ROOT/usr/local/cctbx/%{version} %clean [ $RPM_BUILD_ROOT != "/" ] && rm -rf $RPM_BUILD_ROOT %post #ln -fs /usr/local/rasmol/rasmol_32BIT /usr/local/bin/rasmol %preun #%{__rm} -f /usr/local/bin/rasmol %files %defattr(-,root,root) #%doc ChangeLog.txt ChangeLog.html INSTALL* NOTICE README* /usr/local/cctbx/2009_02_15_2320 /etc/profile.d/%{name}.sh /etc/profile.d/%{name}.csh %changelog * Sun Jul 26 2009 Cameron Mura <cmura(a)virginia.edu) 2009_02_15_2320-1 - Initial RPM build on Fedora 10...

1 1

draw_symops <-- cctbx(/phenix) broken due to boost segfault handling ?
by Cameron Mura 06 Jul '09

06 Jul '09

Hi all, Apologies for the message to both mailing lists, but perhaps someone on either cctbxbb or pymol-users has had some recent experience with this... The subject line more-or-less says it all: I'm wondering if anyone has recently used Rob Campbell's draw_symops (or draw_symops_cctbx) Python modules with a relatively recent (within past year) version of the cctbx (particularly the Boost C++ parts) ? I succeeded in doing this a few years ago to generate images such as the ones near the top of http://muralab.org/~cmura/PyMOL/Sandbox/, but I can't get it to work on my present system (64-bit Fedora 10; Python 2.5.2 with Python 2.6 coexisting). Below is what happens when I... 1) Initiate a new PyMOL session (using a PyMOL v1.1 that I built in-place from source, so should be compatible with local python, system libs, etc) 2) load PDB files (tried with a few diff't files; verified that they contain proper CRYST specifications) 3) do the usual 'run draw_symops_cctbx.py' 4) call draw_symops() -- either implicitly, or by explicit specification of the u.c. params (i.e., using something like "draw_symops_param((45.2,45.2,70.8,90.,90.,120.),'p3121',0.5,0.1)") ...which crashes the pymol session and yields the following trace: > PyMOL>draw_symops_param((45.2,45.2,70.8,90.,90.,120.),'p3121',0.5,0.1) > libc backtrace (48 frames, most recent call last): > /usr/bin/python [0x400649] > /lib64/libc.so.6(__libc_start_main+0xe6) [0x344081e576] > /usr/lib64/libpython2.5.so.1.0(Py_Main+0xb11) [0x34584e7ef1] > *<snip> > ... > ... > </snip>* > /usr/lib64/libpython2.5.so.1.0(PyObject_Call+0x13) [0x345843d493] > /root/cctbx/working/v4/cctbx_build/lib/libboost_python.so > [0x7ff19691db16] > > /root/cctbx/working/v4/cctbx_build/lib/libboost_python.so(boost::python::handle_exception_impl(boost::function0<void>)+0x9f) > [0x7ff19692a2ef] > /root/cctbx/working/v4/cctbx_build/lib/libboost_python.so > [0x7ff196920b08] > > /root/cctbx/working/v4/cctbx_build/lib/libboost_python.so(boost::python::objects::function::call(_object*, > _object*) const+0x125) [0x7ff196920635] > > /root/cctbx/working/v4/cctbx_build/lib/cctbx_sgtbx_ext.so(boost::python::detail::caller_arity<2u>::impl<cctbx::sgtbx::rt_mx > (*)(cctbx::sgtbx::space_group const&, unsigned long), > boost::python::default_call_policies, > boost::mpl::vector3<cctbx::sgtbx::rt_mx, cctbx::sgtbx::space_group > const&, unsigned long> >::operator()(_object*, _object*)+0x12a) > [0x7ff18f820a3a] > /root/cctbx/working/v4/cctbx_build/lib/cctbx_sgtbx_ext.so > [0x7ff18f8034cc] > > /root/cctbx/working/v4/cctbx_build/lib/libboost_python.so(boost::python::throw_error_already_set()+0xe) > [0x7ff196929dde] > /usr/lib64/libstdc++.so.6(__cxa_allocate_exception+0x2f) [0x3444cc2cdf] > /lib64/libc.so.6 [0x3440832f90] > Segmentation fault (libc call stack above) > This crash may be due to a problem in any imported > Python module, including modules which are not part > of the cctbx project. To disable the traps leading > to this message, define these environment variables > (e.g. assign the value 1): > BOOST_ADAPTBX_FPE_DEFAULT > BOOST_ADAPTBX_SIGNALS_DEFAULT > This will NOT solve the problem, just mask it, but > may allow you to proceed in case it is not critical. > PyMOL: abrupt program termination. I tried the recommendation of setting the BOOST_ADAPTBX* env vars, but this doesn't remedy the problem. Also, I tried accomplishing this with many combinations of workflows: 1) old cctbx (e.g., 2005_04_29_1615), latest stable release of cctbx (2009_02_15_2320), bleeding-edge cctbx (2009_07_04_0143) 2) pre-built cctbx binaries versus building myself from src (for each of the above), insuring, again, optimal match between local system libs... 3) cctbx alone, cctbx in the context of Phenix, cctbx +/- its own bundled Python (2.5 and 2.6) Also, I tried this with different combinations of Rob's scripts -- e.g., the 'symop_axes.dat'-dependent "draw_symops.py", versus the 'all_axes_new.py'-dependent "draw_symops_cctbx.py". Interestingly, I found that get_all_axes() from 'all_axes_new.py' works fine for different PDB files. Investigating the above draw_symops stacktrace led me to this post from Ralf Grosse-Kunstleve to python-dev: http://mail.python.org/pipermail/python-dev/2008-September/082639.html It looks like exactly the same problem that's causing the pymol/draw_symops scripts to choke. Also, I can verify that issuing commands such as in Ralf's post (boost_adaptbx.segmentation_fault, boost_adaptbx.divide_by_zero) result in the same segfaults and errors as in his post, so... I believe this is the crux of it. Before delving into Boost codebase, cctbx's exception-handling in C++/Python Boost, etc., I'm wondering if anyone else has seen anything like this...? or has advice on how best to proceed? Any tips or suggestions would be most greatly appreciated. Thanks, Cam -- ------------------------------------------------------- Cameron Mura Assistant Professor University of Virginia Department of Chemistry McCormick Road, Box 400319 Charlottesville, VA 22904-4319 Email: cmura(a)virginia.edu Web: http://muralab.org Tel: 434.924.7824 Fax: 434.924.3710 -------------------------------------------------------

2 3

Reindexing matrix -> transform orientation matrix
by Winter, G (Graeme) 05 May '09

05 May '09

Hi Folks, I'm just trying to figure out how to do the following, given a triclinic orientation matrix (real or reciprocal space, differ by .inverse() -> triclinic unit cell) (like iotbx.lattice_symmetry) stuff -> set of unit cell, Bravais lattice options with reindex operations apply reindex operation to orientation matrix -> new basis a, b, c The former I have historically done using iotbx.lattice_symmetry but given the move from python to cctbx.python, it seems wasteful. I looked at the source code but got a little confused. The latter seems to be more problematic - even reconstructing the change of basis from a known reindexing operation I don't seem to be able to get the right answer :o( Has anyone done anything like this with the CCTBX tools? The context is autoindexing with XDS using a full sweep's worth of images, then inspecting the indexes in SPOT.XDS based on the contents of XPARM.XDS. Thanks, Graeme Graeme Winter Software and MX Support Scientist Diamond Light Source +44 1235 778091 (work) +44 7786 662784 (work mobile) 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

2 2

cctbx geometry modules?
by Artem Evdokimov 31 Dec '08

31 Dec '08

Hi, Sorry to ask a possibly obvious question - but I am very new to cctbx and an absolute novice of python :-) Apologies in advance. Quesiton 1 - is there a web site or some other form of a source for a complete cctbx documentation, most importantly with simple working examples? I've been using some of the examples published in the newsletters but they're scattered here and there. Maybe there is a cctbx wiki somewhere? Question 2 - is there anything inside cctbx that can compute elements of protein geometry - such as dihedral angles, planes, etc. (given selections of atoms)? Thanks in advance! Happy New Year! Cheers, Artem

2 1

Re: [cctbxbb] structure factor interpolation
by celine Besnard 19 Nov '08

19 Nov '08

Hi Ralf and Jon, Thanks a lot for your help. I will try to follow your advices. Celine -- Dr. Céline Besnard EPFL - SB - IPMC - LCR2 CH-1015 Lausanne Suisse e-mail : Celine.Besnard(a)epfl.ch Tel : +41(0)21.693.05.41 Fax : +41(0)21.693.05.04 http://lcr.epfl.ch

1 0

structure factor interpolation
by celine Besnard 07 Nov '08

07 Nov '08

Hi, I am looking for a function that would allow me to interpolate the structure factors values in between the reciprocal lattice points (h,k,l not integer) so that I can get the reciprocal map not of the whole crystal but of "only one" unit-cell. Is there a simple way to do that in the cctbx ? Thanks Celine -- Dr. Céline Besnard EPFL - SB - IPMC - LCR2 CH-1015 Lausanne Suisse e-mail : Celine.Besnard(a)epfl.ch Tel : +41(0)21.693.05.41 Fax : +41(0)21.693.05.04 http://lcr.epfl.ch

3 3

postdoc position in LBNL
by Peter Zwart 21 Oct '08

21 Oct '08

Dear All, We would like to bring the following postdoctoral research position to your attention. Within Lawrence Berkeley Laboratory, we have an opening for a postdoc to work on the development of new software for the analysis of small angle X-ray scattering (SAXS) data. The availability of user friendly, robust software analysis tools for SAXS data is of growing importance as the biological community is increasingly using SAXS as a standard method to analyze complex biological systems. We are developing a new set of tools, including a software library, for SAXS data analysis. The successful candidate will have the opportunity to work with experienced software developers in the field of biophysics and interact with users and beamline staff on SAXS beamline 7.3.3 at the Advanced Light Source (ALS) in Berkeley California. Interested candidates can contact me (Peter Zwart, PHZwart(a)lbl.gov) for more details, or visit the LBL job site at: http://jobs.lbl.gov/ Sincerely, Peter Zwart -- ----------------------------------------------------------------- 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 -----------------------------------------------------------------

1 0

Re: [cctbxbb] interpreting auto generated python interface documentation
by Ralf W. Grosse-Kunstleve 28 Aug '08

28 Aug '08

> Say , for a particular class say http://cctbx.sourceforge.net/current/python/iotbx.pdb.atom_selection.html, > the documentation gives me a "cache" class __init__ method that says > > __init__(self, root, wildcard_escape_char=r'\') > > How do I find out what "root" is in this example. Should I be looking at the c++ code to find out the "types" of these arguments. It is a iotbx.pdb.hierarchy.root instance in this case, mostly referred to as "pdb_hierarchy" or just "hierarchy". Yes, for some things you have to go into the code. The best approach is to insert print statements to see what the objects are, then run the associated tests (e.g. iotbx/pdb/tst_*.py). > In addition a few of the auto-generated documentation pages such as http://cctbx.sourceforge.net/current_cvs/python/iotbx.pdb.atom.html > return an error page. Thanks for pointing this out. I'll look into this next week, when I'm back home. > I just want to know a good way of understanding the varrious interfaces when writing code using cctbx in Python Don't hesitate to ask. You may have many questions at first, but then they'll become fewer and far in between. I'll respond more promptly after the conference trip. Ralf

1 0

Re: [cctbxbb] input_atom_labels_list() not available in Python?
by Ralf W. Grosse-Kunstleve 27 Aug '08

27 Aug '08

Hi Hari, > input_atom_labels_list() > This function seems to be available from both c++ and python. The Nov 2006 letter has a python example that uses > this function > However I dont see that in the python interface with the current version? Is that true or am I using a different build > that does not make this available in python That method is gone, indeed, ever since about April/May this year. Look for "atoms_with_labels": cd cctbx_sources libtbx.find_files -g atoms_with_labels -f Ralf P.S.: Have to run now... will look at your other email later.

1 0

input_atom_labels_list() not available in Python?
by hari jayaram 26 Aug '08

26 Aug '08

Thanks Nathaniel.. I too was using the approach of looking at the phenix scripts to figure things out. But I was hoping that there was a "google hidden" source of documentation for the python interfaces that I was not aware of. For the present moment the c++ documentation is much more useful .After I sent the email out , I did notice the section in the manual which says.. "Most documented C++ interfaces are also available at the Python layer. Unfortunately the documentation tools available are not capable of merging the documentations. Therefore Python users need to also consult the C++ documention." On an unrelated note The Nov 2006 newsletter and the c++ and python documentation lists a function in the class iotbx.pdb.input called input_atom_labels_list() This function seems to be available from both c++ and python. The Nov 2006 letter has a python example that uses this function However I dont see that in the python interface with the current version? Is that true or am I using a different build that does not make this available in python Just curious Hari On Tue, Aug 26, 2008 at 12:50 PM, Nathaniel Echols < nathaniel.echols(a)gmail.com> wrote: > Having spent the last few months trying to learn the python side of cctbx > for my own work, I recommend looking at some of the actual programs that > call these functions (e.g. the various phenix scripts, or testing scripts) > for clues on how to use them. In particular, > iotbx/pdb/tst_atom_selection.py has examples of selection caches. (I'm not > sure how exactly the "cache" class fits in, since it isn't invoked directly > in that script, but hopefully the test script has some approximation of what > you're trying to do.) I haven't found it necessary to look at the C++ code > - boost.python obscures many details anyway, and it's usually pretty obvious > to me what's going on in the python sources based on the > variable/class/function names. > The various CCP4 newsletter articles may have some useful suggestions, but > I haven't looked at these in much detail yet. > > On Tue, Aug 26, 2008 at 9:18 AM, hari jayaram <harijay(a)gmail.com> wrote: > >> Hi , >> I am a few days into learning to use the cctbx libraries. My interests >> are mostly to be able to write scripts that plot various statistics from mtz >> files and use things like iotbx.pdb to manipulate coordinate files. While I >> have managed to write a few scripts and have read the very descriptive >> examples in the newsletters and the few routines inside of phenix. I have >> a problem with getting much out of the python interface documentation at the >> sourceforge site. >> >> Say , for a particular class say >> http://cctbx.sourceforge.net/current/python/iotbx.pdb.atom_selection.html >> , >> the documentation gives me a "cache" class __init__ method that says >> >> __init__(self, root, wildcard_escape_char=r'\') >> >> How do I find out what "root" is in this example. Should I be looking at >> the c++ code to find out the "types" of these arguments. >> >> In addition a few of the auto-generated documentation pages such as >> http://cctbx.sourceforge.net/current_cvs/python/iotbx.pdb.atom.html >> return an error page. >> >> I just want to know a good way of understanding the varrious interfaces >> when writing code using cctbx in Python >> Thanks in advance >> hari >> >> >> _______________________________________________ >> cctbxbb mailing list >> cctbxbb(a)phenix-online.org >> http://www.phenix-online.org/mailman/listinfo/cctbxbb >> >> > > _______________________________________________ > cctbxbb mailing list > cctbxbb(a)phenix-online.org > http://www.phenix-online.org/mailman/listinfo/cctbxbb > >

1 0

interpreting auto generated python interface documentation
by hari jayaram 26 Aug '08

26 Aug '08

Hi , I am a few days into learning to use the cctbx libraries. My interests are mostly to be able to write scripts that plot various statistics from mtz files and use things like iotbx.pdb to manipulate coordinate files. While I have managed to write a few scripts and have read the very descriptive examples in the newsletters and the few routines inside of phenix.I have a problem with getting much out of the python interface documentation at the sourceforge site. Say , for a particular class say http://cctbx.sourceforge.net/current/python/iotbx.pdb.atom_selection.html, the documentation gives me a "cache" class __init__ method that says __init__(self, root, wildcard_escape_char=r'\') How do I find out what "root" is in this example. Should I be looking at the c++ code to find out the "types" of these arguments. In addition a few of the auto-generated documentation pages such as http://cctbx.sourceforge.net/current_cvs/python/iotbx.pdb.atom.html return an error page. I just want to know a good way of understanding the varrious interfaces when writing code using cctbx in Python Thanks in advance hari

2 1

Re: [cctbxbb] getting resolution for any "h" "k" "l" and cctbx
by Ralf W. Grosse-Kunstleve 05 Aug '08

05 Aug '08

Hi Hari, > 1) I was wondering if someone can point me to routines available from > cctbx modules that can give me the "resolution" for a given index "h" > "k" "l". in any space group To get the resolution for a given hkl, you just need the unit cell: from cctbx import uctbx uc = uctbx.unit_cell((11,12,13,93,98,103)) hkl = (2,3,4) print uc.d(miller_index=hkl) print uc.stol(miller_index=hkl) wavelength = 1.3 print uc.two_theta(miller_index=hkl, wavelength=wavelength) Save as "resolution_hkl.py" and run with: cctbx.python resolution_hkl.py The .d(), .stol(), .two_theta() methods are documented here: http://cctbx.sourceforge.net/current/c_plus_plus/classcctbx_1_1uctbx_1_1uni… > 2) In addition , if someone can point me in the direction of a well > illustrated depiction of the derivation of the reciprocal lattice > for the triclinic case , that will be very much appreciated bonus . I'm not sure which textbook to recommend. I liked the old Buerger books a lot, but I don't remember how he introduces the reciprocal lattice. Ralf

1 0

getting resolution for any "h" "k" "l" and cctbx
by hari jayaram 05 Aug '08

05 Aug '08

Hi , I have been programming in python for a little bit, I am new to cctbx . My crystallography knowledge is as good as solving a few structures gives you. I have a few questions. 1) I was wondering if someone can point me to routines available from cctbx modules that can give me the "resolution" for a given index "h" "k" "l". in any space group 2) In addition , if someone can point me in the direction of a well illustrated depiction of the derivation of the reciprocal lattice for the triclinic case , that will be very much appreciated bonus . Thank you in advance Hari

1 0

Re: [cctbxbb] (no subject)
by Ralf W. Grosse-Kunstleve 26 May '08

26 May '08

Hi Xie, > Are cctbx and phenix compatible with a 64 bit operating system? Yes. > I have Ubuntu 8.04 (64 bit) installed on an Intel Core2 duo system. Please try the source bundle: cctbx_python_252_bundle.selfx http://cci.lbl.gov/cctbx_build/ Let me know if there are any problems. Ralf

1 0

(no subject)
by Xie Jiabao 26 May '08

26 May '08

Hi, Are cctbx and phenix compatible with a 64 bit operating system? I have Ubuntu 8.04 (64 bit) installed on an Intel Core2 duo system. Thanks, Xie

1 0

Re: [cctbxbb] problem in "space_group_type"
by Ralf W. Grosse-Kunstleve 16 Apr '08

16 Apr '08

Hi Koen, > The error only occurs when running in Debug mode. In Release mode, > everything works fine, suggesting a memory initialization problem. This is a known feature of Windows: libraries compiled in release and debug mode are not binary compatible. Under Linux I check with valgrind all the time, therefore it is almost certain that the Windows debug/release mix is the only problem. > I'm unable to debug the problem myself or get more details, as I can't > get the library to compile with debugging information in it (Is there an > easy way to do this with the Scons build system? I haven't used Scons > before.) The *one* script to change is libtbx/SConscript Currently, in the if (sys.platform == "win32"): ... if (env_etc.compiler == "win32_cl"): branch you'll find: if (libtbx.env.build_options.debug_symbols): raise RuntimeError("Debug build not supported.") This is the place to start. On Unix platforms it simply works like this: if (libtbx.env.build_options.debug_symbols): opts.insert(0, "-g") I forgot the details of why it is more complicated under Windows. I vaguely recall that there were issues linking with the correct runtime libraries, in particular linking with the correct Python runtime. Idea: replace the RuntimeError above with opts.insert(0, "???") and initialize a new build directory with libtbx\configure.py --build-boost-python-extensions=False --build=debug Since you are not building Python extensions, the Python runtime library problem goes away. Ralf

1 0

problem in "space_group_type"
by Koen De Keyser 16 Apr '08

16 Apr '08

Hi, I'm currently in the process of writing some texture analysis/visualization software (XRD pole figures, EBSD) and would like to use cctbx as a basis for the crystallography calculations, as I don't want to reinvent the wheel. Development is in C++ on Windows, and trying one of the examples from the code, I end up being able to crash the program, when running: space_group_type sg_type("C 2"); This occurs using the latest release and still happens in the (unreleased) code from 13th of April, both times, I compiled the library from source, using Visual Studio 2005 and 2008. It looks to me that the problem is situated in the Hall symbol part of the space group (those structures are indicated as bad_ptr in the debugger, when crashing) The error only occurs when running in Debug mode. In Release mode, everything works fine, suggesting a memory initialization problem. I'm unable to debug the problem myself or get more details, as I can't get the library to compile with debugging information in it (Is there an easy way to do this with the Scons build system? I haven't used Scons before.) thanks, Koen De Keyser -- ir. Koen De Keyser Solid State Sciences Ghent University Krijgslaan 281, S1 9000 Gent, Belgium Tel: +32 9 264 43 54 Fax: +32 9 264 49 96 Koen.DeKeyser(a)UGent.be

1 0

Centos rpm
by Pascal 27 Mar '08

27 Mar '08

Hi, I didn't see any rpms for centos. Do you mind if I'll make one ? Is the spec file for fedora available somewhere ? Regards, Pascal -- Pascal Parois Glasgow university http://centos.parois.net

1 0

Re: [cctbxbb] solvent mask
by Ralf W. Grosse-Kunstleve 18 Oct '07

18 Oct '07

Hi Celine, > Thank you for all those explanation. In fact the problem for me is not > the mean of the recalculated map obtained but the range of it, which is > no longer one. Did you consider using the raw mask with the ones and zeros directly? That's the .data attribute of the bulk_solvent object. It is an integer array, but you can convert it with .data.as_double(). See also the .mask_as_xplor_map() method. The only thing the FFT back and forth does is to introduce the Ewald-sphere truncation, which leads to a smoothing of the original mask. If you want that, a quick-and-dirty approach is to copy the .structure_factors() method, rename it, then replace miller_set.unit_cell().volume() with 1.0 in the computation of scale. I think you still need to divide by the number of grid points. Use the .fft_map() of the structure factors object (miller.array) returned by your modified method to get back a map with truncation effects, but on the same scale is the original mask. > To be more clear, what I am trying to do is to "calculate" data at > different solvent contrasts from a pdb file. So I have a series of > solvent electron densities and I would like to use the complex formula F > = Fprotein in vacuum + solvent density*Fmask where Fmask are the > structure factors from an indicator function that is 1 in the solvent > region and 0 in the protein region (= my cctbx mask). That sounds exactly like what's done by the f_model class, which will refine the "solvent density" for you. I think you can also specify a value. See mmtbx/mmtbx/examples/f_model_manager.py and the Computing Commission Newsletter No. 5 (http://cci.lbl.gov/publications/download/iucrcompcomm_jan2005.pdf). The code in the newsletter article is out of date, but the math is still accurate and examples/f_model_manager.py is the updated code. Fcalc in the newsletter formula (page 10 of the .pdf document) is your "Fprotein in vacuum". ksol is your "solvent density". You can set Uaniso and Bsol to 0 if you want. Ralf

2 1

Re: [cctbxbb] solvent mask
by Ralf W. Grosse-Kunstleve 15 Oct '07

15 Oct '07

Hi Celine, > I have a question about solvent mask in cctbx. I generated a mask using > the mask manager. Then I am calculating the structure factors from the > mask (using mask_manager.structure_factors() ). If I am trying to redo a > map from the calculated structure factors (outside the cctbx and I hope > without any scaling), I get a map with a mean density of zero, That's because the F000 is missing. Without the F000, the mean of an FFT map will always be exactly zero. -- This is a general feature of discrete Fourier transforms. > So I wonder what is exactly this > matrix.col(fft_manager.n_real()).product() doing? That's the number of grid points in the map. The FFT introduces this as a scale factor which we are taking out. -- Again, the number of grid point issue is a general feature of discrete Fourier transforms. > And how do I generate > structure factors for a mask without any scaling? I'm guessing the scaling is not the problem, but just the absence of F000, which you can get as map_of_coeff[0]. I think it will also work if you give an array of Miller indices including (0,0,0). We could add suppress_scaling to the signature of the .structure_factors() method, but the scale is rather arbitrary anyway, so I don't think this will be very useful. Ralf

1 0

coordinates -> space group + matrices
by Kris Andersen 15 Oct '07

15 Oct '07

This is probably a simple (maybe stupid) question, but that never stopped me before: I want to give cctbx a set of primitive vectors and coordinates (cartesian or fractional), and have it return (1) the space group and (2) all the symmetry matrices. For example, given R = [4 0 0 % primitive vectors 0 4 0 0 0 4] x = [0.00 0.0 0.00 % fractional coordinates 0.51 0.51 0.51 0.50 0.50 0.00 0.00 0.50 0.50 0.50 0.00 0.50] I want cctbx to tell me this is space group #160 R3m, that there are 6 symmetry operations, and give me the corresponding matrices. This sounds easy, given everything cctbx can do, so it must be possible. Is there a simple example anywhere? Thanks. -- Kristopher E. Andersen http://www.physics.nau.edu/~andersen/ Northern Arizona University kris.andersen(a)nau.edu Department of Physics and Astronomy (928) 523-7202

3 3

Re: [cctbxbb] coordinates -> space group + matrices
by Ralf W. Grosse-Kunstleve 18 Sep '07

18 Sep '07

Hi Kris, your questions isn't simple and the solution isn't straightforward. The functionality you are looking for isn't implemented in the cctbx. You are right in your assumption that most of what you need is already available, but there are missing pieces. Given a unit cell (your primitive vectors), the cctbx can give you the corresponding highest point-group symmetry (lattice symmetry) in a very robust way, but there is nothing to loop over the compatible space groups and to check if the symmetry operations are compatible with the coordinates. There are some programs out there which do this. I think Ton Spek's PLATON for example. The Superflip program (google) determines the symmetry from a P1 map. There is also J. Appl. Cryst. (2005). 38, 237-238 and J. Appl. Cryst. (1998). 31, 922-928. There may be more. I haven't followed the developments very closely. An idea that has been floating around in my mind for a long time is to use the fast translation function. Make a map given your coordinates, loop over all space group, use the fast translation function to find the origin. This is implemented in phenix.hyss. It should be relatively easy to adapt this implementation for your purpose. You'd get scores (correlation coefficients) instead of yes/no answers. This could be useful to uncover approximate symmetry, but I'm just guessing. If you want to try implementations yourself, I'd be happy to answer specific questions. Ralf ____________________________________________________________________________________ Be a better Globetrotter. Get better travel answers from someone who knows. Yahoo! Answers - Check it out. http://answers.yahoo.com/dir/?link=list&sid=396545469

2 1

Re: [cctbxbb] structure_factors_from_map
by Ralf W. Grosse-Kunstleve 09 Jul '07

09 Jul '07

Hi Philipp, The two FFTs introduce a scale factor equivalent to the number of grid points. Please try adding this code: x = abs(new_array).data() y = abs(very_new_array).data() lc = flex.linear_correlation(x, y) lc.show_summary() lr = flex.linear_regression(x, y) lr.show_summary() The correlation coefficient should be 1.0 exactly, and the linear-regression slope is the same as the "total number of grid points" shown in the output: total number of grid points: 345600 ... is_well_defined: True mean_x: 43.9292797741 mean_y: 15181959.0899 coefficient: 1.0 is_well_defined: True y_intercept: 5.1687812285e-07 slope: 345600.0 Ralf ____________________________________________________________________________________Ready for the edge of your seat? Check out tonight's top picks on Yahoo! TV. http://tv.yahoo.com/

1 0

Re: [cctbxbb] pdb output
by Ralf W. Grosse-Kunstleve 09 Jul '07

09 Jul '07

Hi Celine, There isn't something like a pdb.output object, but there are a few helper functions for writing certain records. E.g. see section 5.6 of this newsletter: http://cci.lbl.gov/publications/download/iucrcompcomm_nov2006.pdf Look for "def format_" in iotbx/iotbx/pdb/__init__.py. Ralf ----- Original Message ---- From: celine Besnard <celine.besnard(a)epfl.ch> To: cctbxbb(a)phenix-online.org Sent: Monday, July 9, 2007 9:02:05 AM Subject: [cctbxbb] pdb output Hello I have a structure object which I want to write a pdb file from. Is there any function that does the contrary of the pdb.input function in the iotbx? I have not been able to find one. Thanks Celine -- Dr. Céline Besnard EPFL - SB - IPMC - LCR2 CH-1015 Lausanne Suisse e-mail : Celine.Besnard(a)epfl.ch Tel : +41(0)21.693.05.41 Fax : +41(0)21.693.05.04 http://lcr.epfl.ch _______________________________________________ cctbxbb mailing list cctbxbb(a)phenix-online.org http://www.phenix-online.org/mailman/listinfo/cctbxbb ____________________________________________________________________________________ Need Mail bonding? Go to the Yahoo! Mail Q&A for great tips from Yahoo! Answers users. http://answers.yahoo.com/dir/?link=list&sid=396546091

1 0

pdb output
by celine Besnard 09 Jul '07

09 Jul '07

Hello I have a structure object which I want to write a pdb file from. Is there any function that does the contrary of the pdb.input function in the iotbx? I have not been able to find one. Thanks Celine -- Dr. Céline Besnard EPFL - SB - IPMC - LCR2 CH-1015 Lausanne Suisse e-mail : Celine.Besnard(a)epfl.ch Tel : +41(0)21.693.05.41 Fax : +41(0)21.693.05.04 http://lcr.epfl.ch

1 0