mailr6040 - /1.3/docs/latex/model-free.tex


Others Months | Index by Date | Thread Index
>>   [Date Prev] [Date Next] [Thread Prev] [Thread Next]

Header


Content

Posted by edward on May 01, 2008 - 16:57:
Author: bugman
Date: Thu May  1 16:49:23 2008
New Revision: 6040

URL: http://svn.gna.org/viewcvs/relax?rev=6040&view=rev
Log:
Converted the model-free.py sample script presented in the model-free chapter 
to the new design.


Modified:
    1.3/docs/latex/model-free.tex

Modified: 1.3/docs/latex/model-free.tex
URL: 
http://svn.gna.org/viewcvs/relax/1.3/docs/latex/model-free.tex?rev=6040&r1=6039&r2=6040&view=diff
==============================================================================
--- 1.3/docs/latex/model-free.tex (original)
+++ 1.3/docs/latex/model-free.tex Thu May  1 16:49:23 2008
@@ -733,45 +733,44 @@
 name = `m4' \\
 run.create(name, `mf') \\
  \\
-\# Nuclei type \\
-nuclei(`N') \\
- \\
 \# Load the sequence. \\
-sequence.read(name, `noe.500.out') \\
+sequence.read(`noe.500.out') \\
  \\
 \# Load the relaxation data. \\
-relax\_data.read(name, `R1', `600', 600.0 * 1e6, `r1.600.out') \\
-relax\_data.read(name, `R2', `600', 600.0 * 1e6, `r2.600.out') \\
-relax\_data.read(name, `NOE', `600', 600.0 * 1e6, `noe.600.out') \\
-relax\_data.read(name, `R1', `500', 500.0 * 1e6, `r1.500.out') \\
-relax\_data.read(name, `R2', `500', 500.0 * 1e6, `r2.500.out') \\
-relax\_data.read(name, `NOE', `500', 500.0 * 1e6, `noe.500.out') \\
+relax\_data.read(`R1', `600', 600.0 * 1e6, `r1.600.out') \\
+relax\_data.read(`R2', `600', 600.0 * 1e6, `r2.600.out') \\
+relax\_data.read(`NOE', `600', 600.0 * 1e6, `noe.600.out') \\
+relax\_data.read(`R1', `500', 500.0 * 1e6, `r1.500.out') \\
+relax\_data.read(`R2', `500', 500.0 * 1e6, `r2.500.out') \\
+relax\_data.read(`NOE', `500', 500.0 * 1e6, `noe.500.out') \\
  \\
 \# Setup other values. \\
-diffusion\_tensor.init(name, 10e-9, fixed=1) \\
-value.set(name, 1.02 * 1e-10, `bond\_length') \\
-value.set(name, -160 * 1e-6, `csa') \\
+diffusion\_tensor.init(10e-9, fixed=True) \\
+value.set(1.02 * 1e-10, `bond\_length') \\
+value.set(-160 * 1e-6, `csa') \\
+value.set(`15N', `heteronucleus') \\
+value.set(`1H', `proton') \\
  \\
 \# Select the model-free model. \\
-model\_free.select\_model(run=name, model=name) \\
+model\_free.select\_model(model=name) \\
  \\
 \# Grid search. \\
-grid\_search(name, inc=11) \\
+grid\_search(inc=11) \\
  \\
 \# Minimise. \\
-minimise(`newton', run=name) \\
+minimise(`newton') \\
  \\
 \# Monte Carlo simulations. \\
-monte\_carlo.setup(name, number=100) \\
-monte\_carlo.create\_data(name) \\
-monte\_carlo.initial\_values(name) \\
-minimise(`newton', run=name) \\
-eliminate(run=name) \\
-monte\_carlo.error\_analysis(name) \\
+monte\_carlo.setup(number=100) \\
+monte\_carlo.create\_data() \\
+monte\_carlo.initial\_values() \\
+minimise(`newton') \\
+eliminate() \\
+monte\_carlo.error\_analysis() \\
  \\
 \# Finish. \\
-results.write(run=name, file=`results', force=1) \\
-state.save(`save', force=1)
+results.write(file=`results', force=True) \\
+state.save(`save', force=True)
 \end{exampleenv}
 
 




Related Messages


Powered by MHonArc, Updated Thu May 01 17:40:10 2008