Re: Is it possible to analyse CPMG experiments with relax?
Hi Troels, Welcome to the relax mailing lists. For now the answer to your question is, unfortunately, no - relax does not officially support relaxation dispersion. The analysis you are running is simple two parameter exponential curve-fitting (http://www.nmr-relax.com/manual/Relaxation_curve_fitting.html). This can be used to find the R2eff or R1rho values if you have measured the full exponential curves, but otherwise you cannot perform a dispersion analysis with this. This may not be of use for you at the moment, but note that relax has unofficial and incomplete support for dispersion analyses (both CPMG-type and R1rho-type data sets). As relax is open source, there are many NMR spectroscopists who have added code to relax (for example see http://gna.org/project/memberlist.php?group=relax). An initial implementation of the relaxation dispersion analysis was added to a relax branch back in 2009 by Sebastian Morin (http://thread.gmane.org/gmane.science.nmr.relax.devel/1728). But as this was not completed at the time, it was never merged back into the relax main line (the source code where official relax releases come from). I have recently restored the branch to a partially working state and added a graphical interface for the analysis - mainly for my own purposes (http://svn.gna.org/viewcvs/relax/branches/). So at some point in the near future relax will be able to perform the analyses you are interested in. As relax is open source, if you are interested and adventurous enough you are most welcome to help in the development. Even if you do not know how to code, there are many other things which can be done. For example calculating the partial derivatives of the analytic solutions to obtain the gradients and Hessians so that with relax you can have access to far more powerful optimisation algorithms than any of the other dispersion software has access to. Or to create test data whereby the solution is know, or to collect the input and output test data from published results. If you have the subversion version control software installed, you can obtain the code by typing either: $ svn co svn://svn.gna.org/svn/relax/branches/relax-disp or: $ svn co http://svn.gna.org/svn/relax/branches/relax-disp If you are more interested in quickly performing the analysis, I would point you to Dr. Flemming Hansen's CATIA program: http://www.biochem.ucl.ac.uk/hansen/catia/ (the old page is http://pound.med.utoronto.ca/~flemming/catia/). This performs numerically integration of the Bloch-McConnell equations, so not the optimisation of the analytic solutions of Meiboom, Richard-Carver, etc. It is also only for CPMG-type data rather than R1rho, whereas the relax branch will handle both. I hope this information helps. Regards, Edward On 30 April 2013 18:40, Troels Emtekær Linnet tlin...@gmail.com wrote: Dear relax users. I am looking into different NMR programs to fit relaxation data for CPMG relaxation dispersion experiments and T1rho. Essentially, I am looking for programs for which can fit functions, which for example nessy provide: http://home.gna.org/nessy/reference.html The Meiboom equation or Richard-Carver equation Nessy is very buggy, and I am looking for a replacement. I should be able to: R2eff = -1.0/time_T2*log(Intensity/averageZero) ncyc_arr=[28, 0, 4, 32, 60, 2, 10, 16, 8, 20, 50, 18, 40, 6, 12, 0, 24] time_T2 = 0.06 second nu = ncyc_arr[i]/time_T2 R2cpmg_slow: tau_cpmg = 1.0/(4*nu) R2eff = R2+ka*(1.0-sin(Domega*tau_cpmg)/(Domega*tau_cpmg)) I have followed the tutorial in the homepage manual: Can relax analyse these kinds of experiments? Should i provide the: relax_fit.relax_time(time to be equal tau_cpmg ? I put in time_T2, even though its wrong. I just wanted to try the program. :-) Script for relaxation curve fitting. # Create the 'rx' data pipe. pipe.create('rx', 'relax_fit') ## Load the backbone amide 15N spins from a PDB file. pdbfile=False if pdbfile: structure.read_pdb(pdbfile) structure.load_spins(spin_id='@N') else: molecule.create(mol_name='protein', mol_type='protein') residue.create(res_num=2, res_name='VAL') spin.create(res_num=2, spin_name='N') residue.create(res_num=3, res_name='PHE') spin.create(res_num=3, spin_name='N') residue.create(res_num=4, res_name='GLY') spin.create(res_num=4, spin_name='N') residue.create(res_num=5, res_name='ARG') spin.create(res_num=5, spin_name='N') residue.create(res_num=6, res_name='CYS') and so on ## Loop over the spectra intensities. Relaxation times should be in seconds. readint=True if readint: spectrum.read_intensities(dir='relax', file='proc_list.txt.0int', spectrum_id='0_0.0', int_method='point sum', heteronuc='N', proton='HN', int_col=3) relax_fit.relax_time(time=0.06, spectrum_id='0_0.0') spectrum.read_intensities(dir='relax',
Re: Missing data for spin system and compatibility in GUI
Hi All, yes there is temperature difference on two different spectrometers. i will re calibrate and possibly re-do the experiments. Many thanks for all the suggestions. Best, Nav On Tue, Apr 30, 2013 at 4:32 PM, Edward d'Auvergne edw...@nmr-relax.comwrote: Hi Nav, Martin is spot on here. The temperature control and temperature calibration has been a topic much discussed on the relax mailing lists. For example, here are some threads where you can find out more information, if you wish: http://thread.gmane.org/gmane.science.nmr.relax.user/83 http://thread.gmane.org/gmane.science.nmr.relax.user/273/focus=274 (this is Chris MacRaild's response to a message by Seb Morin) http://thread.gmane.org/gmane.science.nmr.relax.user/1121 http://thread.gmane.org/gmane.science.nmr.relax.user/1419 http://thread.gmane.org/gmane.science.nmr.relax.user/1368 (this one is from Martin) http://thread.gmane.org/gmane.science.nmr.relax.user/1397 (a continuation of the previous thread) You need to click on all the messages in the thread to follow them. There are many more threads on this area, but I don't have the time to find them all right now. The following message might be of significant help for you: http://thread.gmane.org/gmane.science.nmr.relax.user/1419/focus=1423 Firstly note that if you have a temperature problem, this will not be solved by looking at a single field strength as it is a problem between experiments on the same spectrometer. So even if you perform an analysis with data from a single field strength, the R1 data might be up to 2 degrees warmer or colder than the R2, and the same problem will occur to a different degree (or direction) on the second spectrometer. Any analysis using such data will be meaningless, as this will have a large effect on the diffusion tensor. Considering that the diffusion tensor is the major contributor to liquid state relaxation, the internal dynamics can contribute 20% or less (see the original model-free papers for these numbers), then any internal dynamics will be severely distorted, possibly hidden, and artificial motions will appear. I would recommend you look at this section of the relax manual for more temperature related details: http://www.nmr-relax.com/manual/Temperature_control_calibration.html As for relax supporting an analysis at a single field strength, there is nothing stopping you from performing such analysis. relax is designed with flexibility in mind, so you can perform your model-free analysis any way you can imagine. With the minimisation settings, you can replicate the exact results from Art Palmer's Modelfree, from Dasha, or from Tensor2. However note that I have not written any scripts or GUI to handle this situation, as I have no interest to. So you would need to write the script yourself. I would then recommend looking at the file 'auto_analyses/dauvergne_protocol.py' to get an idea of how to implement a full model-free protocol. It is quite big because of the iterative optimisation of the model-free parameters, then model elimination, then model selection, and finally diffusion tensor optimisation, with convergence testing. These steps, even for single field strength data, need to be iterated until convergence. This can take up to 15 iterations. You will also need to decide how to determine your initial diffusion tensor estimate - and relax can perform this as well. But note that you should be aware of of the problems discovered by Schurr et al., 1994 and Tjandra et al., 1996 of the artificial motions. For example see: http://thread.gmane.org/gmane.science.nmr.relax.user/326/focus=332 Korzhnev's review (I don't have the reference at hand) and my paper (http://dx.doi.org/10.1039/b702202f, this has the Korzhnev reference in it) cover all of the problems you will encounter and hopefully convince you that an analysis of single field strength data would only be useful for perfectly isotropic systems (which is never possible due to water shell differences around the system) which have no significant internal motion (but note that if you see no motion with single field strength data, that does not mean that there is no motion). I hope this helps, Edward On 26 April 2013 14:48, Martin Ballaschk ballas...@fmp-berlin.de wrote: Hi Nav, I had similar problems in the past. The inconsistetncy has to come from somewhere. For me, it was the temperature difference between different magnets. After I figured out that we calibrated our spectrometers the wrong way, I finally got consistent data. Just try to superimpose two high-resolution spectra from your two fields. They should be identical. If you can see peaks that are not 99% on top of each other, I would recommend take a long hard stare at temperature control. How do you control your temperature? I found neat methanol did not work with our spectrometers with cryoprobes, see also [1]. We use
Re: Missing data for spin system and compatibility in GUI
Hi Nav, I wouldn't worry about this. I don't think there is a single person in history who has collected perfect relaxation data the first time. If someone says they have, I definitely would not believe them! Just make sure you handle both temperature compensation (single-scan interleaving and temperature compensation blocks) and temperature control (checking every experiment on every spectrometer with MeOH/Ethylene glycol). Unfortunately the paper from the Gooley lab describing this is not published yet. But all the details are described in the relax manual, so maybe the relax papers might be a possible citation ;) Regards, Edward On 2 May 2013 11:13, Navratna Vajpai navratna.vaj...@gmail.com wrote: Hi All, yes there is temperature difference on two different spectrometers. i will re calibrate and possibly re-do the experiments. Many thanks for all the suggestions. Best, Nav On Tue, Apr 30, 2013 at 4:32 PM, Edward d'Auvergne edw...@nmr-relax.com wrote: Hi Nav, Martin is spot on here. The temperature control and temperature calibration has been a topic much discussed on the relax mailing lists. For example, here are some threads where you can find out more information, if you wish: http://thread.gmane.org/gmane.science.nmr.relax.user/83 http://thread.gmane.org/gmane.science.nmr.relax.user/273/focus=274 (this is Chris MacRaild's response to a message by Seb Morin) http://thread.gmane.org/gmane.science.nmr.relax.user/1121 http://thread.gmane.org/gmane.science.nmr.relax.user/1419 http://thread.gmane.org/gmane.science.nmr.relax.user/1368 (this one is from Martin) http://thread.gmane.org/gmane.science.nmr.relax.user/1397 (a continuation of the previous thread) You need to click on all the messages in the thread to follow them. There are many more threads on this area, but I don't have the time to find them all right now. The following message might be of significant help for you: http://thread.gmane.org/gmane.science.nmr.relax.user/1419/focus=1423 Firstly note that if you have a temperature problem, this will not be solved by looking at a single field strength as it is a problem between experiments on the same spectrometer. So even if you perform an analysis with data from a single field strength, the R1 data might be up to 2 degrees warmer or colder than the R2, and the same problem will occur to a different degree (or direction) on the second spectrometer. Any analysis using such data will be meaningless, as this will have a large effect on the diffusion tensor. Considering that the diffusion tensor is the major contributor to liquid state relaxation, the internal dynamics can contribute 20% or less (see the original model-free papers for these numbers), then any internal dynamics will be severely distorted, possibly hidden, and artificial motions will appear. I would recommend you look at this section of the relax manual for more temperature related details: http://www.nmr-relax.com/manual/Temperature_control_calibration.html As for relax supporting an analysis at a single field strength, there is nothing stopping you from performing such analysis. relax is designed with flexibility in mind, so you can perform your model-free analysis any way you can imagine. With the minimisation settings, you can replicate the exact results from Art Palmer's Modelfree, from Dasha, or from Tensor2. However note that I have not written any scripts or GUI to handle this situation, as I have no interest to. So you would need to write the script yourself. I would then recommend looking at the file 'auto_analyses/dauvergne_protocol.py' to get an idea of how to implement a full model-free protocol. It is quite big because of the iterative optimisation of the model-free parameters, then model elimination, then model selection, and finally diffusion tensor optimisation, with convergence testing. These steps, even for single field strength data, need to be iterated until convergence. This can take up to 15 iterations. You will also need to decide how to determine your initial diffusion tensor estimate - and relax can perform this as well. But note that you should be aware of of the problems discovered by Schurr et al., 1994 and Tjandra et al., 1996 of the artificial motions. For example see: http://thread.gmane.org/gmane.science.nmr.relax.user/326/focus=332 Korzhnev's review (I don't have the reference at hand) and my paper (http://dx.doi.org/10.1039/b702202f, this has the Korzhnev reference in it) cover all of the problems you will encounter and hopefully convince you that an analysis of single field strength data would only be useful for perfectly isotropic systems (which is never possible due to water shell differences around the system) which have no significant internal motion (but note that if you see no motion with single field strength data, that does not mean that there is no motion). I hope this helps,
Re: Is it possible to analyse CPMG experiments with relax?
Dear Edward. Thank you for your generous email, which helped a-lot. I am happy to see the active development, and I would be more than happy to join in. I am quite good in python programming, and are confident i revision programs as svn and git. And I have courses in scientific computing, so I think i get along quite good. My reason for my interest, is that I think I should change my working habits, to something more effect full. My work-flow at the moment, is this. 1) CPMG/T1-rho experiment acquisition with NLS, through VnmrJ. 2) Data reconstruction in qMDD. (3) Main peak positioning in CcpNmr Analysis.) 4) Small peak adjustment, control in SPARKY. 5) Point sum integration in with: seriesTabhttp://spin.niddk.nih.gov/NMRPipe/ref/prog/seriestab.htmlwith: -dx 1 -dy 1 6) Integration analysis in gnuplot/IgorPro,Originlab. The use of IgorPro,Originlab have been used because of easy use of the global fitting routine, but pose a problem, since we only have a very few licenses. And I weigh open-source very high. :-) The last weeks, I have fiddling around my workflow to try something like this: 1) CPMG/T1-rho experiment acquisition with NLS, through VnmrJ. 2) Data reconstruction in qMDD. (3) Main peak positioning in CcpNmr Analysis.) 4) Small peak adjustment, control in SPARKY. 5) Point sum integration in nmrglue,http://springer.r.delivery.net/r/r?2.1.Ee.2Tp.1jiDvx.C5t4Ay..H.a%2ace.44Wy.bW89MQ%5f%5fDIJWFRd0and easy visualization of each integrationhttps://github.com/jjhelmus/nmrglue/wiki/Plot-and-integrate-peaks-from-nmrpipe-format-and-sparky-list. Preparation of data for fitting. 6) Global fitting. Either with nessy, relax or with python scipy leastsq. Here I tried to make a nessy databasehttps://github.com/jjhelmus/nmrglue/wiki/Make-a-nessy-database-for-CPMG-eksperiment, but nessy came out very buggy. And I was about to set out for some python scipy fitting after nmrglue. But I had a hard time imagining that NMR software were not already developed for this, and I was very pleased to see the development of relax, which have not come to my attention before. And especially the inclusion of the python interpreter, and possibility to write scripts, is genius. Which is similar to the where the power of pymol is shining through. An optimal workflow would be this. 1) CPMG/T1-rho experiment acquisition with NLS, through VnmrJ. 2) Data reconstruction in qMDD. (3) Main peak positioning in CcpNmr Analysis.) 4) Small peak adjustment, control in SPARKY. 5) Point sum integration in nmrglue, and easy visualization of each integrationhttps://github.com/jjhelmus/nmrglue/wiki/Plot-and-integrate-peaks-from-nmrpipe-format-and-sparky-list. Preparation of data for fitting. 6) nmrgluehttp://springer.r.delivery.net/r/r?2.1.Ee.2Tp.1jiDvx.C5t4Ay..H.a%2ace.44Wy.bW89MQ%5f%5fDIJWFRd0script preparation for relax, execution, save of result. 7) Graphical inspection of result through relax. So, I think I will try out the current state of the relax-disp branch. Are there any others developing on this branch? Best Troels -- Troels Emtekær Linnet PhD student Copenhagen University SBiNLab, 3-0-41 Ole Maaloes Vej 5 2200 Copenhagen N Tlf: +45 353-22083 2013/5/2 Edward d'Auvergne edw...@nmr-relax.com Hi Troels, Welcome to the relax mailing lists. For now the answer to your question is, unfortunately, no - relax does not officially support relaxation dispersion. The analysis you are running is simple two parameter exponential curve-fitting (http://www.nmr-relax.com/manual/Relaxation_curve_fitting.html). This can be used to find the R2eff or R1rho values if you have measured the full exponential curves, but otherwise you cannot perform a dispersion analysis with this. This may not be of use for you at the moment, but note that relax has unofficial and incomplete support for dispersion analyses (both CPMG-type and R1rho-type data sets). As relax is open source, there are many NMR spectroscopists who have added code to relax (for example see http://gna.org/project/memberlist.php?group=relax). An initial implementation of the relaxation dispersion analysis was added to a relax branch back in 2009 by Sebastian Morin (http://thread.gmane.org/gmane.science.nmr.relax.devel/1728). But as this was not completed at the time, it was never merged back into the relax main line (the source code where official relax releases come from). I have recently restored the branch to a partially working state and added a graphical interface for the analysis - mainly for my own purposes (http://svn.gna.org/viewcvs/relax/branches/). So at some point in the near future relax will be able to perform the analyses you are interested in. As relax is open source, if you are interested and adventurous enough you are most welcome to help in the development. Even if you do not know how to code, there are many other things which can be done. For example calculating the partial derivatives of the analytic solutions to obtain the
Building relax from source using scons and a different python interpretor
Dear Relax users. I am trying to build relax, from source, using scons. We have the python enthought distribution installed, to solve most of our dependencies, and we install packages here, so each computer can use this shared installation. How I can I tell scons, to use python2.7 to build relax? I think it tries to use my local python installation, but here I have not installed minfx, ,wxPython. [tlinnet@tomat ~/software]$ cd relax-trunk/ [tlinnet@tomat relax-trunk]$ scons scons: Reading SConscript files ... The dependency 'minfx' has not been installed (see https://gna.org/projects/minfx/). [tlinnet@tomat relax-trunk]$ python2.7 Enthought Python Distribution -- www.enthought.com Version: 7.3-2 (64-bit) Python 2.7.3 |EPD 7.3-2 (64-bit)| (default, Apr 11 2012, 17:52:16) [GCC 4.1.2 20080704 (Red Hat 4.1.2-44)] on linux2 Type credits, demo or enthought for more information. import minfx import readline import wx import numpy import scipy Best Troels Troels Emtekær Linnet Ved kløvermarken 9, 1.th 2300 København S Mobil: +45 60210234 ___ relax (http://www.nmr-relax.com) This is the relax-users mailing list relax-users@gna.org To unsubscribe from this list, get a password reminder, or change your subscription options, visit the list information page at https://mail.gna.org/listinfo/relax-users
Re: Is it possible to analyse CPMG experiments with relax?
Hi, Its great that you have an interest and know Python - you are in a perfect position to join as a relax developer! See below for more responses: Thank you for your generous email, which helped a-lot. You're welcome! I am happy to see the active development, and I would be more than happy to join in. I am quite good in python programming, and are confident i revision programs as svn and git. And I have courses in scientific computing, so I think i get along quite good. I would recommend you have a read of the relax open source infrastructure chapter of the relax manual (http://www.nmr-relax.com/manual/Open_source_infrastructure.html) and, more importantly, the development chapter of the relax manual (http://www.nmr-relax.com/manual/relax_development.html). The PDF version of the manual is much easier to read (http://download.gna.org/relax/manual/relax.pdf). These chapters describe in full detail everything you would ever need as a relax developer. Note that relax is a very mature project, so learning how to code in such an environment to avoid breaking the rest of the program will give you quite a different skill set. You might also be interested in learning about the minfx project that relax uses for optimisation (https://gna.org/projects/minfx/). This originated as a relax package as the scipy optimisers all contained fatal bugs back in 2003 (I'm not sure they have been fixed as the original developers were MIA even back then and I think have never returned). But it was spun out into its own software distribution. My reason for my interest, is that I think I should change my working habits, to something more effect full. My work-flow at the moment, is this. 1) CPMG/T1-rho experiment acquisition with NLS, through VnmrJ. I have to warn you here that non-linear sampling is notoriously bad for measuring high precision NMR parameters such as relaxation data. I would recommend avoiding this technique if you can. It is great for low precision data required for assignment, for example, but not so good for the high precision data measurements. 2) Data reconstruction in qMDD. (3) Main peak positioning in CcpNmr Analysis.) 4) Small peak adjustment, control in SPARKY. 5) Point sum integration in with: seriesTab with: -dx 1 -dy 1 6) Integration analysis in gnuplot/IgorPro,Originlab. The use of IgorPro,Originlab have been used because of easy use of the global fitting routine, but pose a problem, since we only have a very few licenses. And I weigh open-source very high. :-) The way I perform this is a bit different in that I use peak heights directly from Sparky. The last weeks, I have fiddling around my workflow to try something like this: 1) CPMG/T1-rho experiment acquisition with NLS, through VnmrJ. 2) Data reconstruction in qMDD. (3) Main peak positioning in CcpNmr Analysis.) 4) Small peak adjustment, control in SPARKY. 5) Point sum integration in nmrglue, and easy visualization of each integration. Preparation of data for fitting. 6) Global fitting. Either with nessy, relax or with python scipy leastsq. Here I tried to make a nessy database, but nessy came out very buggy. And I was about to set out for some python scipy fitting after nmrglue. I am currently the maintainer of the NESSY project, but the types of bugs reported require significant amounts of coding to solve the problems. Unfortunately I don't have the time for this - it could be a few months of work. As for nmrglue (http://code.google.com/p/nmrglue/), this appears quite new and this is the first time I have heard of it. It looks like an interesting project. I wonder if they use 3-point quadratic integration for determining the maximum peak height? But I had a hard time imagining that NMR software were not already developed for this, and I was very pleased to see the development of relax, which have not come to my attention before. And especially the inclusion of the python interpreter, and possibility to write scripts, is genius. Which is similar to the where the power of pymol is shining through. Thanks! I am directing the development of relax to have the maximum amount of flexibility. For the basic users who want quick results, there are the auto-analyses which can be used as blackboxes, giving the user the best practice analysis. These are used in the GUI. For the medium level users, the user functions (which are special Python functions which perform a lot of checking of the user input) allow for advanced scripting. For the advanced users, the relax API can be used to build complete new analyses (http://www.nmr-relax.com/api/). I have been developing relax so that in the future it can be used by NMR users as a replacement for Matlab/Mathematica for numerical operations. The relax library - the 'lib' package - is a large collection of NMR specific functions. For example, have a look at the rotation matrix module 'lib.geometry.rotations'
Re: Building relax from source using scons and a different python interpretor
Hi Troels, In the future, for such questions I would recommend asking on the relax development mailing list (relax-devel att gna.org). The answer is complicated as SCons has problems with multiple Python installations. You have to have scons installed for the Python version you are targeting and no other SCons versions present. I have managed to get it to work, but the scons must be run with the correct Python version and also then separately find the correct Python version for compilations (these are unfortunately not coupled in SCons). There is one cheating way and that is to use the 'devel_scripts/manual_c_module.py' script, just modifying it to point your Python version. This script is mainly to get around these bad design issues in SCons for when multiple Python versions are present, as I have lots of different Python versions compiled into a special directory for testing relax (every version from 1.0 to 3.3) and SCons really cannot handle this. As for using Enthought Python, I don't know how well relax will run in such an environment. It'll be interesting to see what happens. Regards, Edward On 2 May 2013 17:26, Troels Emtekær Linnet tlin...@gmail.com wrote: Dear Relax users. I am trying to build relax, from source, using scons. We have the python enthought distribution installed, to solve most of our dependencies, and we install packages here, so each computer can use this shared installation. How I can I tell scons, to use python2.7 to build relax? I think it tries to use my local python installation, but here I have not installed minfx, ,wxPython. [tlinnet@tomat ~/software]$ cd relax-trunk/ [tlinnet@tomat relax-trunk]$ scons scons: Reading SConscript files ... The dependency 'minfx' has not been installed (see https://gna.org/projects/minfx/). [tlinnet@tomat relax-trunk]$ python2.7 Enthought Python Distribution -- www.enthought.com Version: 7.3-2 (64-bit) Python 2.7.3 |EPD 7.3-2 (64-bit)| (default, Apr 11 2012, 17:52:16) [GCC 4.1.2 20080704 (Red Hat 4.1.2-44)] on linux2 Type credits, demo or enthought for more information. import minfx import readline import wx import numpy import scipy Best Troels Troels Emtekær Linnet Ved kløvermarken 9, 1.th 2300 København S Mobil: +45 60210234 ___ relax (http://www.nmr-relax.com) This is the relax-users mailing list relax-users@gna.org To unsubscribe from this list, get a password reminder, or change your subscription options, visit the list information page at https://mail.gna.org/listinfo/relax-users ___ relax (http://www.nmr-relax.com) This is the relax-users mailing list relax-users@gna.org To unsubscribe from this list, get a password reminder, or change your subscription options, visit the list information page at https://mail.gna.org/listinfo/relax-users
Re: Building relax from source using scons and a different python interpretor
Heh, that was actually easier than thought. The Enthought distribution came along with scons installed. /software/python-enthought-dis/epd-7.3-2-rh5-x86_64/bin/scons So no problem building. And then making this shell file, to use the enthought: relax2 #!/software/python-enthought-dis/epd-7.3-2-rh5-x86_64/bin/python2.7 import relax relax.start(profile_flag=False) Best Troels Troels Emtekær Linnet Ved kløvermarken 9, 1.th 2300 København S Mobil: +45 60210234 2013/5/2 Edward d'Auvergne edw...@nmr-relax.com Hi Troels, In the future, for such questions I would recommend asking on the relax development mailing list (relax-devel att gna.org). The answer is complicated as SCons has problems with multiple Python installations. You have to have scons installed for the Python version you are targeting and no other SCons versions present. I have managed to get it to work, but the scons must be run with the correct Python version and also then separately find the correct Python version for compilations (these are unfortunately not coupled in SCons). There is one cheating way and that is to use the 'devel_scripts/manual_c_module.py' script, just modifying it to point your Python version. This script is mainly to get around these bad design issues in SCons for when multiple Python versions are present, as I have lots of different Python versions compiled into a special directory for testing relax (every version from 1.0 to 3.3) and SCons really cannot handle this. As for using Enthought Python, I don't know how well relax will run in such an environment. It'll be interesting to see what happens. Regards, Edward On 2 May 2013 17:26, Troels Emtekær Linnet tlin...@gmail.com wrote: Dear Relax users. I am trying to build relax, from source, using scons. We have the python enthought distribution installed, to solve most of our dependencies, and we install packages here, so each computer can use this shared installation. How I can I tell scons, to use python2.7 to build relax? I think it tries to use my local python installation, but here I have not installed minfx, ,wxPython. [tlinnet@tomat ~/software]$ cd relax-trunk/ [tlinnet@tomat relax-trunk]$ scons scons: Reading SConscript files ... The dependency 'minfx' has not been installed (see https://gna.org/projects/minfx/). [tlinnet@tomat relax-trunk]$ python2.7 Enthought Python Distribution -- www.enthought.com Version: 7.3-2 (64-bit) Python 2.7.3 |EPD 7.3-2 (64-bit)| (default, Apr 11 2012, 17:52:16) [GCC 4.1.2 20080704 (Red Hat 4.1.2-44)] on linux2 Type credits, demo or enthought for more information. import minfx import readline import wx import numpy import scipy Best Troels Troels Emtekær Linnet Ved kløvermarken 9, 1.th 2300 København S Mobil: +45 60210234 ___ relax (http://www.nmr-relax.com) This is the relax-users mailing list relax-users@gna.org To unsubscribe from this list, get a password reminder, or change your subscription options, visit the list information page at https://mail.gna.org/listinfo/relax-users ___ relax (http://www.nmr-relax.com) This is the relax-users mailing list relax-users@gna.org To unsubscribe from this list, get a password reminder, or change your subscription options, visit the list information page at https://mail.gna.org/listinfo/relax-users