RcppExports.cpp 6.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153
  1. // Generated by using Rcpp::compileAttributes() -> do not edit by hand
  2. // Generator token: 10BE3573-1514-4C36-9D1C-5A225CD40393
  3. #include <RcppArmadillo.h>
  4. #include <Rcpp.h>
  5. using namespace Rcpp;
  6. #ifdef RCPP_USE_GLOBAL_ROSTREAM
  7. Rcpp::Rostream<true>& Rcpp::Rcout = Rcpp::Rcpp_cout_get();
  8. Rcpp::Rostream<false>& Rcpp::Rcerr = Rcpp::Rcpp_cerr_get();
  9. #endif
  10. // matrix_multiplication_cpp
  11. arma::mat matrix_multiplication_cpp(arma::mat A, arma::mat B);
  12. RcppExport SEXP _CALDER_matrix_multiplication_cpp(SEXP ASEXP, SEXP BSEXP) {
  13. BEGIN_RCPP
  14. Rcpp::RObject rcpp_result_gen;
  15. Rcpp::RNGScope rcpp_rngScope_gen;
  16. Rcpp::traits::input_parameter< arma::mat >::type A(ASEXP);
  17. Rcpp::traits::input_parameter< arma::mat >::type B(BSEXP);
  18. rcpp_result_gen = Rcpp::wrap(matrix_multiplication_cpp(A, B));
  19. return rcpp_result_gen;
  20. END_RCPP
  21. }
  22. // matrix_multiplication_sym_cpp
  23. arma::mat matrix_multiplication_sym_cpp(arma::mat A);
  24. RcppExport SEXP _CALDER_matrix_multiplication_sym_cpp(SEXP ASEXP) {
  25. BEGIN_RCPP
  26. Rcpp::RObject rcpp_result_gen;
  27. Rcpp::RNGScope rcpp_rngScope_gen;
  28. Rcpp::traits::input_parameter< arma::mat >::type A(ASEXP);
  29. rcpp_result_gen = Rcpp::wrap(matrix_multiplication_sym_cpp(A));
  30. return rcpp_result_gen;
  31. END_RCPP
  32. }
  33. // loglik_lnorm_cpp
  34. double loglik_lnorm_cpp(double sum_ln1, double sum_ln2, double p, double q);
  35. RcppExport SEXP _CALDER_loglik_lnorm_cpp(SEXP sum_ln1SEXP, SEXP sum_ln2SEXP, SEXP pSEXP, SEXP qSEXP) {
  36. BEGIN_RCPP
  37. Rcpp::RObject rcpp_result_gen;
  38. Rcpp::RNGScope rcpp_rngScope_gen;
  39. Rcpp::traits::input_parameter< double >::type sum_ln1(sum_ln1SEXP);
  40. Rcpp::traits::input_parameter< double >::type sum_ln2(sum_ln2SEXP);
  41. Rcpp::traits::input_parameter< double >::type p(pSEXP);
  42. Rcpp::traits::input_parameter< double >::type q(qSEXP);
  43. rcpp_result_gen = Rcpp::wrap(loglik_lnorm_cpp(sum_ln1, sum_ln2, p, q));
  44. return rcpp_result_gen;
  45. END_RCPP
  46. }
  47. // loglik_lnorm_cpp_vec
  48. double loglik_lnorm_cpp_vec(arma::vec vec_values);
  49. RcppExport SEXP _CALDER_loglik_lnorm_cpp_vec(SEXP vec_valuesSEXP) {
  50. BEGIN_RCPP
  51. Rcpp::RObject rcpp_result_gen;
  52. Rcpp::RNGScope rcpp_rngScope_gen;
  53. Rcpp::traits::input_parameter< arma::vec >::type vec_values(vec_valuesSEXP);
  54. rcpp_result_gen = Rcpp::wrap(loglik_lnorm_cpp_vec(vec_values));
  55. return rcpp_result_gen;
  56. END_RCPP
  57. }
  58. // get_A_len
  59. arma::mat get_A_len(arma::mat A);
  60. RcppExport SEXP _CALDER_get_A_len(SEXP ASEXP) {
  61. BEGIN_RCPP
  62. Rcpp::RObject rcpp_result_gen;
  63. Rcpp::RNGScope rcpp_rngScope_gen;
  64. Rcpp::traits::input_parameter< arma::mat >::type A(ASEXP);
  65. rcpp_result_gen = Rcpp::wrap(get_A_len(A));
  66. return rcpp_result_gen;
  67. END_RCPP
  68. }
  69. // get_A_ln1
  70. arma::mat get_A_ln1(arma::mat A);
  71. RcppExport SEXP _CALDER_get_A_ln1(SEXP ASEXP) {
  72. BEGIN_RCPP
  73. Rcpp::RObject rcpp_result_gen;
  74. Rcpp::RNGScope rcpp_rngScope_gen;
  75. Rcpp::traits::input_parameter< arma::mat >::type A(ASEXP);
  76. rcpp_result_gen = Rcpp::wrap(get_A_ln1(A));
  77. return rcpp_result_gen;
  78. END_RCPP
  79. }
  80. // get_A_ln2
  81. arma::mat get_A_ln2(arma::mat A);
  82. RcppExport SEXP _CALDER_get_A_ln2(SEXP ASEXP) {
  83. BEGIN_RCPP
  84. Rcpp::RObject rcpp_result_gen;
  85. Rcpp::RNGScope rcpp_rngScope_gen;
  86. Rcpp::traits::input_parameter< arma::mat >::type A(ASEXP);
  87. rcpp_result_gen = Rcpp::wrap(get_A_ln2(A));
  88. return rcpp_result_gen;
  89. END_RCPP
  90. }
  91. // loglik_lnorm_cpp_mat
  92. arma::mat loglik_lnorm_cpp_mat(arma::mat sum_ln1, arma::mat sum_ln2, arma::mat ps, arma::mat qs);
  93. RcppExport SEXP _CALDER_loglik_lnorm_cpp_mat(SEXP sum_ln1SEXP, SEXP sum_ln2SEXP, SEXP psSEXP, SEXP qsSEXP) {
  94. BEGIN_RCPP
  95. Rcpp::RObject rcpp_result_gen;
  96. Rcpp::RNGScope rcpp_rngScope_gen;
  97. Rcpp::traits::input_parameter< arma::mat >::type sum_ln1(sum_ln1SEXP);
  98. Rcpp::traits::input_parameter< arma::mat >::type sum_ln2(sum_ln2SEXP);
  99. Rcpp::traits::input_parameter< arma::mat >::type ps(psSEXP);
  100. Rcpp::traits::input_parameter< arma::mat >::type qs(qsSEXP);
  101. rcpp_result_gen = Rcpp::wrap(loglik_lnorm_cpp_mat(sum_ln1, sum_ln2, ps, qs));
  102. return rcpp_result_gen;
  103. END_RCPP
  104. }
  105. // zigzag_loglik_ancestors_v4_5
  106. List zigzag_loglik_ancestors_v4_5(arma::mat A, int k, int min_n_bins);
  107. RcppExport SEXP _CALDER_zigzag_loglik_ancestors_v4_5(SEXP ASEXP, SEXP kSEXP, SEXP min_n_binsSEXP) {
  108. BEGIN_RCPP
  109. Rcpp::RObject rcpp_result_gen;
  110. Rcpp::RNGScope rcpp_rngScope_gen;
  111. Rcpp::traits::input_parameter< arma::mat >::type A(ASEXP);
  112. Rcpp::traits::input_parameter< int >::type k(kSEXP);
  113. Rcpp::traits::input_parameter< int >::type min_n_bins(min_n_binsSEXP);
  114. rcpp_result_gen = Rcpp::wrap(zigzag_loglik_ancestors_v4_5(A, k, min_n_bins));
  115. return rcpp_result_gen;
  116. END_RCPP
  117. }
  118. // compute_L
  119. List compute_L(arma::mat A, arma::mat L, int k);
  120. RcppExport SEXP _CALDER_compute_L(SEXP ASEXP, SEXP LSEXP, SEXP kSEXP) {
  121. BEGIN_RCPP
  122. Rcpp::RObject rcpp_result_gen;
  123. Rcpp::RNGScope rcpp_rngScope_gen;
  124. Rcpp::traits::input_parameter< arma::mat >::type A(ASEXP);
  125. Rcpp::traits::input_parameter< arma::mat >::type L(LSEXP);
  126. Rcpp::traits::input_parameter< int >::type k(kSEXP);
  127. rcpp_result_gen = Rcpp::wrap(compute_L(A, L, k));
  128. return rcpp_result_gen;
  129. END_RCPP
  130. }
  131. static const R_CallMethodDef CallEntries[] = {
  132. {"_CALDER_matrix_multiplication_cpp", (DL_FUNC) &_CALDER_matrix_multiplication_cpp, 2},
  133. {"_CALDER_matrix_multiplication_sym_cpp", (DL_FUNC) &_CALDER_matrix_multiplication_sym_cpp, 1},
  134. {"_CALDER_loglik_lnorm_cpp", (DL_FUNC) &_CALDER_loglik_lnorm_cpp, 4},
  135. {"_CALDER_loglik_lnorm_cpp_vec", (DL_FUNC) &_CALDER_loglik_lnorm_cpp_vec, 1},
  136. {"_CALDER_get_A_len", (DL_FUNC) &_CALDER_get_A_len, 1},
  137. {"_CALDER_get_A_ln1", (DL_FUNC) &_CALDER_get_A_ln1, 1},
  138. {"_CALDER_get_A_ln2", (DL_FUNC) &_CALDER_get_A_ln2, 1},
  139. {"_CALDER_loglik_lnorm_cpp_mat", (DL_FUNC) &_CALDER_loglik_lnorm_cpp_mat, 4},
  140. {"_CALDER_zigzag_loglik_ancestors_v4_5", (DL_FUNC) &_CALDER_zigzag_loglik_ancestors_v4_5, 3},
  141. {"_CALDER_compute_L", (DL_FUNC) &_CALDER_compute_L, 3},
  142. {NULL, NULL, 0}
  143. };
  144. RcppExport void R_init_CALDER(DllInfo *dll) {
  145. R_registerRoutines(dll, NULL, CallEntries, NULL, NULL);
  146. R_useDynamicSymbols(dll, FALSE);
  147. }