phynexis-deps/lapack
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

remove a TEMP = that is good to make the real arithmetic code looks like the complex arithmetic one but weird otherwise

Browse Source
This commit is contained in:
Julien Langou
2026-07-12 10:08:38 -06:00
parent e8785e4e1e
commit 301097afd8
2 changed files with 4 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files
BLAS/SRC/dtrsm.f
+2 -4
View File
@@ -377,9 +377,8 @@
270 CONTINUE
END IF
DO 290 J = 1,K - 1
TEMP = A(J,K)
DO 280 I = 1,M
B(I,J) = B(I,J) - TEMP*B(I,K)
B(I,J) = B(I,J) - A(J,K)*B(I,K)
280 CONTINUE
290 CONTINUE
DO 300 I = 1,M
@@ -394,9 +393,8 @@
320 CONTINUE
END IF
DO 340 J = K + 1,N
TEMP = A(J,K)
DO 330 I = 1,M
B(I,J) = B(I,J) - TEMP*B(I,K)
B(I,J) = B(I,J) - A(J,K)*B(I,K)
330 CONTINUE
340 CONTINUE
DO 350 I = 1,M
BLAS/SRC/strsm.f
+2 -4
View File
@@ -377,9 +377,8 @@
270 CONTINUE
END IF
DO 290 J = 1,K - 1
TEMP = A(J,K)
DO 280 I = 1,M
B(I,J) = B(I,J) - TEMP*B(I,K)
B(I,J) = B(I,J) - A(J,K)*B(I,K)
280 CONTINUE
290 CONTINUE
DO 300 I = 1,M
@@ -394,9 +393,8 @@
320 CONTINUE
END IF
DO 340 J = K + 1,N
TEMP = A(J,K)
DO 330 I = 1,M
B(I,J) = B(I,J) - TEMP*B(I,K)
B(I,J) = B(I,J) - A(J,K)*B(I,K)
330 CONTINUE
340 CONTINUE
DO 350 I = 1,M