Replace WAXPY and WSCAL in routines WGEFA and WGESL

util/blas.f90
- Replaced WAXPY and WSCAL with explicit loops (generated by AI)
  in the WGEFA and WGESL routines.  These are needed for some
  integrators.

Signed-off-by: Bob Yantosca <[email protected]>

Author: yantosca

util/blas.f90

@@ -124,13 +124,19 @@ SUBROUTINE WGEFA(N,A,Ipvt,info)
                t = A(l,k); A(l,k) = A(k,k); A(k,k) = t
             END IF
             t = -ONE/A(k,k)
-            CALL WSCAL(n-k,t,A(k+1,k),1)
+!            CALL WSCAL(n-k,t,A(k+1,k),1)
+            DO i = k+1, n
+               A(i,k) = t * A(i,k)
+            END DO
             DO j = k+1, n
                t = A(l,j)
                IF (l /= k) THEN
                   A(l,j) = A(k,j); A(k,j) = t
                END IF
-               CALL WAXPY(n-k,t,A(k+1,k),1,A(k+1,j),1)
+               !CALL WAXPY(n-k,t,A(k+1,k),1,A(k+1,j),1)
+               DO i = k+1, n
+                  A(i,j) = A(i,j) + t * A(i,k)
+               END DO
             END DO         
          END IF
 
@@ -176,15 +182,21 @@ SUBROUTINE WGESL(Trans,N,A,Ipvt,b)
                b(l) = b(k)
                b(k) = t
             END IF
-            CALL WAXPY(n-k,t,a(k+1,k),1,b(k+1),1)
+            !CALL WAXPY(n-k,t,a(k+1,k),1,b(k+1),1)
+            DO i = k+1, n
+               b(i) = b(i) + t * a(i,k)
+            END DO
           END DO
          END IF
 !        now solve  U*x = y
          DO kb = 1, n
             k = n + 1 - kb
             b(k) = b(k)/a(k,k)
             t = -b(k)
-            CALL WAXPY(k-1,t,a(1,k),1,b(1),1)
+            !CALL WAXPY(k-1,t,a(1,k),1,b(1),1)
+            DO i = 1, k-1
+               b(i) = b(i) + t * a(i,k)
+            END DO
          END DO
 
       CASE ('t','T')  !  Solve transpose(A) * x = b