Displaying 3 results from an estimated 3 matches for "dgetrf".
2022 Nov 09
1
det(diag(c(NaN, 1))) should be NaN, not 0
Hello,
Currently, determinant(A) calculates the determinant of 'A' by factorizing
A=LU and computing prod(diag(U)) [or the logarithm of the absolute value].
The factorization is done by LAPACK routine DGETRF, which gives a status
code INFO, documented [1] as follows:
*> INFO is INTEGER
*> = 0: successful exit
*> < 0: if INFO = -i, the i-th argument had an illegal value
*> > 0: if INFO = i, U(i,i) is exactly zero. The factorization
*>...
2003 Jun 23
3
FW: S4 classes, creating in C
..., info;
if (!isMatrix(aa) || !isReal(aa)) {
error("A must be a double precision matrix");
}
adims = GET_DIM(aa); m = INTEGER(adims)[0]; n = INTEGER(adims)[1];
pivot = PROTECT(NEW_INTEGER(m < n ? m : n));
F77_CALL(dgetrf)(&m, &n, REAL(aa), &m, INTEGER(pivot), &info);
check_Lapack_error(info, "dtrtrf");
val = PROTECT(NEW_OBJECT(MAKE_CLASS("LUdecomposition")));
SET_SLOT(val, install("a"), aa);
SET_SLOT(val, install("pi...
2015 Mar 25
4
F77_CALL/NAME problem
...int*, int*, double*, int*, int*, double*, int*, int*);
void solve( int *p, double *A, double *Ainv)
{
...
F77_CALL(DGESV)(p, p, Ain, p, ipiv, Bin, p, &info);
}
In order to create the ".so" file to load in R I downloaded the dgesv.f file
as well as the dependencies (dgetf2.f dgetrf.f dgetrs.f dlaswp.f). As I was
used to I ran in a terminal
R CMD COMPILE *.f
R CMD SHLIB MY_C_FILE.c *.o
to get the MY_C_FILE.so file. However, when I try to load it in a R session
I get the following error message "undefined symbol: DGESV_". Similar errors
are occurring with some other...