Skip to content

Commit f6a2c88

Browse files
committed
Fix conjugate-transpose spmv.
1 parent a42928f commit f6a2c88

File tree

2 files changed

+2
-0
lines changed

2 files changed

+2
-0
lines changed

src/stdlib_specialmatrices_symtridiagonal.fypp

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -156,6 +156,7 @@ submodule (stdlib_specialmatrices) symtridiagonal_matrices
156156
alpha_ = 1.0_${k1}$ ; if (present(alpha)) alpha_ = alpha
157157
beta_ = 0.0_${k1}$ ; if (present(beta)) beta_ = beta
158158
op_ = "N" ; if (present(op)) op_ = op
159+
if (op_ == "H") op_ = "C"
159160

160161
! Prepare Lapack arguments.
161162
n = A%n ; ldx = n ; ldy = n ; y = 0.0_${k1}$

src/stdlib_specialmatrices_tridiagonal.fypp

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -166,6 +166,7 @@ submodule (stdlib_specialmatrices) tridiagonal_matrices
166166
alpha_ = 1.0_${k1}$ ; if (present(alpha)) alpha_ = alpha
167167
beta_ = 0.0_${k1}$ ; if (present(beta)) beta_ = beta
168168
op_ = "N" ; if (present(op)) op_ = op
169+
if (op_ == "H") op_ = "C"
169170

170171
! Prepare Lapack arguments.
171172
n = A%n ; ldx = n ; ldy = n ; y = 0.0_${k1}$

0 commit comments

Comments
 (0)