corr_dev: clarification

Author: jvdp1

src/stdlib_experimental_stats_corr.fypp

@@ -200,8 +200,8 @@ contains
 
       select case(dim)
         case(1)
-          do i = 1, size(x, 2)
-            do j = 1, size(x, 2)
+          do i = 1, size(res, 2)
+            do j = 1, size(res, 1)
              mask_ = merge(.true., .false., mask(:, i) .and. mask(:, j))
              centeri_ = merge( x(:, i) - mean(x(:, i), mask = mask_),&
               #:if t1[0] == 'r'
@@ -225,8 +225,8 @@ contains
             end do
           end do
         case(2)
-          do i = 1, size(x, 1)
-            do j = 1, size(x, 1)
+          do i = 1, size(res, 2)
+            do j = 1, size(res, 1)
              mask_ = merge(.true., .false., mask(i, :) .and. mask(j, :))
              centeri_ = merge( x(i, :) - mean(x(i, :), mask = mask_),&
               #:if t1[0] == 'r'
@@ -272,15 +272,13 @@ contains
 
       select case(dim)
         case(1)
-          do i = 1, size(x, 2)
-            do j = 1, size(x, 2)
+          do i = 1, size(res, 2)
+            do j = 1, size(res, 1)
              mask_ = merge(.true., .false., mask(:, i) .and. mask(:, j))
              centeri_ = merge( x(:, i) - mean(x(:, i), mask = mask_),&
-                0._dp,&
-                mask_)
+                0._dp, mask_)
              centerj_ = merge( x(:, j) - mean(x(:, j), mask = mask_),&
-                0._dp,&
-                mask_)
+                0._dp, mask_)
 
              res(j, i) = dot_product( centerj_, centeri_)&
                  /sqrt(dot_product( centeri_, centeri_)*&
@@ -289,15 +287,13 @@ contains
             end do
           end do
         case(2)
-          do i = 1, size(x, 1)
-            do j = 1, size(x, 1)
+          do i = 1, size(res, 2)
+            do j = 1, size(res, 1)
              mask_ = merge(.true., .false., mask(i, :) .and. mask(j, :))
              centeri_ = merge( x(i, :) - mean(x(i, :), mask = mask_),&
-                0._dp,&
-                mask_)
+                0._dp, mask_)
              centerj_ = merge( x(j, :) - mean(x(j, :), mask = mask_),&
-                0._dp,&
-                mask_)
+                0._dp, mask_)
 
              res(j, i) = dot_product( centerj_, centeri_)&
                  /sqrt(dot_product( centeri_, centeri_)*&