yMat (実際には xMat)の第一ブロックだけ Cholesky 分解して出力する部分を SDPARA に追加しました。。。
void IO::printSolutionyMat(BlockStruct& bs, Solutions& currentPt,
Parameter& param, FILE* fpout)
{
if (fpout != NULL) {
SDPA_INT info;
fprintf(fpout,"yMat_Cholesky = \n");
dpotrf_f77((char *)"Lower",¤tPt.xMat.SDP_block[0].nRow,currentPt.xMat.SDP_block[0].de_ele,
¤tPt.xMat.SDP_block[0].nRow,&info,strlen("Lower"));
currentPt.xMat.SDP_block[0].display(fpout,param.YPrint);
}
}
void IO::printSolutionyMat(BlockStruct& bs, Solutions& currentPt,
Parameter& param, FILE* fpout)
{
if (fpout != NULL) {
SDPA_INT info;
fprintf(fpout,"yMat_Cholesky = \n");
dpotrf_f77((char *)"Lower",¤tPt.xMat.SDP_block[0].nRow,currentPt.xMat.SDP_block[0].de_ele,
¤tPt.xMat.SDP_block[0].nRow,&info,strlen("Lower"));
currentPt.xMat.SDP_block[0].display(fpout,param.YPrint);
}
}