<p><br>
On May 18, 2011 2:53 AM, <<a href="mailto:soeren.vogel@uzh.ch">soeren.vogel@uzh.ch</a>> wrote:<br>
><br>
> Hello Dirk<br>
><br>
> That was exactly the trick! Thanks a lot! Our class and code work so far, but it may take more time to debug or find bad design. Perhaps we will be back with more detailed questions. Anyway, one question already emerged. When librarying the class, assigning an object, do some calculations, whatever -- everything goes fine. Yet, when using R CMD BATCH on a file listing the same commands as by hand, R produces some nice information at the end, however, here it produces an segmentation fault, which does -- according to my small knowledge -- say something bad:<br>
><br>
> > proc.time()<br>
> user system elapsed<br>
> 0.702 0.017 0.714<br>
> R(72992) malloc: *** error for object 0x2c9b604: incorrect checksum for freed object - object was probably modified after being freed.<br>
> *** set a breakpoint in malloc_error_break to debug<br>
><br>
> *** caught segfault ***<br>
> address 0x4909d58, cause 'memory not mapped'<br>
><br>
> Traceback:<br>
> 1: save(list = ls(envir = .GlobalEnv, all.names = TRUE), file = outfile, version = version, ascii = ascii, compress = compress, envir = .GlobalEnv, precheck = FALSE)<br>
> 2: save.image(name)<br>
> 3: sys.save.image(".RData")<br>
> aborting ...<br>
><br>
> Any idea what goes wrong there?<br>
><br>
> Regards<br>
> Sören + Carlo<br>
><br>
As you can see, the problem occurs when R is saving the worksheet. At present an instance of a Rcpp module class cannot be saved. It is likely that this will need to be addressed by the person designing each C++ class that will be expressed in a module. R can't serialize the contents of memory that it doesn't "own".</p>
<p>The current fix is "don't do that". Add the --no-save or --vanilla in your call to R CMD BATCH to suppress saving the worksheet.</p>
<p>> On 14.05.2011, at 18:17, Dirk Eddelbuettel wrote:<br>
><br>
> Soeren,<br>
><br>
> On 13 May 2011 at 20:07, <a href="mailto:soeren.vogel@uzh.ch">soeren.vogel@uzh.ch</a> wrote:<br>
> | Compilation with R CMD CHECK FOO fails with the following error in 00install.out:<br>
> |<br>
> | Error in dyn.load(file, DLLpath = DLLpath, ...) :<br>
> | unable to load shared object '/Users/sovo/GUTS/FOO.Rcheck/FOO/libs/i386/FOO.so':<br>
> | dlopen(/Users/sovo/GUTS/FOO.Rcheck/FOO/libs/i386/FOO.so, 6): Symbol not found: __ZN3FOOC1Ev<br>
> | Referenced from: /Users/sovo/GUTS/FOO.Rcheck/FOO/libs/i386/FOO.so<br>
> | Expected in: flat namespace<br>
> |<br>
> | The (original) class and its functions compile fine with R CMD SHLIB. So we guess that this error has something to do with the Rcpp modules implementation.<br>
><br>
> It took me a few minutes with trial and error, but essentially the mistake<br>
> seems to have been this:<br>
><br>
> | > | FOO();<br>
> | > | ~FOO();<br>
><br>
> You declared the constructor and deconstructor but there was no code<br>
> anywhere. A pair of empty braces ' {}; ' is all it takes.<br>
><br>
> In the process, I renamed things from 'foo' (just to be plain) to 'Baz'. So<br>
> the file is now<br>
><br>
> -----------------------------------------------------------------------------<br>
> // baz_module.cpp<br>
> #include <Rcpp.h><br>
><br>
> // from Baz.h<br>
> class Baz<br>
> {<br>
> private:<br>
> double dtau;<br>
> std::vector<double> D, ee, ff, S;<br>
><br>
> public:<br>
> int M;<br>
> std::vector<double> C, t, s, par;<br>
> std::vector<int> y;<br>
><br>
> Baz() {};<br>
> ~Baz() {};<br>
><br>
> double do_bar(std::vector<double> z);<br>
> };<br>
><br>
> // from Baz.cpp<br>
> double Baz::do_bar(std::vector<double> z) {<br>
> // whatever it does<br>
> return 42.0; // need to return something here<br>
> }<br>
><br>
> RCPP_MODULE(baz){<br>
> using namespace Rcpp ;<br>
> class_<Baz>( "Baz" )<br>
> .constructor()<br>
> .field( "M" , &Baz::M )<br>
> .field( "C" , &Baz::C )<br>
> .method( "do_bar", &Baz::do_bar )<br>
> ;<br>
> }<br>
> -----------------------------------------------------------------------------<br>
><br>
> and with the line<br>
><br>
> RcppModules: yada, baz<br>
><br>
> we get this module loaded:<br>
><br>
><br>
> R> library(mypackage)<br>
> R> Baz<br>
> C++ class 'Baz' <0x1c37db0><br>
> Constructors:<br>
> Baz()<br>
><br>
> Fields:<br>
> std::vector<double, std::allocator<double> > C<br>
> int M<br>
><br>
> Methods:<br>
> double do_bar(std::vector<double, std::allocator<double> >)<br>
><br>
> R><br>
> R> bb <- new(Baz)<br>
> R> bb<br>
> C++ object <0x18c8f90> of class 'Baz' <0x1c37db0><br>
> R> bb$M <- 12<br>
> R> print(bb$M)<br>
> [1] 12<br>
> R> print(bb$C)<br>
> numeric(0)<br>
> R> bb$C <- seq(1.1, 5.5, by=1.1)<br>
> R> bb$C<br>
> [1] 1.1 2.2 3.3 4.4 5.5<br>
> R> R> print(bb$do_bar(1:4))<br>
> [1] 42<br>
> R><br>
><br>
><br>
> I hope you can take it from here.<br>
><br>
> Cheers, Dirk<br>
><br>
> --<br>
> Gauss once played himself in a zero-sum game and won $50.<br>
> -- #11 at <a href="http://www.gaussfacts.com">http://www.gaussfacts.com</a><br>
><br>
> _______________________________________________<br>
> Rcpp-devel mailing list<br>
> <a href="mailto:Rcpp-devel@lists.r-forge.r-project.org">Rcpp-devel@lists.r-forge.r-project.org</a><br>
> <a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel</a><br>
</p>