explicit cast

Author: jvdp1

src/stdlib_stats_moment.fypp

@@ -28,7 +28,7 @@ contains
           return
         end if
 
-        n = size(x, kind = int64)
+        n = real(size(x, kind = int64), ${k1}$)
 
         if (present(center)) then
          res = sum((x - center)**order) / n
@@ -58,7 +58,7 @@ contains
           return
         end if
 
-        n = size(x, kind = int64)
+        n = real(size(x, kind = int64), dp)
 
         if (present(center)) then
          res = sum((real(x, dp) - center)**order) / n
@@ -118,7 +118,7 @@ contains
           return
         end if
 
-        n = size(x, dim)
+        n = real(size(x, dim), ${k1}$)
 
         res = 0
         select case(dim)
@@ -193,7 +193,7 @@ contains
           return
         end if
 
-        n = size(x, dim)
+        n = real(size(x, dim), dp)
 
         res = 0
         select case(dim)
@@ -234,7 +234,7 @@ contains
 
         real(${k1}$) :: n
 
-        n = count(mask, kind = int64)
+        n = real(count(mask, kind = int64), ${k1}$)
 
         if (present(center)) then
          res = sum((x - center)**order, mask) / n
@@ -259,7 +259,7 @@ contains
 
         real(dp) :: n
 
-        n = count(mask, kind = int64)
+        n = real(count(mask, kind = int64), dp)
 
         if (present(center)) then
          res = sum((real(x, dp) - center)**order, mask) / n
@@ -309,7 +309,7 @@ contains
         real(${k1}$) :: n${reduced_shape('x', rank, 'dim')}$
         ${t1}$, allocatable :: mean_${ranksuffix(rank-1)}$
 
-        n = count(mask, dim)
+        n = real(count(mask, dim), ${k1}$)
 
         res = 0
         select case(dim)
@@ -387,7 +387,7 @@ contains
         real(dp) :: n${reduced_shape('x', rank, 'dim')}$
         real(dp), allocatable :: mean_${ranksuffix(rank-1)}$
 
-        n = count(mask, dim)
+        n = real(count(mask, dim), dp)
 
         res = 0
         select case(dim)

src/stdlib_stats_var.fypp

@@ -27,7 +27,7 @@ contains
           return
         end if
 
-        n = size(x, kind = int64)
+        n = real(size(x, kind = int64), ${k1}$)
         mean = sum(x) / n
 
         #:if t1[0] == 'r'
@@ -57,7 +57,7 @@ contains
           return
         end if
 
-        n = size(x, kind = int64)
+        n = real(size(x, kind = int64), dp)
         mean = sum(real(x, dp)) / n
 
         res = sum((real(x, dp) - mean)**2) / (n - merge(1, 0, optval(corrected, .true.)))
@@ -161,7 +161,7 @@ contains
         real(${k1}$) :: n
         ${t1}$ :: mean
 
-        n = count(mask, kind = int64)
+        n = real(count(mask, kind = int64), ${k1}$)
         mean = sum(x, mask) / n
 
         #:if t1[0] == 'r'
@@ -187,7 +187,7 @@ contains
 
         real(dp) :: n, mean
 
-        n = count(mask, kind = int64)
+        n = real(count(mask, kind = int64), dp)
         mean = sum(real(x, dp), mask) / n
 
         res = sum((real(x, dp) - mean)**2, mask) / (n -&

src/tests/stats/test_corr.f90

@@ -233,7 +233,7 @@ subroutine test_int64(x, x2)
             , 'int64 check 1')
         call check( ieee_is_nan(corr(x, 1, .false.))&
             , 'int64 check 2')
-        call check( ieee_is_nan(corr(x, 1, x == 1.)), 'int64 check 3')
+        call check( ieee_is_nan(corr(x, 1, x == 1)), 'int64 check 3')
         call check( abs(corr(x, 1, x < 5) - 1._dp) < dptol, 'int64 check 4')
         call check( ieee_is_nan(corr(x, 1, x < -999)), 'int64 check 5')