R Under development (unstable) (2024-07-20 r86909 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2024 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> if (requireNamespace("RUnit", quietly=TRUE) && requireNamespace("RcppXPtrUtils", quietly=TRUE) && requireNamespace("r2sundials", quietly=TRUE) && requireNamespace("RcppArmadillo", quietly=TRUE)) {
+    library(RUnit)
+    library(RcppXPtrUtils)
+    library(r2sundials)
+    library(RcppArmadillo)
+ 
+    testSuite <- defineTestSuite(
+       name = "r2sundials unit tests",
+       dirs = system.file("unitTests", package = "r2sundials"),
+       testFuncRegexp = "^[Tt]est.+"
+    )
+    Sys.setenv("R_TESTS"="")
+    tests <- runTestSuite(testSuite)
+ 
+    printTextProtocol(tests)
+ 
+    if (getErrors(tests)$nFail > 0) stop("RUnit test failure")
+    if (getErrors(tests)$nErr > 0) stop("Errors in RUnit tests")
+ }
Loading required package: rmumps

Generated code for function definition: 
--------------------------------------------------------

// [[Rcpp::depends(RcppArmadillo, r2sundials, rmumps)]]

#include <RcppArmadillo.h>
#include <Rcpp.h>
#include <r2sundials.h>
#include <rmumps.h>

using namespace Rcpp;

// [[Rcpp::plugins(cpp14)]]
using namespace arma;
#include <r2sundials.h>
// [[Rcpp::export]]
SEXP getXPtr();

int d_robertson(double t, const vec& y, vec& ydot, RObject& param, NumericVector& psens) {
  NumericVector p(param);
  ydot[0] = -p["k1"]*y[0] + p["k3"]*y[1]*y[2];
  ydot[2] = p["k2"]*y[1]*y[1];
  ydot[1] = -ydot[0] - ydot[2];
  return(CV_SUCCESS);
}

SEXP getXPtr() {
  typedef int (*funcPtr)( double t, const vec& y, vec& ydot, RObject& param, NumericVector& psens );
  return XPtr<funcPtr>(new funcPtr(& d_robertson ));
}

Generated extern "C" functions 
--------------------------------------------------------


#include <Rcpp.h>
#ifdef RCPP_USE_GLOBAL_ROSTREAM
Rcpp::Rostream<true>&  Rcpp::Rcout = Rcpp::Rcpp_cout_get();
Rcpp::Rostream<false>& Rcpp::Rcerr = Rcpp::Rcpp_cerr_get();
#endif

// getXPtr
SEXP getXPtr();
RcppExport SEXP sourceCpp_1_getXPtr() {
BEGIN_RCPP
    Rcpp::RObject rcpp_result_gen;
    Rcpp::RNGScope rcpp_rngScope_gen;
    rcpp_result_gen = Rcpp::wrap(getXPtr());
    return rcpp_result_gen;
END_RCPP
}

Generated R functions 
-------------------------------------------------------

`.sourceCpp_1_DLLInfo` <- dyn.load('D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b439af5ba9/sourceCpp_2.dll')

getXPtr <- Rcpp:::sourceCppFunction(function() {}, FALSE, `.sourceCpp_1_DLLInfo`, 'sourceCpp_1_getXPtr')

rm(`.sourceCpp_1_DLLInfo`)

Building shared library
--------------------------------------------------------

DIR: D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b439af5ba9

D:/RCompile/recent/R/bin/x64/R CMD SHLIB -o "sourceCpp_2.dll" "file2f9b434e83842.cpp" 
using C++ compiler: 'g++.exe (GCC) 13.2.0'
using C++14
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b439af5ba9'
g++  -std=gnu++14 -I"D:/RCompile/recent/R/include" -DNDEBUG -I../inst/include -fopenmp  -I"D:/RCompile/CRANpkg/lib/4.5/Rcpp/include" -I"D:/RCompile/CRANpkg/lib/4.5/RcppArmadillo/include" -I"D:/RCompile/CRANincoming/R-devel/lib/r2sundials/include" -I"D:/RCompile/CRANpkg/lib/4.5/rmumps/include" -I"D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13"   -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include"     -O2 -Wall  -mfpmath=sse -msse2 -mstackrealign  -c file2f9b434e83842.cpp -o file2f9b434e83842.o
g++ -shared -s -static-libgcc -o sourceCpp_2.dll tmp.def file2f9b434e83842.o -fopenmp -LD:/RCompile/recent/R/bin/x64 -lRlapack -LD:/RCompile/recent/R/bin/x64 -lRblas -lgfortran -lm -lquadmath -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib/x64 -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib -LD:/RCompile/recent/R/bin/x64 -lR
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b439af5ba9'
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b439af5ba9'
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b439af5ba9'

Generated code for function definition: 
--------------------------------------------------------

// [[Rcpp::depends(RcppArmadillo, r2sundials, rmumps)]]

#include <RcppArmadillo.h>
#include <Rcpp.h>
#include <r2sundials.h>
#include <rmumps.h>

using namespace Rcpp;

// [[Rcpp::plugins(cpp14)]]
using namespace arma;
#include <r2sundials.h>
// [[Rcpp::export]]
SEXP getXPtr();

int jac_robertson(double t, const vec& y, const vec& ydot, mat& J, RObject& param, NumericVector& psens, vec& tmp1, vec& tmp2, vec& tmp3) {
  NumericVector p(param);
  J(0, 0) = -p["k1"];
  J(1, 0) = p["k1"]; 
  J(2, 0) = 0.;

  J(0, 1) = p["k3"]*y[2];
  J(1, 1) = -p["k3"]*y[2]-2*p["k2"]*y[1];
  J(2, 1) = 2*p["k2"]*y[1];

  J(0, 2) = p["k3"]*y[1];
  J(1, 2) = -p["k3"]*y[1];
  J(2, 2) = 0.;
  return(CV_SUCCESS);
}

SEXP getXPtr() {
  typedef int (*funcPtr)( double t, const vec& y, const vec& ydot, mat& J, RObject& param, NumericVector& psens, vec& tmp1, vec& tmp2, vec& tmp3 );
  return XPtr<funcPtr>(new funcPtr(& jac_robertson ));
}

Generated extern "C" functions 
--------------------------------------------------------


#include <Rcpp.h>
#ifdef RCPP_USE_GLOBAL_ROSTREAM
Rcpp::Rostream<true>&  Rcpp::Rcout = Rcpp::Rcpp_cout_get();
Rcpp::Rostream<false>& Rcpp::Rcerr = Rcpp::Rcpp_cerr_get();
#endif

// getXPtr
SEXP getXPtr();
RcppExport SEXP sourceCpp_3_getXPtr() {
BEGIN_RCPP
    Rcpp::RObject rcpp_result_gen;
    Rcpp::RNGScope rcpp_rngScope_gen;
    rcpp_result_gen = Rcpp::wrap(getXPtr());
    return rcpp_result_gen;
END_RCPP
}

Generated R functions 
-------------------------------------------------------

`.sourceCpp_3_DLLInfo` <- dyn.load('D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4468d6d9b/sourceCpp_4.dll')

getXPtr <- Rcpp:::sourceCppFunction(function() {}, FALSE, `.sourceCpp_3_DLLInfo`, 'sourceCpp_3_getXPtr')

rm(`.sourceCpp_3_DLLInfo`)

Building shared library
--------------------------------------------------------

DIR: D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4468d6d9b

D:/RCompile/recent/R/bin/x64/R CMD SHLIB -o "sourceCpp_4.dll" "file2f9b468e875d3.cpp" 
using C++ compiler: 'g++.exe (GCC) 13.2.0'
using C++14
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4468d6d9b'
g++  -std=gnu++14 -I"D:/RCompile/recent/R/include" -DNDEBUG -I../inst/include -fopenmp  -I"D:/RCompile/CRANpkg/lib/4.5/Rcpp/include" -I"D:/RCompile/CRANpkg/lib/4.5/RcppArmadillo/include" -I"D:/RCompile/CRANincoming/R-devel/lib/r2sundials/include" -I"D:/RCompile/CRANpkg/lib/4.5/rmumps/include" -I"D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13"   -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include"     -O2 -Wall  -mfpmath=sse -msse2 -mstackrealign  -c file2f9b468e875d3.cpp -o file2f9b468e875d3.o
g++ -shared -s -static-libgcc -o sourceCpp_4.dll tmp.def file2f9b468e875d3.o -fopenmp -LD:/RCompile/recent/R/bin/x64 -lRlapack -LD:/RCompile/recent/R/bin/x64 -lRblas -lgfortran -lm -lquadmath -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib/x64 -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib -LD:/RCompile/recent/R/bin/x64 -lR
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4468d6d9b'
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4468d6d9b'
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4468d6d9b'

Generated code for function definition: 
--------------------------------------------------------

// [[Rcpp::depends(RcppArmadillo, r2sundials, rmumps)]]

#include <RcppArmadillo.h>
#include <Rcpp.h>
#include <r2sundials.h>
#include <rmumps.h>

using namespace Rcpp;

// [[Rcpp::plugins(cpp14)]]
using namespace arma;
#include <r2sundials.h>
// [[Rcpp::export]]
SEXP getXPtr();

int spjac_robertson(double t, const vec& y, const vec& ydot, uvec& ir, uvec& pj, vec& v, int n, int nz, RObject& param, NumericVector& psens, vec& tmp1, vec& tmp2, vec& tmp3) {
  if (nz < 8)
    stop("spjac_robertson: not enough room for non zeros, must have at least 8, instead got %d", nz);
  NumericVector prm(param);
  int i=0;
  pj[0] = 0; // init pj
  // first column
  ir[i] = 0;
  v[i++] = -prm["k1"];
  ir[i] = 1;
  v[i++] = prm["k1"]; 
  pj[1] = i;
  // second column
  ir[i] = 0;
  v[i++] = prm["k3"]*y[2];
  ir[i] = 1;
  v[i++] = -prm["k3"]*y[2]-2*prm["k2"]*y[1];
  ir[i] = 2;
  v[i++] = 2*prm["k2"]*y[1];
  pj[2] = i;
  // third column
  ir[i] = 0;
  v[i++] = prm["k3"]*y[1];
  ir[i] = 1;
  v[i++] = -prm["k3"]*y[1];
  ir[i] = 2;
  v[i++] = 0; // just to hold the place for a full main diagonal
  pj[3] = i;
  return(CV_SUCCESS);
}

SEXP getXPtr() {
  typedef int (*funcPtr)( double t, const vec& y, const vec& ydot, uvec& ir, uvec& pj, vec& v, int n, int nz, RObject& param, NumericVector& psens, vec& tmp1, vec& tmp2, vec& tmp3 );
  return XPtr<funcPtr>(new funcPtr(& spjac_robertson ));
}

Generated extern "C" functions 
--------------------------------------------------------


#include <Rcpp.h>
#ifdef RCPP_USE_GLOBAL_ROSTREAM
Rcpp::Rostream<true>&  Rcpp::Rcout = Rcpp::Rcpp_cout_get();
Rcpp::Rostream<false>& Rcpp::Rcerr = Rcpp::Rcpp_cerr_get();
#endif

// getXPtr
SEXP getXPtr();
RcppExport SEXP sourceCpp_5_getXPtr() {
BEGIN_RCPP
    Rcpp::RObject rcpp_result_gen;
    Rcpp::RNGScope rcpp_rngScope_gen;
    rcpp_result_gen = Rcpp::wrap(getXPtr());
    return rcpp_result_gen;
END_RCPP
}

Generated R functions 
-------------------------------------------------------

`.sourceCpp_5_DLLInfo` <- dyn.load('D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4233b70e4/sourceCpp_6.dll')

getXPtr <- Rcpp:::sourceCppFunction(function() {}, FALSE, `.sourceCpp_5_DLLInfo`, 'sourceCpp_5_getXPtr')

rm(`.sourceCpp_5_DLLInfo`)

Building shared library
--------------------------------------------------------

DIR: D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4233b70e4

D:/RCompile/recent/R/bin/x64/R CMD SHLIB -o "sourceCpp_6.dll" "file2f9b46daeaed.cpp" 
using C++ compiler: 'g++.exe (GCC) 13.2.0'
using C++14
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4233b70e4'
g++  -std=gnu++14 -I"D:/RCompile/recent/R/include" -DNDEBUG -I../inst/include -fopenmp  -I"D:/RCompile/CRANpkg/lib/4.5/Rcpp/include" -I"D:/RCompile/CRANpkg/lib/4.5/RcppArmadillo/include" -I"D:/RCompile/CRANincoming/R-devel/lib/r2sundials/include" -I"D:/RCompile/CRANpkg/lib/4.5/rmumps/include" -I"D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13"   -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include"     -O2 -Wall  -mfpmath=sse -msse2 -mstackrealign  -c file2f9b46daeaed.cpp -o file2f9b46daeaed.o
g++ -shared -s -static-libgcc -o sourceCpp_6.dll tmp.def file2f9b46daeaed.o -fopenmp -LD:/RCompile/recent/R/bin/x64 -lRlapack -LD:/RCompile/recent/R/bin/x64 -lRblas -lgfortran -lm -lquadmath -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib/x64 -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib -LD:/RCompile/recent/R/bin/x64 -lR
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4233b70e4'
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4233b70e4'
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4233b70e4'

Generated code for function definition: 
--------------------------------------------------------

// [[Rcpp::depends(RcppArmadillo, r2sundials, rmumps)]]

#include <RcppArmadillo.h>
#include <Rcpp.h>
#include <r2sundials.h>
#include <rmumps.h>

using namespace Rcpp;

// [[Rcpp::plugins(cpp14)]]
using namespace arma;
#include <r2sundials.h>
// [[Rcpp::export]]
SEXP getXPtr();

int sens_robertson1(int Ns, double t, const vec& y, const vec& ydot, int iS, const vec& yS, vec& ySdot, RObject& param, NumericVector& p, vec& tmp1, vec& tmp2) {
  // calculate (∂f /∂y)s_i(t) + (∂f /∂p_i) for i = iS
  // (∂f /∂y)s_i(t)
//print(p);
//stop("print");
  ySdot[0] = -p["k1"]*yS[0] + p["k3"]*y[2]*yS[1] + p["k3"]*y[1]*yS[2];
  ySdot[1] = p["k1"]*yS[0] - (p["k3"]*y[2]+2*p["k2"]*y[1])*yS[1] - p["k3"]*y[1]*yS[2]; 
  ySdot[2] = 2*p["k2"]*y[1]*yS[1];
  // + (∂f /∂p_i)
  switch(iS) {
    case 0:
      ySdot[0] -= y[0];
      ySdot[1] += y[0];
      break;
    case 1:
      ySdot[1] -= y[1]*y[1];
      ySdot[2] += y[1]*y[1];
      break;
    case 2:
      ySdot[0] += y[1]*y[2];
      ySdot[1] -= y[1]*y[2];
  }
  return(CV_SUCCESS);
}

SEXP getXPtr() {
  typedef int (*funcPtr)( int Ns, double t, const vec& y, const vec& ydot, int iS, const vec& yS, vec& ySdot, RObject& param, NumericVector& p, vec& tmp1, vec& tmp2 );
  return XPtr<funcPtr>(new funcPtr(& sens_robertson1 ));
}

Generated extern "C" functions 
--------------------------------------------------------


#include <Rcpp.h>
#ifdef RCPP_USE_GLOBAL_ROSTREAM
Rcpp::Rostream<true>&  Rcpp::Rcout = Rcpp::Rcpp_cout_get();
Rcpp::Rostream<false>& Rcpp::Rcerr = Rcpp::Rcpp_cerr_get();
#endif

// getXPtr
SEXP getXPtr();
RcppExport SEXP sourceCpp_7_getXPtr() {
BEGIN_RCPP
    Rcpp::RObject rcpp_result_gen;
    Rcpp::RNGScope rcpp_rngScope_gen;
    rcpp_result_gen = Rcpp::wrap(getXPtr());
    return rcpp_result_gen;
END_RCPP
}

Generated R functions 
-------------------------------------------------------

`.sourceCpp_7_DLLInfo` <- dyn.load('D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4256e31f1/sourceCpp_8.dll')

getXPtr <- Rcpp:::sourceCppFunction(function() {}, FALSE, `.sourceCpp_7_DLLInfo`, 'sourceCpp_7_getXPtr')

rm(`.sourceCpp_7_DLLInfo`)

Building shared library
--------------------------------------------------------

DIR: D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4256e31f1

D:/RCompile/recent/R/bin/x64/R CMD SHLIB -o "sourceCpp_8.dll" "file2f9b4502175ed.cpp" 
using C++ compiler: 'g++.exe (GCC) 13.2.0'
using C++14
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4256e31f1'
g++  -std=gnu++14 -I"D:/RCompile/recent/R/include" -DNDEBUG -I../inst/include -fopenmp  -I"D:/RCompile/CRANpkg/lib/4.5/Rcpp/include" -I"D:/RCompile/CRANpkg/lib/4.5/RcppArmadillo/include" -I"D:/RCompile/CRANincoming/R-devel/lib/r2sundials/include" -I"D:/RCompile/CRANpkg/lib/4.5/rmumps/include" -I"D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13"   -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include"     -O2 -Wall  -mfpmath=sse -msse2 -mstackrealign  -c file2f9b4502175ed.cpp -o file2f9b4502175ed.o
g++ -shared -s -static-libgcc -o sourceCpp_8.dll tmp.def file2f9b4502175ed.o -fopenmp -LD:/RCompile/recent/R/bin/x64 -lRlapack -LD:/RCompile/recent/R/bin/x64 -lRblas -lgfortran -lm -lquadmath -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib/x64 -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib -LD:/RCompile/recent/R/bin/x64 -lR
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4256e31f1'
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4256e31f1'
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4256e31f1'

Generated code for function definition: 
--------------------------------------------------------

// [[Rcpp::depends(RcppArmadillo, r2sundials, rmumps)]]

#include <RcppArmadillo.h>
#include <Rcpp.h>
#include <r2sundials.h>
#include <rmumps.h>

using namespace Rcpp;

// [[Rcpp::plugins(cpp14)]]
using namespace arma;
#include <r2sundials.h>
// [[Rcpp::export]]
SEXP getXPtr();

int d_ball(double t, const vec& y, vec& ydot, RObject& param, NumericVector& psens) {
  NumericVector p(param);
  ydot[0] = y[2];
  ydot[1] = y[3];
  ydot[2] = y[1] > 0 ? 0. : -y[2]; // falling till y=0 then damping
  ydot[3] = y[1] > 0 ? -p["g"] : -y[3]; // falling till y=0 then damping
  return(CV_SUCCESS);
}

SEXP getXPtr() {
  typedef int (*funcPtr)( double t, const vec& y, vec& ydot, RObject& param, NumericVector& psens );
  return XPtr<funcPtr>(new funcPtr(& d_ball ));
}

Generated extern "C" functions 
--------------------------------------------------------


#include <Rcpp.h>
#ifdef RCPP_USE_GLOBAL_ROSTREAM
Rcpp::Rostream<true>&  Rcpp::Rcout = Rcpp::Rcpp_cout_get();
Rcpp::Rostream<false>& Rcpp::Rcerr = Rcpp::Rcpp_cerr_get();
#endif

// getXPtr
SEXP getXPtr();
RcppExport SEXP sourceCpp_9_getXPtr() {
BEGIN_RCPP
    Rcpp::RObject rcpp_result_gen;
    Rcpp::RNGScope rcpp_rngScope_gen;
    rcpp_result_gen = Rcpp::wrap(getXPtr());
    return rcpp_result_gen;
END_RCPP
}

Generated R functions 
-------------------------------------------------------

`.sourceCpp_9_DLLInfo` <- dyn.load('D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b433089ae/sourceCpp_10.dll')

getXPtr <- Rcpp:::sourceCppFunction(function() {}, FALSE, `.sourceCpp_9_DLLInfo`, 'sourceCpp_9_getXPtr')

rm(`.sourceCpp_9_DLLInfo`)

Building shared library
--------------------------------------------------------

DIR: D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b433089ae

D:/RCompile/recent/R/bin/x64/R CMD SHLIB -o "sourceCpp_10.dll" "file2f9b42491ace.cpp" 
using C++ compiler: 'g++.exe (GCC) 13.2.0'
using C++14
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b433089ae'
g++  -std=gnu++14 -I"D:/RCompile/recent/R/include" -DNDEBUG -I../inst/include -fopenmp  -I"D:/RCompile/CRANpkg/lib/4.5/Rcpp/include" -I"D:/RCompile/CRANpkg/lib/4.5/RcppArmadillo/include" -I"D:/RCompile/CRANincoming/R-devel/lib/r2sundials/include" -I"D:/RCompile/CRANpkg/lib/4.5/rmumps/include" -I"D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13"   -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include"     -O2 -Wall  -mfpmath=sse -msse2 -mstackrealign  -c file2f9b42491ace.cpp -o file2f9b42491ace.o
g++ -shared -s -static-libgcc -o sourceCpp_10.dll tmp.def file2f9b42491ace.o -fopenmp -LD:/RCompile/recent/R/bin/x64 -lRlapack -LD:/RCompile/recent/R/bin/x64 -lRblas -lgfortran -lm -lquadmath -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib/x64 -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib -LD:/RCompile/recent/R/bin/x64 -lR
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b433089ae'
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b433089ae'
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b433089ae'

Generated code for function definition: 
--------------------------------------------------------

// [[Rcpp::depends(RcppArmadillo, r2sundials, rmumps)]]

#include <RcppArmadillo.h>
#include <Rcpp.h>
#include <r2sundials.h>
#include <rmumps.h>

using namespace Rcpp;

// [[Rcpp::plugins(cpp14)]]
using namespace arma;
#include <r2sundials.h>
// [[Rcpp::export]]
SEXP getXPtr();

int root_ball(double t, const vec& y, vec& vroot, RObject& param, NumericVector& psens) {
  vroot[0] = y[1]; // y==0
  return(CV_SUCCESS);
}

SEXP getXPtr() {
  typedef int (*funcPtr)( double t, const vec& y, vec& vroot, RObject& param, NumericVector& psens );
  return XPtr<funcPtr>(new funcPtr(& root_ball ));
}

Generated extern "C" functions 
--------------------------------------------------------


#include <Rcpp.h>
#ifdef RCPP_USE_GLOBAL_ROSTREAM
Rcpp::Rostream<true>&  Rcpp::Rcout = Rcpp::Rcpp_cout_get();
Rcpp::Rostream<false>& Rcpp::Rcerr = Rcpp::Rcpp_cerr_get();
#endif

// getXPtr
SEXP getXPtr();
RcppExport SEXP sourceCpp_11_getXPtr() {
BEGIN_RCPP
    Rcpp::RObject rcpp_result_gen;
    Rcpp::RNGScope rcpp_rngScope_gen;
    rcpp_result_gen = Rcpp::wrap(getXPtr());
    return rcpp_result_gen;
END_RCPP
}

Generated R functions 
-------------------------------------------------------

`.sourceCpp_11_DLLInfo` <- dyn.load('D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4139978f9/sourceCpp_12.dll')

getXPtr <- Rcpp:::sourceCppFunction(function() {}, FALSE, `.sourceCpp_11_DLLInfo`, 'sourceCpp_11_getXPtr')

rm(`.sourceCpp_11_DLLInfo`)

Building shared library
--------------------------------------------------------

DIR: D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4139978f9

D:/RCompile/recent/R/bin/x64/R CMD SHLIB -o "sourceCpp_12.dll" "file2f9b419b017fa.cpp" 
using C++ compiler: 'g++.exe (GCC) 13.2.0'
using C++14
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4139978f9'
g++  -std=gnu++14 -I"D:/RCompile/recent/R/include" -DNDEBUG -I../inst/include -fopenmp  -I"D:/RCompile/CRANpkg/lib/4.5/Rcpp/include" -I"D:/RCompile/CRANpkg/lib/4.5/RcppArmadillo/include" -I"D:/RCompile/CRANincoming/R-devel/lib/r2sundials/include" -I"D:/RCompile/CRANpkg/lib/4.5/rmumps/include" -I"D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13"   -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include"     -O2 -Wall  -mfpmath=sse -msse2 -mstackrealign  -c file2f9b419b017fa.cpp -o file2f9b419b017fa.o
g++ -shared -s -static-libgcc -o sourceCpp_12.dll tmp.def file2f9b419b017fa.o -fopenmp -LD:/RCompile/recent/R/bin/x64 -lRlapack -LD:/RCompile/recent/R/bin/x64 -lRblas -lgfortran -lm -lquadmath -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib/x64 -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib -LD:/RCompile/recent/R/bin/x64 -lR
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4139978f9'
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4139978f9'
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4139978f9'

Generated code for function definition: 
--------------------------------------------------------

// [[Rcpp::depends(RcppArmadillo, r2sundials, rmumps)]]

#include <RcppArmadillo.h>
#include <Rcpp.h>
#include <r2sundials.h>
#include <rmumps.h>

using namespace Rcpp;

// [[Rcpp::plugins(cpp14)]]
using namespace arma;
#include <r2sundials.h>
// [[Rcpp::export]]
SEXP getXPtr();

int event_ball(double t, const vec& y, vec& ynew, int Ns, std::vector<vec>& ySv, const ivec& rootsfound, RObject& param, NumericVector& psens) {
  NumericVector p(param);
  static int nbounce=0;
  if (y[3] > 0) // we cross 0 in ascending trajectory, it can happen when y < 0 in limits of abstol
    return(R2SUNDIALS_EVENT_IGNORE);
  ynew=y;
  if (++nbounce < p["nbounce"]) {
    // here nbounce=1:4
    ynew[2] *= 1.-p["kx"]; // horizontal speed is lowered
    ynew[3] *= -(1.-p["ky"]); // vertical speed is lowered and reflected
    return(R2SUNDIALS_EVENT_HOLD);
  } else {
    // here nbounce=5
    nbounce=0; // reinit counter for possible next calls to cvode
    return(R2SUNDIALS_EVENT_STOP);
  }
}

SEXP getXPtr() {
  typedef int (*funcPtr)( double t, const vec& y, vec& ynew, int Ns, std::vector<vec>& ySv, const ivec& rootsfound, RObject& param, NumericVector& psens );
  return XPtr<funcPtr>(new funcPtr(& event_ball ));
}

Generated extern "C" functions 
--------------------------------------------------------


#include <Rcpp.h>
#ifdef RCPP_USE_GLOBAL_ROSTREAM
Rcpp::Rostream<true>&  Rcpp::Rcout = Rcpp::Rcpp_cout_get();
Rcpp::Rostream<false>& Rcpp::Rcerr = Rcpp::Rcpp_cerr_get();
#endif

// getXPtr
SEXP getXPtr();
RcppExport SEXP sourceCpp_13_getXPtr() {
BEGIN_RCPP
    Rcpp::RObject rcpp_result_gen;
    Rcpp::RNGScope rcpp_rngScope_gen;
    rcpp_result_gen = Rcpp::wrap(getXPtr());
    return rcpp_result_gen;
END_RCPP
}

Generated R functions 
-------------------------------------------------------

`.sourceCpp_13_DLLInfo` <- dyn.load('D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b47b261857/sourceCpp_14.dll')

getXPtr <- Rcpp:::sourceCppFunction(function() {}, FALSE, `.sourceCpp_13_DLLInfo`, 'sourceCpp_13_getXPtr')

rm(`.sourceCpp_13_DLLInfo`)

Building shared library
--------------------------------------------------------

DIR: D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b47b261857

D:/RCompile/recent/R/bin/x64/R CMD SHLIB -o "sourceCpp_14.dll" "file2f9b44d0f5d5f.cpp" 
using C++ compiler: 'g++.exe (GCC) 13.2.0'
using C++14
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b47b261857'
g++  -std=gnu++14 -I"D:/RCompile/recent/R/include" -DNDEBUG -I../inst/include -fopenmp  -I"D:/RCompile/CRANpkg/lib/4.5/Rcpp/include" -I"D:/RCompile/CRANpkg/lib/4.5/RcppArmadillo/include" -I"D:/RCompile/CRANincoming/R-devel/lib/r2sundials/include" -I"D:/RCompile/CRANpkg/lib/4.5/rmumps/include" -I"D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13"   -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include"     -O2 -Wall  -mfpmath=sse -msse2 -mstackrealign  -c file2f9b44d0f5d5f.cpp -o file2f9b44d0f5d5f.o
g++ -shared -s -static-libgcc -o sourceCpp_14.dll tmp.def file2f9b44d0f5d5f.o -fopenmp -LD:/RCompile/recent/R/bin/x64 -lRlapack -LD:/RCompile/recent/R/bin/x64 -lRblas -lgfortran -lm -lquadmath -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib/x64 -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib -LD:/RCompile/recent/R/bin/x64 -lR
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b47b261857'
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b47b261857'
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b47b261857'

Generated code for function definition: 
--------------------------------------------------------

// [[Rcpp::depends(RcppArmadillo, r2sundials, rmumps)]]

#include <RcppArmadillo.h>
#include <Rcpp.h>
#include <r2sundials.h>
#include <rmumps.h>

using namespace Rcpp;

// [[Rcpp::plugins(cpp14)]]
using namespace arma;
#include <r2sundials.h>
// [[Rcpp::export]]
SEXP getXPtr();

int d_exp(double t, const vec& y, vec& ydot, RObject& param, NumericVector& psens) {
  ydot[0] = -psens["nu"]*(y[0]-psens["lim"]);
  return(CV_SUCCESS);
}

SEXP getXPtr() {
  typedef int (*funcPtr)( double t, const vec& y, vec& ydot, RObject& param, NumericVector& psens );
  return XPtr<funcPtr>(new funcPtr(& d_exp ));
}

Generated extern "C" functions 
--------------------------------------------------------


#include <Rcpp.h>
#ifdef RCPP_USE_GLOBAL_ROSTREAM
Rcpp::Rostream<true>&  Rcpp::Rcout = Rcpp::Rcpp_cout_get();
Rcpp::Rostream<false>& Rcpp::Rcerr = Rcpp::Rcpp_cerr_get();
#endif

// getXPtr
SEXP getXPtr();
RcppExport SEXP sourceCpp_15_getXPtr() {
BEGIN_RCPP
    Rcpp::RObject rcpp_result_gen;
    Rcpp::RNGScope rcpp_rngScope_gen;
    rcpp_result_gen = Rcpp::wrap(getXPtr());
    return rcpp_result_gen;
END_RCPP
}

Generated R functions 
-------------------------------------------------------

`.sourceCpp_15_DLLInfo` <- dyn.load('D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4322a2fe5/sourceCpp_16.dll')

getXPtr <- Rcpp:::sourceCppFunction(function() {}, FALSE, `.sourceCpp_15_DLLInfo`, 'sourceCpp_15_getXPtr')

rm(`.sourceCpp_15_DLLInfo`)

Building shared library
--------------------------------------------------------

DIR: D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4322a2fe5

D:/RCompile/recent/R/bin/x64/R CMD SHLIB -o "sourceCpp_16.dll" "file2f9b4752639c6.cpp" 
using C++ compiler: 'g++.exe (GCC) 13.2.0'
using C++14
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4322a2fe5'
g++  -std=gnu++14 -I"D:/RCompile/recent/R/include" -DNDEBUG -I../inst/include -fopenmp  -I"D:/RCompile/CRANpkg/lib/4.5/Rcpp/include" -I"D:/RCompile/CRANpkg/lib/4.5/RcppArmadillo/include" -I"D:/RCompile/CRANincoming/R-devel/lib/r2sundials/include" -I"D:/RCompile/CRANpkg/lib/4.5/rmumps/include" -I"D:/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13"   -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include"     -O2 -Wall  -mfpmath=sse -msse2 -mstackrealign  -c file2f9b4752639c6.cpp -o file2f9b4752639c6.o
g++ -shared -s -static-libgcc -o sourceCpp_16.dll tmp.def file2f9b4752639c6.o -fopenmp -LD:/RCompile/recent/R/bin/x64 -lRlapack -LD:/RCompile/recent/R/bin/x64 -lRblas -lgfortran -lm -lquadmath -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib/x64 -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib -LD:/RCompile/recent/R/bin/x64 -lR
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4322a2fe5'
make[1]: Entering directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4322a2fe5'
make[1]: Leaving directory '/d/RCompile/CRANincoming/R-devel/r2sundials.Rcheck/tests/lib/sourceCpp-x86_64-w64-mingw32-1.0.13/sourcecpp_2f9b4322a2fe5'


Executing test function test.dense  ...  done successfully.



Executing test function test.expdecay  ...  done successfully.



Executing test function test.r_vs_cpp  ...  done successfully.



Executing test function test.root.cpp  ...  done successfully.



Executing test function test.root.r  ...  done successfully.



Executing test function test.sensitivity  ...  done successfully.



Executing test function test.sparse  ...  done successfully.

RUNIT TEST PROTOCOL -- Mon Jul 22 17:12:20 2024 
*********************************************** 
Number of test functions: 7 
Number of errors: 0 
Number of failures: 0 

 
1 Test Suite : 
r2sundials unit tests - 7 test functions, 0 errors, 0 failures



Details 
*************************** 
Test Suite: r2sundials unit tests 
Test function regexp: ^[Tt]est.+ 
Test file regexp: ^runit.+\.[rR]$ 
Involved directory: 
D:/RCompile/CRANincoming/R-devel/lib/r2sundials/unitTests 
--------------------------- 
Test file: D:/RCompile/CRANincoming/R-devel/lib/r2sundials/unitTests/runit.r2cvodes.R 
test.dense: (1 checks) ... OK (0 seconds)
test.expdecay: (1 checks) ... OK (0 seconds)
test.r_vs_cpp: (1 checks) ... OK (0 seconds)
test.root.cpp: (1 checks) ... OK (0 seconds)
test.root.r: (2 checks) ... OK (0 seconds)
test.sensitivity: (1 checks) ... OK (0 seconds)
test.sparse: (1 checks) ... OK (0 seconds)
Warning messages:
1: In RNGkind(kind = testSuite$rngKind, normal.kind = testSuite$rngNormalKind) :
  RNGkind: Marsaglia-Multicarry has poor statistical properties
2: In RNGkind(kind = testSuite$rngKind, normal.kind = testSuite$rngNormalKind) :
  RNGkind: severe deviations from normality for Kinderman-Ramage + Marsaglia-Multicarry
> 
> proc.time()
   user  system elapsed 
   0.85    0.25   55.03