parameters over values + misc change

Author: Calluumm

test/stats/test_distribution_beta.fypp

@@ -105,11 +105,12 @@ contains
     #:for k1, t1 in RC_KINDS_TYPES
     subroutine test_beta_pdf_${t1[0]}$${k1}$
         ${t1}$ :: x1, x2(3,4), ba, bb, loc
+                integer, parameter :: n = 15
         integer :: i
-        real(${k1}$) :: res(15)
+                real(${k1}$) :: res(n)
     #:if t1[0] == "r"
       #! for real type
-        real(${k1}$), parameter :: ans(15) =                                   &
+                real(${k1}$), parameter :: ans(*) =                                    &
                  [2.45759999999999978E+00_${k1}$,            &
                   2.45759999999999978E+00_${k1}$,            &
                   2.45759999999999978E+00_${k1}$,            &
@@ -127,7 +128,7 @@ contains
                   3.83999999999999897E-02_${k1}$]
     #:else
       #! for complex type
-        real(${k1}$), parameter :: ans(15) =                                   &
+                real(${k1}$), parameter :: ans(*) =                                    &
                  [2.77620563793055197E+00_${k1}$,            &
                   2.77620563793055197E+00_${k1}$,            &
                   2.77620563793055197E+00_${k1}$,            &
@@ -147,13 +148,17 @@ contains
 
         print *, "Test beta_distribution_pdf_${t1[0]}$${k1}$"
     #:if t1[0] == "r"
-        ba = 2.0_${k1}$; bb = 5.0_${k1}$; loc = 0._${k1}$
+        ba = 2.0_${k1}$
+        bb = 5.0_${k1}$
+        loc = 0._${k1}$
         x1 = 0.2_${k1}$
         x2 = reshape([0.05_${k1}$, 0.1_${k1}$, 0.15_${k1}$, 0.25_${k1}$,      &
                       0.3_${k1}$, 0.35_${k1}$, 0.4_${k1}$, 0.45_${k1}$,       &
                       0.5_${k1}$, 0.6_${k1}$, 0.7_${k1}$, 0.8_${k1}$], [3,4])
     #:else
-        ba = (2.0_${k1}$, 0.7_${k1}$); bb = (5.0_${k1}$, 3.0_${k1}$); loc = (0._${k1}$, 0._${k1}$)
+        ba = (2.0_${k1}$, 0.7_${k1}$)
+        bb = (5.0_${k1}$, 3.0_${k1}$)
+        loc = (0._${k1}$, 0._${k1}$)
         x1 = (0.2_${k1}$, 0.3_${k1}$)
         x2 = reshape([(0.05_${k1}$, 0.05_${k1}$), (0.1_${k1}$, 0.1_${k1}$),   &
                       (0.15_${k1}$, 0.2_${k1}$), (0.25_${k1}$, 0.25_${k1}$),   &
@@ -166,7 +171,7 @@ contains
         res(1:3) = beta_pdf(x1, ba, bb, loc)
         res(4:15) = reshape(beta_pdf(x2, ba, bb, loc), [12])
 
-        do i = 1, 15
+        do i = 1, n
             call check(abs(res(i) - ans(i)) < max(${k1}$tol, 1.0e-11_${k1}$), &
                        msg="beta_distribution_pdf_${t1[0]}$${k1}$ failed", &
                        warn=warn)