SourceForge Ticket System
norep****@sourc*****
2012年 11月 29日 (木) 03:36:05 JST
#30208: Estimate of "Block size" is needed each time that DSYSV is called. Open Date: 2012-11-28 13:47 Last Update: 2012-11-29 03:36 URL for this Ticket: http://sourceforge.jp/ticket/browse.php?group_id=6231&tid=30208 RSS feed for this Ticket: http://sourceforge.jp/ticket/ticket_rss.php?group_id=6231&tid=30208 --------------------------------------------------------------------- Last Changes/Comment for this Ticket: 2012-11-29 03:36 Update by: mikiya_fujii * Ticket Close date set to 2012-11-29 03:36 * 状況 Update from オープン to 完了. * 解決法 Update from なし to 修正済み. Comment: fixed at r1146. --------------------------------------------------------------------- Ticket Status: Reporter: mikiya_fujii Owner: mikiya_fujii Type: バグ Status: 完了 Priority: 5 - 中 MileStone: ver0.2 Component: (未割り当て) Severity: 4 Resolution: 修正済み --------------------------------------------------------------------- Ticket details: "Block size" is needed in each time that DSYSV is called.[[BR]] In Lapack. cc, following sections should be called in each time. {{{ 218 #pragma omp critical 219 { 220 if(!this->calculatedDsysvBlockSize){ 221 lwork = -1; 222 double tempWork[3]={0.0, 0.0, 0.0}; 223 dsysv(&uplo, &size, &nrhs, convertedMatrix, &lda, ipiv, tempB, &ldb, tempWork, &lwork, &info); 224 this->calculatedDsysvBlockSize = true; 225 this->dsysvBlockSize = tempWork[0]/size; 226 } 227 } }}} In, Lapack_GNU.cpp following sections should be called in each time. {{{ 218 #pragma omp critical 219 { 220 if(!this->calculatedDsysvBlockSize){ 221 lwork = -1; 222 double tempWork[3]={0.0, 0.0, 0.0}; 223 info = LAPACKE_dsysv_work(LAPACK_COL_MAJOR, uplo, size, nrhs, convertedMatrix, lda, ipiv, tempB, ldb, tempWork, lwork); 224 this->calculatedDsysvBlockSize = true; 225 this->dsysvBlockSize = tempWork[0]/size; 226 } 227 } }}} -- This is Ticket on Project MolDS. MolDS Project hosted on SourceForge.JP. Project URL: http://sourceforge.jp/projects/molds SourceForge.JP: http://sourceforge.jp URL for this Ticket: http://sourceforge.jp/ticket/browse.php?group_id=6231&tid=30208 RSS feed for this Ticket: http://sourceforge.jp/ticket/ticket_rss.php?group_id=6231&tid=30208