search for: nrhs

...tions is often of interest in statistical computations, i'd like to get a working example of using Lapack for regression. However, i run into an error. My matrix-lapack-example.c file: #include <R_ext/Lapack.h> void reg(const char* trans, const int* m, const int* n, const int* nrhs, double* a, const int* lda, double* b, const int* ldb, double* work, const int* lwork, int* info) { F77_CALL(dgels)(trans, m, n, nrhs, a, lda, b, ldb, work, lwork, info); } My matrix-lapack-example.R file: dyn.load( "matrix-lapack-example.so" ) regression <- funct...

...const int* n, and I believe 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 * DOU...

...= 5, size = 7, vis = TRUE )" crashes due to corruption of C stack by Fortran function DPOSV, when compiled with the new gfortran and with -O2. To see the problem, one can just look at the disassembly of DPOSV (LAPACK), neither the package nor R is not necessary: SUBROUTINE DPOSV( UPLO, N, NRHS, A, LDA, B, LDB, INFO ) CHARACTER????????? UPLO In one case, DPOSV calls DPOTRS before returning. The new gfortran with -O2 performs tail-call optimization, jumping to DPOTRS. In the annotated disassembly snippet, at 11747f1, DPOSV tries to ensure that there is constant 1 as string length of UP...

Hi, I have tried to pinpoint potential problems which could lead to the LAPACK issues that are currently seen in R. I built the current R trunk using AR=gcc-ar RANLIB=gcc-ranlib ./configure --prefix=$HOME --enable-lto --enable-BLAS-shlib=no --without-recommended-packages and used this to find problem areas. There are quite a few warnings that were flagged, due to mismatches in function

...with positive definite A. My code looks like the following in solve.c ========================== #include<stdio.h> #include <R_ext/BLAS.h> #include <R_ext/Lapack.h> int main(){ double A[4]={1,0.5,0.5,1}; double B[2]={3,4}; char uplo='U'; int n = 2, nrhs=1, lda=2, ldb=2, info, i; F77_CALL(dposv)(&uplo,&n, &nrhs, A, &lda, B, &ldb, &info); for(i=0; i<2; i++){ printf("%f\n", B[i]); } return info; } ========================== When I am trying to link to BLAS/LAPACK using gcc -std=gnu99 solve.c...

...ts.data.cpp -o ts.data.o g++ -Ic:/src/R/R-2.5.1/include -Wall -O2 -c tseries.cpp -o tseries.o tseries.cpp: In function `tseries linear_model(const tseries&, const tseries&)': tseries.cpp:1058: warning: unused variable 'trans' tseries.cpp:1059: warning: unused variable 'nrhs' tseries.cpp:1060: warning: unused variable 'info' c:/src/R/Rtools/MinGW/bin/../lib/gcc/mingw32/3.4.5/../../../../include/c++/ 3.4.5/bits/stl_uninitiali zed.h: In member function `std::vector<std::string, std::allocator<std::string> > tseries::getColNam es() const': c:/s...

Dear R-devel, I am trying to use Fortran DGESV subroutine into C. Here it is the relevant part 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,

...Bf4Wvmmo+{aI-t!ilr;qI zsRj76++x2ZaDng+U+17kt at eb54lZRL)zGmeD{ZnXETI^vOeOH9u+-iQ=;a!^v}6`# zOTaw8N<T`O&-E4_)y`V%tv2j(0lWO{u-gJq`_MqFAgTq3eIz<tirtqb<RWCE=cwC! zx?vf}Zt8SnHOq@%TG-LbGP7E&V3);Z5qGt_seT+UE)n65%K7E-<cheI)Z}ZgafQRf zjA>6&kj8bHJpRW&#dd*vWgND&duj2eNvmjS=a#nrhs)sl50~M2?pUGY at W-xVVo6hC zMoar at w=2C`!(rt0EkBeCt65WV0Eb-TN1>P|A{FSa2tYvw_J#Yk_I^$3`<dv6odGC% zpfQ<Vs-^bU_O<P;l?zLeVi)(&x3E_z2&rfiQZe(02HGiR7Hc5mqDeAPSP at VulQngf zR+F#MxWRPc3H(})UlYVn4F`wd4CAL3KXv%2$4>)(#^7fx1CHb0 at xXHJH4_j%5$*%p zHI4Yb5WeO-e7}FqBz&F;Uo%-4R`1c;H+az}...