[Rcpp-devel] Printing intermediate values in C++ code used in Rcpp

Christofer Bogaso bogaso.christofer at gmail.com
Sat Dec 17 18:45:13 CET 2022


Nevermind. Please ignore my previous email. I followed your remaining
instructions and can print the statement now.

Thanks and regards,

On Sat, Dec 17, 2022 at 11:04 PM Christofer Bogaso
<bogaso.christofer at gmail.com> wrote:
>
> Hi,
>
> Thanks for your reply.
>
> Can you please help which setup I need to change?
>
> My Rcpp appears to be latest Rcpp_1.0.9, as I get from session
> information as below,
>
> > sessionInfo()
>
> R version 4.2.1 (2022-06-23)
>
> Platform: x86_64-apple-darwin17.0 (64-bit)
>
> Running under: macOS Big Sur ... 10.16
>
>
> Matrix products: default
>
> BLAS:   /Library/Frameworks/R.framework/Versions/4.2/Resources/lib/libRblas.0.dylib
>
> LAPACK: /Library/Frameworks/R.framework/Versions/4.2/Resources/lib/libRlapack.dylib
>
>
> locale:
>
> [1] C/UTF-8/C/C/C/C
>
>
> attached base packages:
>
> [1] stats     graphics  grDevices utils     datasets  methods   base
>
>
> other attached packages:
>
> [1] GCPM_1.2.2     devtools_2.4.5 usethis_2.1.6  copula_1.1-1   dplyr_1.0.9
>
>
> loaded via a namespace (and not attached):
>
>  [1] tidyselect_1.1.2    remotes_2.4.2       purrr_0.3.4
>
>  [4] lattice_0.20-45     pcaPP_2.0-3         vctrs_0.4.1
>
>  [7] generics_0.1.3      miniUI_0.1.1.1      htmltools_0.5.3
>
> [10] stats4_4.2.1        utf8_1.2.2          rlang_1.0.4
>
> [13] pkgbuild_1.3.1      urlchecker_1.0.1    pillar_1.8.1
>
> [16] later_1.3.0         glue_1.6.2          DBI_1.1.3
>
> [19] sessioninfo_1.2.2   lifecycle_1.0.1     stringr_1.4.1
>
> [22] pspline_1.0-19      htmlwidgets_1.5.4   mvtnorm_1.1-3
>
> [25] memoise_2.0.1       callr_3.7.1         fastmap_1.1.0
>
> [28] httpuv_1.6.6        ps_1.7.1            ADGofTest_0.3
>
> [31] parallel_4.2.1      fansi_1.0.3         Rcpp_1.0.9
>
> [34] xtable_1.8-4        promises_1.2.0.1    cachem_1.0.6
>
> [37] pkgload_1.3.0       mime_0.12           fs_1.5.2
>
> [40] RcppProgress_0.4.2  digest_0.6.29       stringi_1.7.8
>
> [43] processx_3.7.0      shiny_1.7.2         gsl_2.1-7.1
>
> [46] numDeriv_2016.8-1.1 grid_4.2.1          stabledist_0.7-1
>
> [49] cli_3.4.1           tools_4.2.1         magrittr_2.0.3
>
> [52] tibble_3.1.8        profvis_0.3.7       crayon_1.5.1
>
> [55] pkgconfig_2.0.3     ellipsis_0.3.2      Matrix_1.5-3
>
> [58] prettyunits_1.1.1   assertthat_0.2.1    R6_2.5.1
>
> [61] compiler_4.2.1
>
> On Sat, Dec 17, 2022 at 7:45 PM Dirk Eddelbuettel <edd at debian.org> wrote:
> >
> >
> > On 17 December 2022 at 17:29, Christofer Bogaso wrote:
> > | Hi,
> > |
> > | I am working with a package called
> > |
> > | https://cran.r-project.org/src/contrib/GCPM_1.2.2.tar.gz
> > |
> > | The source code contains C++ codes which are available in the src folder.
> > |
> > | In this folder, the C++ codes are available in cpploss.cpp
> > |
> > | In this cpp file, there is a function like below,
> > |
> > | #ifdef _OPENMP
> > | #include <omp.h>
> > | #endif
> > | // [[Rcpp::depends(RcppProgress)]]
> > | #include <progress.hpp>
> > | #include "cpploss.h"
> > | #include <Rcpp.h>
> > | #include <Rmath.h>
> > | #include <iostream>
> > |
> > | using namespace Rcpp;
> > |
> > | // [[Rcpp::export]]
> > | SEXP  GCPM_cpploss(SEXP default_distr_a,SEXP link_function_a, SEXP
> > | S_a,SEXP Sigma_a, SEXP W_a, SEXP PD_a, SEXP PL_a, SEXP calc_rc_a, SEXP
> > | loss_thr_a, SEXP max_entries_a){
> > |   NumericMatrix S(S_a), W(W_a),Sigma(Sigma_a);
> > |   NumericVector PD(PD_a),
> > | PL(PL_a),max_entries(max_entries_a),default_distr(default_distr_a),link_function(link_function_a),calc_rc(calc_rc_a),loss_thr(loss_thr_a);
> > |   List ret;
> > |
> > | etc.
> > |
> > | However as I run this C++ codes, I want to print some intermediate
> > | values generated by underlying C++ codes for some debugging purposes.
> > | So I added below line (just an example)
> > |
> > | Rcpp::Rcout << "Some Value" << std::endl << 1.23 << std::endl;
> > |
> > | After adding this line, then after re-building the package and
> > | re-installng it, when I run the R code, I dont get above line printed
> > | in the R console.
> >
> > Then you are running a different build of the package. The statement, if
> > unconditional, _will_ print.  Check your setup.
> >
> >   > Rcpp::cppFunction('void foo() { Rcpp::Rcout << "boo" << std::endl; }')
> >   > foo()
> >   boo
> >   >
> >
> > I am currently working a lot with a new package spdl for just this: logging.
> >
> >   > Rcpp::cppFunction('void spdlfoo() { spdl::debug("Something {} or other {}", 42, "zing"); }', depends="RcppSpdlog", include="#include <spdl.h>")
> >   > spdl::setup("demo", "warn")  # new logger 'demo' at level 'warn'
> >   > spdlfoo()   # not shown as 'debug' < 'warn'
> >   > spdl::set_level("debug")  # lower the level
> >   > spdlfoo()   # now shown
> >   [2022-12-17 08:13:50.491] [demo] [Process: 3817047] [debug] Something 42 or other zing
> >   >
> >
> > Of course you normally want that in a package. Import spdl, and LinkingTo
> > RcppSpdlog (which came first) does that, more in the vignette of RcppSpdlog.
> >
> > | Could you please help how can I get designated intermadiate values
> > | printed in my R console when I run the R/C++ code?
> > |
> > | Additionally, what is the meaning of the statements like NumericVector
> > | PD(PD_a) etc? I understand that PD_a is the function argument. But
> > | what is the meaning of PD(PD_a)?
> >
> > Instantiate a new object named 'PD' from an object named 'PD_a'.
> >
> > Dirk
> >
> > --
> > dirk.eddelbuettel.com | @eddelbuettel | edd at debian.org


More information about the Rcpp-devel mailing list