Member Site › Forums › Rosetta 3 › Rosetta 3 – Applications › RosettaVIP: ERROR: seqpos <= size()
- This topic has 5 replies, 2 voices, and was last updated 12 years, 5 months ago by attesor.
-
AuthorPosts
-
-
July 26, 2012 at 2:30 pm #1356Anonymous
Hi, I tried RosettaVIP shipped in Rosetta3.4. In some cases it breaks after a few iterations with this error:
ERROR: seqpos <= size()
ERROR:: Exit from: src/core/conformation/Conformation.hh line: 289
I can not figure out the root for the error from the source file. The command line I used is:
vip.linuxgccrelease -s mypdb -cp:cutoff 6.0 -sasa_calculator_probe_radius 1.0 -cp:output myoutput -cp:print_reports -run:silent
I have not observed other anomalies in the screen printout before the ERROR, neither in the reports.txt file.
-
July 26, 2012 at 3:22 pm #7467Anonymous
Turn off -run:silent and see if there are more messages relating to the error? That message itself it just a generic message saying that the code tried to access a residue that doesn’t exist (say, residue 101 of a 100 residue protein).
-
July 26, 2012 at 3:56 pm #7470attesor
I turned off the silent switch and ran VIP again. Sometimes I could reproduce the ERROR (I guess because RosettaHoles used in VIP is stochastic (see RosettaVIP manual)). Whenever the ERROR happens, I see a very strange energy value (no other WARNINGs or ERRORs observed):
Total weighted score: -255.319
VIP: Comparing new energy 1.73326e-316 with old energy -255.193
VIP: Rejecting attempted mutation - finished!
-
July 26, 2012 at 5:39 pm #7473Anonymous
Stochastic errors like this are generally due to bugs in the code, not errors in your input. Can you share your input data so that I can see if I can reproduce it? I’ll file a bug report in our tracking system either way and bring it to the attention of the code’s author.
-
July 26, 2012 at 11:32 pm #7477Anonymous
This is the author’s comment:
“Yeah, I’m aware of the error… it occurs when RosettaHoles doesn’t find any voids. Holes fills voids with cavity balls, which I essentially tacked on to the pose as extra ‘residues’, so if there are none, vip tries to read the first cavity ball and crashes, hence the seqpos error. Because the cavity finding is stochastic, the error is also, however, it shouldn’t cause any problems aside from being a bit annoying.
I’ll implement a more elegant exit at some point and it will be fixed for the next release.”
-
July 27, 2012 at 10:29 am #7486attesor
Thanks a lot for the help!
As a matter of fact, it does cause another problem: the final PDB file is empty if this error is reported.
-
-
AuthorPosts
- You must be logged in to reply to this topic.