diff --git a/SRC/pclarz.f b/SRC/pclarz.f index d79bd29d..f031385e 100644 --- a/SRC/pclarz.f +++ b/SRC/pclarz.f @@ -394,7 +394,7 @@ SUBROUTINE PCLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL CAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL CGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), RDEST, @@ -404,7 +404,7 @@ SUBROUTINE PCLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL CAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), LDC ) + $ 1, C( IOFFC1 ), LDC ) CALL CGERC( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) END IF @@ -437,7 +437,7 @@ SUBROUTINE PCLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL CAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK, MAX( 1, NQC2 ) ) + $ WORK, 1 ) * CALL CGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK, MAX( 1, NQC2 ), RDEST, @@ -447,7 +447,7 @@ SUBROUTINE PCLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL CAXPY( NQC2, -TAULOC( 1 ), WORK, - $ MAX( 1, NQC2 ), C( IOFFC1 ), + $ 1, C( IOFFC1 ), $ LDC ) CALL CGERC( MPV, NQC2, -TAULOC( 1 ), V( IOFFV ), $ 1, WORK, 1, C( IOFFC2 ), LDC ) @@ -488,7 +488,7 @@ SUBROUTINE PCLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL CAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL CGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), @@ -498,7 +498,7 @@ SUBROUTINE PCLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL CAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), + $ 1, C( IOFFC1 ), $ LDC ) CALL CGERC( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) @@ -554,7 +554,7 @@ SUBROUTINE PCLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL CAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL CGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), RDEST, @@ -564,7 +564,7 @@ SUBROUTINE PCLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL CAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), LDC ) + $ 1, C( IOFFC1 ), LDC ) CALL CGERC( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) END IF @@ -606,7 +606,7 @@ SUBROUTINE PCLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL CAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL CGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), RDEST, @@ -616,7 +616,7 @@ SUBROUTINE PCLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL CAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), LDC ) + $ 1, C( IOFFC1 ), LDC ) CALL CGERC( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) END IF diff --git a/SRC/pdlarz.f b/SRC/pdlarz.f index 9dcf6864..e89003a9 100644 --- a/SRC/pdlarz.f +++ b/SRC/pdlarz.f @@ -393,7 +393,7 @@ SUBROUTINE PDLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL DAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL DGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), RDEST, @@ -403,7 +403,7 @@ SUBROUTINE PDLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL DAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), LDC ) + $ 1, C( IOFFC1 ), LDC ) CALL DGER( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) END IF @@ -436,7 +436,7 @@ SUBROUTINE PDLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL DAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK, MAX( 1, NQC2 ) ) + $ WORK, 1 ) * CALL DGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK, MAX( 1, NQC2 ), RDEST, @@ -446,7 +446,7 @@ SUBROUTINE PDLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL DAXPY( NQC2, -TAULOC( 1 ), WORK, - $ MAX( 1, NQC2 ), C( IOFFC1 ), + $ 1, C( IOFFC1 ), $ LDC ) CALL DGER( MPV, NQC2, -TAULOC( 1 ), V( IOFFV ), $ 1, WORK, 1, C( IOFFC2 ), LDC ) @@ -487,7 +487,7 @@ SUBROUTINE PDLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL DAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL DGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), @@ -497,7 +497,7 @@ SUBROUTINE PDLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL DAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), + $ 1, C( IOFFC1 ), $ LDC ) CALL DGER( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) @@ -553,7 +553,7 @@ SUBROUTINE PDLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL DAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL DGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), RDEST, @@ -563,7 +563,7 @@ SUBROUTINE PDLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL DAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), LDC ) + $ 1, C( IOFFC1 ), LDC ) CALL DGER( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) END IF @@ -605,7 +605,7 @@ SUBROUTINE PDLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL DAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL DGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), RDEST, @@ -615,7 +615,7 @@ SUBROUTINE PDLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL DAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), LDC ) + $ 1, C( IOFFC1 ), LDC ) CALL DGER( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) END IF diff --git a/SRC/pslarz.f b/SRC/pslarz.f index 5617c803..a4d20f6a 100644 --- a/SRC/pslarz.f +++ b/SRC/pslarz.f @@ -393,7 +393,7 @@ SUBROUTINE PSLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL SAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL SGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), RDEST, @@ -403,7 +403,7 @@ SUBROUTINE PSLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL SAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), LDC ) + $ 1, C( IOFFC1 ), LDC ) CALL SGER( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) END IF @@ -436,7 +436,7 @@ SUBROUTINE PSLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL SAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK, MAX( 1, NQC2 ) ) + $ WORK, 1 ) * CALL SGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK, MAX( 1, NQC2 ), RDEST, @@ -446,7 +446,7 @@ SUBROUTINE PSLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL SAXPY( NQC2, -TAULOC( 1 ), WORK, - $ MAX( 1, NQC2 ), C( IOFFC1 ), + $ 1, C( IOFFC1 ), $ LDC ) CALL SGER( MPV, NQC2, -TAULOC( 1 ), V( IOFFV ), $ 1, WORK, 1, C( IOFFC2 ), LDC ) @@ -487,7 +487,7 @@ SUBROUTINE PSLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL SAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL SGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), @@ -497,7 +497,7 @@ SUBROUTINE PSLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL SAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), + $ 1, C( IOFFC1 ), $ LDC ) CALL SGER( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) @@ -553,7 +553,7 @@ SUBROUTINE PSLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL SAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL SGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), RDEST, @@ -563,7 +563,7 @@ SUBROUTINE PSLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL SAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), LDC ) + $ 1, C( IOFFC1 ), LDC ) CALL SGER( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) END IF @@ -605,7 +605,7 @@ SUBROUTINE PSLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL SAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL SGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), RDEST, @@ -615,7 +615,7 @@ SUBROUTINE PSLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL SAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), LDC ) + $ 1, C( IOFFC1 ), LDC ) CALL SGER( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) END IF diff --git a/SRC/pzlarz.f b/SRC/pzlarz.f index b77316d1..71818cb3 100644 --- a/SRC/pzlarz.f +++ b/SRC/pzlarz.f @@ -394,7 +394,7 @@ SUBROUTINE PZLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL ZAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL ZGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), RDEST, @@ -404,7 +404,7 @@ SUBROUTINE PZLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL ZAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), LDC ) + $ 1, C( IOFFC1 ), LDC ) CALL ZGERC( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) END IF @@ -437,7 +437,7 @@ SUBROUTINE PZLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL ZAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK, MAX( 1, NQC2 ) ) + $ WORK, 1 ) * CALL ZGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK, MAX( 1, NQC2 ), RDEST, @@ -447,7 +447,7 @@ SUBROUTINE PZLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL ZAXPY( NQC2, -TAULOC( 1 ), WORK, - $ MAX( 1, NQC2 ), C( IOFFC1 ), + $ 1, C( IOFFC1 ), $ LDC ) CALL ZGERC( MPV, NQC2, -TAULOC( 1 ), V( IOFFV ), $ 1, WORK, 1, C( IOFFC2 ), LDC ) @@ -488,7 +488,7 @@ SUBROUTINE PZLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL ZAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL ZGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), @@ -498,7 +498,7 @@ SUBROUTINE PZLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL ZAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), + $ 1, C( IOFFC1 ), $ LDC ) CALL ZGERC( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) @@ -554,7 +554,7 @@ SUBROUTINE PZLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL ZAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL ZGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), RDEST, @@ -564,7 +564,7 @@ SUBROUTINE PZLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL ZAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), LDC ) + $ 1, C( IOFFC1 ), LDC ) CALL ZGERC( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) END IF @@ -606,7 +606,7 @@ SUBROUTINE PZLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, END IF IF( MYROW.EQ.ICROW1 ) $ CALL ZAXPY( NQC2, ONE, C( IOFFC1 ), LDC, - $ WORK( IPW ), MAX( 1, NQC2 ) ) + $ WORK( IPW ), 1 ) * CALL ZGSUM2D( ICTXT, 'Columnwise', ' ', NQC2, 1, $ WORK( IPW ), MAX( 1, NQC2 ), RDEST, @@ -616,7 +616,7 @@ SUBROUTINE PZLARZ( SIDE, M, N, L, V, IV, JV, DESCV, INCV, TAU, C, * IF( MYROW.EQ.ICROW1 ) $ CALL ZAXPY( NQC2, -TAULOC( 1 ), WORK( IPW ), - $ MAX( 1, NQC2 ), C( IOFFC1 ), LDC ) + $ 1, C( IOFFC1 ), LDC ) CALL ZGERC( MPV, NQC2, -TAULOC( 1 ), WORK, 1, $ WORK( IPW ), 1, C( IOFFC2 ), LDC ) END IF