replace loop by array operation

HugoMVale · Aug 3, 2024 · a01b3e9 · a01b3e9
1 parent d681db8
commit a01b3e9
Showing 1 changed file with 1 addition and 3 deletions.
diff --git a/src/odrpack_core.f90 b/src/odrpack_core.f90
@@ -4740,9 +4740,7 @@ subroutine dodstp &
       ! Do QR factorization (with column pivoting of TFJACB if ALPHA = 0)
       if (alpha == zero) then
          ipvt = 1
-         do k = 1, np
-            kpvt(k) = 0
-         end do
+         kpvt(1:np) = 0
       else
          ipvt = 0
       end if