modify: add new lines between inline assignments

Author: Mahmood-Sinan

src/specialmatrices/stdlib_specialmatrices_sym_tridiagonal.fypp

@@ -156,9 +156,12 @@ submodule (stdlib_specialmatrices) sym_tridiagonal_matrices
         end if
 
         ! Deal with optional arguments.
-        alpha_ = 1.0_${k1}$ ; if (present(alpha)) alpha_ = alpha
-        beta_  = 0.0_${k1}$ ; if (present(beta))  beta_  = beta
-        op_    = "N"        ; if (present(op))    op_    = op
+        alpha_ = 1.0_${k1}$
+        if (present(alpha)) alpha_ = alpha
+        beta_ = 0.0_${k1}$
+        if (present(beta)) beta_ = beta
+        op_ = "N"
+        if (present(op)) op_ = op
         if (op_ == "H") op_ = "C"
 
         #:if t1.startswith('real')
@@ -167,11 +170,14 @@ submodule (stdlib_specialmatrices) sym_tridiagonal_matrices
         #:endif
 
         ! Prepare Lapack arguments.
-        n = A%n ; ldx = n ; ldy = n ;
+        n = A%n
+        ldx = n
+        ldy = n
         nrhs = #{if rank==1}# 1 #{else}# size(x, dim=2, kind=ilp) #{endif}#
 
         ! Pointer trick.
-        xmat(1:n, 1:nrhs) => x ; ymat(1:n, 1:nrhs) => y
+        xmat(1:n, 1:nrhs) => x
+        ymat(1:n, 1:nrhs) => y
         #:if t1.startswith('complex')
         call glagtm(op_, n, nrhs, alpha_, A%du, A%dv, A%du, xmat, ldx, beta_, ymat, ldy)
         #:else
@@ -207,13 +213,15 @@ submodule (stdlib_specialmatrices) sym_tridiagonal_matrices
             #:else
             allocate(B(n,n), source=zero_${k1}$)
             #:endif
-            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%du(i - 1)
                 B(i, i) = A%dv(i)
                 B(i, i + 1) = A%du(i)
             enddo
-            B(n , n -1) = A%du(n - 1); B(n, n) = A%dv(n)
+            B(n , n -1) = A%du(n - 1)
+            B(n, n) = A%dv(n)
         end associate
     end function
     #:endfor
@@ -246,15 +254,17 @@ submodule (stdlib_specialmatrices) sym_tridiagonal_matrices
         type(sym_tridiagonal_${s1}$_type), intent(in) :: A
         type(sym_tridiagonal_${s1}$_type) :: B
         B = sym_tridiagonal(A%du, A%dv)
-        B%du = alpha*B%du; B%dv = alpha*B%dv;
+        B%du = alpha*B%du
+        B%dv = alpha*B%dv
     end function
 
     pure module function scalar_multiplication_bis_sym_tridiagonal_${s1}$(A, alpha) result(B)
         type(sym_tridiagonal_${s1}$_type), intent(in) :: A
         ${t1}$, intent(in) :: alpha
         type(sym_tridiagonal_${s1}$_type) :: B
         B = sym_tridiagonal(A%du, A%dv)
-        B%du = alpha*B%du; B%dv = alpha*B%dv;
+        B%du = alpha*B%du
+        B%dv = alpha*B%dv
     end function
     #:endfor
 
@@ -264,15 +274,17 @@ submodule (stdlib_specialmatrices) sym_tridiagonal_matrices
         type(sym_tridiagonal_${s1}$_type), intent(in) :: B
         type(sym_tridiagonal_${s1}$_type) :: C
         C = sym_tridiagonal(A%du, A%dv)
-        C%du = C%du + B%du; C%dv = C%dv + B%dv
+        C%du = C%du + B%du
+        C%dv = C%dv + B%dv
     end function
 
     pure module function matrix_sub_sym_tridiagonal_${s1}$(A, B) result(C)
         type(sym_tridiagonal_${s1}$_type), intent(in) :: A
         type(sym_tridiagonal_${s1}$_type), intent(in) :: B
         type(sym_tridiagonal_${s1}$_type) :: C
         C = sym_tridiagonal(A%du, A%dv)
-        C%du = C%du - B%du; C%dv = C%dv - B%dv
+        C%du = C%du - B%du
+        C%dv = C%dv - B%dv
     end function
     #:endfor
 

src/specialmatrices/stdlib_specialmatrices_tridiagonal.fypp

@@ -164,9 +164,12 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
         end if
 
         ! Deal with optional arguments.
-        alpha_ = 1.0_${k1}$ ; if (present(alpha)) alpha_ = alpha
-        beta_  = 0.0_${k1}$ ; if (present(beta))  beta_  = beta
-        op_    = "N"        ; if (present(op))    op_    = op
+        alpha_ = 1.0_${k1}$
+        if (present(alpha)) alpha_ = alpha
+        beta_ = 0.0_${k1}$
+        if (present(beta)) beta_ = beta
+        op_ = "N"
+        if (present(op)) op_ = op
         if (op_ == "H") op_ = "C"
 
         #:if t1.startswith('real')
@@ -175,11 +178,14 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
         #:endif
 
         ! Prepare Lapack arguments.
-        n = A%n ; ldx = n ; ldy = n ;
+        n = A%n
+        ldx = n
+        ldy = n
         nrhs = #{if rank==1}# 1 #{else}# size(x, dim=2, kind=ilp) #{endif}#
 
         ! Pointer trick.
-        xmat(1:n, 1:nrhs) => x ; ymat(1:n, 1:nrhs) => y
+        xmat(1:n, 1:nrhs) => x
+        ymat(1:n, 1:nrhs) => y
         #:if t1.startswith('complex')
         call glagtm(op_, n, nrhs, alpha_, A%dl, A%dv, A%du, xmat, ldx, beta_, ymat, ldy)
         #:else
@@ -216,13 +222,15 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
         #:else
         allocate(B(n, n), source=zero_${k1}$)
         #:endif
-        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 associate
     end function
     #:endfor
@@ -255,15 +263,19 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
         type(tridiagonal_${s1}$_type), intent(in) :: A
         type(tridiagonal_${s1}$_type) :: B
         B = tridiagonal(A%dl, A%dv, A%du)
-        B%dl = alpha*B%dl; B%dv = alpha*B%dv; B%du = alpha*B%du
+        B%dl = alpha*B%dl
+        B%dv = alpha*B%dv
+        B%du = alpha*B%du
     end function
 
     pure module function scalar_multiplication_bis_tridiagonal_${s1}$(A, alpha) result(B)
         type(tridiagonal_${s1}$_type), intent(in) :: A
         ${t1}$, intent(in) :: alpha
         type(tridiagonal_${s1}$_type) :: B
         B = tridiagonal(A%dl, A%dv, A%du)
-        B%dl = alpha*B%dl; B%dv = alpha*B%dv; B%du = alpha*B%du
+        B%dl = alpha*B%dl
+        B%dv = alpha*B%dv
+        B%du = alpha*B%du
     end function
     #:endfor
 
@@ -273,15 +285,19 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
         type(tridiagonal_${s1}$_type), intent(in) :: B
         type(tridiagonal_${s1}$_type) :: C
         C = tridiagonal(A%dl, A%dv, A%du)
-        C%dl = C%dl + B%dl; C%dv = C%dv + B%dv; C%du = C%du + B%du
+        C%dl = C%dl + B%dl
+        C%dv = C%dv + B%dv
+        C%du = C%du + B%du
     end function
 
     pure module function matrix_sub_tridiagonal_${s1}$(A, B) result(C)
         type(tridiagonal_${s1}$_type), intent(in) :: A
         type(tridiagonal_${s1}$_type), intent(in) :: B
         type(tridiagonal_${s1}$_type) :: C
         C = tridiagonal(A%dl, A%dv, A%du)
-        C%dl = C%dl - B%dl; C%dv = C%dv - B%dv; C%du = C%du - B%du
+        C%dl = C%dl - B%dl
+        C%dv = C%dv - B%dv
+        C%du = C%du - B%du
     end function
     #:endfor