new2

Author: srinjoy933

test/linalg/test_linalg_specialmatrices.fypp

@@ -154,10 +154,16 @@ contains
             integer, parameter :: wp = ${k1}$
             type(tridiagonal_${s1}$_type) :: A
             ${t1}$, allocatable :: Amat(:,:)
-            ! Maintainer suggested using parameter for constant values
-            ${t1}$, parameter :: dl(0) = [${t1}$ ::]
-            ${t1}$, parameter :: du(0) = [${t1}$ ::]
+            
+            #:if t1.startswith('complex')
+            ${t1}$, parameter :: dl(0) = [(cmplx(0.0_wp, 0.0_wp, kind=wp), i=1,0)]
+            ${t1}$, parameter :: du(0) = [(cmplx(0.0_wp, 0.0_wp, kind=wp), i=1,0)]
+            ${t1}$, parameter :: dv(1) = [cmplx(5.0_wp, 0.0_wp, kind=wp)]
+            #:else
+            ${t1}$, parameter :: dl(0) = [(0.0_wp, i=1,0)]
+            ${t1}$, parameter :: du(0) = [(0.0_wp, i=1,0)]
             ${t1}$, parameter :: dv(1) = [5.0_wp]
+            #:endif
 
             A = tridiagonal(dl, dv, du) 
             Amat = dense(A)
@@ -182,14 +188,24 @@ contains
             integer, parameter :: n = 3
             type(tridiagonal_${s1}$_type) :: A, B, C
 
-            ! Define components as parameters as suggested by maintainer
-            ${t1}$, parameter :: dl1(n-1) = [1.0_wp, 1.0_wp]
-            ${t1}$, parameter :: dv1(n)   = [2.0_wp, 2.0_wp, 2.0_wp]
-            ${t1}$, parameter :: du1(n-1) = [3.0_wp, 3.0_wp]
+            ! Explicitly separate real vs complex to satisfy strict compile-time type checks
+            #:if t1.startswith('complex')
+            ${t1}$, parameter :: dl1(n-1) = [(cmplx(1.0_wp, 0.0_wp, kind=wp), i=1,n-1)]
+            ${t1}$, parameter :: dv1(n)   = [(cmplx(2.0_wp, 0.0_wp, kind=wp), i=1,n)]
+            ${t1}$, parameter :: du1(n-1) = [(cmplx(3.0_wp, 0.0_wp, kind=wp), i=1,n-1)]
+            
+            ${t1}$, parameter :: dl2(n-1) = [(cmplx(4.0_wp, 0.0_wp, kind=wp), i=1,n-1)]
+            ${t1}$, parameter :: dv2(n)   = [(cmplx(5.0_wp, 0.0_wp, kind=wp), i=1,n)]
+            ${t1}$, parameter :: du2(n-1) = [(cmplx(6.0_wp, 0.0_wp, kind=wp), i=1,n-1)]
+            #:else
+            ${t1}$, parameter :: dl1(n-1) = [(1.0_wp, i=1,n-1)]
+            ${t1}$, parameter :: dv1(n)   = [(2.0_wp, i=1,n)]
+            ${t1}$, parameter :: du1(n-1) = [(3.0_wp, i=1,n-1)]
 
-            ${t1}$, parameter :: dl2(n-1) = [4.0_wp, 4.0_wp]
-            ${t1}$, parameter :: dv2(n)   = [5.0_wp, 5.0_wp, 5.0_wp]
-            ${t1}$, parameter :: du2(n-1) = [6.0_wp, 6.0_wp]
+            ${t1}$, parameter :: dl2(n-1) = [(4.0_wp, i=1,n-1)]
+            ${t1}$, parameter :: dv2(n)   = [(5.0_wp, i=1,n)]
+            ${t1}$, parameter :: du2(n-1) = [(6.0_wp, i=1,n-1)]
+            #:endif
 
             A = tridiagonal(dl1, dv1, du1)
             B = tridiagonal(dl2, dv2, du2)