Update src/stdlib_linalg_pinv.fypp

Author: perazz

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