changed

Author: srinjoy933

src/specialmatrices/stdlib_specialmatrices_tridiagonal.fypp

@@ -15,7 +15,7 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
     !-----                      -----
     !--------------------------------
 
-    #:for k1, t1, s1 in (KINDS_TYPES)
+    #:for k1, t1, s1 in KINDS_TYPES
     pure module function initialize_tridiagonal_pure_${s1}$(dl, dv, du) result(A)
         !! Construct a `tridiagonal` matrix from the rank-1 arrays
         !! `dl`, `dv` and `du`.
@@ -147,7 +147,7 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
     !-----------------------------------------
 
     !! spmv_tridiag
-    #:for k1, t1, s1 in (KINDS_TYPES)
+    #:for k1, t1, s1 in KINDS_TYPES
     #:for rank in RANKS
     module subroutine spmv_tridiag_${rank}$d_${s1}$(A, x, y, alpha, beta, op)
         type(tridiagonal_${s1}$_type), intent(in) :: A
@@ -226,7 +226,7 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
     !-----                           -----
     !-------------------------------------
 
-    #:for k1, t1, s1 in (KINDS_TYPES)
+    #:for k1, t1, s1 in KINDS_TYPES
     pure module function tridiagonal_to_dense_${s1}$(A) result(B)
         !! Convert a `tridiagonal` matrix to its dense representation.
         type(tridiagonal_${s1}$_type), intent(in) :: A
@@ -247,19 +247,21 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
         if (n == 1) then
             B(1, 1) = A%dv(1)
         else
-            B(1, 1) = A%dv(1) ; B(1, 2) = A%du(1)
+            B(1, 1) = A%dv(1)
+            B(1, 2) = A%du(1)
             do concurrent (i=2:n-1)
                 B(i, i-1) = A%dl(i-1)
                 B(i, i) = A%dv(i)
                 B(i, i+1) = A%du(i)
             enddo
-            B(n, n-1) = A%dl(n-1) ; B(n, n) = A%dv(n)
+            B(n, n-1) = A%dl(n-1)
+            B(n, n) = A%dv(n)
         end if
         end associate
     end function
     #:endfor
 
-    #:for k1, t1, s1 in (KINDS_TYPES)
+    #:for k1, t1, s1 in KINDS_TYPES
     pure module function transpose_tridiagonal_${s1}$(A) result(B)
         type(tridiagonal_${s1}$_type), intent(in) :: A
         !! Input matrix.
@@ -268,7 +270,7 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
     end function
     #:endfor
 
-    #:for k1, t1, s1 in (KINDS_TYPES)
+    #:for k1, t1, s1 in KINDS_TYPES
     pure module function hermitian_tridiagonal_${s1}$(A) result(B)
         type(tridiagonal_${s1}$_type), intent(in) :: A
         !! Input matrix.
@@ -281,7 +283,7 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
     end function
     #:endfor
 
-    #:for k1, t1, s1 in (KINDS_TYPES)
+    #:for k1, t1, s1 in KINDS_TYPES
     pure module function scalar_multiplication_tridiagonal_${s1}$(alpha, A) result(B)
         ${t1}$, intent(in) :: alpha
         type(tridiagonal_${s1}$_type), intent(in) :: A
@@ -299,7 +301,7 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
     end function
     #:endfor
 
-    #:for k1, t1, s1 in (KINDS_TYPES)
+    #:for k1, t1, s1 in KINDS_TYPES
     pure module function matrix_add_tridiagonal_${s1}$(A, B) result(C)
         type(tridiagonal_${s1}$_type), intent(in) :: A
         type(tridiagonal_${s1}$_type), intent(in) :: B

test/linalg/test_linalg_specialmatrices.fypp

@@ -155,7 +155,7 @@ contains
     subroutine test_tridiagonal_1x1(error)
         !> Test 1x1 matrix edge case for dense conversion
         type(error_type), allocatable, intent(out) :: error
-        #:for k1, t1, s1 in (KINDS_TYPES)
+        #:for k1, t1, s1 in KINDS_TYPES
         block
             integer, parameter :: wp = ${k1}$
             type(tridiagonal_${s1}$_type) :: A
@@ -174,7 +174,7 @@ contains
             A = tridiagonal(dl, dv, du) 
             Amat = dense(A)
 
-            ! Check if it compiled and converted properly without segfaulting
+           ! Check if the 1x1 matrix converted properly at runtime without segfaulting
             call check(error, size(Amat, 1) == 1, .true.)
             if (allocated(error)) return
             call check(error, size(Amat, 2) == 1, .true.)
@@ -188,7 +188,7 @@ contains
     subroutine test_tridiagonal_arithmetic(error)
         !> Test arithmetic operations and optimization
         type(error_type), allocatable, intent(out) :: error
-        #:for k1, t1, s1 in (KINDS_TYPES)
+        #:for k1, t1, s1 in KINDS_TYPES
         block
             integer, parameter :: wp = ${k1}$
             integer, parameter :: n = 3