use std::isfinite instead of arma::is_finite

Author: dcgerard

DESCRIPTION

@@ -26,7 +26,7 @@ BugReports: https://github.com/dcgerard/updog/issues
 URL: https://github.com/dcgerard/updog/
 Encoding: UTF-8
 LazyData: true
-RoxygenNote: 7.2.3
+RoxygenNote: 7.3.3
 LinkingTo: Rcpp,
     RcppArmadillo
 Imports: 

man/updog-package.Rd

@@ -125,7 +125,7 @@ NumericVector grad_for_mu_sigma2_wrapper(arma::Col<double> muSigma2, NumericMatr
 // [[Rcpp::export]]
 double dpen_dh(double h, double mu_h, double sigma2_h) {
   double deriv;
-  if (arma::is_finite(sigma2_h)) {
+  if (std::isfinite(sigma2_h)) {
     deriv = -1.0 * (1.0 + (std::log(h) - mu_h) / sigma2_h) / h;
   } else {
     deriv = 0.0;
@@ -153,7 +153,7 @@ double dpen_dh(double h, double mu_h, double sigma2_h) {
 double dpen_deps(double eps, double mu_eps, double sigma2_eps) {
   double deriv;
 
-  if (arma::is_finite(sigma2_eps)) {
+  if (std::isfinite(sigma2_eps)) {
     deriv = -1.0 * (1.0 - 2.0 * eps + (logit(eps) - mu_eps) / sigma2_eps) / (eps * (1.0 - eps));
   } else {
     deriv = 0.0;

src/gradients.cpp

@@ -221,7 +221,7 @@ double pen_bias(double h, double mu_h, double sigma2_h) {
   }
 
   double pen;
-  if (arma::is_finite(sigma2_h)) {
+  if (std::isfinite(sigma2_h)) {
     pen = -std::log(h) - std::pow(std::log(h) - mu_h, 2) / (2.0 * sigma2_h);
   } else {
     pen = 0.0;
@@ -256,7 +256,7 @@ double pen_seq_error(double eps, double mu_eps, double sigma2_eps) {
   }
 
   double pen;
-  if (arma::is_finite(sigma2_eps)) {
+  if (std::isfinite(sigma2_eps)) {
     pen = -log(eps * (1.0 - eps)) - std::pow(logit(eps) - mu_eps, 2) / (2.0 * sigma2_eps);
   } else {
     pen = 0.0;