search for: isreal

Hi, I am looking for C/C++ codes for computing generalized inverse of a matrix. Can anyone help me in this regard? Thanks, Mahbub.

...it is called with a dataframe from the R side with: .Call("myfunc", somedataframe) On the C side (actually C++ side) I use code like this: SEXP colnames = getAttrib(df, R_NamesSymbol) cname = string(CHAR(STRING_ELT(colnames,i)) SEXP coldata = VECTOR_ELT(df,i) (data for i-th column) if(isReal(colData)) x = REAL(colData)[j]; else if(isInteger(colData)) i = INTEGER(colData)[j]; else if(isString(colData)) s = CHAR(STRING_ELT(colData,j)) The problem is that the last test (isString) never passes, even when I pass in a frame for which one or more cols contain character strings. W...

...lem should be something with the dll generation/linkage in windows but I can't figure out. As a matter of test I did the following C code: #include <R.h> #include <Rinternals.h> SEXP thisisatest(SEXP); SEXP thisisatest(SEXP a) { long int i; if (!isReal(a)) printf("Vector should be double.\n"); for (i=LENGTH(a)-1; i >=0; i--) { REAL(a)[i] = REAL(a)[i] + 1; } return (a); } Linux output: R CMD SHLIB thisisthetest.c gcc -I/usr/lib/R/include -I/usr/local/include...

...2003/Tutorials/RExtensions/slide s.pdf> http://www.ci.tuwien.ac.at/Conferences/DSC-2003/Tutorials/RExtensions/slides .pdf e.g. SEXP La_DGE_dc(SEXP A) { SEXP aa = PROTECT(duplicate(A)); SEXP adims, pivot, val; int m, n, 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(pi...

...LT(pars, i)); j++, start += LENGTH(ans)) { SEXP ans_del; double origPar, xx, delta; origPar = REAL(VECTOR_ELT(pars, i))[j]; xx = fabs(origPar); delta = (xx == 0) ? eps : xx*eps; REAL(VECTOR_ELT(pars, i))[j] += rDir[i] * delta; PROTECT(ans_del = eval(expr, rho)); if(!isReal(ans_del)) ans_del = coerceVector(ans_del, REALSXP); UNPROTECT(1); for(k = 0; k < LENGTH(ans); k++) { if (!R_FINITE(REAL(ans_del)[k])) *error(_("Missing value or an infinity produced when evaluating the model"));* REAL(gradient)[start + k] = rDir[i] * (REAL(a...

Hello, I would like some advice on a SATA RAID controller that is able to work in a 32 bit/33mhz PCI slot. I have been looking for a while but all I can find is 64 bit cards that cost over $300. I just need it to mirror 2 drives, and of course work in Centos 4.2. I am not looking for RAID 5 or anything like that. All the cards I find that would do this are fakeraid. Also, has anyone used

Dear Samba List, My Windows XP machine is unable to reconnect my Linux samba network drive when I login. I installed: Win2000_PlainPassword.reg which allows me to connect manually, but it is a requirement that this be automatic on each login. Has anyone found a way to achieve this? Thanks, Dick

Dear R core team Using the following code produces an empty plot (similar to col = NA): > plot(1:9, col = factor(rep(1:3,3), labels = c("red", "blue", "black"))) My question: Shouldn't one get at least a warning (or an error) if one tries to use a factor as col argument? Thanks for an answer. Regards, Christoph Buser

Dear R developers, I've run into a weird behavior of the numericDeriv function (from the stats package) which I also posted on StackOverflow (question has same title as this email, except for the version of R). Running the code bellow we can see that the numericDeriv function gives an error as the derivative of x^a wrt a is x^a * log(x) and log is not defined for negative numbers. However,

...filename = STRING_ELT(CAR(args), 0); - threshold = REAL(CADDR(args))[0]; - if (strlen(CHAR(filename))) - R_InitMemReporting(filename, append_mode, threshold); + filename = STRING_ELT(CAR(ap), 0); + + ap = CDR(ap); + append_mode = asLogical(CAR(ap)); + + ap = CDR(ap); + if (!isReal(CAR(ap)) || (LENGTH(CAR(ap))) != 1) + error(_("invalid '%s' argument"), "threshold"); + R_MemReportingThreshold = REAL(CAR(ap))[0]; + + ap = CDR(ap); + if (!isReal(CAR(ap)) || (LENGTH(CAR(ap))) != 1) + error(_("invalid '%s' argument"), "ne...

...y have to add new data structure return true; } In above function we are calling a new function TII.isReMaterializablePossible(*DefMI, aa) which will act as early heuristic and return false by checking if instruction is definitely not rematerialize. We have found some cases from TargetInstrInfo::isReallyTriviallyReMaterializableGeneric and code for same is here: bool TargetInstrInfo::isReMaterializablePossible( const MachineInstr &MI, AliasAnalysis *AA) const { const MachineFunction &MF = *MI.getParent()->getParent(); const MachineRegisterInfo &MRI = MF.getRegInfo(); //...

...l make no difference on most platforms: -fPIC is needed on Sparc (and still used there), but -fpic should give slightly better performance on PowerPC (although -fPIC is used on PPC64 as it is needed to build libR.so there). o More use is made of inlining for small utility functions such as isReal. Because this can only be done portably with C99 constructs (and we know of no actual implementation), this is only done for the GNU C compiler. o There is an experimental feature to allow shared installations of sub-architectures. See the R-admin manual. o All platforms now use R...

...l make no difference on most platforms: -fPIC is needed on Sparc (and still used there), but -fpic should give slightly better performance on PowerPC (although -fPIC is used on PPC64 as it is needed to build libR.so there). o More use is made of inlining for small utility functions such as isReal. Because this can only be done portably with C99 constructs (and we know of no actual implementation), this is only done for the GNU C compiler. o There is an experimental feature to allow shared installations of sub-architectures. See the R-admin manual. o All platforms now use R...