raw_moment_dev: correction of test_rawmoment

Author: jvdp1

src/tests/stats/test_moment.f90

@@ -627,7 +627,6 @@ subroutine test_csp(x1, x2)
         complex(sp), intent(in) :: x1(:), x2(:, :)
 
         integer :: order
-        complex(sp), allocatable :: x3(:, :, :)
 
         order = 1
 

src/tests/stats/test_rawmoment.f90

@@ -45,13 +45,6 @@ subroutine test_sp(x1, x2)
         zero2_1 = 0
         zero2_2 = 0
 
-        allocate(zero3_1(size(x3, 2), size(x3, 3))&
-                 ,zero3_2(size(x3, 1), size(x3, 3))&
-                 ,zero3_3(size(x3, 1), size(x3, 2)))
-        zero3_1 = 0
-        zero3_2 = 0
-        zero3_3 = 0
-
         order = 1
 
         !1dim
@@ -110,6 +103,13 @@ subroutine test_sp(x1, x2)
         x3(:,:,2)=x2*2;
         x3(:,:,3)=x2*4;
 
+        allocate(zero3_1(size(x3, 2), size(x3, 3))&
+                 ,zero3_2(size(x3, 1), size(x3, 3))&
+                 ,zero3_3(size(x3, 1), size(x3, 2)))
+        zero3_1 = 0
+        zero3_2 = 0
+        zero3_3 = 0
+
         print*,' test_sp_3dim', order
         call assert( abs(moment(x3, order, center = 0.) - mean(x3)) < sptol)
         call assert( all( abs( moment(x3, order, dim = 1, center = 0._sp) -&
@@ -270,14 +270,6 @@ subroutine test_int32(x1, x2)
         zero2_1 = 0
         zero2_2 = 0
 
-        allocate(zero3_1(size(x3, 2), size(x3, 3))&
-                 ,zero3_2(size(x3, 1), size(x3, 3))&
-                 ,zero3_3(size(x3, 1), size(x3, 2)))
-        zero3_1 = 0
-        zero3_2 = 0
-        zero3_3 = 0
-
-
         order = 1
 
         !1dim
@@ -337,6 +329,13 @@ subroutine test_int32(x1, x2)
         x3(:,:,2)=x2*2;
         x3(:,:,3)=x2*4;
 
+        allocate(zero3_1(size(x3, 2), size(x3, 3))&
+                 ,zero3_2(size(x3, 1), size(x3, 3))&
+                 ,zero3_3(size(x3, 1), size(x3, 2)))
+        zero3_1 = 0
+        zero3_2 = 0
+        zero3_3 = 0
+
         print*,' test_sp_3dim', order
         call assert( abs(moment(x3, order, center = 0._dp) - mean(x3)) < sptol)
         call assert( all( abs( moment(x3, order, dim = 1, center = 0._dp) -&
@@ -500,7 +499,6 @@ subroutine test_csp(x1, x2)
         complex(sp), intent(in) :: x1(:), x2(:, :)
 
         integer :: order
-        complex(sp), allocatable :: x3(:, :, :)
         complex(sp), allocatable :: zero2_1(:), zero2_2(:)
 
         allocate(zero2_1(size(x2, 2)), zero2_2(size(x2, 1)))