Re: [gmx-users] Re:g_clustsize problem

2006-03-13 Thread David van der Spoel

Wang Zhun wrote:

Hi, David,
The command I use is
g_cluster -f traj.trr -s topol.tpr -cutoff 0.2 -minstruct 2 -o rmsd-clust.xpm -g cluster.log 
 choose group 1  Protein

This is fine, results look OK. So what's the problem?


and
g_clustsize -f traj.trr -o csize.xpm -cut 0.2  -n index.ndx
 choose group 1  Protein

This is not meaningful.


the following is part of cluster log which says 12 clusters found rather than 1 
cluster:
Using linkage method for clustering
Using RMSD cutoff 0.2 nm
The RMSD ranges from 0.0738896 to 2.25219 nm
Average RMSD is 0.61356
Number of structures for matrix 4001
Energy of the matrix is 1004.6 nm

Found 12 clusters


cl. | #st rmsd | middle rmsd | cluster members
  1 |   1  |  0  |  0
  2 |   8 .612 |  5 .271 |  5 10 15 20 25 30 35
|  | | 40
  3 |   21.084 | 45 .169 | 45 50
  4 |  141.636 | 55 .292 | 55 60 65 70 75 80 85
|  | | 90 95100105110115120
  5 |   32.092 |125 .197 |125130135
  6 |  132.640 |140 .292 |140145150155160165170
|  | |175180185190195200
  7 |  313.199 |205 .433 |205210215220225230235
|  | |240245250255260265270
|  | |275280285290295300305
|  | |310315320325330335340
|  | |345350355
  8 |   23.572 |360 .143 |360365
  9 |   74.119 |370 .264 |370375380385390395400
 10 |  124.633 |405 .283 |405410415420425430435
|  | |440445450455460
 11 |  325.164 |465 .356 |465470475480485490495
|  | |500505510515520525530
|  | |535540545550555560565
|  | |570575580585590595600
|  | |605610615620
 12 | 38765.784 |810 .806 |625630635640645650655
|  | |660665670675680685690
|  | |695700705710715720725
|  | |730735740745750755760
|  | |765770775780785790795
|  | |800805810815820825830
|  | |835840845850855860865
|  | |870875880885890895900
|  | |905910915920925930935
|  | |940945950955960965970
//



___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php



--
David.

David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,  75124 Uppsala, Sweden
phone:  46 18 471 4205  fax: 46 18 511 755
[EMAIL PROTECTED]   [EMAIL PROTECTED]   http://folding.bmc.uu.se

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: [gmx-users] Range checking error

2006-03-13 Thread David van der Spoel

Gaurav Chopra wrote:

Hi

I am doing a steepest descent for the initial mdrun for a peptide and it 
gives me the following error. Please advice.


Thanks
Gaurav
---
Program mdrun, VERSION 3.3
Source code file: nsgrid.c, line: 226

Range checking error:
Explanation: During neighborsearching, we assign each particle to a grid
based on its coordinates. If your system contains collisions or parameter
errors that give particles very high velocities you might end up with some
coordinates being +-Infinity or NaN (not-a-number). Obviously, we cannot
put these on a grid, so this is usually where we detect those errors.
Make sure your system is properly energy-minimized and that the potential
energy seems reasonable before trying again.

Please follow the advice in the above two sentences. Check your structure.



Variable ci has value -2147483648. It should have been within [ 0 .. 2000 ]
Please report this to the mailing list (gmx-users@gromacs.org)
---

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the www 
interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php



--
David.

David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,  75124 Uppsala, Sweden
phone:  46 18 471 4205  fax: 46 18 511 755
[EMAIL PROTECTED]   [EMAIL PROTECTED]   http://folding.bmc.uu.se

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: [gmx-users] xtc files

2006-03-13 Thread David van der Spoel

YOLANDA SMALL wrote:

Hi All,

How do I convert xtc files exported for a subset of atoms in my system to pdb
files?

When I run trjconv, I get the following error message:

Fatal error: Index[0] 3611 is larger than the number of atoms in the trajectory
file (3152)

In my mdp file I specify the option 'xtc_grps' to export a subset of solvent and
protein atoms to an xtc file. This subset of atoms is defined in my index file
and contains 3152 atoms. Then I discard the .trr files to save space.  It looks
like the atom index gets renumbered from 1 in the new .xtc file so that the new
atom indeces no longer correspond to the atoms in the original coordinate file.
 How do I recover the coordinates with matching atom names from the .xtc file?


Sorry I misread your question.
You make a pdb file from your original tpr using the same index group 
you used for the xtc-grps, then you use this pdb file as reference for 
trjconv on the xtc to pdb conversion.

Thanks,
Yolanda

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php



--
David.

David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,  75124 Uppsala, Sweden
phone:  46 18 471 4205  fax: 46 18 511 755
[EMAIL PROTECTED]   [EMAIL PROTECTED]   http://folding.bmc.uu.se

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: [gmx-users] xtc files

2006-03-13 Thread David van der Spoel

YOLANDA SMALL wrote:

Hi All,

How do I convert xtc files exported for a subset of atoms in my system to pdb
files?

When I run trjconv, I get the following error message:

Fatal error: Index[0] 3611 is larger than the number of atoms in the trajectory
file (3152)

In my mdp file I specify the option 'xtc_grps' to export a subset of solvent and
protein atoms to an xtc file. This subset of atoms is defined in my index file
and contains 3152 atoms. Then I discard the .trr files to save space.  It looks
like the atom index gets renumbered from 1 in the new .xtc file so that the new
atom indeces no longer correspond to the atoms in the original coordinate file.
 How do I recover the coordinates with matching atom names from the .xtc file?


you need to make an inverse index files.
First you had a selection
1   --> 1
3   --> 2
5   --> 3
7   --> 4
11  --> 5
13  --> 6

Now you have to generate the inverse of your xtc-grp



Thanks,
Yolanda

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php



--
David.

David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,  75124 Uppsala, Sweden
phone:  46 18 471 4205  fax: 46 18 511 755
[EMAIL PROTECTED]   [EMAIL PROTECTED]   http://folding.bmc.uu.se

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] keep atoms planar

2006-03-13 Thread Rongliang Wu
hello, 

i found in ffoplsaa.itp :

[ dihedraltypes ]
; Improper OPLS dihedrals to keep groups planar.
; (OPLS doesnt use impropers for chiral atoms).
; Since these functions are periodic of the form 1-cos(2*x), the are actually
; implemented as proper dihedrals [1+cos(2*x+180)] for the moment, 
; to keep things compatible.
; The defines are used in ffoplsaa.rtp or directly in your .top ile.

; O?-C -X -Y improper torsion. C can be C_2 or C_3 too.
#define improper_O_C_X_Y180.0 43.93200   2

; X-NO-ON-NO improper torsion.
#define improper_X_NO_ON_NO 180.0 43.93200   2

; N2-X-N2-N2 improper torsion.
#define improper_N2_X_N2_N2 180.0 43.93200   2

; Z -N?-X -Y improper torsion
#define improper_Z_N_X_Y180.0  4.18400   2

; Z -CM-X -Y improper torsion. CM can be C= too.
#define improper_Z_CM_X_Y   180.0 62.76000   2

; Z -CA-X -Y improper torsion. CA is any ring carbon 
(CA,CB,CN,CV,CW,CR,CK,CQ,CS,C*)
#define improper_Z_CA_X_Y   180.0  4.60240   2


are X Y Z designate any atom types? 
--
Rongliang Wu
2006-03-14


___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] Re:g_clustsize problem

2006-03-13 Thread Wang Zhun
Hi, David,
The command I use is
g_cluster -f traj.trr -s topol.tpr -cutoff 0.2 -minstruct 2 -o rmsd-clust.xpm 
-g cluster.log 
 choose group 1  Protein
and
g_clustsize -f traj.trr -o csize.xpm -cut 0.2  -n index.ndx
 choose group 1  Protein
the following is part of cluster log which says 12 clusters found rather than 1 
cluster:
Using linkage method for clustering
Using RMSD cutoff 0.2 nm
The RMSD ranges from 0.0738896 to 2.25219 nm
Average RMSD is 0.61356
Number of structures for matrix 4001
Energy of the matrix is 1004.6 nm

Found 12 clusters


cl. | #st rmsd | middle rmsd | cluster members
  1 |   1  |  0  |  0
  2 |   8 .612 |  5 .271 |  5 10 15 20 25 30 35
|  | | 40
  3 |   21.084 | 45 .169 | 45 50
  4 |  141.636 | 55 .292 | 55 60 65 70 75 80 85
|  | | 90 95100105110115120
  5 |   32.092 |125 .197 |125130135
  6 |  132.640 |140 .292 |140145150155160165170
|  | |175180185190195200
  7 |  313.199 |205 .433 |205210215220225230235
|  | |240245250255260265270
|  | |275280285290295300305
|  | |310315320325330335340
|  | |345350355
  8 |   23.572 |360 .143 |360365
  9 |   74.119 |370 .264 |370375380385390395400
 10 |  124.633 |405 .283 |405410415420425430435
|  | |440445450455460
 11 |  325.164 |465 .356 |465470475480485490495
|  | |500505510515520525530
|  | |535540545550555560565
|  | |570575580585590595600
|  | |605610615620
 12 | 38765.784 |810 .806 |625630635640645650655
|  | |660665670675680685690
|  | |695700705710715720725
|  | |730735740745750755760
|  | |765770775780785790795
|  | |800805810815820825830
|  | |835840845850855860865
|  | |870875880885890895900
|  | |905910915920925930935
|  | |940945950955960965970
//



___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: [gmx-users] xtc files

2006-03-13 Thread Yang Ye

YOLANDA SMALL wrote:

Hi All,

How do I convert xtc files exported for a subset of atoms in my system to pdb
files?

When I run trjconv, I get the following error message:

Fatal error: Index[0] 3611 is larger than the number of atoms in the trajectory
file (3152)

In my mdp file I specify the option 'xtc_grps' to export a subset of solvent and
protein atoms to an xtc file. This subset of atoms is defined in my index file
and contains 3152 atoms. Then I discard the .trr files to save space.  It looks
like the atom index gets renumbered from 1 in the new .xtc file so that the new
atom indeces no longer correspond to the atoms in the original coordinate file.
  

Did you use -shuffle -sort? Even in such case, the renumbering won't
take number larger  than the total number of atoms. You may have
selected the group in the index file which includes all the atoms.
More information from you is appreciated.

P.S.: Use gmxcheck -f to check how many atoms inside the xtc.

 How do I recover the coordinates with matching atom names from the .xtc file?

Thanks,
Yolanda

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php

  




___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] xtc files

2006-03-13 Thread YOLANDA SMALL
Hi All,

How do I convert xtc files exported for a subset of atoms in my system to pdb
files?

When I run trjconv, I get the following error message:

Fatal error: Index[0] 3611 is larger than the number of atoms in the trajectory
file (3152)

In my mdp file I specify the option 'xtc_grps' to export a subset of solvent and
protein atoms to an xtc file. This subset of atoms is defined in my index file
and contains 3152 atoms. Then I discard the .trr files to save space.  It looks
like the atom index gets renumbered from 1 in the new .xtc file so that the new
atom indeces no longer correspond to the atoms in the original coordinate file.
 How do I recover the coordinates with matching atom names from the .xtc file?

Thanks,
Yolanda

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] Rychaert-Belleman's dihedral

2006-03-13 Thread Rongliang Wu
hello all,
 i'm using the opls force field to simulate a system of molecules with benzene 
ring, i found the default dihedral fun is 3(Rychaert-Belleman's dihedral), and 
do i have to add those dihedrals in the benzene to the [ impropers ] so as to 
insure the atoms are in a plane? or can the "fun 3" keep the stoms of the 
benzene within the plane? is it necessary to use Rychaert-Belleman's dihedral 
for the oplsaa force filed?

Regards! 

--
Rongliang Wu
2006-03-14


___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] Range checking error

2006-03-13 Thread Gaurav Chopra

Hi

I am doing a steepest descent for the initial mdrun for a peptide and it 
gives me the following error. Please advice.


Thanks
Gaurav
---
Program mdrun, VERSION 3.3
Source code file: nsgrid.c, line: 226

Range checking error:
Explanation: During neighborsearching, we assign each particle to a grid
based on its coordinates. If your system contains collisions or parameter
errors that give particles very high velocities you might end up with some
coordinates being +-Infinity or NaN (not-a-number). Obviously, we cannot
put these on a grid, so this is usually where we detect those errors.
Make sure your system is properly energy-minimized and that the potential
energy seems reasonable before trying again.

Variable ci has value -2147483648. It should have been within [ 0 .. 2000 ]
Please report this to the mailing list (gmx-users@gromacs.org)
---

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: [gmx-users] Why my protein always falls out of the water box?

2006-03-13 Thread Jochen Hub

Shanjie Huang wrote:


Hi dear friends,
I performed a molecular dynamics of a ligand-receptor complex for 5ns;
but I found that my complex always fall out from the water box quickly
in the first 1ns. What may be the problem? 
Can I constrain some residues which are far away from the active site to

aviod this problem? But I am worry about the extra constrain would
interfere my MD result, for in the normal situation there should be no
constrains. 


Thanks.


When using

comm-grps= Protein

in your mdp file the center of mass of you protein will stay at its 
starting position.


good luck,
Jochen

--

Jochen Hub
Max Planck Institute for Biophysical Chemistry
Computational biomolecular dynamics group
Am Fassberg 11
D-37077 Goettingen, Germany
Email: jhub[at]gwdg.de



___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: [gmx-users] Why my protein always falls out of the water box?

2006-03-13 Thread MGiò
Hi you!it may not be as bad as it looks like. Usually the one happening to your protein is only a visualization problem, if an atom of your protein crosses the periodic simulation box, it might appear as if the whole molecule jumps out of the box, because the coordinates of one of its periodic images are shown instead of those referring to the real box. You can simply assess if this is the case by running a single step dynamics, at a low temperature, if the protein jumps out all the same, it can't be because of the real protein movement! In this case, you can just convert you trajectory by 
trjconv -center -pbc wholethis way you should have a trajectory in which your protein quietly stays in its box.hope it helps!greetings,MGiòOn 3/13/06, 
Shanjie Huang <[EMAIL PROTECTED]> wrote:
Hi dear friends,I performed a molecular dynamics of a ligand-receptor complex for 5ns;but I found that my complex always fall out from the water box quicklyin the first 1ns. What may be the problem?Can I constrain some residues which are far away from the active site to
aviod this problem? But I am worry about the extra constrain wouldinterfere my MD result, for in the normal situation there should be noconstrains.Thanks.--  Shanjie Huang  
[EMAIL PROTECTED]--http://www.fastmail.fm - A no graphics, no pop-ups email service___gmx-users mailing list
gmx-users@gromacs.orghttp://www.gromacs.org/mailman/listinfo/gmx-usersPlease don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED].Can't post? Read http://www.gromacs.org/mailing_lists/users.php

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] Why my protein always falls out of the water box?

2006-03-13 Thread Shanjie Huang
Hi dear friends,
I performed a molecular dynamics of a ligand-receptor complex for 5ns;
but I found that my complex always fall out from the water box quickly
in the first 1ns. What may be the problem? 
Can I constrain some residues which are far away from the active site to
aviod this problem? But I am worry about the extra constrain would
interfere my MD result, for in the normal situation there should be no
constrains. 

Thanks.

 
-- 
  Shanjie Huang
  [EMAIL PROTECTED]

-- 
http://www.fastmail.fm - A no graphics, no pop-ups email service

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: [gmx-users] can not write a gro file without atom names

2006-03-13 Thread Jochen Hub
Do you want to create a really short trajectory with only two frames, or 
do you want to create a structure with both structures given in your gro 
files in it?


In the latter case, turn you gro files into pdb (editconf), use "cat 
a.pdb b.pdb > c.pdb" and remove the ENDMDL and TER lines by hand.


Jochen

Viswanadham Sridhara wrote:



Hi everyone,
I am trying to concatenate two .gro files obtained from two separate 
trajcetories.

It says "can not write a gro file without atom names" with trjcat command.
Any suggestions or any more info you required??
Thanks,
-Viswam.
--
Viswanadham Sridhara,
Graduate Research Assistant,
Old Dominion University, "VIRGINIA".



___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php




--

Jochen Hub
Max Planck Institute for Biophysical Chemistry
Computational biomolecular dynamics group
Am Fassberg 11
D-37077 Goettingen, Germany
Email: jhub[at]gwdg.de



___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] Re: gmx-users Digest, Vol 23, Issue 34

2006-03-13 Thread Maxim Fedorov

> 
> > 1. Does anybody have a 'ready-to-run' topology
> > for SW water (the polarisable water model from J.Phys. Chem. 
> > B (2001), 105, 2618? 
> > As I see the sw.itp in the standart distribution
> > has improper order of directives. Am I right or it is just 
> > the Friday night and I should take a rest for a while...? 
> 
> http://www.gromacs.org/pipermail/gmx-users/2003-December/008346.html
> 
Great!
Thanks a lot, Dallas.
I was almost sure that it should be something like this:  a simple but
not-so-obvious mistake :-).

All the best,

Maxim.


> Catch ya,
> 
> Dr. Dallas Warren
> Lecturer
> Department of Pharmaceutical Biology and Pharmacology
> Victorian College of Pharmacy, Monash University
> 381 Royal Parade, Parkville VIC 3010
> [EMAIL PROTECTED]
> +61 3 9903 9524
> -
> When the only tool you own is a hammer, every problem begins to 
> resemblea nail.
> 

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] Re: Gromacs3.3 Compilation on IBM powerpc

2006-03-13 Thread Yinghong



Dear Erik:
 
The last error is just "configure:6705: error: cannot compute sizeof 
(int), 77"
 
 
Xie Yinghong
 
 
> Find the _last_ error in config.log, and check what the compiler 
said  > about it.> Cheers,> ErikOn 
Mar 14, 2006, at 2:11 AM, Yinghong wrote:> Dear 
Erik:>> I have carefully read through config.log 
file, and I found out some  > possible problems for 
compilation failure, could you please help me  > to locate 
the real reason and solution?>> As I think, maybe 
this error is from xlc, but it seemed much  > related to 
"conftest.c" file, which was reported for several times  > as 
errors.>> How should I deal it with? By the way, I 
have passed through the  > compilation of single version. The 
following errors were from MPI  > compilation. Of course, I 
only added a "--enable-mpi" flag when  > compling MPI 
version, as described in the last 
email.>>> config.log:> 
...> ...> ...> ## --- 
##> ## Core tests. ##> ## --- 
##> ...> ...> configure:2586: checking 
for C compiler version> configure:2589: xlc --version 
&5> xlc: 1501-216 command option --version 
is not recognized - passed  > to ld> 
...> ...> ...> configure:2806: 
checking whether we are using the GNU C compiler> configure:2830: 
xlc -c -O3 -qarch=pwr4 -qtune=pwr4 -qansialias -w - > 
I/afs/ipp-garching.mpg.de/home/x/xieyh/fftw_213/include conftest.c 
>&5> "conftest.c", line 20.14: 1506-275 (S) Unexpected text 
me encountered.> "conftest.c", line 20.8: 1506-045 (S) Undeclared 
identifier choke.> configure:2836: $? = 1> 
configure: failed program was:> | /* confdefs.h.  
*/> |> | #define PACKAGE_NAME 
"GROMACS"> | #define PACKAGE_TARNAME "gromacs"> | 
#define PACKAGE_VERSION "3.3"> | #define PACKAGE_STRING "GROMACS 
3.3"> | #define PACKAGE_BUGREPORT "gmx-users at 
gromacs.org"> | #define PACKAGE "gromacs"> | 
#define VERSION "3.3"> | #define GMX_SOFTWARE_SQRT> 
| #define BUILD_TIME "Mon Mar 13 08:13:33 NFT 2006"> | #define 
BUILD_USER "xieyh at 
psi24"> | #define BUILD_MACHINE "AIX 2 
0010594A4C00"> | /* end confdefs.h.  */> 
|> | int> | main ()> | 
{> | #ifndef __GNUC__> 
|    choke me> | 
#endif> |> |   ;> 
|   return 0;> | }>> 
configure:2889: xlc -c -g -I/afs/ipp-garching.mpg.de/home/x/xieyh/ 
> fftw_213/include conftest.c >&5> 
configure:2895: $? = 0> configure:2899: test 
-z>  
|| test ! -s conftest.err> ...> ...> 
configure:3061: xlc -c -O3 -qarch=pwr4 -qtune=pwr4 -qansialias -w - 
> I/afs/ipp-garching.mpg.de/home/x/xieyh/fftw_213/include 
conftest.c >&5> "conftest.c", line 2.3: 1506-166 (S) 
Definition of function choke  > requires 
parentheses.> "conftest.c", line 2.9: 1506-046 (S) Syntax 
error.> configure:3067: $? = 1> configure: failed 
program was:> | #ifndef __cplusplus> |   
choke me> | #endif> ...> 
...> configure:4775: mpcc -E 
-I/afs/ipp-garching.mpg.de/home/x/xieyh/ > fftw_213/include 
conftest.c> "conftest.c", line 18.10: 1506-296 (S) #include 
file  >  not found.> 
...> ...> | #include > 
| #include > | #ifdef 
F77_DUMMY_MAIN> |> | #  ifdef 
__cplusplus> |  extern 
"C"> | #  endif> |    int 
F77_DUMMY_MAIN() { return 1; }> |> | 
#endif> | int> | main ()> | 
{> |> |   FILE *f = fopen ("conftest.val", 
"w");> |   if (! f)> 
| exit (1);> |   if (((long) 
(sizeof (int))) < 0)> | 
{> |   long i = longval 
();> |   if (i != ((long) 
(sizeof (int> |   exit 
(1);> |   fprintf (f, "%ld\n", 
i);> | }> |   
else> | {> 
|   unsigned long i = ulongval 
();> |   if (i != ((long) 
(sizeof (int> |   exit 
(1);> |   fprintf (f, "%lu\n", 
i);> | }> |   exit 
(ferror (f) || fclose (f) != 0);> |> |   
;> |   return 0;> | }> 
configure:6705: error: cannot compute sizeof (int), 77> See 
`config.log' for more 
details. Xie 
Yinghong 
Hi, Check config.log for the actual error. 
You might also want to make>> sure that you can compile a small 
MPI test program. 
Cheers, 
Erik On Mar 13, 2006, at 8:52 PM, Yinghong 
wrote:> Dear 
users:>> I have successfully 
installed fftw2.1.3 and single version>>> Gromacs3.3 on IBM 
powerpc Aix5.2.0.0 with the following>>> environmental 
variables:>> setenv CC 
xlc>>> setenv MPICC mpcc>>> setenv F77 
xlf>>> setenv CXX xlC>>> setenv CPPFLAGS 
-I.../fftw_213/include>>> setenv LDFLAGS 
-L.../fftw_213/lib>>> setenv CFLAGS "-O3 -qarch=pwr4 
-qtune=pwr4 -qansialias -w">>> setenv FFLAGS "-O3 
-qarch=pwr4 -qtune=pwr4 -qansialias -qstrict -w">>> Now, I 
would like to compile the parallel version of Gromacs3.3,>>> 
and I did it as the follows:>>> 1. make 
disclean;>>> 2. ./configure --prefix=$HOME/gromacs33 
--exec-prefix=$HOME/>>> gromacs33 --enable-float 
--disable-fortran --disable-largefile -->>> without-x 
--disable->>> vectorized-sqrt --with-fft=fftw2 --enable-mpi 
--program-suffix=_mpi>> It is worthy 
noting that the all parameters for ./configure are the>>> 
same as perviously successful single version, except that 
"--enable->>> mpi --program-suffix=_mpi" is

Re: [gmx-users] Re: Gromacs3.3 Compilation on IBM powerpc

2006-03-13 Thread Erik Lindahl
Find the _last_ error in config.log, and check what the compiler said  
about it.



Cheers,

Erik

On Mar 14, 2006, at 2:11 AM, Yinghong wrote:


Dear Erik:

I have carefully read through config.log file, and I found out some  
possible problems for compilation failure, could you please help me  
to locate the real reason and solution?


As I think, maybe this error is from xlc, but it seemed much  
related to "conftest.c" file, which was reported for several times  
as errors.


How should I deal it with? By the way, I have passed through the  
compilation of single version. The following errors were from MPI  
compilation. Of course, I only added a "--enable-mpi" flag when  
compling MPI version, as described in the last email.



config.log:
...
...
...
## --- ##
## Core tests. ##
## --- ##
...
...
configure:2586: checking for C compiler version
configure:2589: xlc --version &5
xlc: 1501-216 command option --version is not recognized - passed  
to ld

...
...
...
configure:2806: checking whether we are using the GNU C compiler
configure:2830: xlc -c -O3 -qarch=pwr4 -qtune=pwr4 -qansialias -w - 
I/afs/ipp-garching.mpg.de/home/x/xieyh/fftw_213/include conftest.c >&5

"conftest.c", line 20.14: 1506-275 (S) Unexpected text me encountered.
"conftest.c", line 20.8: 1506-045 (S) Undeclared identifier choke.
configure:2836: $? = 1
configure: failed program was:
| /* confdefs.h.  */
|
| #define PACKAGE_NAME "GROMACS"
| #define PACKAGE_TARNAME "gromacs"
| #define PACKAGE_VERSION "3.3"
| #define PACKAGE_STRING "GROMACS 3.3"
| #define PACKAGE_BUGREPORT "gmx-users@gromacs.org"
| #define PACKAGE "gromacs"
| #define VERSION "3.3"
| #define GMX_SOFTWARE_SQRT
| #define BUILD_TIME "Mon Mar 13 08:13:33 NFT 2006"
| #define BUILD_USER "[EMAIL PROTECTED]"
| #define BUILD_MACHINE "AIX 2 0010594A4C00"
| /* end confdefs.h.  */
|
| int
| main ()
| {
| #ifndef __GNUC__
|choke me
| #endif
|
|   ;
|   return 0;
| }

configure:2889: xlc -c -g -I/afs/ipp-garching.mpg.de/home/x/xieyh/ 
fftw_213/include conftest.c >&5

configure:2895: $? = 0
configure:2899: test -z
 || test ! -s conftest.err
...
...
configure:3061: xlc -c -O3 -qarch=pwr4 -qtune=pwr4 -qansialias -w - 
I/afs/ipp-garching.mpg.de/home/x/xieyh/fftw_213/include conftest.c >&5
"conftest.c", line 2.3: 1506-166 (S) Definition of function choke  
requires parentheses.

"conftest.c", line 2.9: 1506-046 (S) Syntax error.
configure:3067: $? = 1
configure: failed program was:
| #ifndef __cplusplus
|   choke me
| #endif
...
...
configure:4775: mpcc -E -I/afs/ipp-garching.mpg.de/home/x/xieyh/ 
fftw_213/include conftest.c
"conftest.c", line 18.10: 1506-296 (S) #include file  
 not found.

...
...
| #include 
| #include 
| #ifdef F77_DUMMY_MAIN
|
| #  ifdef __cplusplus
|  extern "C"
| #  endif
|int F77_DUMMY_MAIN() { return 1; }
|
| #endif
| int
| main ()
| {
|
|   FILE *f = fopen ("conftest.val", "w");
|   if (! f)
| exit (1);
|   if (((long) (sizeof (int))) < 0)
| {
|   long i = longval ();
|   if (i != ((long) (sizeof (int
|   exit (1);
|   fprintf (f, "%ld\n", i);
| }
|   else
| {
|   unsigned long i = ulongval ();
|   if (i != ((long) (sizeof (int
|   exit (1);
|   fprintf (f, "%lu\n", i);
| }
|   exit (ferror (f) || fclose (f) != 0);
|
|   ;
|   return 0;
| }
configure:6705: error: cannot compute sizeof (int), 77
See `config.log' for more details.



Xie Yinghong



Hi,

Check config.log for the actual error. You might also want to make
sure that you can compile a small MPI test program.

Cheers,

Erik

On Mar 13, 2006, at 8:52 PM, Yinghong wrote:


Dear users:

I have successfully installed fftw2.1.3 and single version
Gromacs3.3 on IBM powerpc Aix5.2.0.0 with the following
environmental variables:

setenv CC xlc
setenv MPICC mpcc
setenv F77 xlf
setenv CXX xlC
setenv CPPFLAGS -I.../fftw_213/include
setenv LDFLAGS -L.../fftw_213/lib
setenv CFLAGS "-O3 -qarch=pwr4 -qtune=pwr4 -qansialias -w"
setenv FFLAGS "-O3 -qarch=pwr4 -qtune=pwr4 -qansialias -qstrict -w"
Now, I would like to compile the parallel version of Gromacs3.3,
and I did it as the follows:
1. make disclean;
2. ./configure --prefix=$HOME/gromacs33 --exec-prefix=$HOME/
gromacs33 --enable-float --disable-fortran --disable-largefile --
without-x --disable-
vectorized-sqrt --with-fft=fftw2 --enable-mpi --program-suffix=_mpi

It is worthy noting that the all parameters for ./configure are the
same as perviously successful single version, except that "--enable-
mpi --program-suffix=_mpi" is added in compiling MPI version.
However, compilation can not pass because of the following reason.
If anybody encountered such an error, could you give me some
suggestion? How to deal with it? Any help will be extermely
appreciately.

...
...
checking for style of include used by make... GNU
checking dependency style of xlc... none
checking for mpxlc... mpcc
checking whether the MPI cc command works... yes
checking how to run 

[gmx-users] Re: Gromacs3.3 Compilation on IBM powerpc

2006-03-13 Thread Yinghong
Dear Erik:

I have carefully read through config.log file, and I found out some possible 
problems for compilation failure, could you please help me to locate the real 
reason and solution?

As I think, maybe this error is from xlc, but it seemed much related to 
"conftest.c" file, which was reported for several times as errors. 

How should I deal it with? By the way, I have passed through the compilation of 
single version. The following errors were from MPI compilation. Of course, I 
only added a "--enable-mpi" flag when compling MPI version, as described in the 
last email.


config.log:
...
...
...
## --- ##
## Core tests. ##
## --- ##
...
...
configure:2586: checking for C compiler version
configure:2589: xlc --version &5
xlc: 1501-216 command option --version is not recognized - passed to ld
...
...
...
configure:2806: checking whether we are using the GNU C compiler
configure:2830: xlc -c -O3 -qarch=pwr4 -qtune=pwr4 -qansialias -w 
-I/afs/ipp-garching.mpg.de/home/x/xieyh/fftw_213/include conftest.c >&5
"conftest.c", line 20.14: 1506-275 (S) Unexpected text me encountered.
"conftest.c", line 20.8: 1506-045 (S) Undeclared identifier choke.
configure:2836: $? = 1
configure: failed program was:
| /* confdefs.h.  */
| 
| #define PACKAGE_NAME "GROMACS"
| #define PACKAGE_TARNAME "gromacs"
| #define PACKAGE_VERSION "3.3"
| #define PACKAGE_STRING "GROMACS 3.3"
| #define PACKAGE_BUGREPORT "gmx-users@gromacs.org"
| #define PACKAGE "gromacs"
| #define VERSION "3.3"
| #define GMX_SOFTWARE_SQRT
| #define BUILD_TIME "Mon Mar 13 08:13:33 NFT 2006"
| #define BUILD_USER "[EMAIL PROTECTED]"
| #define BUILD_MACHINE "AIX 2 0010594A4C00"
| /* end confdefs.h.  */
| 
| int
| main ()
| {
| #ifndef __GNUC__
|choke me
| #endif
| 
|   ;
|   return 0;
| }

configure:2889: xlc -c -g 
-I/afs/ipp-garching.mpg.de/home/x/xieyh/fftw_213/include conftest.c >&5
configure:2895: $? = 0
configure:2899: test -z 
 || test ! -s conftest.err
...
...
configure:3061: xlc -c -O3 -qarch=pwr4 -qtune=pwr4 -qansialias -w 
-I/afs/ipp-garching.mpg.de/home/x/xieyh/fftw_213/include conftest.c >&5
"conftest.c", line 2.3: 1506-166 (S) Definition of function choke requires 
parentheses.
"conftest.c", line 2.9: 1506-046 (S) Syntax error.
configure:3067: $? = 1
configure: failed program was:
| #ifndef __cplusplus
|   choke me
| #endif
...
...
configure:4775: mpcc -E 
-I/afs/ipp-garching.mpg.de/home/x/xieyh/fftw_213/include conftest.c
"conftest.c", line 18.10: 1506-296 (S) #include file  not 
found.
...
...
| #include 
| #include 
| #ifdef F77_DUMMY_MAIN
| 
| #  ifdef __cplusplus
|  extern "C"
| #  endif
|int F77_DUMMY_MAIN() { return 1; }
| 
| #endif
| int
| main ()
| {
| 
|   FILE *f = fopen ("conftest.val", "w");
|   if (! f)
| exit (1);
|   if (((long) (sizeof (int))) < 0)
| {
|   long i = longval ();
|   if (i != ((long) (sizeof (int
|   exit (1);
|   fprintf (f, "%ld\n", i);
| }
|   else
| {
|   unsigned long i = ulongval ();
|   if (i != ((long) (sizeof (int
|   exit (1);
|   fprintf (f, "%lu\n", i);
| }
|   exit (ferror (f) || fclose (f) != 0);
| 
|   ;
|   return 0;
| }
configure:6705: error: cannot compute sizeof (int), 77
See `config.log' for more details.



Xie Yinghong


>Hi,
>
>Check config.log for the actual error. You might also want to make  
>sure that you can compile a small MPI test program.
>
>Cheers,
>
>Erik
>
>On Mar 13, 2006, at 8:52 PM, Yinghong wrote:
>
>> Dear users:
>>
>> I have successfully installed fftw2.1.3 and single version  
>> Gromacs3.3 on IBM powerpc Aix5.2.0.0 with the following  
>> environmental variables:
>>
>> setenv CC xlc
>> setenv MPICC mpcc
>> setenv F77 xlf
>> setenv CXX xlC
>> setenv CPPFLAGS -I.../fftw_213/include
>> setenv LDFLAGS -L.../fftw_213/lib
>> setenv CFLAGS "-O3 -qarch=pwr4 -qtune=pwr4 -qansialias -w"
>> setenv FFLAGS "-O3 -qarch=pwr4 -qtune=pwr4 -qansialias -qstrict -w"
>> Now, I would like to compile the parallel version of Gromacs3.3,  
>> and I did it as the follows:
>> 1. make disclean;
>> 2. ./configure --prefix=$HOME/gromacs33 --exec-prefix=$HOME/ 
>> gromacs33 --enable-float --disable-fortran --disable-largefile -- 
>> without-x --disable-
>> vectorized-sqrt --with-fft=fftw2 --enable-mpi --program-suffix=_mpi
>>
>> It is worthy noting that the all parameters for ./configure are the  
>> same as perviously successful single version, except that "--enable- 
>> mpi --program-suffix=_mpi" is added in compiling MPI version.  
>> However, compilation can not pass because of the following reason.  
>> If anybody encountered such an error, could you give me some  
>> suggestion? How to deal with it? Any help will be extermely  
>> appreciately.
>>
>> ...
>> ...
>> checking for style of include used by make... GNU
>> checking dependency style of xlc... none
>> checking for mpxlc... mpcc
>> checking whether the MPI cc command works... yes
>> checking how to run the C preprocessor... mpc

Re: [gmx-users] Gromacs3.3 Compilation on IBM powerpc

2006-03-13 Thread Erik Lindahl

Hi,

Check config.log for the actual error. You might also want to make  
sure that you can compile a small MPI test program.


Cheers,

Erik

On Mar 13, 2006, at 8:52 PM, Yinghong wrote:


Dear users:

I have successfully installed fftw2.1.3 and single version  
Gromacs3.3 on IBM powerpc Aix5.2.0.0 with the following  
environmental variables:


setenv CC xlc
setenv MPICC mpcc
setenv F77 xlf
setenv CXX xlC
setenv CPPFLAGS -I.../fftw_213/include
setenv LDFLAGS -L.../fftw_213/lib
setenv CFLAGS "-O3 -qarch=pwr4 -qtune=pwr4 -qansialias -w"
setenv FFLAGS "-O3 -qarch=pwr4 -qtune=pwr4 -qansialias -qstrict -w"
Now, I would like to compile the parallel version of Gromacs3.3,  
and I did it as the follows:

1. make disclean;
2. ./configure --prefix=$HOME/gromacs33 --exec-prefix=$HOME/ 
gromacs33 --enable-float --disable-fortran --disable-largefile -- 
without-x --disable-

vectorized-sqrt --with-fft=fftw2 --enable-mpi --program-suffix=_mpi

It is worthy noting that the all parameters for ./configure are the  
same as perviously successful single version, except that "--enable- 
mpi --program-suffix=_mpi" is added in compiling MPI version.  
However, compilation can not pass because of the following reason.  
If anybody encountered such an error, could you give me some  
suggestion? How to deal with it? Any help will be extermely  
appreciately.


...
...
checking for style of include used by make... GNU
checking dependency style of xlc... none
checking for mpxlc... mpcc
checking whether the MPI cc command works... yes
checking how to run the C preprocessor... mpcc -E
**
* Using CFLAGS from environment variable *
**
checking whether byte ordering is bigendian... yes
checking for egrep... grep -E
checking for ANSI C header files... no
checking for sys/types.h... yes
checking for sys/stat.h... yes
checking for stdlib.h... yes
checking for string.h... yes
checking for memory.h... yes
checking for strings.h... yes
checking for inttypes.h... yes
checking for stdint.h... yes
checking for unistd.h... yes
checking for int... yes
checking size of int... configure: error: cannot compute sizeof  
(int), 77

See `config.log' for more details.


Xie Yinghong
Hong Kong University
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: [gmx-users] g_clustsize problem

2006-03-13 Thread David van der Spoel

Wang Zhun wrote:

Hi,

I'm working on peptide folding with GMX3.2.1 now. When I carried out 
g_clustsize, it says:
Group 1 ( Protein) has   211 elements
Group 2 (   Protein-H) has   167 elements
Group 3 ( C-alpha) has22 elements
Group 4 (Backbone) has66 elements
Group 5 (   MainChain) has89 elements
Group 6 (MainChain+Cb) has   108 elements
Group 7 ( MainChain+H) has   112 elements
Group 8 (   SideChain) has99 elements
Group 9 ( SideChain-H) has78 elements
Group10 ( Prot-Masses) has   211 elements
Group11 ( Non-Protein) has 20439 elements
Group12 ( SOL) has 20439 elements
Group13 (   Other) has 20439 elements
Select a group: 1
Selected 1: 'Protein'
Last frame   4000 time 2.000
cmid: 1, cmax: 1, max_size: 211
Fatal error: Lo: 0.00, Mid: 1.00, Hi: 1.00
What's the meaning of the last two lines?


there is only one cluster, which is logical since you have one molecule. 
 you probably want to use g_cluster

Thank you in advance!


___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php



--
David.

David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,  75124 Uppsala, Sweden
phone:  46 18 471 4205  fax: 46 18 511 755
[EMAIL PROTECTED]   [EMAIL PROTECTED]   http://folding.bmc.uu.se

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php