@@ -161,15 +161,21 @@ SUBROUTINE ZLAHILB( N, NRHS, A, LDA, X, LDX, B, LDB, WORK,
161161      INTEGER  NMAX_EXACT, NMAX_APPROX, SIZE_D
162162      PARAMETER  (NMAX_EXACT =  6 , NMAX_APPROX =  11 , SIZE_D =  8 )
163163* 
164- *      d's are generated from random permutation of those eight elements.
165-       COMPLEX * 16  d1(8 ), d2(8 ), invd1(8 ), invd2(8 )
166-       DATA  D1 / (- 1 ,0 ),(0 ,1 ),(- 1 ,- 1 ),(0 ,- 1 ),(1 ,0 ),(- 1 ,1 ),(1 ,1 ),(1 ,- 1 )/ 
167-       DATA  D2 / (- 1 ,0 ),(0 ,- 1 ),(- 1 ,1 ),(0 ,1 ),(1 ,0 ),(- 1 ,- 1 ),(1 ,- 1 ),(1 ,1 )/ 
164+ *      D's are generated from random permutation of those eight elements.
165+       COMPLEX * 16  D1(8 ), D2(8 ), INVD1(8 ), INVD2(8 )
166+       DATA  D1 / (- 1.0D0 ,0.0D0 ),(0.0D0 ,1.0D0 ),(- 1.0D0 ,- 1.0D0 ),
167+      $     (0.0D0 ,- 1.0D0 ),(1.0D0 ,0.0D0 ),(- 1.0D0 ,1.0D0 ),(1.0D0 ,1.0D0 ),
168+      $     (1.0D0 ,- 1.0D0 )/ 
169+       DATA  D2 / (- 1.0D0 ,0.0D0 ),(0.0D0 ,- 1.0D0 ),(- 1.0D0 ,1.0D0 ),
170+      $     (0.0D0 ,1.0D0 ),(1.0D0 ,0.0D0 ),(- 1.0D0 ,- 1.0D0 ),(1.0D0 ,- 1.0D0 ),
171+      $     (1.0D0 ,1.0D0 )/ 
168172
169-       DATA  INVD1 / (- 1 ,0 ),(0 ,- 1 ),(- .5 ,.5 ),(0 ,1 ),(1 ,0 ),
170-      $     (- .5 ,- .5 ),(.5 ,- .5 ),(.5 ,.5 )/ 
171-       DATA  INVD2 / (- 1 ,0 ),(0 ,1 ),(- .5 ,- .5 ),(0 ,- 1 ),(1 ,0 ),
172-      $     (- .5 ,.5 ),(.5 ,.5 ),(.5 ,- .5 )/ 
173+       DATA  INVD1 / (- 1.0D0 ,0.0D0 ),(0.0D0 ,- 1.0D0 ),(- 0.5D0 ,0.5D0 ),
174+      $     (0.0D0 ,1.0D0 ),(1.0D0 ,0.0D0 ),(- 0.5D0 ,- 0.5D0 ),(0.5D0 ,- 0.5D0 ),
175+      $     (0.5D0 ,0.5D0 )/ 
176+       DATA  INVD2 / (- 1.0D0 ,0.0D0 ),(0.0D0 ,1.0D0 ),(- 0.5D0 ,- 0.5D0 ),
177+      $     (0.0D0 ,- 1.0D0 ),(1.0D0 ,0.0D0 ),(- 0.5D0 ,0.5D0 ),(0.5D0 ,0.5D0 ),
178+      $     (0.5D0 ,- 0.5D0 )/ 
173179*      ..
174180*      .. External Subroutines ..
175181      EXTERNAL  XERBLA
0 commit comments