Merge pull request #158 from kyungminlee/master

langou · web-flow · commit e280b846b096 · 2026-04-08T11:53:22.000-06:00
fix bug in topology restoration in pzungql and pzunml2
diff --git a/SRC/pzungql.f b/SRC/pzungql.f
@@ -289,8 +289,8 @@ SUBROUTINE PZUNGQL( M, N, K, A, IA, JA, DESCA, TAU, WORK, LWORK,
 *
    10 CONTINUE
 *
-      CALL PB_TOPGET( ICTXT, 'Broadcast', 'Rowwise', ROWBTOP )
-      CALL PB_TOPGET( ICTXT, 'Broadcast', 'Columnwise', COLBTOP )
+      CALL PB_TOPSET( ICTXT, 'Broadcast', 'Rowwise', ROWBTOP )
+      CALL PB_TOPSET( ICTXT, 'Broadcast', 'Columnwise', COLBTOP )
 *
       WORK( 1 ) = DCMPLX( DBLE( LWMIN ) )
 *
diff --git a/SRC/pzunml2.f b/SRC/pzunml2.f
@@ -391,8 +391,8 @@ SUBROUTINE PZUNML2( SIDE, TRANS, M, N, K, A, IA, JA, DESCA, TAU,
 *
    10 CONTINUE
 *
-      CALL PB_TOPGET( ICTXT, 'Broadcast', 'Rowwise', ROWBTOP )
-      CALL PB_TOPGET( ICTXT, 'Broadcast', 'Columnwise', COLBTOP )
+      CALL PB_TOPSET( ICTXT, 'Broadcast', 'Rowwise', ROWBTOP )
+      CALL PB_TOPSET( ICTXT, 'Broadcast', 'Columnwise', COLBTOP )
 *
       WORK( 1 ) = DCMPLX( DBLE( LWMIN ) )
 *

Original file line number	Diff line number	Diff line change
`@@ -289,8 +289,8 @@ SUBROUTINE PZUNGQL( M, N, K, A, IA, JA, DESCA, TAU, WORK, LWORK,`
`289`	`289`	`*`
`290`	`290`	`10 CONTINUE`
`291`	`291`	`*`
`292`		`- CALL PB_TOPGET( ICTXT, 'Broadcast', 'Rowwise', ROWBTOP )`
`293`		`- CALL PB_TOPGET( ICTXT, 'Broadcast', 'Columnwise', COLBTOP )`
	`292`	`+ CALL PB_TOPSET( ICTXT, 'Broadcast', 'Rowwise', ROWBTOP )`
	`293`	`+ CALL PB_TOPSET( ICTXT, 'Broadcast', 'Columnwise', COLBTOP )`
`294`	`294`	`*`
`295`	`295`	`WORK( 1 ) = DCMPLX( DBLE( LWMIN ) )`
`296`	`296`	`*`
Original file line number	Diff line number	Diff line change
`@@ -391,8 +391,8 @@ SUBROUTINE PZUNML2( SIDE, TRANS, M, N, K, A, IA, JA, DESCA, TAU,`
`391`	`391`	`*`
`392`	`392`	`10 CONTINUE`
`393`	`393`	`*`
`394`		`- CALL PB_TOPGET( ICTXT, 'Broadcast', 'Rowwise', ROWBTOP )`
`395`		`- CALL PB_TOPGET( ICTXT, 'Broadcast', 'Columnwise', COLBTOP )`
	`394`	`+ CALL PB_TOPSET( ICTXT, 'Broadcast', 'Rowwise', ROWBTOP )`
	`395`	`+ CALL PB_TOPSET( ICTXT, 'Broadcast', 'Columnwise', COLBTOP )`
`396`	`396`	`*`
`397`	`397`	`WORK( 1 ) = DCMPLX( DBLE( LWMIN ) )`
`398`	`398`	`*`