IUBio Biosequences .. Software .. Molbio soft .. Network News .. FTP

Problem with profilesearch (with a solution)

Heikki Lehvaslaiho Heikki.Lehvaslaiho at csc.fi
Fri Mar 31 07:15:09 EST 1995


Friends,

Some time ago, our GCG program profilesearch started crashing middle of
searches with an error message:

    Sequence information overflows string buffer!

After some head scratching and a mail message to gcg.com, things were
cleared out.

The problem was there only when the whole SWISS-PROT database was
searched. Profilesearch does something very stupid. It buffers all of the
sequence
descriptions as it searches the database and it happens that descriptions
in SWISS-PROT  are longer than in PIR.

To get going again, one has to modify the critical parameter and
recompile:

profilesearch.f:

!*      The names and description of each sequence is saved in an
!*      array called bigstring.  This is the major item using memory.
!*      If your memory is limited change this so that only the sequence 
!*      entry name is saved, and decrease the size of MaxStrBuff.  

       Integer MAXSTRBUFF
        Parameter (MAXSTRBUFF=4 000 000)


But to what value one should give to MAXSTRBUFF?  The default value
crashed the program at 36 k sequences. I have it now at its maximmum 7
000 000 which should be enough until SWISS-PROT reaches 60 k entries. (It
is currently 43470+483). Luckily, a note from GCG said that they "will
change the way the program handles the documentation for individual
sequences and try to remove the limit of 60,000".

        -Heikki

______ _/      _/____________________________________________________
      _/      _/
     _/  _/  _/  Heikki Lehvaslaiho  <Heikki.Lehvaslaiho at CSC.FI>
    _/_/_/_/_/  X.400: /G=Heikki/S=Lehvaslaiho/O=CSC/A=fumail/C=fi/
   _/  _/  _/  Center for Scientific Computing
  _/  _/  _/  Tietotie 6, P.O. Box 405, FIN-02101 Espoo FINLAND
 _/  _/  _/   Phone: +358 0 457 2076       FAX: +358 0 457 2302
    _/         URL:  http://www.csc.fi/molbio/
__ _/_/_/_/_/________________________________________________________



More information about the Info-gcg mailing list

Send comments to us at biosci-help [At] net.bio.net