[ADMB Users] Possible overflow in DF_dvlog

[Mollie Brooks]

Checking out a new copy and rebuilding got rid of the singular matrix error and the out of bounds error. So that's an improvement. 

On closer inspection, I realized that both computers were giving the "Possible overflow in DF_dvlog" warning this whole time. I went into the src code and found that this error could have been caused by accidentally taking the log of 0. I'm trying different starting values to see if I can avoid this. If that doesn't work I'll try posfun().

There is one strange thing though...
svn said that I updated to a revision in the 600s, but when I use the -version command, it still says I'm on revision 591.
Does the -version command have to be updated manually?

thanks,
Mollie 

On 12 Sep 2012, at 5:23 PM, John Sibert wrote:

> It looks like the LU decomposition code was reverted in r602.
> http://admb-project.org/redmine/projects/issues/repository/revisions
> 
> John Sibert
> Emeritus Researcher, SOEST
> University of Hawaii at Manoa
> 
> Visit the ADMB project http://admb-project.org/
> 
> On 09/12/2012 09:23 AM, Johnoel Ancheta wrote:
>> Try and recheckout a fresh copy first and rebuild all, then rerun...
>> 
>> On Wed, Sep 12, 2012 at 8:11 AM, Mollie Brooks <mbrooks at ufl.edu <mailto:mbrooks at ufl.edu>> wrote:
>> 
>>    I just noticed that I'm only getting this behaviour on one
>>    computer but not the other.
>> 
>>    The one it works on gives this version info
>> 
>>    bash-3.2$ ./pond -version
>>    ADMB Program: ./pond
>> 
>>    ADMB-11.X safe libraries compiled with Gnu C++ 4.2.1(64bit)
>>    Build date: Aug 31 2012
>>    Revision: 591:595M
>> 
>>    Copyright (c) 2008-2012 ADMB Foundation and Regents of the
>>    University of California
>> 
>>    The computer where it doesn't work gives this version info
>>    [mbrooks at r11a-s18 Bernoulli]$ ./pond -version
>>    ADMB Program: ./pond
>> 
>>    ADMB-trunk safe libraries compiled with Gnu C++ 4.4.6(64bit)
>>    Build date: Aug 31 2012
>>    Revision: 591
>> 
>>    Could this be a problem with the ludecomp reversion? I'm using
>>    choleski_decomp() in this program.
>> 
>>    thanks,
>>    Mollie
>> 
>>    Mollie Brooks
>>    Postdoctoral Researcher, Ponciano Lab
>>    Biology Department, University of Florida
>>    http://people.biology.ufl.edu/mbrooks
>> 
>> 
>>    On 12 Sep 2012, at 1:44 PM, Mollie Brooks wrote:
>> 
>>    > Hi,
>>    > I've been getting a lot of these warnings.
>>    > Possible overflow in DF_dvlog
>>    >
>>    > Then a lot of these
>>    > Error in matrix inverse -- matrix singular in inv(dmatrix)
>>    >
>>    > and then one of these that crashes the program
>>    > matrix bound exceeded -- row index too low in
>>    dmatrix::operator()value was0
>>    >
>>    > I'm almost certain that I'm accessing the matrix elements
>>    correctly and the matrix shouldn't be singular at the initial values.
>>    > Is there anything else that could lead to this series of errors?
>>    >
>>    > The other output is
>>    >  Inner f = -inf
>>    > Newton raphson 1
>>    > f = nan max g = nan
>>    > Newton raphson 2
>>    > f = nan max g = nan
>>    >
>>    > thanks,
>>    > Mollie
>>    > _______________________________________________
>>    > Users mailing list
>>    > Users at admb-project.org <mailto:Users at admb-project.org>
>>    > http://lists.admb-project.org/mailman/listinfo/users
>>    >
>> 
>>    _______________________________________________
>>    Users mailing list
>>    Users at admb-project.org <mailto:Users at admb-project.org>
>>    http://lists.admb-project.org/mailman/listinfo/users
>> 
>> 
>> 
>> 
>> _______________________________________________
>> Users mailing list
>> Users at admb-project.org
>> http://lists.admb-project.org/mailman/listinfo/users
> 
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.admb-project.org/pipermail/users/attachments/20120912/2f9f13c6/attachment.html>