Skip to content

Commit 9384776

Browse files
martin-frbgmartin-frbg
committed
Fix MSVC versions of the inline c/zdot function
1 parent 68ff451 commit 9384776

File tree

998 files changed

+33351
-16386
lines changed

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

998 files changed

+33351
-16386
lines changed

‎lapack-netlib/SRC/cbbcsd.c‎

Lines changed: 33 additions & 16 deletions
Original file line numberDiff line numberDiff line change
@@ -378,19 +378,24 @@ static integer smaxloc_(float *w, integer s, integer e, integer *n)
378378
if (w[i-1]>m) mi=i ,m=w[i-1];
379379
return mi-s+1;
380380
}
381+
381382
static inline void cdotc_(complex *z, integer *n_, complex *x, integer *incx_, complex *y, integer *incy_) {
382383
integer n = *n_, incx = *incx_, incy = *incy_, i;
383384
#ifdef _MSC_VER
384385
_Fcomplex zdotc = {0.0, 0.0};
385386
if (incx == 1 && incy == 1) {
386387
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
387-
zdotc._Val[0] += conjf(Cf(&x[i]))._Val[0] * Cf(&y[i])._Val[0];
388-
zdotc._Val[1] += conjf(Cf(&x[i]))._Val[1] * Cf(&y[i])._Val[1];
388+
zdotc._Val[0] += Cf(&x[i])._Val[0] * Cf(&y[i])._Val[0]
389+
+ Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1];
390+
zdotc._Val[1] -= Cf(&x[i])._Val[1] * Cf(&y[i])._Val[0]
391+
- Cf(&x[i])._Val[0] * Cf(&y[i])._Val[1];
389392
}
390393
} else {
391394
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
392-
zdotc._Val[0] += conjf(Cf(&x[i*incx]))._Val[0] * Cf(&y[i*incy])._Val[0];
393-
zdotc._Val[1] += conjf(Cf(&x[i*incx]))._Val[1] * Cf(&y[i*incy])._Val[1];
395+
zdotc._Val[0] += Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[0]
396+
+ Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1];
397+
zdotc._Val[1] -= Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1]
398+
- Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[1];
394399
}
395400
}
396401
pCf(z) = zdotc;
@@ -415,13 +420,17 @@ static inline void zdotc_(doublecomplex *z, integer *n_, doublecomplex *x, integ
415420
_Dcomplex zdotc = {0.0, 0.0};
416421
if (incx == 1 && incy == 1) {
417422
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
418-
zdotc._Val[0] += conj(Cd(&x[i]))._Val[0] * Cd(&y[i])._Val[0];
419-
zdotc._Val[1] += conj(Cd(&x[i]))._Val[1] * Cd(&y[i])._Val[1];
423+
zdotc._Val[0] += conj(Cd(&x[i]))._Val[0] * Cd(&y[i])._Val[0]
424+
+ Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
425+
zdotc._Val[1] += conj(Cd(&x[i]))._Val[1] * Cd(&y[i])._Val[1]
426+
- Cd(&x[i])._Val[0] * Cd(&y[i])._Val[1];
420427
}
421428
} else {
422429
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
423-
zdotc._Val[0] += conj(Cd(&x[i*incx]))._Val[0] * Cd(&y[i*incy])._Val[0];
424-
zdotc._Val[1] += conj(Cd(&x[i*incx]))._Val[1] * Cd(&y[i*incy])._Val[1];
430+
zdotc._Val[0] += conj(Cd(&x[i*incx]))._Val[0] * Cd(&y[i*incy])._Val[0]
431+
+ Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
432+
zdotc._Val[1] += conj(Cd(&x[i*incx]))._Val[1] * Cd(&y[i*incy])._Val[1]
433+
- Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[1];
425434
}
426435
}
427436
pCd(z) = zdotc;
@@ -446,13 +455,17 @@ static inline void cdotu_(complex *z, integer *n_, complex *x, integer *incx_, c
446455
_Fcomplex zdotc = {0.0, 0.0};
447456
if (incx == 1 && incy == 1) {
448457
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
449-
zdotc._Val[0] += Cf(&x[i])._Val[0] * Cf(&y[i])._Val[0];
450-
zdotc._Val[1] += Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1];
458+
zdotc._Val[0] += Cf(&x[i])._Val[0] * Cf(&y[i])._Val[0]
459+
- Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1];
460+
zdotc._Val[1] += Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1]
461+
+ Cf(&x[i])._Val[0] * Cf(&y[i])._Val[1];
451462
}
452463
} else {
453464
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
454-
zdotc._Val[0] += Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[0];
455-
zdotc._Val[1] += Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1];
465+
zdotc._Val[0] += Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[0]
466+
- Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1];
467+
zdotc._Val[1] += Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1]
468+
+ Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[1];
456469
}
457470
}
458471
pCf(z) = zdotc;
@@ -477,13 +490,17 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
477490
_Dcomplex zdotc = {0.0, 0.0};
478491
if (incx == 1 && incy == 1) {
479492
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
480-
zdotc._Val[0] += Cd(&x[i])._Val[0] * Cd(&y[i])._Val[0];
481-
zdotc._Val[1] += Cd(&x[i])._Val[1] * Cd(&y[i])._Val[1];
493+
zdotc._Val[0] += Cd(&x[i])._Val[0] * Cd(&y[i])._Val[0]
494+
- Cd(&x[i])._Val[1] * Cd(&y[i])._Val[1];
495+
zdotc._Val[1] += Cd(&x[i])._Val[1] * Cd(&y[i])._Val[1]
496+
+ Cd(&x[i])._Val[0] * Cd(&y[i])._Val[1];
482497
}
483498
} else {
484499
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
485-
zdotc._Val[0] += Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[0];
486-
zdotc._Val[1] += Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
500+
zdotc._Val[0] += Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[0]
501+
- Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
502+
zdotc._Val[1] += Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1]
503+
+ Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[1];
487504
}
488505
}
489506
pCd(z) = zdotc;

‎lapack-netlib/SRC/cbdsqr.c‎

Lines changed: 33 additions & 16 deletions
Original file line numberDiff line numberDiff line change
@@ -378,19 +378,24 @@ static integer smaxloc_(float *w, integer s, integer e, integer *n)
378378
if (w[i-1]>m) mi=i ,m=w[i-1];
379379
return mi-s+1;
380380
}
381+
381382
static inline void cdotc_(complex *z, integer *n_, complex *x, integer *incx_, complex *y, integer *incy_) {
382383
integer n = *n_, incx = *incx_, incy = *incy_, i;
383384
#ifdef _MSC_VER
384385
_Fcomplex zdotc = {0.0, 0.0};
385386
if (incx == 1 && incy == 1) {
386387
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
387-
zdotc._Val[0] += conjf(Cf(&x[i]))._Val[0] * Cf(&y[i])._Val[0];
388-
zdotc._Val[1] += conjf(Cf(&x[i]))._Val[1] * Cf(&y[i])._Val[1];
388+
zdotc._Val[0] += Cf(&x[i])._Val[0] * Cf(&y[i])._Val[0]
389+
+ Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1];
390+
zdotc._Val[1] -= Cf(&x[i])._Val[1] * Cf(&y[i])._Val[0]
391+
- Cf(&x[i])._Val[0] * Cf(&y[i])._Val[1];
389392
}
390393
} else {
391394
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
392-
zdotc._Val[0] += conjf(Cf(&x[i*incx]))._Val[0] * Cf(&y[i*incy])._Val[0];
393-
zdotc._Val[1] += conjf(Cf(&x[i*incx]))._Val[1] * Cf(&y[i*incy])._Val[1];
395+
zdotc._Val[0] += Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[0]
396+
+ Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1];
397+
zdotc._Val[1] -= Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1]
398+
- Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[1];
394399
}
395400
}
396401
pCf(z) = zdotc;
@@ -415,13 +420,17 @@ static inline void zdotc_(doublecomplex *z, integer *n_, doublecomplex *x, integ
415420
_Dcomplex zdotc = {0.0, 0.0};
416421
if (incx == 1 && incy == 1) {
417422
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
418-
zdotc._Val[0] += conj(Cd(&x[i]))._Val[0] * Cd(&y[i])._Val[0];
419-
zdotc._Val[1] += conj(Cd(&x[i]))._Val[1] * Cd(&y[i])._Val[1];
423+
zdotc._Val[0] += conj(Cd(&x[i]))._Val[0] * Cd(&y[i])._Val[0]
424+
+ Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
425+
zdotc._Val[1] += conj(Cd(&x[i]))._Val[1] * Cd(&y[i])._Val[1]
426+
- Cd(&x[i])._Val[0] * Cd(&y[i])._Val[1];
420427
}
421428
} else {
422429
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
423-
zdotc._Val[0] += conj(Cd(&x[i*incx]))._Val[0] * Cd(&y[i*incy])._Val[0];
424-
zdotc._Val[1] += conj(Cd(&x[i*incx]))._Val[1] * Cd(&y[i*incy])._Val[1];
430+
zdotc._Val[0] += conj(Cd(&x[i*incx]))._Val[0] * Cd(&y[i*incy])._Val[0]
431+
+ Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
432+
zdotc._Val[1] += conj(Cd(&x[i*incx]))._Val[1] * Cd(&y[i*incy])._Val[1]
433+
- Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[1];
425434
}
426435
}
427436
pCd(z) = zdotc;
@@ -446,13 +455,17 @@ static inline void cdotu_(complex *z, integer *n_, complex *x, integer *incx_, c
446455
_Fcomplex zdotc = {0.0, 0.0};
447456
if (incx == 1 && incy == 1) {
448457
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
449-
zdotc._Val[0] += Cf(&x[i])._Val[0] * Cf(&y[i])._Val[0];
450-
zdotc._Val[1] += Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1];
458+
zdotc._Val[0] += Cf(&x[i])._Val[0] * Cf(&y[i])._Val[0]
459+
- Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1];
460+
zdotc._Val[1] += Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1]
461+
+ Cf(&x[i])._Val[0] * Cf(&y[i])._Val[1];
451462
}
452463
} else {
453464
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
454-
zdotc._Val[0] += Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[0];
455-
zdotc._Val[1] += Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1];
465+
zdotc._Val[0] += Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[0]
466+
- Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1];
467+
zdotc._Val[1] += Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1]
468+
+ Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[1];
456469
}
457470
}
458471
pCf(z) = zdotc;
@@ -477,13 +490,17 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
477490
_Dcomplex zdotc = {0.0, 0.0};
478491
if (incx == 1 && incy == 1) {
479492
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
480-
zdotc._Val[0] += Cd(&x[i])._Val[0] * Cd(&y[i])._Val[0];
481-
zdotc._Val[1] += Cd(&x[i])._Val[1] * Cd(&y[i])._Val[1];
493+
zdotc._Val[0] += Cd(&x[i])._Val[0] * Cd(&y[i])._Val[0]
494+
- Cd(&x[i])._Val[1] * Cd(&y[i])._Val[1];
495+
zdotc._Val[1] += Cd(&x[i])._Val[1] * Cd(&y[i])._Val[1]
496+
+ Cd(&x[i])._Val[0] * Cd(&y[i])._Val[1];
482497
}
483498
} else {
484499
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
485-
zdotc._Val[0] += Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[0];
486-
zdotc._Val[1] += Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
500+
zdotc._Val[0] += Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[0]
501+
- Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
502+
zdotc._Val[1] += Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1]
503+
+ Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[1];
487504
}
488505
}
489506
pCd(z) = zdotc;

‎lapack-netlib/SRC/cgbbrd.c‎

Lines changed: 34 additions & 17 deletions
Original file line numberDiff line numberDiff line change
@@ -378,19 +378,24 @@ static integer smaxloc_(float *w, integer s, integer e, integer *n)
378378
if (w[i-1]>m) mi=i ,m=w[i-1];
379379
return mi-s+1;
380380
}
381+
381382
static inline void cdotc_(complex *z, integer *n_, complex *x, integer *incx_, complex *y, integer *incy_) {
382383
integer n = *n_, incx = *incx_, incy = *incy_, i;
383384
#ifdef _MSC_VER
384385
_Fcomplex zdotc = {0.0, 0.0};
385386
if (incx == 1 && incy == 1) {
386387
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
387-
zdotc._Val[0] += conjf(Cf(&x[i]))._Val[0] * Cf(&y[i])._Val[0];
388-
zdotc._Val[1] += conjf(Cf(&x[i]))._Val[1] * Cf(&y[i])._Val[1];
388+
zdotc._Val[0] += Cf(&x[i])._Val[0] * Cf(&y[i])._Val[0]
389+
+ Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1];
390+
zdotc._Val[1] -= Cf(&x[i])._Val[1] * Cf(&y[i])._Val[0]
391+
- Cf(&x[i])._Val[0] * Cf(&y[i])._Val[1];
389392
}
390393
} else {
391394
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
392-
zdotc._Val[0] += conjf(Cf(&x[i*incx]))._Val[0] * Cf(&y[i*incy])._Val[0];
393-
zdotc._Val[1] += conjf(Cf(&x[i*incx]))._Val[1] * Cf(&y[i*incy])._Val[1];
395+
zdotc._Val[0] += Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[0]
396+
+ Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1];
397+
zdotc._Val[1] -= Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1]
398+
- Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[1];
394399
}
395400
}
396401
pCf(z) = zdotc;
@@ -415,13 +420,17 @@ static inline void zdotc_(doublecomplex *z, integer *n_, doublecomplex *x, integ
415420
_Dcomplex zdotc = {0.0, 0.0};
416421
if (incx == 1 && incy == 1) {
417422
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
418-
zdotc._Val[0] += conj(Cd(&x[i]))._Val[0] * Cd(&y[i])._Val[0];
419-
zdotc._Val[1] += conj(Cd(&x[i]))._Val[1] * Cd(&y[i])._Val[1];
423+
zdotc._Val[0] += conj(Cd(&x[i]))._Val[0] * Cd(&y[i])._Val[0]
424+
+ Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
425+
zdotc._Val[1] += conj(Cd(&x[i]))._Val[1] * Cd(&y[i])._Val[1]
426+
- Cd(&x[i])._Val[0] * Cd(&y[i])._Val[1];
420427
}
421428
} else {
422429
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
423-
zdotc._Val[0] += conj(Cd(&x[i*incx]))._Val[0] * Cd(&y[i*incy])._Val[0];
424-
zdotc._Val[1] += conj(Cd(&x[i*incx]))._Val[1] * Cd(&y[i*incy])._Val[1];
430+
zdotc._Val[0] += conj(Cd(&x[i*incx]))._Val[0] * Cd(&y[i*incy])._Val[0]
431+
+ Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
432+
zdotc._Val[1] += conj(Cd(&x[i*incx]))._Val[1] * Cd(&y[i*incy])._Val[1]
433+
- Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[1];
425434
}
426435
}
427436
pCd(z) = zdotc;
@@ -446,13 +455,17 @@ static inline void cdotu_(complex *z, integer *n_, complex *x, integer *incx_, c
446455
_Fcomplex zdotc = {0.0, 0.0};
447456
if (incx == 1 && incy == 1) {
448457
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
449-
zdotc._Val[0] += Cf(&x[i])._Val[0] * Cf(&y[i])._Val[0];
450-
zdotc._Val[1] += Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1];
458+
zdotc._Val[0] += Cf(&x[i])._Val[0] * Cf(&y[i])._Val[0]
459+
- Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1];
460+
zdotc._Val[1] += Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1]
461+
+ Cf(&x[i])._Val[0] * Cf(&y[i])._Val[1];
451462
}
452463
} else {
453464
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
454-
zdotc._Val[0] += Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[0];
455-
zdotc._Val[1] += Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1];
465+
zdotc._Val[0] += Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[0]
466+
- Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1];
467+
zdotc._Val[1] += Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1]
468+
+ Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[1];
456469
}
457470
}
458471
pCf(z) = zdotc;
@@ -477,13 +490,17 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
477490
_Dcomplex zdotc = {0.0, 0.0};
478491
if (incx == 1 && incy == 1) {
479492
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
480-
zdotc._Val[0] += Cd(&x[i])._Val[0] * Cd(&y[i])._Val[0];
481-
zdotc._Val[1] += Cd(&x[i])._Val[1] * Cd(&y[i])._Val[1];
493+
zdotc._Val[0] += Cd(&x[i])._Val[0] * Cd(&y[i])._Val[0]
494+
- Cd(&x[i])._Val[1] * Cd(&y[i])._Val[1];
495+
zdotc._Val[1] += Cd(&x[i])._Val[1] * Cd(&y[i])._Val[1]
496+
+ Cd(&x[i])._Val[0] * Cd(&y[i])._Val[1];
482497
}
483498
} else {
484499
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
485-
zdotc._Val[0] += Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[0];
486-
zdotc._Val[1] += Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
500+
zdotc._Val[0] += Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[0]
501+
- Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
502+
zdotc._Val[1] += Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1]
503+
+ Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[1];
487504
}
488505
}
489506
pCd(z) = zdotc;
@@ -720,7 +737,7 @@ f"> */
720737
/* System generated locals */
721738
integer ab_dim1, ab_offset, c_dim1, c_offset, pt_dim1, pt_offset, q_dim1,
722739
q_offset, i__1, i__2, i__3, i__4, i__5, i__6, i__7;
723-
complex q__1, q__2, q__3;
740+
complex q__1={0.,0.}, q__2={0.,0.}, q__3={0.,0.};
724741

725742
/* Local variables */
726743
integer inca;

‎lapack-netlib/SRC/cgbcon.c‎

Lines changed: 34 additions & 17 deletions
Original file line numberDiff line numberDiff line change
@@ -378,19 +378,24 @@ static integer smaxloc_(float *w, integer s, integer e, integer *n)
378378
if (w[i-1]>m) mi=i ,m=w[i-1];
379379
return mi-s+1;
380380
}
381+
381382
static inline void cdotc_(complex *z, integer *n_, complex *x, integer *incx_, complex *y, integer *incy_) {
382383
integer n = *n_, incx = *incx_, incy = *incy_, i;
383384
#ifdef _MSC_VER
384385
_Fcomplex zdotc = {0.0, 0.0};
385386
if (incx == 1 && incy == 1) {
386387
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
387-
zdotc._Val[0] += conjf(Cf(&x[i]))._Val[0] * Cf(&y[i])._Val[0];
388-
zdotc._Val[1] += conjf(Cf(&x[i]))._Val[1] * Cf(&y[i])._Val[1];
388+
zdotc._Val[0] += Cf(&x[i])._Val[0] * Cf(&y[i])._Val[0]
389+
+ Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1];
390+
zdotc._Val[1] -= Cf(&x[i])._Val[1] * Cf(&y[i])._Val[0]
391+
- Cf(&x[i])._Val[0] * Cf(&y[i])._Val[1];
389392
}
390393
} else {
391394
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
392-
zdotc._Val[0] += conjf(Cf(&x[i*incx]))._Val[0] * Cf(&y[i*incy])._Val[0];
393-
zdotc._Val[1] += conjf(Cf(&x[i*incx]))._Val[1] * Cf(&y[i*incy])._Val[1];
395+
zdotc._Val[0] += Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[0]
396+
+ Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1];
397+
zdotc._Val[1] -= Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1]
398+
- Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[1];
394399
}
395400
}
396401
pCf(z) = zdotc;
@@ -415,13 +420,17 @@ static inline void zdotc_(doublecomplex *z, integer *n_, doublecomplex *x, integ
415420
_Dcomplex zdotc = {0.0, 0.0};
416421
if (incx == 1 && incy == 1) {
417422
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
418-
zdotc._Val[0] += conj(Cd(&x[i]))._Val[0] * Cd(&y[i])._Val[0];
419-
zdotc._Val[1] += conj(Cd(&x[i]))._Val[1] * Cd(&y[i])._Val[1];
423+
zdotc._Val[0] += conj(Cd(&x[i]))._Val[0] * Cd(&y[i])._Val[0]
424+
+ Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
425+
zdotc._Val[1] += conj(Cd(&x[i]))._Val[1] * Cd(&y[i])._Val[1]
426+
- Cd(&x[i])._Val[0] * Cd(&y[i])._Val[1];
420427
}
421428
} else {
422429
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
423-
zdotc._Val[0] += conj(Cd(&x[i*incx]))._Val[0] * Cd(&y[i*incy])._Val[0];
424-
zdotc._Val[1] += conj(Cd(&x[i*incx]))._Val[1] * Cd(&y[i*incy])._Val[1];
430+
zdotc._Val[0] += conj(Cd(&x[i*incx]))._Val[0] * Cd(&y[i*incy])._Val[0]
431+
+ Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
432+
zdotc._Val[1] += conj(Cd(&x[i*incx]))._Val[1] * Cd(&y[i*incy])._Val[1]
433+
- Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[1];
425434
}
426435
}
427436
pCd(z) = zdotc;
@@ -446,13 +455,17 @@ static inline void cdotu_(complex *z, integer *n_, complex *x, integer *incx_, c
446455
_Fcomplex zdotc = {0.0, 0.0};
447456
if (incx == 1 && incy == 1) {
448457
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
449-
zdotc._Val[0] += Cf(&x[i])._Val[0] * Cf(&y[i])._Val[0];
450-
zdotc._Val[1] += Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1];
458+
zdotc._Val[0] += Cf(&x[i])._Val[0] * Cf(&y[i])._Val[0]
459+
- Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1];
460+
zdotc._Val[1] += Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1]
461+
+ Cf(&x[i])._Val[0] * Cf(&y[i])._Val[1];
451462
}
452463
} else {
453464
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
454-
zdotc._Val[0] += Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[0];
455-
zdotc._Val[1] += Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1];
465+
zdotc._Val[0] += Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[0]
466+
- Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1];
467+
zdotc._Val[1] += Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1]
468+
+ Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[1];
456469
}
457470
}
458471
pCf(z) = zdotc;
@@ -477,13 +490,17 @@ static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integ
477490
_Dcomplex zdotc = {0.0, 0.0};
478491
if (incx == 1 && incy == 1) {
479492
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
480-
zdotc._Val[0] += Cd(&x[i])._Val[0] * Cd(&y[i])._Val[0];
481-
zdotc._Val[1] += Cd(&x[i])._Val[1] * Cd(&y[i])._Val[1];
493+
zdotc._Val[0] += Cd(&x[i])._Val[0] * Cd(&y[i])._Val[0]
494+
- Cd(&x[i])._Val[1] * Cd(&y[i])._Val[1];
495+
zdotc._Val[1] += Cd(&x[i])._Val[1] * Cd(&y[i])._Val[1]
496+
+ Cd(&x[i])._Val[0] * Cd(&y[i])._Val[1];
482497
}
483498
} else {
484499
for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
485-
zdotc._Val[0] += Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[0];
486-
zdotc._Val[1] += Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
500+
zdotc._Val[0] += Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[0]
501+
- Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
502+
zdotc._Val[1] += Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1]
503+
+ Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[1];
487504
}
488505
}
489506
pCd(z) = zdotc;
@@ -671,7 +688,7 @@ f"> */
671688
/* System generated locals */
672689
integer ab_dim1, ab_offset, i__1, i__2, i__3;
673690
real r__1, r__2;
674-
complex q__1, q__2;
691+
complex q__1={0.,0.}, q__2={0.,0.};
675692

676693
/* Local variables */
677694
integer kase, kase1, j;

0 commit comments

Comments
 (0)