Fix pow

Author: shy86

DESCRIPTION

@@ -9,10 +9,12 @@ Description: The traditional linear regression trend, Modified Mann-Kendall (MK)
     Hamed, K. H., & Rao, A. R. (1998). A modified Mann-Kendall trend test for 
     autocorrelated data. Journal of hydrology, 204(1-4), 182-196. 
     <doi:10.1016/S0022-1694(97)00125-X>.
-Version: 0.1.1
+Version: 0.1.2
 Authors@R: c(
-    person("Dongdong", "Kong", role = c("aut", "cre"), 
-        email = "kongdd.sysu@gmail.com", comment = c(ORCID = "0000-0003-1836-8172")))
+  person("Dongdong", "Kong", role = c("aut", "cre"), 
+         email = "kongdd.sysu@gmail.com", comment = c(ORCID = "0000-0003-1836-8172")), 
+  person('Heyang', 'Song', role = "aut",
+         email = "7800556@gmail.com", comment = c(ORCID = "0000-0002-4192-5603")))
 License: MIT + file LICENSE
 Encoding: UTF-8
 Roxygen: list(markdown = TRUE)

src/movmean.cpp

@@ -26,7 +26,7 @@ arma::mat sgmat_S(int halfwin = 5, int d = 2) {
 
     for (int i = 0; i < frame; i++) {
         for (int j =0; j <= d; j++) {
-            mat(i, j) = pow(i - halfwin, j);
+            mat(i, j) = pow(i * 1.0 - halfwin, j * 1.0);
         }
     }
     // Rcout << mat << std::endl;
@@ -78,7 +78,7 @@ arma::mat sgmat_wB(const arma::mat S, const arma::colvec w) {
 //' @examples
 //' y <- c(1, 3, 2, 5, 6, 8, 10, 1)
 //' w <- seq_along(y)/length(y)
-//' 
+//'
 //' halfwin = 2
 //' d = 2
 //' s1 <- smooth_wSG(y, halfwin, d, w)
@@ -180,7 +180,7 @@ NumericVector movmean(
     const arma::colvec y,
     int halfwin = 1,
     bool SG_style = false,
-    Nullable<NumericVector> w = R_NilValue) 
+    Nullable<NumericVector> w = R_NilValue)
 {
     int n = y.size();
     arma::colvec yy(y);