Re: [Rd] Ancient C /Fortran code linpack error

2017-02-10 Thread Göran Broström
Thanks Berend, I will make that change and submit to CRAN. Best, Göran On 2017-02-10 16:13, Berend Hasselman wrote: On 10 Feb 2017, at 14:53, Göran Broström wrote: Thanks to all who answered my third question. I learned something, but: On 2017-02-09 17:44, Martin

Re: [Rd] Ancient C /Fortran code linpack error

2017-02-10 Thread Berend Hasselman
> On 10 Feb 2017, at 14:53, Göran Broström wrote: > > Thanks to all who answered my third question. I learned something, but: > > On 2017-02-09 17:44, Martin Maechler wrote: >> On 9 Feb 2017, at 16:00, Göran Broström wrote: In my

Re: [Rd] Ancient C /Fortran code linpack error

2017-02-10 Thread Göran Broström
Thanks to all who answered my third question. I learned something, but: On 2017-02-09 17:44, Martin Maechler wrote: On 9 Feb 2017, at 16:00, Göran Broström wrote: In my package 'glmmML' I'm using old C code and linpack in the optimizing procedure. Specifically, one

Re: [Rd] Ancient C /Fortran code linpack error

2017-02-09 Thread Berend Hasselman
> On 9 Feb 2017, at 17:44, Martin Maechler wrote: > [snip] ... >> You should/could use macro R_FINITE to test each entry of the hessian. >> In package nleqslv I test for a "correct" jacobian like this in file >> nleqslv.c in function fcnjac: > >>for (j = 0; j

Re: [Rd] Ancient C /Fortran code linpack error

2017-02-09 Thread Martin Maechler
> > On 9 Feb 2017, at 16:00, Göran Broström wrote: > > > > In my package 'glmmML' I'm using old C code and linpack in the optimizing > > procedure. Specifically, one part of the code looks like this: > > > >F77_CALL(dpoco)(*hessian, , , , work, info); > >if

Re: [Rd] Ancient C /Fortran code linpack error

2017-02-09 Thread Tomas Kalibera
On 02/09/2017 05:05 PM, Berend Hasselman wrote: On 9 Feb 2017, at 16:00, Göran Broström wrote: In my package 'glmmML' I'm using old C code and linpack in the optimizing procedure. Specifically, one part of the code looks like this: F77_CALL(dpoco)(*hessian, , , ,

Re: [Rd] Ancient C /Fortran code linpack error

2017-02-09 Thread Berend Hasselman
> On 9 Feb 2017, at 16:00, Göran Broström wrote: > > In my package 'glmmML' I'm using old C code and linpack in the optimizing > procedure. Specifically, one part of the code looks like this: > >F77_CALL(dpoco)(*hessian, , , , work, info); >if (*info == 0){ >

[Rd] Ancient C /Fortran code linpack error

2017-02-09 Thread Göran Broström
In my package 'glmmML' I'm using old C code and linpack in the optimizing procedure. Specifically, one part of the code looks like this: F77_CALL(dpoco)(*hessian, , , , work, info); if (*info == 0){ F77_CALL(dpodi)(*hessian, , , det, ); This usually works OK,