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 Maechler wrote:
On 9 F
> 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 package 'glmmML' I'm using old C code and linpa
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 part of the code looks li
> 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 < *n; j++)
>>for (i
> > 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, &bdim, &bdim, &rcond, work, info);
> >if (*info == 0)
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, &bdim, &bdim, &rcond, work, inf
> 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, &bdim, &bdim, &rcond, work, info);
>if (*info == 0){
>F77
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, &bdim, &bdim, &rcond, work, info);
if (*info == 0){
F77_CALL(dpodi)(*hessian, &bdim, &bdim, det, &job);
.