Hello,
This is not the appropriate mailing list. Use R-devel for questions
about C, etc ...
One thing that might help you is the inline package.
require( inline )
fx <- cfunction( signature( s = "numeric" ), '
SEXP result;
PROTECT(result = NEW_NUMERIC(1));
double* ptr=NUMERIC_POINTER(result);
double t = *REAL(s);
double u = t-floor(t)-0.5;
if(u>0) *ptr=-1+4*u; else *ptr=-1-4*u;
Rprintf("The value is %f", *ptr);
UNPROTECT(1);
return result;
', verbose = TRUE )
fx( 10 )
The verbose = TRUE argument will show you how inline runs the show.
Romain
Le 18/06/10 16:18, michael meyer a écrit :
Greetings,
I am trying to call simple C-code from R.
I am on Windows XP with RTools installed.
The C-function is
#include<R.h>
#include<Rinternals.h>
#include<Rmath.h>
#include<Rdefines.h>
// prevent name mangling
extern "C" {
SEXP __cdecl test(SEXP s){
SEXP result;
PROTECT(result = NEW_NUMERIC(1));
double* ptr=NUMERIC_POINTER(result);
double t = *REAL(s);
double u = t-floor(t)-0.5;
if(u>0) *ptr=-1+4*u; else *ptr=-1-4*u;
Rprintf("The value is %f", *ptr);
UNPROTECT(1);
return result;
}
};
It is compiled with
R CMD SHLIB source.c
If you want C++, then name your file source.cpp
SHLIB compiles files with .c extensions with a C compiler, which will
not be happy about extern "C"
with flag
MAKEFLAGS="CC=g++"
If I compile with the default flags I get an error message about an
undefined reference to "__gxx_personality_v0".
However when I call this code from R with
test<- function(t){
.Call("test",t)
}
dyn.load("./source.dll")
test(0)
dyn.unload("./source.dll")
then R crashes.
I have a vague idea of the issue of calling conventions and was hoping that
the __cdecl
specifier would force the appropriate convention.
I also have Cygwin installed as part of the Python(x,y) distribution but I
am assuming that
R CMD SHLIB source.c
calls the right compiler.
What could the problem be?
Many thanks,
Michael
--
Romain Francois
Professional R Enthusiast
+33(0) 6 28 91 30 30
http://romainfrancois.blog.free.fr
|- http://bit.ly/98Uf7u : Rcpp 0.8.1
|- http://bit.ly/c6YnCi : graph gallery collage
`- http://bit.ly/bZ7ltC : inline 0.3.5
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.