@@ -584,7 +584,7 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
584
584
585
585
#:for k1, t1, s1 in (KINDS_TYPES)
586
586
! Tridiag + Tridiag = Tridiag
587
- pure function matrix_add_tridiag_tridiag_${s1}$(A, B) result(C)
587
+ pure module function matrix_add_tridiag_tridiag_${s1}$(A, B) result(C)
588
588
type(tridiagonal_${s1}$_type), intent(in) :: A, B
589
589
type(tridiagonal_${s1}$_type) :: C
590
590
C = tridiagonal(A%dl, A%dv, A%du)
@@ -668,7 +668,7 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
668
668
669
669
#:for k1, t1, s1 in (KINDS_TYPES)
670
670
! Tridiag - Tridiag = Tridiag
671
- pure function matrix_sub_tridiag_tridiag_${s1}$(A, B) result(C)
671
+ pure module function matrix_sub_tridiag_tridiag_${s1}$(A, B) result(C)
672
672
type(tridiagonal_${s1}$_type), intent(in) :: A, B
673
673
type(tridiagonal_${s1}$_type) :: C
674
674
C = tridiagonal(A%dl, A%dv, A%du)
@@ -746,7 +746,6 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
746
746
C%dl = C%dl - B%dl ; C%dv = C%dv - B%dv ; C%du = conjg(C%dl)
747
747
end function
748
748
#:endif
749
-
750
749
#:endfor
751
750
752
751
end submodule
0 commit comments