Displaying 6 results from an estimated 6 matches for "bdim".
Did you mean:
bdi
2017 Feb 09
3
Ancient C /Fortran code linpack error
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);
........
This usually works OK, but with an ill-conditioned data set (from a user
of glmmML) it happened that the hessian was all nan. However, dpoco
retur...
2017 Feb 09
3
Ancient C /Fortran code linpack error
...Feb 2017, at 16:00, G?ran Brostr?m <goran.brostrom at umu.se> 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_CALL(dpodi)(*hessian, &bdim, &bdim, det, &job);
> > ........
> >
> > This usually works OK, but with an ill-conditioned data set (from a user of glmmML) it happened that the...
2017 Feb 10
1
Ancient C /Fortran code linpack error
...rostr?m <goran.brostrom at umu.se> 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_CALL(dpodi)(*hessian, &bdim, &bdim, det, &job);
>>>> ........
>>>>
>>>> This usually works OK, but with an ill-conditioned data set (from a user...
2008 Apr 10
1
Computing time when calling C functions - why does an extra function call induce such an overhead?
...all to mysolve2) then the computations take only 0.03 secs. Can anyone tell me why there is such a "time-overhead" in introducing the extra function call to mysolve2?
I am on windows XP using R 2.6.1
Best regards
S?ren
SEXP mysolve1(SEXP Ain, SEXP Bin, SEXP tolin)
{
int *Adims, *Bdims;
double tol = asReal(tolin);
SEXP A, B;
PROTECT(A = duplicate(Ain));
PROTECT(B = duplicate(Bin));
Adims = INTEGER(coerceVector(getAttrib(A, R_DimSymbol), INTSXP));
Bdims = INTEGER(coerceVector(getAttrib(B, R_DimSymbol), INTSXP));
int nrA, ncB;
double *Ap, *Bp;
nrA = Adims[0...
2017 Feb 09
0
Ancient C /Fortran code linpack error
> On 9 Feb 2017, at 16:00, G?ran Brostr?m <goran.brostrom at umu.se> 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_CALL(dpodi)(*hessian, &bdim, &bdim, det, &job);
> ........
>
> This usually works OK, but with an ill-conditioned data set (from a user of glmmML) it happened that the hessian was all nan. How...
2017 Feb 10
0
Ancient C /Fortran code linpack error
...at 16:00, G?ran Brostr?m <goran.brostrom at umu.se> 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_CALL(dpodi)(*hessian, &bdim, &bdim, det, &job);
>>> ........
>>>
>>> This usually works OK, but with an ill-conditioned data set (from a user of glmmML) it hap...