Update src/stdlib_linalg_pinv.fypp

fortran-lang · Dec 18, 2024 · 38bae4d · 38bae4d
1 parent 3d10d39
commit 38bae4d
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/stdlib_linalg_pinv.fypp b/src/stdlib_linalg_pinv.fypp
@@ -72,7 +72,7 @@ submodule(stdlib_linalg) stdlib_linalg_pseudoinverse
          ! Get pseudo-inverse: A_pinv = V * (diag(1/s) * U^H) = V * (U * diag(1/s))^H
 
          ! 1) compute (U * diag(1/s)) in-place
-         forall (i=1:m,j=1:k) u(i,j) = s(j)*u(i,j)
+         do concurrent (i=1:m,j=1:k); u(i,j) = s(j)*u(i,j); end do
 
          ! 2) commutate matmul: A_pinv = V^H * (U * diag(1/s))^H = ((U * diag(1/s)) * V^H)^H. 
          !    This avoids one matrix transpose