search for: dgemm

Hi, I'm trying to follow the R-admin instructions for using MKL10 as the external BLAS compiling R-2.6.2 under Linux on a RH EL head node of a cluster. The configure process seems to have problems when it checks for dgemm in the BLAS. I'm using configure as: ./configure CC=icc F77=ifort --with-lapack="$MKL" --with-blas="$MKL" where $MKL is defined as in R-admin section A.3.1.4. checking for cblas_cdotu_sub in vecLib framework... no checking for dgemm_ in -L/usr1/util/Intel/mkl/10.0.1.0...

Hi, I am trying to compile a single module (dgemm.c) using clang and generate dgemm.ll. Command: clang -emit-llvm dgemm.c -S -o dgemm.ll Some of the functions are declared as static and clang ignores these functions: One of the functions is: static void innerloop(double a, const double* b, double* c, double beta) { *b = a*beta; return; } Is...

I am using R to multiply some large (30k x 30k double) matrices on a 64 core machine (xeon phi). I added some timers to src/main/array.c to see where the time is going. All of the time is being spent in the matprod function, most of that time is spent in dgemm. 15 seconds is in matprod in some code that is checking if there are NaNs. > system.time (C <- B %*% A) nancheck: wall time 15.240282s ? dgemm: wall time 43.111064s matprod: wall time 58.351572s ??? user?? system? elapsed 2710.154?? 20.999?? 58.398 The NaN checking code is not being vector...

...fore, but in this case I need to do matrix operations in my C code. As I have never done that before, I'm trying to write some simple examples to make sure I understand the basics. I am stuck on the first one. I'm trying to write a function to multiply two matrices using the blas routine dgemm. The name of my example package is CMATRIX. My code is as follows. I have a file matrix.c in my src directory: #include <R.h> #include <R_ext/Utils.h> #include <R_ext/Lapack.h> #include <R_ext/BLAS.h> //Computes C = A*B void R_matrix_multiply(double * A, double * B, int...

Hi I am new to writing C code and am trying to write an R extension in C. I have hit a wall with F77_CALL(dgemm) in that it produces wrong results. The code below is a simplified example that multiplies the matrices Ab and Bm to give Cm. The results below show clearly that Cm is wrong. Am= 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 Bm= 1 1 1 1 1 1 1 1 1 1...

...> define internal i32 @add(i32 %a, i32 %b) nounwind { > [...] > } > > Olivier. > > > On Mon, Dec 14, 2009 at 8:49 PM, Arvind Sudarsanam > <arvind.sudarsanam at aggiemail.usu.edu> wrote: >> >> Hi, >> >> I am trying to compile a single module (dgemm.c) using clang and >> generate dgemm.ll. >> Command: clang -emit-llvm dgemm.c -S -o dgemm.ll >> Some of the functions are declared as static and clang ignores these >> functions: >> >> One of the functions is: >> static void innerloop(double a, const double...

If I configure in R-patched and R-devel, with the same options to configure, then R-patched says checking for ATL_xerbla in -latlas... yes checking for cblas_dgemm in -lcblas... yes checking for dgemm_ in -lf77blas... yes which R-devel says checking for ATL_xerbla in -latlas... yes checking for cblas_dgemm in -lcblas... yes checking for dgemm in -lf77blas... no This is on the same computer and using the same libraries. Observe in R-patched we check for dge...

...386-pc-linux-gnu" define i32 @use() nounwind { [...] } define internal i32 @add(i32 %a, i32 %b) nounwind { [...] } Olivier. On Mon, Dec 14, 2009 at 8:49 PM, Arvind Sudarsanam < arvind.sudarsanam at aggiemail.usu.edu> wrote: > Hi, > > I am trying to compile a single module (dgemm.c) using clang and > generate dgemm.ll. > Command: clang -emit-llvm dgemm.c -S -o dgemm.ll > Some of the functions are declared as static and clang ignores these > functions: > > One of the functions is: > static void innerloop(double a, const double* b, double* c, double beta)...

...t; > [...] > > } > > > > Olivier. > > > > > > On Mon, Dec 14, 2009 at 8:49 PM, Arvind Sudarsanam > > <arvind.sudarsanam at aggiemail.usu.edu> wrote: > >> > >> Hi, > >> > >> I am trying to compile a single module (dgemm.c) using clang and > >> generate dgemm.ll. > >> Command: clang -emit-llvm dgemm.c -S -o dgemm.ll > >> Some of the functions are declared as static and clang ignores these > >> functions: > >> > >> One of the functions is: > >> static vo...

...41:42 +0000 writes: > I am using R to multiply some large (30k x 30k double) matrices on a > 64 core machine (xeon phi). I added some timers to src/main/array.c > to see where the time is going. All of the time is being spent in the > matprod function, most of that time is spent in dgemm. 15 seconds is > in matprod in some code that is checking if there are NaNs. > > system.time (C <- B %*% A) > nancheck: wall time 15.240282s > dgemm: wall time 43.111064s > matprod: wall time 58.351572s > user system elapsed > 2710.154 20.999 58.398 >...

...:42 +0000 writes: > I am using R to multiply some large (30k x 30k double) > matrices on a 64 core machine (xeon phi). I added some timers > to src/main/array.c to see where the time is going. All of the > time is being spent in the matprod function, most of that time > is spent in dgemm. 15 seconds is in matprod in some code that > is checking if there are NaNs. > > system.time (C <- B %*% A) > nancheck: wall time 15.240282s > dgemm: wall time 43.111064s > matprod: wall time 58.351572s > user system elapsed > 2710.154 20.999 58.398 >...

Dear all, I would like to submit an issue that we are facing. Indeed, in our environment, we are optimizing the R code to speed up some mathematical calculations as matrix products using the INTEL libraries ( MKL) ( https://software.intel.com/en-us/mkl ) With the last version of the MKL libraries Intel 2018, we are facing to an issue with *all INTERNAL command* that are executing in R.

...ites: >> I am using R to multiply some large (30k x 30k double) matrices on a >> 64 core machine (xeon phi). I added some timers to src/main/array.c >> to see where the time is going. All of the time is being spent in the >> matprod function, most of that time is spent in dgemm. 15 seconds is >> in matprod in some code that is checking if there are NaNs. >>> system.time (C <- B %*% A) >> nancheck: wall time 15.240282s >> dgemm: wall time 43.111064s >> matprod: wall time 58.351572s >> user system elapsed >> 271...

Hi, I was wondering if I could use the matprod function from array.c in my own C routine. I tried the following as a test /* my_matprod.c */ # include <Rinternals.h> /* for REAL, SEXP etc */ # include <R_ext/Applic.h> /* array.c says need for dgemm */ /* following copied from array.c */ static void matprod(double *x, int nrx, int ncx, double *y, int nry, int ncy, double *z) { char *transa = "N", *transb = "N"; int i, j, k; double one = 1.0, zero = 0.0, sum; Rboolean have_na = FALSE; if (nrx &gt...

Hi, I am new to LLVM and was trying out various tool of LLVM. I compiled my source file containing a implementation of matrix-matrix multiplication, using llvm-g++ -emit-llvm -O3 -c dgemm.cpp -o dgemm.bc and tried running the source file on the llvm interpreter using lli --force-interpreter=true dgemm.bc I get the error Could not resolve external global address: __dso_handle followed by some dump. Am I missing something? Do I need to include anything else while compiling the...

...writes: >> I am using R to multiply some large (30k x 30k double) matrices on a >> 64 core machine (xeon phi). I added some timers to src/main/array.c >> to see where the time is going. All of the time is being spent in the >> matprod function, most of that time is spent in dgemm. 15 seconds is >> in matprod in some code that is checking if there are NaNs. >>> system.time (C <- B %*% A) >> nancheck: wall time 15.240282s >> dgemm: wall time 43.111064s >> matprod: wall time 58.351572s >> user system elapsed >> 2710...

This is not the correct group to post this question. Sorry, but I subscribe nowhere else. What is the BLAS routine to multiply a matrix by a scalar? Thanks. -- Mehmet Balcilar, PhD Assistant Professor Manas University College of Economics & Administrative Sciences Department of Economics Prospect Tinctik 56 Bishkek Kyrgyzstan Tel: +996 (312) 54 19 42 +996 (312) 54 19 43 +996

The original report was that compiling with gcc2.95.2 gcc and the related g77 on Irix 6.,5.13 produced an R which wailed its tests with: >> On entry to DGEMM parameter number 1 had an illegal value To some extend this confused me, as the source only seems to mention DGEMM in lowercase. Eventually, after the debug code I added did nothing, I realized that the R which was build was not using the src/appl/blas.f source file: it was using the SGI libb...

...to load shared library '/home/mhitczen/Cstuff/testmore.so': /home/mhitczen/Cstuff/testmore.so: undefined symbol: dpotrf_ This error goes away and everything works when I simply call the BLAS routine (In the testmore.c file I simply remove funct2, leaving funct which calls the routine "dgemm" found in the BLAS.h file). I read the "Writing R Extensions" and created a file: /R-2.9.0/src/Makevars that contains the line: PKG_LIBS = $(LAPACK_LIBS) $(BLAS_LIBS) $(FLIBS) What am I doing wrong? Thanks for your time and effort. My testmore.c file is: #include <Rinternals.h...

Hi all, I have a C function, say Cfun, that calls Lapack's DGEMM routine and I need to create a shared library to use Cfun inside R. The C file is the following #include<stdio.h> #include<R.h> #include<R_ext/Lapack.h> void Cfun(double *res, double *X, int *n, int *q) { char *ptr_TRANSA, TRANSA='T', *ptr_TRANSB, TRANSB='N';...