Add symmetric tridiagonal matrices module by Mahmood-Sinan · Pull Request #1091 · fortran-lang/stdlib

Mahmood-Sinan · 2026-01-21T06:15:21Z

This draft PR adds support for symmetric tridiagonal matrices inside specialmatrices in stdlib.
It includes implementation, tests, documentation and examples.
The PR implements the same arithmetic and matrix operations as the existing tridiagonal module.

Remaining tasks:

Add validation for the op argument in spmv, following PR Add validation for op argument in spmv routine for special matrices #1087
(this will be added here once that implementation is merged).

Feedback is welcome.

codecov · 2026-01-21T06:51:02Z

Codecov Report

❌ Patch coverage is 0% with 7 lines in your changes missing coverage. Please review.
✅ Project coverage is 68.00%. Comparing base (dae9d4b) to head (d6b8572).
⚠️ Report is 49 commits behind head on master.

Files with missing lines	Patch %	Lines
...pecialmatrices/example_sym_tridiagonal_dp_type.f90	0.00%	7 Missing ⚠️

Additional details and impacted files

@@            Coverage Diff             @@
##           master    #1091      +/-   ##
==========================================
- Coverage   68.04%   68.00%   -0.04%     
==========================================
  Files         404      405       +1     
  Lines       12948    12955       +7     
  Branches     1395     1395              
==========================================
  Hits         8810     8810              
- Misses       4138     4145       +7

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

🚀 New features to boost your workflow:

❄️ Test Analytics: Detect flaky tests, report on failures, and find test suite problems.

loiseaujc · 2026-01-22T18:22:39Z

It looks pretty nice. One little thing though (even if I overlooked it when implementing the tridiagonal_type ) : you could define a single pure subroutine to construct the symtridiagonal matrices (make it two if you include the constant diagonal cases) and then have it called by the pure and impure constructors. It's a little thing but that would reduce a bit code duplication I guess.

I'll go through the code a bit more thoroughly tomorrow but it looks pretty good to me already.

…tants

Mahmood-Sinan · 2026-01-25T10:44:30Z

One little thing though (even if I overlooked it when implementing the tridiagonal_type ) : you could define a single pure subroutine to construct the symtridiagonal matrices (make it two if you include the constant diagonal cases) and then have it called by the pure and impure constructors. It's a little thing but that would reduce a bit code duplication I guess.

Thanks for the suggestion! You’re right, I implemented so. Now the tridiagonal and symtridiagonal construction is now handled by two pure helper routines(one from arrays and one from constants), which is called by both the pure and impure constructors.

…es directory

loiseaujc

I'll take a deeper look at it by the end of the week, but it looks quite good to me so far.

loiseaujc

LGTM.

…rogram name

loiseaujc · 2026-02-12T09:06:41Z

@perazz, @jalvesz and @jvdp1 : This PR is ready for merging. The only question I have left is a triviality though. What convention should we use for naming types : sym_tridiagonal_dp_type or symtridiagonal_dp_type ? My preference goes to the first one personally.

Mahmood-Sinan · 2026-02-20T16:35:40Z

@perazz, @jalvesz and @jvdp1 : This PR is ready for merging. The only question I have left is a triviality though. What convention should we use for naming types : sym_tridiagonal_dp_type or symtridiagonal_dp_type ? My preference goes to the first one personally.

@perazz @jalvesz @jvdp1 Please look into this.

jvdp1

Thank you @Mahmood-Sinan for this PR. Please find below some comments and suggestions.

jvdp1 · 2026-02-20T16:39:57Z

+
+`A = ` [[stdlib_specialmatrices(module):sym_tridiagonal(interface)]] `(du, dv)`
+
+- To construct a symmetric tridiagonal matrix of size `n x n` with constant diagonal elements `du` and `dv`:


by "constant" do you mean "scalar"

Yes, I changed it so scalar now.

jvdp1 · 2026-02-20T17:14:19Z

+            allocate(y1(n), source = 0.0_wp)  ; allocate(y2(n), source=0.0_wp)
+
+            ! Test y = A @ x
+            y1 = matmul(Amat, x) ; call spmv(A, x, y2)


Suggested change

y1 = matmul(Amat, x) ; call spmv(A, x, y2)

y1 = matmul(Amat, x)

call spmv(A, x, y2)

Mahmood-Sinan · 2026-02-21T08:55:19Z

Thank you @Mahmood-Sinan for this PR. Please find below some comments and suggestions.

Thank you so much for reviewing. I have made those changes. Please take a look at them.

Mahmood-Sinan · 2026-03-09T03:17:21Z

@jvdp1 I have made the changes you asked for, please check

jvdp1

Thank you @Mahmood-Sinan for this PR (sorry for my delayed review). LGTM, and this PR is ready to be merged IMO.

loiseaujc · 2026-03-30T06:28:02Z

I'll bother you for one last set of changes and then we can merge. Following #1162, can you make sure thesym_tridiagonal implementation handles the case of 1 x 1 matrices gracefully ?

loiseaujc · 2026-03-31T06:27:45Z

@Mahmood-Sinan : LGTM. Thanks for reacting rapidly to the 1x1 edge case.

@jvdp1 : Do you mind if I merge this today so that @Mahmood-Sinan can include it in their GSoC application ?

Mahmood-Sinan · 2026-04-03T12:40:09Z

@jvdp1 : Do you mind if I merge this today so that @Mahmood-Sinan can include it in their GSoC application ?

@jvdp1 please check.

jvdp1 · 2026-04-03T12:59:50Z

Sorry @Mahmood-Sinan and @loiseaujc for this late answer. I will merge it.

Mahmood-Sinan added 4 commits January 21, 2026 10:50

added function definitions

492c378

added function implementations and examples

ce3f60f

added tests for sym_tridiag and renamed the test suite name

624bdfa

docs: add symmetric tridiagonal matrices and cdp example

af1df17

add build functions for construction of matrices from arrays and cons…

eecc723

…tants

Mahmood-Sinan added 3 commits January 25, 2026 20:05

merge

bf53d95

merged the latest stdlib and moved sym_tridiagonal into specialmatric…

ab6b3d1

…es directory

move build constructors under interfaces

f692605

loiseaujc reviewed Jan 27, 2026

View reviewed changes

Comment thread src/specialmatrices/stdlib_specialmatrices.fypp

Comment thread src/specialmatrices/stdlib_specialmatrices_sym_tridiagonal.fypp Outdated

Comment thread src/specialmatrices/stdlib_specialmatrices_sym_tridiagonal.fypp Outdated

Mahmood-Sinan added 3 commits January 28, 2026 14:26

remove if condition for linalg_error_handling inside constructors

97464f4

modify the fypp logic inside spmv for rank based preprocessing

a17f62c

Merge remote-tracking branch 'upstream/master' into symtridiag

6e4291b

loiseaujc marked this pull request as ready for review February 5, 2026 09:52

loiseaujc approved these changes Feb 5, 2026

View reviewed changes

Mahmood-Sinan added 2 commits February 6, 2026 14:50

Merge remote-tracking branch 'upstream/master' into symtridiag

173283c

add: spmv error check for operator for sym tridiag, modify: example p…

fc95858

…rogram name

loiseaujc approved these changes Feb 9, 2026

View reviewed changes

jvdp1 reviewed Feb 20, 2026

View reviewed changes

Mahmood-Sinan added 7 commits February 20, 2026 23:21

Merge remote-tracking branch 'upstream/master' into symtridiag

70796c9

update: docs to add arguments and missing sentence

867cf4e

modify: example and src

7054d66

modify: add new lines between inline assignments

3dff7c4

add: new lines in test

5081221

add debug statements in test

2121952

change constant to scalar inside comments

154f312

Mahmood-Sinan added 3 commits February 21, 2026 11:54

modify: if(present()) to optval inside spmv

3dc8641

change: constant to scalar inside docs

dd191d7

modify: docs to add err argument

296c35a

jvdp1 approved these changes Mar 14, 2026

View reviewed changes

Mahmood-Sinan added 2 commits March 30, 2026 12:37

Merge remote-tracking branch 'upstream/master' into symtridiag

32f31dc

add: handle 1x1 case, test for 1x1 and test for arithmetic operators

d6b8572

loiseaujc approved these changes Mar 31, 2026

View reviewed changes

jvdp1 merged commit 3b447b9 into fortran-lang:master Apr 3, 2026
75 of 77 checks passed


		`A = ` [[stdlib_specialmatrices(module):sym_tridiagonal(interface)]] `(du, dv)`

		- To construct a symmetric tridiagonal matrix of size `n x n` with constant diagonal elements `du` and `dv`:

	y1 = matmul(Amat, x) ; call spmv(A, x, y2)
	y1 = matmul(Amat, x)
	call spmv(A, x, y2)

Conversation

Mahmood-Sinan commented Jan 21, 2026 • edited Loading Uh oh! There was an error while loading. Please reload this page.

Uh oh!

Uh oh!

codecov Bot commented Jan 21, 2026 • edited Loading Uh oh! There was an error while loading. Please reload this page.

Uh oh!

Codecov Report

Uh oh!

loiseaujc commented Jan 22, 2026

Uh oh!

Mahmood-Sinan commented Jan 25, 2026

Uh oh!

loiseaujc left a comment

Choose a reason for hiding this comment

Uh oh!

Uh oh!

Uh oh!

Uh oh!

loiseaujc left a comment

Choose a reason for hiding this comment

Uh oh!

loiseaujc commented Feb 12, 2026

Uh oh!

Mahmood-Sinan commented Feb 20, 2026

Uh oh!

jvdp1 left a comment

Choose a reason for hiding this comment

Uh oh!

Uh oh!

jvdp1 Feb 20, 2026

Choose a reason for hiding this comment

Uh oh!

Mahmood-Sinan Feb 21, 2026

Choose a reason for hiding this comment

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Uh oh!

jvdp1 Feb 20, 2026

Choose a reason for hiding this comment

Uh oh!

Uh oh!

Mahmood-Sinan commented Feb 21, 2026

Uh oh!

Mahmood-Sinan commented Mar 9, 2026

Uh oh!

jvdp1 left a comment

Choose a reason for hiding this comment

Uh oh!

loiseaujc commented Mar 30, 2026

Uh oh!

loiseaujc commented Mar 31, 2026

Uh oh!

Mahmood-Sinan commented Apr 3, 2026

Uh oh!

jvdp1 commented Apr 3, 2026

Uh oh!

Uh oh!

Reviewers

Assignees

Labels

Projects

Milestone

Development

Uh oh!

3 participants

Mahmood-Sinan commented Jan 21, 2026 •

edited

Loading

codecov Bot commented Jan 21, 2026 •

edited

Loading