]> git.pld-linux.org Git - packages/lapack.git/blame - blas-nan.patch
projects / packages / lapack.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
- updated to 3.11.0
[packages/lapack.git] / blas-nan.patch
CommitLineData
8a7beb04
JB		 1Disable optimizations based on (0*x = 0) assumption;
2they make the code non IEC-60559/IEEE-754 compliant wrt. NaN (0*x=0 _unless_ x=NaN)
8a7beb04
JB		 3--- lapack-3.2.2/BLAS/SRC/zgerc.f.orig 2009-01-02 22:57:02.000000000 +0100
4+++ lapack-3.2.2/BLAS/SRC/zgerc.f 2010-10-19 09:27:26.711744396 +0200
5@@ -128,12 +128,12 @@
6 END IF
7 IF (INCX.EQ.1) THEN
8 DO 20 J = 1,N
9- IF (Y(JY).NE.ZERO) THEN
10+c IF (Y(JY).NE.ZERO) THEN
11 TEMP = ALPHA*DCONJG(Y(JY))
12 DO 10 I = 1,M
13 A(I,J) = A(I,J) + X(I)*TEMP
14 10 CONTINUE
15- END IF
16+c END IF
17 JY = JY + INCY
18 20 CONTINUE
19 ELSE
20@@ -143,14 +143,14 @@
21 KX = 1 - (M-1)*INCX
22 END IF
23 DO 40 J = 1,N
24- IF (Y(JY).NE.ZERO) THEN
25+c IF (Y(JY).NE.ZERO) THEN
26 TEMP = ALPHA*DCONJG(Y(JY))
27 IX = KX
28 DO 30 I = 1,M
29 A(I,J) = A(I,J) + X(IX)*TEMP
30 IX = IX + INCX
31 30 CONTINUE
32- END IF
33+c END IF
34 JY = JY + INCY
35 40 CONTINUE
36 END IF
37--- lapack-3.2.2/BLAS/SRC/zher2.f.orig 2009-01-02 22:57:02.000000000 +0100
38+++ lapack-3.2.2/BLAS/SRC/zher2.f 2010-10-19 09:29:15.131742720 +0200
39@@ -171,7 +171,7 @@
40 *
41 IF ((INCX.EQ.1) .AND. (INCY.EQ.1)) THEN
42 DO 20 J = 1,N
43- IF ((X(J).NE.ZERO) .OR. (Y(J).NE.ZERO)) THEN
44+c IF ((X(J).NE.ZERO) .OR. (Y(J).NE.ZERO)) THEN
45 TEMP1 = ALPHA*DCONJG(Y(J))
46 TEMP2 = DCONJG(ALPHA*X(J))
47 DO 10 I = 1,J - 1
48@@ -179,13 +179,13 @@
49 10 CONTINUE
50 A(J,J) = DBLE(A(J,J)) +
51 + DBLE(X(J)*TEMP1+Y(J)*TEMP2)
52- ELSE
53- A(J,J) = DBLE(A(J,J))
54- END IF
55+c ELSE
56+c A(J,J) = DBLE(A(J,J))
57+c END IF
58 20 CONTINUE
59 ELSE
60 DO 40 J = 1,N
61- IF ((X(JX).NE.ZERO) .OR. (Y(JY).NE.ZERO)) THEN
62+c IF ((X(JX).NE.ZERO) .OR. (Y(JY).NE.ZERO)) THEN
63 TEMP1 = ALPHA*DCONJG(Y(JY))
64 TEMP2 = DCONJG(ALPHA*X(JX))
65 IX = KX
66@@ -197,9 +197,9 @@
67 30 CONTINUE
68 A(J,J) = DBLE(A(J,J)) +
69 + DBLE(X(JX)*TEMP1+Y(JY)*TEMP2)
70- ELSE
71- A(J,J) = DBLE(A(J,J))
72- END IF
73+c ELSE
74+c A(J,J) = DBLE(A(J,J))
75+c END IF
76 JX = JX + INCX
77 JY = JY + INCY
78 40 CONTINUE
79@@ -210,7 +210,7 @@
80 *
81 IF ((INCX.EQ.1) .AND. (INCY.EQ.1)) THEN
82 DO 60 J = 1,N
83- IF ((X(J).NE.ZERO) .OR. (Y(J).NE.ZERO)) THEN
84+c IF ((X(J).NE.ZERO) .OR. (Y(J).NE.ZERO)) THEN
85 TEMP1 = ALPHA*DCONJG(Y(J))
86 TEMP2 = DCONJG(ALPHA*X(J))
87 A(J,J) = DBLE(A(J,J)) +
88@@ -218,13 +218,13 @@
89 DO 50 I = J + 1,N
90 A(I,J) = A(I,J) + X(I)*TEMP1 + Y(I)*TEMP2
91 50 CONTINUE
92- ELSE
93- A(J,J) = DBLE(A(J,J))
94- END IF
95+c ELSE
96+c A(J,J) = DBLE(A(J,J))
97+c END IF
98 60 CONTINUE
99 ELSE
100 DO 80 J = 1,N
101- IF ((X(JX).NE.ZERO) .OR. (Y(JY).NE.ZERO)) THEN
102+c IF ((X(JX).NE.ZERO) .OR. (Y(JY).NE.ZERO)) THEN
103 TEMP1 = ALPHA*DCONJG(Y(JY))
104 TEMP2 = DCONJG(ALPHA*X(JX))
105 A(J,J) = DBLE(A(J,J)) +
106@@ -236,9 +236,9 @@
107 IY = IY + INCY
108 A(I,J) = A(I,J) + X(IX)*TEMP1 + Y(IY)*TEMP2
109 70 CONTINUE
110- ELSE
111- A(J,J) = DBLE(A(J,J))
112- END IF
113+c ELSE
114+c A(J,J) = DBLE(A(J,J))
115+c END IF
116 JX = JX + INCX
117 JY = JY + INCY
118 80 CONTINUE
119--- lapack-3.2.2/BLAS/SRC/zher2k.f.orig 2009-01-02 22:57:02.000000000 +0100
120+++ lapack-3.2.2/BLAS/SRC/zher2k.f 2010-10-19 09:29:55.943741043 +0200
121@@ -260,7 +260,7 @@
122 C(J,J) = DBLE(C(J,J))
123 END IF
124 DO 120 L = 1,K
125- IF ((A(J,L).NE.ZERO) .OR. (B(J,L).NE.ZERO)) THEN
126+c IF ((A(J,L).NE.ZERO) .OR. (B(J,L).NE.ZERO)) THEN
127 TEMP1 = ALPHA*DCONJG(B(J,L))
128 TEMP2 = DCONJG(ALPHA*A(J,L))
129 DO 110 I = 1,J - 1
130@@ -269,7 +269,7 @@
131 110 CONTINUE
132 C(J,J) = DBLE(C(J,J)) +
133 + DBLE(A(J,L)*TEMP1+B(J,L)*TEMP2)
134- END IF
135+c END IF
136 120 CONTINUE
137 130 CONTINUE
138 ELSE
139@@ -287,7 +287,7 @@
140 C(J,J) = DBLE(C(J,J))
141 END IF
142 DO 170 L = 1,K
143- IF ((A(J,L).NE.ZERO) .OR. (B(J,L).NE.ZERO)) THEN
144+c IF ((A(J,L).NE.ZERO) .OR. (B(J,L).NE.ZERO)) THEN
145 TEMP1 = ALPHA*DCONJG(B(J,L))
146 TEMP2 = DCONJG(ALPHA*A(J,L))
147 DO 160 I = J + 1,N
148@@ -296,7 +296,7 @@
149 160 CONTINUE
150 C(J,J) = DBLE(C(J,J)) +
151 + DBLE(A(J,L)*TEMP1+B(J,L)*TEMP2)
152- END IF
153+c END IF
154 170 CONTINUE
155 180 CONTINUE
156 END IF
157--- lapack-3.2.2/BLAS/SRC/ztrmm.f.orig 2009-01-02 22:57:02.000000000 +0100
158+++ lapack-3.2.2/BLAS/SRC/ztrmm.f 2010-10-19 09:31:16.815739367 +0200
159@@ -211,27 +211,27 @@
160 IF (UPPER) THEN
161 DO 50 J = 1,N
162 DO 40 K = 1,M
163- IF (B(K,J).NE.ZERO) THEN
164+c IF (B(K,J).NE.ZERO) THEN
165 TEMP = ALPHA*B(K,J)
166 DO 30 I = 1,K - 1
167 B(I,J) = B(I,J) + TEMP*A(I,K)
168 30 CONTINUE
169 IF (NOUNIT) TEMP = TEMP*A(K,K)
170 B(K,J) = TEMP
171- END IF
172+c END IF
173 40 CONTINUE
174 50 CONTINUE
175 ELSE
176 DO 80 J = 1,N
177 DO 70 K = M,1,-1
178- IF (B(K,J).NE.ZERO) THEN
179+c IF (B(K,J).NE.ZERO) THEN
180 TEMP = ALPHA*B(K,J)
181 B(K,J) = TEMP
182 IF (NOUNIT) B(K,J) = B(K,J)*A(K,K)
183 DO 60 I = K + 1,M
184 B(I,J) = B(I,J) + TEMP*A(I,K)
185 60 CONTINUE
186- END IF
187+c END IF
188 70 CONTINUE
189 80 CONTINUE
190 END IF
191@@ -290,12 +290,12 @@
192 B(I,J) = TEMP*B(I,J)
193 170 CONTINUE
194 DO 190 K = 1,J - 1
195- IF (A(K,J).NE.ZERO) THEN
196+c IF (A(K,J).NE.ZERO) THEN
197 TEMP = ALPHA*A(K,J)
198 DO 180 I = 1,M
199 B(I,J) = B(I,J) + TEMP*B(I,K)
200 180 CONTINUE
201- END IF
202+c END IF
203 190 CONTINUE
204 200 CONTINUE
205 ELSE
206@@ -306,12 +306,12 @@
207 B(I,J) = TEMP*B(I,J)
208 210 CONTINUE
209 DO 230 K = J + 1,N
210- IF (A(K,J).NE.ZERO) THEN
211+c IF (A(K,J).NE.ZERO) THEN
212 TEMP = ALPHA*A(K,J)
213 DO 220 I = 1,M
214 B(I,J) = B(I,J) + TEMP*B(I,K)
215 220 CONTINUE
216- END IF
217+c END IF
218 230 CONTINUE
219 240 CONTINUE
220 END IF
221@@ -350,7 +350,7 @@
222 ELSE
223 DO 320 K = N,1,-1
224 DO 300 J = K + 1,N
225- IF (A(J,K).NE.ZERO) THEN
226+c IF (A(J,K).NE.ZERO) THEN
227 IF (NOCONJ) THEN
228 TEMP = ALPHA*A(J,K)
229 ELSE
230@@ -359,7 +359,7 @@
231 DO 290 I = 1,M
232 B(I,J) = B(I,J) + TEMP*B(I,K)
233 290 CONTINUE
234- END IF
235+c END IF
236 300 CONTINUE
237 TEMP = ALPHA
238 IF (NOUNIT) THEN
239--- lapack-3.2.2/BLAS/SRC/ztrmv.f.orig 2009-01-02 22:57:02.000000000 +0100
240+++ lapack-3.2.2/BLAS/SRC/ztrmv.f 2010-10-19 09:32:09.067739367 +0200
241@@ -169,18 +169,18 @@
242 IF (LSAME(UPLO,'U')) THEN
243 IF (INCX.EQ.1) THEN
244 DO 20 J = 1,N
245- IF (X(J).NE.ZERO) THEN
246+c IF (X(J).NE.ZERO) THEN
247 TEMP = X(J)
248 DO 10 I = 1,J - 1
249 X(I) = X(I) + TEMP*A(I,J)
250 10 CONTINUE
251 IF (NOUNIT) X(J) = X(J)*A(J,J)
252- END IF
253+c END IF
254 20 CONTINUE
255 ELSE
256 JX = KX
257 DO 40 J = 1,N
258- IF (X(JX).NE.ZERO) THEN
259+c IF (X(JX).NE.ZERO) THEN
260 TEMP = X(JX)
261 IX = KX
262 DO 30 I = 1,J - 1
263@@ -188,26 +188,26 @@
264 IX = IX + INCX
265 30 CONTINUE
266 IF (NOUNIT) X(JX) = X(JX)*A(J,J)
267- END IF
268+c END IF
269 JX = JX + INCX
270 40 CONTINUE
271 END IF
272 ELSE
273 IF (INCX.EQ.1) THEN
274 DO 60 J = N,1,-1
275- IF (X(J).NE.ZERO) THEN
276+c IF (X(J).NE.ZERO) THEN
277 TEMP = X(J)
278 DO 50 I = N,J + 1,-1
279 X(I) = X(I) + TEMP*A(I,J)
280 50 CONTINUE
281 IF (NOUNIT) X(J) = X(J)*A(J,J)
282- END IF
283+c END IF
284 60 CONTINUE
285 ELSE
286 KX = KX + (N-1)*INCX
287 JX = KX
288 DO 80 J = N,1,-1
289- IF (X(JX).NE.ZERO) THEN
290+c IF (X(JX).NE.ZERO) THEN
291 TEMP = X(JX)
292 IX = KX
293 DO 70 I = N,J + 1,-1
294@@ -215,7 +215,7 @@
295 IX = IX - INCX
296 70 CONTINUE
297 IF (NOUNIT) X(JX) = X(JX)*A(J,J)
298- END IF
299+c END IF
300 JX = JX - INCX
301 80 CONTINUE
302 END IF
303--- lapack-3.2.2/BLAS/SRC/ztrsm.f.orig 2009-01-02 22:57:02.000000000 +0100
304+++ lapack-3.2.2/BLAS/SRC/ztrsm.f 2010-10-19 09:37:56.767745234 +0200
305@@ -218,12 +218,12 @@
306 30 CONTINUE
307 END IF
308 DO 50 K = M,1,-1
309- IF (B(K,J).NE.ZERO) THEN
310+c IF (B(K,J).NE.ZERO) THEN
311 IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
312 DO 40 I = 1,K - 1
313 B(I,J) = B(I,J) - B(K,J)*A(I,K)
314 40 CONTINUE
315- END IF
316+c END IF
317 50 CONTINUE
318 60 CONTINUE
319 ELSE
320@@ -234,12 +234,12 @@
321 70 CONTINUE
322 END IF
323 DO 90 K = 1,M
324- IF (B(K,J).NE.ZERO) THEN
325+c IF (B(K,J).NE.ZERO) THEN
326 IF (NOUNIT) B(K,J) = B(K,J)/A(K,K)
327 DO 80 I = K + 1,M
328 B(I,J) = B(I,J) - B(K,J)*A(I,K)
329 80 CONTINUE
330- END IF
331+c END IF
332 90 CONTINUE
333 100 CONTINUE
334 END IF
335@@ -299,11 +299,11 @@
336 190 CONTINUE
337 END IF
338 DO 210 K = 1,J - 1
339- IF (A(K,J).NE.ZERO) THEN
340+c IF (A(K,J).NE.ZERO) THEN
341 DO 200 I = 1,M
342 B(I,J) = B(I,J) - A(K,J)*B(I,K)
343 200 CONTINUE
344- END IF
345+c END IF
346 210 CONTINUE
347 IF (NOUNIT) THEN
348 TEMP = ONE/A(J,J)
349@@ -320,11 +320,11 @@
350 240 CONTINUE
351 END IF
352 DO 260 K = J + 1,N
353- IF (A(K,J).NE.ZERO) THEN
354+c IF (A(K,J).NE.ZERO) THEN
355 DO 250 I = 1,M
356 B(I,J) = B(I,J) - A(K,J)*B(I,K)
357 250 CONTINUE
358- END IF
359+c END IF
360 260 CONTINUE
361 IF (NOUNIT) THEN
362 TEMP = ONE/A(J,J)
363@@ -352,7 +352,7 @@
364 290 CONTINUE
365 END IF
366 DO 310 J = 1,K - 1
367- IF (A(J,K).NE.ZERO) THEN
368+c IF (A(J,K).NE.ZERO) THEN
369 IF (NOCONJ) THEN
370 TEMP = A(J,K)
371 ELSE
372@@ -361,7 +361,7 @@
373 DO 300 I = 1,M
374 B(I,J) = B(I,J) - TEMP*B(I,K)
375 300 CONTINUE
376- END IF
377+c END IF
378 310 CONTINUE
379 IF (ALPHA.NE.ONE) THEN
380 DO 320 I = 1,M
381@@ -382,7 +382,7 @@
382 340 CONTINUE
383 END IF
384 DO 360 J = K + 1,N
385- IF (A(J,K).NE.ZERO) THEN
386+c IF (A(J,K).NE.ZERO) THEN
387 IF (NOCONJ) THEN
388 TEMP = A(J,K)
389 ELSE
390@@ -391,7 +391,7 @@
391 DO 350 I = 1,M
392 B(I,J) = B(I,J) - TEMP*B(I,K)
393 350 CONTINUE
394- END IF
395+c END IF
396 360 CONTINUE
397 IF (ALPHA.NE.ONE) THEN
398 DO 370 I = 1,M
399--- lapack-3.2.2/BLAS/SRC/ztrsv.f.orig 2009-01-02 22:57:02.000000000 +0100
400+++ lapack-3.2.2/BLAS/SRC/ztrsv.f 2010-10-19 09:40:07.839740205 +0200
401@@ -172,18 +172,18 @@
402 IF (LSAME(UPLO,'U')) THEN
403 IF (INCX.EQ.1) THEN
404 DO 20 J = N,1,-1
405- IF (X(J).NE.ZERO) THEN
406+c IF (X(J).NE.ZERO) THEN
407 IF (NOUNIT) X(J) = X(J)/A(J,J)
408 TEMP = X(J)
409 DO 10 I = J - 1,1,-1
410 X(I) = X(I) - TEMP*A(I,J)
411 10 CONTINUE
412- END IF
413+c END IF
414 20 CONTINUE
415 ELSE
416 JX = KX + (N-1)*INCX
417 DO 40 J = N,1,-1
418- IF (X(JX).NE.ZERO) THEN
419+c IF (X(JX).NE.ZERO) THEN
420 IF (NOUNIT) X(JX) = X(JX)/A(J,J)
421 TEMP = X(JX)
422 IX = JX
423@@ -191,25 +191,25 @@
424 IX = IX - INCX
425 X(IX) = X(IX) - TEMP*A(I,J)
426 30 CONTINUE
427- END IF
428+c END IF
429 JX = JX - INCX
430 40 CONTINUE
431 END IF
432 ELSE
433 IF (INCX.EQ.1) THEN
434 DO 60 J = 1,N
435- IF (X(J).NE.ZERO) THEN
436+c IF (X(J).NE.ZERO) THEN
437 IF (NOUNIT) X(J) = X(J)/A(J,J)
438 TEMP = X(J)
439 DO 50 I = J + 1,N
440 X(I) = X(I) - TEMP*A(I,J)
441 50 CONTINUE
442- END IF
443+c END IF
444 60 CONTINUE
445 ELSE
446 JX = KX
447 DO 80 J = 1,N
448- IF (X(JX).NE.ZERO) THEN
449+c IF (X(JX).NE.ZERO) THEN
450 IF (NOUNIT) X(JX) = X(JX)/A(J,J)
451 TEMP = X(JX)
452 IX = JX
453@@ -217,7 +217,7 @@
454 IX = IX + INCX
455 X(IX) = X(IX) - TEMP*A(I,J)
456 70 CONTINUE
457- END IF
458+c END IF
459 JX = JX + INCX
460 80 CONTINUE
461 END IF
8a7beb04
JB		 462--- lapack-3.2.2/BLAS/SRC/zgeru.f.orig 2009-01-02 22:57:02.000000000 +0100
463+++ lapack-3.2.2/BLAS/SRC/zgeru.f 2010-10-19 09:41:25.835741043 +0200
464@@ -128,12 +128,12 @@
465 END IF
466 IF (INCX.EQ.1) THEN
467 DO 20 J = 1,N
468- IF (Y(JY).NE.ZERO) THEN
469+c IF (Y(JY).NE.ZERO) THEN
470 TEMP = ALPHA*Y(JY)
471 DO 10 I = 1,M
472 A(I,J) = A(I,J) + X(I)*TEMP
473 10 CONTINUE
474- END IF
475+c END IF
476 JY = JY + INCY
477 20 CONTINUE
478 ELSE
479@@ -143,14 +143,14 @@
480 KX = 1 - (M-1)*INCX
481 END IF
482 DO 40 J = 1,N
483- IF (Y(JY).NE.ZERO) THEN
484+c IF (Y(JY).NE.ZERO) THEN
485 TEMP = ALPHA*Y(JY)
486 IX = KX
487 DO 30 I = 1,M
488 A(I,J) = A(I,J) + X(IX)*TEMP
489 IX = IX + INCX
490 30 CONTINUE
491- END IF
492+c END IF
493 JY = JY + INCY
494 40 CONTINUE
495 END IF
496--- lapack-3.2.2/BLAS/SRC/zher.f.orig 2009-01-02 22:57:02.000000000 +0100
497+++ lapack-3.2.2/BLAS/SRC/zher.f 2010-10-19 09:43:12.851741882 +0200
498@@ -148,20 +148,20 @@
499 *
500 IF (INCX.EQ.1) THEN
501 DO 20 J = 1,N
502- IF (X(J).NE.ZERO) THEN
503+c IF (X(J).NE.ZERO) THEN
504 TEMP = ALPHA*DCONJG(X(J))
505 DO 10 I = 1,J - 1
506 A(I,J) = A(I,J) + X(I)*TEMP
507 10 CONTINUE
508 A(J,J) = DBLE(A(J,J)) + DBLE(X(J)*TEMP)
509- ELSE
510- A(J,J) = DBLE(A(J,J))
511- END IF
512+c ELSE
513+c A(J,J) = DBLE(A(J,J))
514+c END IF
515 20 CONTINUE
516 ELSE
517 JX = KX
518 DO 40 J = 1,N
519- IF (X(JX).NE.ZERO) THEN
520+c IF (X(JX).NE.ZERO) THEN
521 TEMP = ALPHA*DCONJG(X(JX))
522 IX = KX
523 DO 30 I = 1,J - 1
524@@ -169,9 +169,9 @@
525 IX = IX + INCX
526 30 CONTINUE
527 A(J,J) = DBLE(A(J,J)) + DBLE(X(JX)*TEMP)
528- ELSE
529- A(J,J) = DBLE(A(J,J))
530- END IF
531+c ELSE
532+c A(J,J) = DBLE(A(J,J))
533+c END IF
534 JX = JX + INCX
535 40 CONTINUE
536 END IF
537@@ -181,20 +181,20 @@
538 *
539 IF (INCX.EQ.1) THEN
540 DO 60 J = 1,N
541- IF (X(J).NE.ZERO) THEN
542+c IF (X(J).NE.ZERO) THEN
543 TEMP = ALPHA*DCONJG(X(J))
544 A(J,J) = DBLE(A(J,J)) + DBLE(TEMP*X(J))
545 DO 50 I = J + 1,N
546 A(I,J) = A(I,J) + X(I)*TEMP
547 50 CONTINUE
548- ELSE
549- A(J,J) = DBLE(A(J,J))
550- END IF
551+c ELSE
552+c A(J,J) = DBLE(A(J,J))
553+c END IF
554 60 CONTINUE
555 ELSE
556 JX = KX
557 DO 80 J = 1,N
558- IF (X(JX).NE.ZERO) THEN
559+c IF (X(JX).NE.ZERO) THEN
560 TEMP = ALPHA*DCONJG(X(JX))
561 A(J,J) = DBLE(A(J,J)) + DBLE(TEMP*X(JX))
562 IX = JX
563@@ -202,9 +202,9 @@
564 IX = IX + INCX
565 A(I,J) = A(I,J) + X(IX)*TEMP
566 70 CONTINUE
567- ELSE
568- A(J,J) = DBLE(A(J,J))
569- END IF
570+c ELSE
571+c A(J,J) = DBLE(A(J,J))
572+c END IF
573 JX = JX + INCX
574 80 CONTINUE
575 END IF
576--- lapack-3.2.2/BLAS/SRC/zhpr.f.orig 2009-01-02 22:57:02.000000000 +0100
577+++ lapack-3.2.2/BLAS/SRC/zhpr.f 2010-10-19 09:44:30.431739367 +0200
578@@ -143,7 +143,7 @@
579 *
580 IF (INCX.EQ.1) THEN
581 DO 20 J = 1,N
582- IF (X(J).NE.ZERO) THEN
583+c IF (X(J).NE.ZERO) THEN
584 TEMP = ALPHA*DCONJG(X(J))
585 K = KK
586 DO 10 I = 1,J - 1
587@@ -151,15 +151,15 @@
588 K = K + 1
589 10 CONTINUE
590 AP(KK+J-1) = DBLE(AP(KK+J-1)) + DBLE(X(J)*TEMP)
591- ELSE
592- AP(KK+J-1) = DBLE(AP(KK+J-1))
593- END IF
594+c ELSE
595+c AP(KK+J-1) = DBLE(AP(KK+J-1))
596+c END IF
597 KK = KK + J
598 20 CONTINUE
599 ELSE
600 JX = KX
601 DO 40 J = 1,N
602- IF (X(JX).NE.ZERO) THEN
603+c IF (X(JX).NE.ZERO) THEN
604 TEMP = ALPHA*DCONJG(X(JX))
605 IX = KX
606 DO 30 K = KK,KK + J - 2
607@@ -167,9 +167,9 @@
608 IX = IX + INCX
609 30 CONTINUE
610 AP(KK+J-1) = DBLE(AP(KK+J-1)) + DBLE(X(JX)*TEMP)
611- ELSE
612- AP(KK+J-1) = DBLE(AP(KK+J-1))
613- END IF
614+c ELSE
615+c AP(KK+J-1) = DBLE(AP(KK+J-1))
616+c END IF
617 JX = JX + INCX
618 KK = KK + J
619 40 CONTINUE
620@@ -180,7 +180,7 @@
621 *
622 IF (INCX.EQ.1) THEN
623 DO 60 J = 1,N
624- IF (X(J).NE.ZERO) THEN
625+c IF (X(J).NE.ZERO) THEN
626 TEMP = ALPHA*DCONJG(X(J))
627 AP(KK) = DBLE(AP(KK)) + DBLE(TEMP*X(J))
628 K = KK + 1
629@@ -188,15 +188,15 @@
630 AP(K) = AP(K) + X(I)*TEMP
631 K = K + 1
632 50 CONTINUE
633- ELSE
634- AP(KK) = DBLE(AP(KK))
635- END IF
636+c ELSE
637+c AP(KK) = DBLE(AP(KK))
638+c END IF
639 KK = KK + N - J + 1
640 60 CONTINUE
641 ELSE
642 JX = KX
643 DO 80 J = 1,N
644- IF (X(JX).NE.ZERO) THEN
645+c IF (X(JX).NE.ZERO) THEN
646 TEMP = ALPHA*DCONJG(X(JX))
647 AP(KK) = DBLE(AP(KK)) + DBLE(TEMP*X(JX))
648 IX = JX
649@@ -204,9 +204,9 @@
650 IX = IX + INCX
651 AP(K) = AP(K) + X(IX)*TEMP
652 70 CONTINUE
653- ELSE
654- AP(KK) = DBLE(AP(KK))
655- END IF
656+c ELSE
657+c AP(KK) = DBLE(AP(KK))
658+c END IF
659 JX = JX + INCX
660 KK = KK + N - J + 1
661 80 CONTINUE
662--- lapack-3.2.2/BLAS/SRC/zhpr2.f.orig 2009-01-02 22:57:02.000000000 +0100
663+++ lapack-3.2.2/BLAS/SRC/zhpr2.f 2010-10-19 09:51:54.631741882 +0200
664@@ -166,7 +166,7 @@
665 *
666 IF ((INCX.EQ.1) .AND. (INCY.EQ.1)) THEN
667 DO 20 J = 1,N
668- IF ((X(J).NE.ZERO) .OR. (Y(J).NE.ZERO)) THEN
669+c IF ((X(J).NE.ZERO) .OR. (Y(J).NE.ZERO)) THEN
670 TEMP1 = ALPHA*DCONJG(Y(J))
671 TEMP2 = DCONJG(ALPHA*X(J))
672 K = KK
673@@ -176,14 +176,14 @@
674 10 CONTINUE
675 AP(KK+J-1) = DBLE(AP(KK+J-1)) +
676 + DBLE(X(J)*TEMP1+Y(J)*TEMP2)
677- ELSE
678- AP(KK+J-1) = DBLE(AP(KK+J-1))
679- END IF
680+c ELSE
681+c AP(KK+J-1) = DBLE(AP(KK+J-1))
682+c END IF
683 KK = KK + J
684 20 CONTINUE
685 ELSE
686 DO 40 J = 1,N
687- IF ((X(JX).NE.ZERO) .OR. (Y(JY).NE.ZERO)) THEN
688+c IF ((X(JX).NE.ZERO) .OR. (Y(JY).NE.ZERO)) THEN
689 TEMP1 = ALPHA*DCONJG(Y(JY))
690 TEMP2 = DCONJG(ALPHA*X(JX))
691 IX = KX
692@@ -195,9 +195,9 @@
693 30 CONTINUE
694 AP(KK+J-1) = DBLE(AP(KK+J-1)) +
695 + DBLE(X(JX)*TEMP1+Y(JY)*TEMP2)
696- ELSE
697- AP(KK+J-1) = DBLE(AP(KK+J-1))
698- END IF
699+c ELSE
700+c AP(KK+J-1) = DBLE(AP(KK+J-1))
701+c END IF
702 JX = JX + INCX
703 JY = JY + INCY
704 KK = KK + J
705@@ -209,7 +209,7 @@
706 *
707 IF ((INCX.EQ.1) .AND. (INCY.EQ.1)) THEN
708 DO 60 J = 1,N
709- IF ((X(J).NE.ZERO) .OR. (Y(J).NE.ZERO)) THEN
710+c IF ((X(J).NE.ZERO) .OR. (Y(J).NE.ZERO)) THEN
711 TEMP1 = ALPHA*DCONJG(Y(J))
712 TEMP2 = DCONJG(ALPHA*X(J))
713 AP(KK) = DBLE(AP(KK)) +
714@@ -219,14 +219,14 @@
715 AP(K) = AP(K) + X(I)*TEMP1 + Y(I)*TEMP2
716 K = K + 1
717 50 CONTINUE
718- ELSE
719- AP(KK) = DBLE(AP(KK))
720- END IF
721+c ELSE
722+c AP(KK) = DBLE(AP(KK))
723+c END IF
724 KK = KK + N - J + 1
725 60 CONTINUE
726 ELSE
727 DO 80 J = 1,N
728- IF ((X(JX).NE.ZERO) .OR. (Y(JY).NE.ZERO)) THEN
729+c IF ((X(JX).NE.ZERO) .OR. (Y(JY).NE.ZERO)) THEN
730 TEMP1 = ALPHA*DCONJG(Y(JY))
731 TEMP2 = DCONJG(ALPHA*X(JX))
732 AP(KK) = DBLE(AP(KK)) +
733@@ -238,9 +238,9 @@
734 IY = IY + INCY
735 AP(K) = AP(K) + X(IX)*TEMP1 + Y(IY)*TEMP2
736 70 CONTINUE
737- ELSE
738- AP(KK) = DBLE(AP(KK))
739- END IF
740+c ELSE
741+c AP(KK) = DBLE(AP(KK))
742+c END IF
743 JX = JX + INCX
744 JY = JY + INCY
745 KK = KK + N - J + 1
746--- lapack-3.2.2/BLAS/SRC/zsyr2k.f.orig 2009-01-02 22:57:02.000000000 +0100
747+++ lapack-3.2.2/BLAS/SRC/zsyr2k.f 2010-10-19 21:48:39.731741882 +0200
748@@ -245,14 +245,14 @@
749 100 CONTINUE
750 END IF
751 DO 120 L = 1,K
752- IF ((A(J,L).NE.ZERO) .OR. (B(J,L).NE.ZERO)) THEN
753+c IF ((A(J,L).NE.ZERO) .OR. (B(J,L).NE.ZERO)) THEN
754 TEMP1 = ALPHA*B(J,L)
755 TEMP2 = ALPHA*A(J,L)
756 DO 110 I = 1,J
757 C(I,J) = C(I,J) + A(I,L)*TEMP1 +
758 + B(I,L)*TEMP2
759 110 CONTINUE
760- END IF
761+c END IF
762 120 CONTINUE
763 130 CONTINUE
764 ELSE
765@@ -267,14 +267,14 @@
766 150 CONTINUE
767 END IF
768 DO 170 L = 1,K
769- IF ((A(J,L).NE.ZERO) .OR. (B(J,L).NE.ZERO)) THEN
770+c IF ((A(J,L).NE.ZERO) .OR. (B(J,L).NE.ZERO)) THEN
771 TEMP1 = ALPHA*B(J,L)
772 TEMP2 = ALPHA*A(J,L)
773 DO 160 I = J,N
774 C(I,J) = C(I,J) + A(I,L)*TEMP1 +
775 + B(I,L)*TEMP2
776 160 CONTINUE
777- END IF
778+c END IF
779 170 CONTINUE
780 180 CONTINUE
781 END IF
782--- lapack-3.2.2/BLAS/SRC/zsyrk.f.orig 2009-01-02 22:57:02.000000000 +0100
783+++ lapack-3.2.2/BLAS/SRC/zsyrk.f 2010-10-19 21:49:12.923740205 +0200
784@@ -226,12 +226,12 @@
785 100 CONTINUE
786 END IF
787 DO 120 L = 1,K
788- IF (A(J,L).NE.ZERO) THEN
789+c IF (A(J,L).NE.ZERO) THEN
790 TEMP = ALPHA*A(J,L)
791 DO 110 I = 1,J
792 C(I,J) = C(I,J) + TEMP*A(I,L)
793 110 CONTINUE
794- END IF
795+c END IF
796 120 CONTINUE
797 130 CONTINUE
798 ELSE
799@@ -246,12 +246,12 @@
800 150 CONTINUE
801 END IF
802 DO 170 L = 1,K
803- IF (A(J,L).NE.ZERO) THEN
804+c IF (A(J,L).NE.ZERO) THEN
805 TEMP = ALPHA*A(J,L)
806 DO 160 I = J,N
807 C(I,J) = C(I,J) + TEMP*A(I,L)
808 160 CONTINUE
809- END IF
810+c END IF
811 170 CONTINUE
812 180 CONTINUE
813 END IF
814--- lapack-3.2.2/BLAS/SRC/ztbmv.f.orig 2009-01-02 22:57:02.000000000 +0100
815+++ lapack-3.2.2/BLAS/SRC/ztbmv.f 2010-10-19 21:50:19.351741882 +0200
816@@ -209,19 +209,19 @@
817 KPLUS1 = K + 1
818 IF (INCX.EQ.1) THEN
819 DO 20 J = 1,N
820- IF (X(J).NE.ZERO) THEN
821+c IF (X(J).NE.ZERO) THEN
822 TEMP = X(J)
823 L = KPLUS1 - J
824 DO 10 I = MAX(1,J-K),J - 1
825 X(I) = X(I) + TEMP*A(L+I,J)
826 10 CONTINUE
827 IF (NOUNIT) X(J) = X(J)*A(KPLUS1,J)
828- END IF
829+c END IF
830 20 CONTINUE
831 ELSE
832 JX = KX
833 DO 40 J = 1,N
834- IF (X(JX).NE.ZERO) THEN
835+c IF (X(JX).NE.ZERO) THEN
836 TEMP = X(JX)
837 IX = KX
838 L = KPLUS1 - J
839@@ -230,7 +230,7 @@
840 IX = IX + INCX
841 30 CONTINUE
842 IF (NOUNIT) X(JX) = X(JX)*A(KPLUS1,J)
843- END IF
844+c END IF
845 JX = JX + INCX
846 IF (J.GT.K) KX = KX + INCX
847 40 CONTINUE
848@@ -238,20 +238,20 @@
849 ELSE
850 IF (INCX.EQ.1) THEN
851 DO 60 J = N,1,-1
852- IF (X(J).NE.ZERO) THEN
853+c IF (X(J).NE.ZERO) THEN
854 TEMP = X(J)
855 L = 1 - J
856 DO 50 I = MIN(N,J+K),J + 1,-1
857 X(I) = X(I) + TEMP*A(L+I,J)
858 50 CONTINUE
859 IF (NOUNIT) X(J) = X(J)*A(1,J)
860- END IF
861+c END IF
862 60 CONTINUE
863 ELSE
864 KX = KX + (N-1)*INCX
865 JX = KX
866 DO 80 J = N,1,-1
867- IF (X(JX).NE.ZERO) THEN
868+c IF (X(JX).NE.ZERO) THEN
869 TEMP = X(JX)
870 IX = KX
871 L = 1 - J
872@@ -260,7 +260,7 @@
873 IX = IX - INCX
874 70 CONTINUE
875 IF (NOUNIT) X(JX) = X(JX)*A(1,J)
876- END IF
877+c END IF
878 JX = JX - INCX
879 IF ((N-J).GE.K) KX = KX - INCX
880 80 CONTINUE
881--- lapack-3.2.2/BLAS/SRC/ztbsv.f.orig 2009-01-02 22:57:02.000000000 +0100
882+++ lapack-3.2.2/BLAS/SRC/ztbsv.f 2010-10-19 21:51:36.671745234 +0200
883@@ -213,21 +213,21 @@
884 KPLUS1 = K + 1
885 IF (INCX.EQ.1) THEN
886 DO 20 J = N,1,-1
887- IF (X(J).NE.ZERO) THEN
888+c IF (X(J).NE.ZERO) THEN
889 L = KPLUS1 - J
890 IF (NOUNIT) X(J) = X(J)/A(KPLUS1,J)
891 TEMP = X(J)
892 DO 10 I = J - 1,MAX(1,J-K),-1
893 X(I) = X(I) - TEMP*A(L+I,J)
894 10 CONTINUE
895- END IF
896+c END IF
897 20 CONTINUE
898 ELSE
899 KX = KX + (N-1)*INCX
900 JX = KX
901 DO 40 J = N,1,-1
902 KX = KX - INCX
903- IF (X(JX).NE.ZERO) THEN
904+c IF (X(JX).NE.ZERO) THEN
905 IX = KX
906 L = KPLUS1 - J
907 IF (NOUNIT) X(JX) = X(JX)/A(KPLUS1,J)
908@@ -236,27 +236,27 @@
909 X(IX) = X(IX) - TEMP*A(L+I,J)
910 IX = IX - INCX
911 30 CONTINUE
912- END IF
913+c END IF
914 JX = JX - INCX
915 40 CONTINUE
916 END IF
917 ELSE
918 IF (INCX.EQ.1) THEN
919 DO 60 J = 1,N
920- IF (X(J).NE.ZERO) THEN
921+c IF (X(J).NE.ZERO) THEN
922 L = 1 - J
923 IF (NOUNIT) X(J) = X(J)/A(1,J)
924 TEMP = X(J)
925 DO 50 I = J + 1,MIN(N,J+K)
926 X(I) = X(I) - TEMP*A(L+I,J)
927 50 CONTINUE
928- END IF
929+c END IF
930 60 CONTINUE
931 ELSE
932 JX = KX
933 DO 80 J = 1,N
934 KX = KX + INCX
935- IF (X(JX).NE.ZERO) THEN
936+c IF (X(JX).NE.ZERO) THEN
937 IX = KX
938 L = 1 - J
939 IF (NOUNIT) X(JX) = X(JX)/A(1,J)
940@@ -265,7 +265,7 @@
941 X(IX) = X(IX) - TEMP*A(L+I,J)
942 IX = IX + INCX
943 70 CONTINUE
944- END IF
945+c END IF
946 JX = JX + INCX
947 80 CONTINUE
948 END IF
949--- lapack-3.2.2/BLAS/SRC/ztpmv.f.orig 2009-01-02 22:57:02.000000000 +0100
950+++ lapack-3.2.2/BLAS/SRC/ztpmv.f 2010-10-19 22:55:50.851739367 +0200
951@@ -165,7 +165,7 @@
952 KK = 1
953 IF (INCX.EQ.1) THEN
954 DO 20 J = 1,N
955- IF (X(J).NE.ZERO) THEN
956+c IF (X(J).NE.ZERO) THEN
957 TEMP = X(J)
958 K = KK
959 DO 10 I = 1,J - 1
960@@ -173,13 +173,13 @@
961 K = K + 1
962 10 CONTINUE
963 IF (NOUNIT) X(J) = X(J)*AP(KK+J-1)
964- END IF
965+c END IF
966 KK = KK + J
967 20 CONTINUE
968 ELSE
969 JX = KX
970 DO 40 J = 1,N
971- IF (X(JX).NE.ZERO) THEN
972+c IF (X(JX).NE.ZERO) THEN
973 TEMP = X(JX)
974 IX = KX
975 DO 30 K = KK,KK + J - 2
976@@ -187,7 +187,7 @@
977 IX = IX + INCX
978 30 CONTINUE
979 IF (NOUNIT) X(JX) = X(JX)*AP(KK+J-1)
980- END IF
981+c END IF
982 JX = JX + INCX
983 KK = KK + J
984 40 CONTINUE
985@@ -196,7 +196,7 @@
986 KK = (N* (N+1))/2
987 IF (INCX.EQ.1) THEN
988 DO 60 J = N,1,-1
989- IF (X(J).NE.ZERO) THEN
990+c IF (X(J).NE.ZERO) THEN
991 TEMP = X(J)
992 K = KK
993 DO 50 I = N,J + 1,-1
994@@ -204,14 +204,14 @@
995 K = K - 1
996 50 CONTINUE
997 IF (NOUNIT) X(J) = X(J)*AP(KK-N+J)
998- END IF
999+c END IF
1000 KK = KK - (N-J+1)
1001 60 CONTINUE
1002 ELSE
1003 KX = KX + (N-1)*INCX
1004 JX = KX
1005 DO 80 J = N,1,-1
1006- IF (X(JX).NE.ZERO) THEN
1007+c IF (X(JX).NE.ZERO) THEN
1008 TEMP = X(JX)
1009 IX = KX
1010 DO 70 K = KK,KK - (N- (J+1)),-1
1011@@ -219,7 +219,7 @@
1012 IX = IX - INCX
1013 70 CONTINUE
1014 IF (NOUNIT) X(JX) = X(JX)*AP(KK-N+J)
1015- END IF
1016+c END IF
1017 JX = JX - INCX
1018 KK = KK - (N-J+1)
1019 80 CONTINUE
1020--- lapack-3.2.2/BLAS/SRC/ztpsv.f.orig 2009-01-02 22:57:02.000000000 +0100
1021+++ lapack-3.2.2/BLAS/SRC/ztpsv.f 2010-10-19 22:57:19.063744396 +0200
1022@@ -168,7 +168,7 @@
1023 KK = (N* (N+1))/2
1024 IF (INCX.EQ.1) THEN
1025 DO 20 J = N,1,-1
1026- IF (X(J).NE.ZERO) THEN
1027+c IF (X(J).NE.ZERO) THEN
1028 IF (NOUNIT) X(J) = X(J)/AP(KK)
1029 TEMP = X(J)
1030 K = KK - 1
1031@@ -176,13 +176,13 @@
1032 X(I) = X(I) - TEMP*AP(K)
1033 K = K - 1
1034 10 CONTINUE
1035- END IF
1036+c END IF
1037 KK = KK - J
1038 20 CONTINUE
1039 ELSE
1040 JX = KX + (N-1)*INCX
1041 DO 40 J = N,1,-1
1042- IF (X(JX).NE.ZERO) THEN
1043+c IF (X(JX).NE.ZERO) THEN
1044 IF (NOUNIT) X(JX) = X(JX)/AP(KK)
1045 TEMP = X(JX)
1046 IX = JX
1047@@ -190,7 +190,7 @@
1048 IX = IX - INCX
1049 X(IX) = X(IX) - TEMP*AP(K)
1050 30 CONTINUE
1051- END IF
1052+c END IF
1053 JX = JX - INCX
1054 KK = KK - J
1055 40 CONTINUE
1056@@ -199,7 +199,7 @@
1057 KK = 1
1058 IF (INCX.EQ.1) THEN
1059 DO 60 J = 1,N
1060- IF (X(J).NE.ZERO) THEN
1061+c IF (X(J).NE.ZERO) THEN
1062 IF (NOUNIT) X(J) = X(J)/AP(KK)
1063 TEMP = X(J)
1064 K = KK + 1
1065@@ -207,13 +207,13 @@
1066 X(I) = X(I) - TEMP*AP(K)
1067 K = K + 1
1068 50 CONTINUE
1069- END IF
1070+c END IF
1071 KK = KK + (N-J+1)
1072 60 CONTINUE
1073 ELSE
1074 JX = KX
1075 DO 80 J = 1,N
1076- IF (X(JX).NE.ZERO) THEN
1077+c IF (X(JX).NE.ZERO) THEN
1078 IF (NOUNIT) X(JX) = X(JX)/AP(KK)
1079 TEMP = X(JX)
1080 IX = JX
1081@@ -221,7 +221,7 @@
1082 IX = IX + INCX
1083 X(IX) = X(IX) - TEMP*AP(K)
1084 70 CONTINUE
1085- END IF
1086+c END IF
1087 JX = JX + INCX
1088 KK = KK + (N-J+1)
1089 80 CONTINUE
The LAPACK libraries for numerical linear algebra
This page took 0.291984 seconds and 4 git commands to generate.