We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 76368a0 commit 3d10d39Copy full SHA for 3d10d39
src/stdlib_linalg_pinv.fypp
@@ -67,7 +67,7 @@ submodule(stdlib_linalg) stdlib_linalg_pseudoinverse
67
68
!> Discard singular values
69
cutoff = tolerance*maxval(s)
70
- s = merge(1/s,0.0_${rk}$,s>cutoff)
+ s = merge(1.0_${rk}$/s,0.0_${rk}$,s>cutoff)
71
72
! Get pseudo-inverse: A_pinv = V * (diag(1/s) * U^H) = V * (U * diag(1/s))^H
73
0 commit comments