search for: transa

...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'; ptr_TRANSA=&TRANSA; ptr_TRANSB=&TRANSB; double *ptr_ALPHA, ALPHA=1.0, *ptr_BETA, BETA=0.0; ptr_ALPHA=&ALPHA; ptr_BETA=&BETA; dgemm_(ptr_TRANSA, ptr_TRANSB, q, q, n, ptr_ALPHA, X, n, X, n, ptr_BETA, res, q); } In Ubuntu I...

...g 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 > 0 && ncx > 0 && nry > 0 && ncy > 0) { /* Don't trust the BLAS to handle NA/NaNs correctly: PR#4582 * The test is...

...(I am *new* to writing external C-functions, so feel free to laugh at my code - or, perhaps, suggest changes): #include <Rinternals.h> #include <R_ext/Applic.h> /* for dgemm */ static void matprod(double *x, int nrx, int ncx, double *y, int nry, int ncy, double *z) { char *transa = "N", *transb = "N"; double one = 1.0, zero = 0.0; F77_CALL(dgemm)(transa, transb, &nrx, &ncy, &ncx, &one, x, &nrx, y, &nry, &zero, z, &nrx); } SEXP trProd2(SEXP x, SEXP y) { int nrx, ncx, nry, ncy, mode, i; SEXP xdims, ydims,...

...ons, so > feel free to laugh at my code - or, perhaps, suggest changes): > > #include <Rinternals.h> > #include <R_ext/Applic.h> /* for dgemm */ > > static void matprod(double *x, int nrx, int ncx, > double *y, int nry, int ncy, double *z) { > char *transa = "N", *transb = "N"; > double one = 1.0, zero = 0.0; > F77_CALL(dgemm)(transa, transb, &nrx, &ncy, &ncx, &one, > x, &nrx, y, &nry, &zero, z, &nrx); } > > SEXP trProd2(SEXP x, SEXP y) > { > int nrx, ncx, nry,...

...#include <Rmath.h> #include <R.h> #include <R_ext/Applic.h> #include <R_ext/Lapack.h> #include <R_ext/BLAS.h> #include <R_ext/RS.h> // FUNCTION CALLING BLAS void funct(double *x, int *dim, double *y, double *result){ double one=1.0; double zero =0.0; char *transa="N"; char *transb="N"; F77_CALL(dgemm) (transa,transb,dim,dim,dim,&one,x,dim,y,dim,&zero,result,dim); } // FUNCTION CALLING LAPACK void funct4(double *x, int *dim, double *result){ char *UorL="L"; int info=0; F77_CALL(dpotrf)(UorL,dim,x,dim, &in...

...main code routine more than 100 times, it crashes R. Or if I call .Call("foobar"....) seperately more than two tims it crashes R. For the most part I am doing matirx multiplies using EXP matrixprod (SEXP x , SEXP y ) { int nrx , ncx , nry , ncy , mode; SEXP xdims , ydims , ans ; char *transa = "N" , *transb = "N" ; double one = 1.0 , zero = 0.0 ; xdims = getAttrib (x , R_DimSymbol ) ; ydims = getAttrib (y , R_DimSymbol ) ; mode = REALSXP; nrx = INTEGER( xdims ) [ 0 ] ; ncx = INTEGER( xdims ) [ 1 ] ; nry = INTEGER( ydims ) [ 0 ] ; ncy = INTEGER( ydims ) [ 1 ] ; PROTE...

...on of the matprod function in src/main/array.c. > > ? Radford Neal > > ------------------------------------------------------------------------ > > 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; > ? ?LDOUBLE sum; > ? ?Rboolean do_it_here; > > ? ?/* NOTE: ncx must be equal to nry. */ > > ? ?if (nrx==1 && ncy==1) { > > ? ? ? ?/* Do dot product quickly. */ >...

...ns regarding compilation. Below is my modified version of the matprod function in src/main/array.c. Radford Neal ------------------------------------------------------------------------ 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; LDOUBLE sum; Rboolean do_it_here; /* NOTE: ncx must be equal to nry. */ if (nrx==1 && ncy==1) { /* Do dot product quickly. */ sum = 0.0; for (i = 0; i&lt...

...rmation is currently useless. In the original version of the package "omxDGEMM" is defined as inline. I changed this in the hope that it would resolve the issue. This was not the case. "omxunsafedgemm" is called by "omxDGEMM" like that F77_CALL(omxunsafedgemm)(&transa, &transb, &(nrow), &(ncol), &(nmid), &alpha, a->data, &(a->leading), b->data, &(b->leading),&beta, result->data, &(result->leading)); Any ideas how to obtain a sen...

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

...;t work. Does anyone have a clue? Thanks Using: Debian 3.0 Kernel 2.4,27 Samba 3.0.22 Ldap 2.1.30-2 -- ___________________________________________________________________________________________________ Atenciosamente, Daniel Felipe Martin GetNet - Tecnologia em Captura e Processamento de Transa??es Infra-estrutura Fone: +55 (51) 3598-9800 Fax: +55 (51) 3598-9801 Ramal 2301 __________________________________________________________________________________________________

I set up a new windows machine and joined it to my existing domain (two Samba DCs). Now it turns out there is a computer in AD, but the associated DNS A/AAAA records of the new machine are missing. I already tried to leave and rejoin, but that didn?t fix the issue. Due to the long ACL I don?t really like the idea of creating the records manually. How can I fix or diagnose that situation?

...08:54 0:00 \_ /usr/sbin/smbd -D user3 32292 0.0 2.5 35092 26056 ? S 09:03 0:02 \_ /usr/sbin/smbd -D root 9788 0.0 0.3 9180 3496 ? S 16:01 0:00 \_ /usr/sbin/smbd -D ---> Excerpt from samba.log [2006/07/12 14:51:21, 3] smbd/process.c:process_smb(1110) Transaction 429419 of length 332 [2006/07/12 14:51:21, 3] smbd/process.c:switch_message(914) switch message SMBtrans (pid 31842) conn 0x84b3a98 [2006/07/12 14:51:21, 3] smbd/sec_ctx.c:set_sec_ctx(241) setting sec ctx (1008, 100) - sec_ctx_stack_ndx = 0 [2006/07/12 14:51:21, 3] smbd/ipc.c:handle_trans(...

Hello, Neven! > i have red your inquiery about slow samba when RO operation is happening on the same file. We have the same problem with SAMBA on SOLARIS. Did you maybe solved the problem somehow ? > Unfortunately, I could not solve this problem. I tweaked many SAMBA's parameters but nothing could help me with this problem. I did not get ANY reaction either from Samba mail list or