Re: AW: update: segmentation faults with AMBER topology and RAMD tcl scripts

From: Balazs Jojart (jojartb_at_gmail.com)
Date: Wed Feb 24 2010 - 01:52:26 CST

Dear Vlad,
I think the topology is OK, because it was used in a former simulation
only using standard MD. In the output file (RAMD), the the structure
information looks ok. I tested the loadcoords part as well, by adding
"puts coordinates(1)" lines in the tcl code, and the coordinates are
also missing - same error "no such element in array".
Our system administrator compiled the newest CVS of namd against tcl8.5,
fortunately I don't obtain any segfault error, but unfortunately, the
same error was occurred with the example files as well, as described
below. The outputs are attached.
Do we have to compile with another version of tcl (for example Tcl8.4.19)?
Cheers,
Balazs

Cojocaru,Vlad wrote:
> Dear Balazs,
>
> This is a more friendly error ...
> It looks as if your definition of the ramd and/or protein atoms is somehow wrong .. You have an atom that somehow does not have a mass .. Can you double check ?
>
> Otherwise, if you could send me a zipped file with the files required to reproduce your problem (topology, coordinates, box info files) I can try to see what is wrong ..
>
> vlad
>
>
> -----Ursprüngliche Nachricht-----
> Von: Balazs Jojart [mailto:jojartb_at_gmail.com]
> Gesendet: Di 23.02.2010 17:05
> An: Cojocaru,Vlad
> Betreff: Re: namd-l: update: segmentation faults with AMBER topology and RAMD tcl scripts
>
> Dear Vlad,
> Using the RAMD method with the compiled namd2.7beta2 against Tcl8.5 and
> charmm topologies I obtain the following error message:
>
> TCL: can't read "masses(5962)": no such element in array
> FATAL ERROR: can't read "masses(5962)": no such element in array
> while executing
> "expr $ligMass + $masses($ramdAtom)"
> ("foreach" body line 2)
>
> the whole error message is attached.
> I'm going to test the charmm topologies, with the namd2.6 + tcl8.5
> version compilation.
> Cheers,
> Balazs
>
>
> Vlad Cojocaru wrote:
>
>> Dear all,
>>
>> To test the errors I've been commenting on recently I compiled today's
>> CVS code of NAMD (23 February 2010) using exactly the same protocol as
>> for my last running compilation (13 August 2009). Tcl version is
>> tcl8.4.19. Architecture: x86_64 + infniband + intel 11.0 compiler +
>> mvapich1.4rc1. The compilation worked absolutely fine.
>>
>> Then I ran the RAMD tests which resulted in "Segmentation violation" .
>> Then I tested the new executable for one of my standard MD simulation
>> using AMBER topology. The same undefined error occurred.
>>
>> To conclude, the latest CVS code produces segmentation faults on AMBER
>> topologies regardless whether RAMD is switched on or not. Therefore,
>> this error apparently has nothing to do with the RAMD scripts. It
>> looks as if the amber topology is read correctly but the program
>> exists with segmentation violation immediately after reading the
>> topology (below last lines of the output)
>>
>> I could not test against CHARMM topologies as I don't have one but
>> another NAMD user mentioned recently that the code runs fine on CHARMM
>> topologies while observing the same errors with AMBER topologies.
>>
>> To test the RAMD scripts until the problem with the AMBER topologies
>> will be fixed, I suggest to compile namd 2.6 against tcl8.4.19 or above.
>> If you still see errors, please let me know.
>>
>> Cheers
>> Vlad
>>
>> -----------------output sample------------------------------
>> Reading parm file (1r9o_1.top) ...
>> PARM file in AMBER 7 format
>> Warning: Encounter 10-12 H-bond term
>> Warning: Found 18886 H-H bonds.
>> Info: SUMMARY OF PARAMETERS:
>> Info: 62 BONDS
>> Info: 131 ANGLES
>> Info: 49 DIHEDRAL
>> Info: 0 IMPROPER
>> Info: 0 CROSSTERM
>> Info: 0 VDW
>> Info: 210 VDW_PAIRS
>> Info: TIME FOR READING PDB FILE: 0.158216
>> Info:
>> Info: Reading from binary file 1r9o_1_eq10_rst.500000.coor
>> Info: ****************************
>> Info: STRUCTURE SUMMARY:
>> Info: 64265 ATOMS
>> Info: 64360 BONDS
>> Info: 13994 ANGLES
>> Info: 29662 DIHEDRALS
>> Info: 0 IMPROPERS
>> Info: 0 CROSSTERMS
>> Info: 98534 EXCLUSIONS
>> Info: 60475 RIGID BONDS
>> Info: 132320 DEGREES OF FREEDOM
>> Info: 22676 HYDROGEN GROUPS
>> Info: 4 ATOMS IN LARGEST HYDROGEN GROUP
>> Info: 22676 MIGRATION GROUPS
>> Info: 4 ATOMS IN LARGEST MIGRATION GROUP
>> Info: TOTAL MASS = 394234 amu
>> Info: TOTAL CHARGE = -1.48128e-05 e
>> Info: MASS DENSITY = 1.00973 g/cm^3
>> Info: ATOM DENSITY = 0.0991208 atoms/A^3
>> Info: *****************************
>> Info:
>> Info: Entering startup at 30.2467 s, 156.332 MB of memory in use
>> Info: Startup phase 0 took 0.00019908 s, 156.332 MB of memory in use
>>
>>
>>
>
>
>


This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 05:22:45 CST