Displaying 3 results from an estimated 3 matches for "lptr".
Did you mean:
ldtr
2007 Dec 31
2
the woes of NA
Joyful.
I'm adapting a FORTRAN 77 package for use with R. Pretty
straightforward.
Except for a glitch it took me some time to figure out. This existing
package has subroutines which have parameters called "NA". So, I called
subroutines like
bnodes <- function(n, lst, lptr, lend, nodes, nb, na, nt)
{
ensure.all.numeric(list(n, lst, lptr, lend, nodes, nb, n.a, nt),
"all arguments to -bnodes- must be numeric")
out <- .Fortran("bnodes", N=as.integer(n), LIST=as.integer(lst),
LPTR=as.integer(lptr), LEND=a...
2012 Dec 15
3
interfacing with .Call
...*modif, int *dimif, double *Lsum)
{ ...}
//* the interface part *//
#define getDim(A) INTEGER(coerceVector(getAttrib(A,R_DimSymbol), INTSXP))
SEXP Projector5(SEXP L, SEXP G, SEXP W, SEXP xymod, SEXP modif)
{
//* digest SEXPs from R *//
int *dimL, *dimG, *dimW, *dimxy, *dimif;
double *lptr, *gptr, *wptr, *ifptr;
int *xyptr;
dimL=getDim(L);
PROTECT(L=coerceVector(L, REALSXP));
lptr=REAL(L);
dimG=getDim(G);
PROTECT(G=coerceVector(G, REALSXP));
gptr=REAL(G);
dimW=getDim(W);
PROTECT(W=coerceVector(W, REALSXP));
wptr=REAL(W);
dimxy=getDim(xymod)...
2012 Dec 06
1
Use .Call interface
...t;Projector.h"
SEXP Projector2(SEXP L, SEXP G, SEXP W, SEXP xymod, SEXP modif){
int nprot=0;
PROTECT(L=AS_NUMERIC(L));nprot++;
PROTECT(G=AS_NUMERIC(G));nprot++;
PROTECT(W=AS_NUMERIC(W));nprot++;
PROTECT(xymod=AS_INTEGER(xymod));nprot++;
PROTECT(modif=AS_NUMERIC(modif));nprot++;
double *lptr; lptr=REAL(L);
double *gptr; gptr=REAL(G);
double *wptr; wptr=REAL(W);
int *xyptr; xyptr=INTEGER(xymod);
double *ifptr; ifptr=REAL(modif);
int *dimL; dimL=INTEGER(GET_DIM(L));
int *dimG; dimG=INTEGER(GET_DIM(G));
int *dimW; dimW=INTEGER(GET_DIM(W));
int *dimxy; dimxy=INTEGER(GET_DIM(xymod))...