search for: ipiv

...k is simply the value of the single element in the block. If the diagonal block is of dimension 2x2 then the determinant of the block may be computed according to the well-known formula b11*b22-b12*b21, where bij is the value in the i'th row and j'th column of the block. int i, q, info, *ipiv, one = 1; double *b, *A, *D, det; /* ** A and D are upper triangular matrices in packed storage ** A[] = a00, a01, a11, a02, a12, a22, a03, a13, a23, a33, ... ** use the following macro to address the element in the ** i'th row and j'th column for i <= j */ #define UMA...

...TECT(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]; ncB = Bdims[1]; Ap = REAL(A); Bp = REAL(B); /* int info, *ipiv = (int *) R_alloc(nrA, sizeof(int)); */ /* F77_CALL(dgesv)(&nrA, &ncB, Ap, &nrA, ipiv, Bp, &nrA, &info); */ mysolve2(Ap, &nrA, Bp, &ncB, &tol); UNPROTECT(2); return B; } void mysolve2(double *A, int *nrA, double *B, int *ncB, double *tolin) { int info;...

...t of the C file I am currently writing #include<stdio.h> #include<R.h> #include<Rmath.h> #include<math.h> void F77_NAME(DGESV)( 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. Howeve...

...*** *** 4827,4832 **** --- 4832,4990 ---- popad ret + ; + ; check_cpuid by <appro at fy.chalmers.se>. + ; + ; This routine attempts to fetch config file based on CPU architecture + ; instead of [partial] IP number. Attempted file name matches following + ; expression: + ; + ; [XY86|ipiv|ia64][-smp] + ; + ; where X is either 'a' for AMD or 'i' for all others, Y varies from 3 + ; to 6 on non-AMD CPUs and from 4 to 7 on AMD with exclusion for AMD + ; x86-64 where Y turns 'w' (as in "wide":-). Intel P4 which is denoted + ; by 'ipiv', while Int...

...ieve that "fact" should've been char instead of int, i.e. F77_NAME(dgesvx)(const char* fact, const char* trans, const int* n, My reasoning: In the R-source: src/modules/lapack/dlapack.f I find the following lines: * SUBROUTINE DGESVX( FACT, TRANS, N, NRHS, A, LDA, AF, LDAF, IPIV, * EQUED, R, C, B, LDB, X, LDX, RCOND, FERR, BERR, * WORK, IWORK, INFO ) * * .. Scalar Arguments .. * CHARACTER EQUED, FACT, TRANS * INTEGER INFO, LDA, LDAF, LDB, LDX, N, NRHS * DOUBLE PRECISION RCOND w...

...iJRTeHHXw4W/vjyc08vr0qIjVjfeDDHNsKMH2X7UDCmGKY1ifeG6XuYGty +izFd9YNS/U+rpehV1DaR0l1v20+2p7z9JEPQdHFfezpYycSA8Ya1u0ueG9GUfsuxyYBP4HS ZyaHJsIia0+tfX5r7IRleE42+XD0rcd79ybgUB8xiX1Gp/T5rnPdDtbHJm/Pq/hOWivfuP0h Z1+C57QvITNPnc29+vgcH66PCxG8LTL6VHPbgt+pT2/0bZxbjOhj0g6Ky2zdgVxPH/2imdHH lczd+ois3bW3bTbHOX0G1+IpIvs/RJ/6nbRfn2B0Na6rT/FpB4XT5H8ox0Yothhl4YWKjfVp o9mKmdenk7KYeJ8+bqr4VVYN2aWP4GH6rGk1nSNnIn06o9dxZfBlYq8+/iAfztH6RCu8XRT9 hrK74w/Ib0y0W/0acUWmNBzpc3tCmw4qnVJMnvfpE9aH34vd+rRft+yLM6GPHr0N7G57EWlP 6HOsQIdP3q4gp0/TYx3koo31Gc8jn1r961AfmU6xh3qXPnF19PRR5TfYmtJOEI4N9YlHVyN7 T6mqYZ8+dwv0L8S3jnfzTn1kuclNGeozM4...