rawmoment: start to add center in test

Author: jvdp1

src/stdlib_experimental_stats.fypp

@@ -213,9 +213,10 @@ module stdlib_experimental_stats
     #:for k1, t1 in RC_KINDS_TYPES
       #:for rank in RANKS
         #:set RName = rname("moment_all",rank, t1, k1)
-        module function ${RName}$(x, order, mask) result(res)
+        module function ${RName}$(x, order, center, mask) result(res)
           ${t1}$, intent(in) :: x${ranksuffix(rank)}$
           integer, intent(in) :: order
+          ${t1}$, intent(in), optional :: center
           logical, intent(in), optional :: mask
           ${t1}$ :: res
         end function ${RName}$
@@ -225,9 +226,10 @@ module stdlib_experimental_stats
     #:for k1, t1 in INT_KINDS_TYPES
       #:for rank in RANKS
         #:set RName = rname("moment_all",rank, t1, k1, 'dp')
-        module function ${RName}$(x, order, mask) result(res)
+        module function ${RName}$(x, order, center, mask) result(res)
           ${t1}$, intent(in) :: x${ranksuffix(rank)}$
           integer, intent(in) :: order
+          real(dp), intent(in), optional :: center
           logical, intent(in), optional :: mask
           real(dp) :: res
         end function ${RName}$
@@ -237,10 +239,11 @@ module stdlib_experimental_stats
     #:for k1, t1 in RC_KINDS_TYPES
       #:for rank in RANKS
         #:set RName = rname("moment",rank, t1, k1)
-        module function ${RName}$(x, order, dim, mask) result(res)
+        module function ${RName}$(x, order, center, dim, mask) result(res)
           ${t1}$, intent(in) :: x${ranksuffix(rank)}$
           integer, intent(in) :: order
           integer, intent(in) :: dim
+          ${t1}$, intent(in), optional :: center${reduced_shape('x', rank, 'dim')}$
           logical, intent(in), optional :: mask
           ${t1}$ :: res${reduced_shape('x', rank, 'dim')}$
         end function ${RName}$
@@ -250,10 +253,11 @@ module stdlib_experimental_stats
     #:for k1, t1 in INT_KINDS_TYPES
       #:for rank in RANKS
         #:set RName = rname("moment",rank, t1, k1, 'dp')
-        module function ${RName}$(x, order, dim, mask) result(res)
+        module function ${RName}$(x, order, center, dim, mask) result(res)
           ${t1}$, intent(in) :: x${ranksuffix(rank)}$
           integer, intent(in) :: order
           integer, intent(in) :: dim
+          real(dp),intent(in), optional :: center${reduced_shape('x', rank, 'dim')}$
           logical, intent(in), optional :: mask
           real(dp) :: res${reduced_shape('x', rank, 'dim')}$
         end function ${RName}$
@@ -263,9 +267,10 @@ module stdlib_experimental_stats
     #:for k1, t1 in RC_KINDS_TYPES
       #:for rank in RANKS
         #:set RName = rname("moment_mask_all",rank, t1, k1)
-        module function ${RName}$(x, order, mask) result(res)
+        module function ${RName}$(x, order, center, mask) result(res)
           ${t1}$, intent(in) :: x${ranksuffix(rank)}$
           integer, intent(in) :: order
+          ${t1}$, intent(in), optional :: center
           logical, intent(in) :: mask${ranksuffix(rank)}$
           ${t1}$ :: res
         end function ${RName}$
@@ -275,9 +280,10 @@ module stdlib_experimental_stats
     #:for k1, t1 in INT_KINDS_TYPES
       #:for rank in RANKS
         #:set RName = rname("moment_mask_all",rank, t1, k1, 'dp')
-        module function ${RName}$(x, order, mask) result(res)
+        module function ${RName}$(x, order, center, mask) result(res)
           ${t1}$, intent(in) :: x${ranksuffix(rank)}$
           integer, intent(in) :: order
+          real(dp),intent(in), optional :: center
           logical, intent(in) :: mask${ranksuffix(rank)}$
           real(dp) :: res
         end function ${RName}$
@@ -287,10 +293,11 @@ module stdlib_experimental_stats
     #:for k1, t1 in RC_KINDS_TYPES
       #:for rank in RANKS
         #:set RName = rname("moment_mask",rank, t1, k1)
-        module function ${RName}$(x, order, dim, mask) result(res)
+        module function ${RName}$(x, order, center, dim, mask) result(res)
           ${t1}$, intent(in) :: x${ranksuffix(rank)}$
           integer, intent(in) :: order
           integer, intent(in) :: dim
+          ${t1}$, intent(in), optional :: center${reduced_shape('x', rank, 'dim')}$
           logical, intent(in) :: mask${ranksuffix(rank)}$
           ${t1}$ :: res${reduced_shape('x', rank, 'dim')}$
         end function ${RName}$
@@ -300,10 +307,11 @@ module stdlib_experimental_stats
     #:for k1, t1 in INT_KINDS_TYPES
       #:for rank in RANKS
         #:set RName = rname("moment_mask",rank, t1, k1, 'dp')
-        module function ${RName}$(x, order, dim, mask) result(res)
+        module function ${RName}$(x, order, center, dim, mask) result(res)
           ${t1}$, intent(in) :: x${ranksuffix(rank)}$
           integer, intent(in) :: order
           integer, intent(in) :: dim
+          real(dp), intent(in), optional :: center${reduced_shape('x', rank, 'dim')}$
           logical, intent(in) :: mask${ranksuffix(rank)}$
           real(dp) :: res${reduced_shape('x', rank, 'dim')}$
         end function ${RName}$