diff --git a/src/maths/sparse/spoutput.c b/src/maths/sparse/spoutput.c index 6464e4fa0..cc5af829a 100644 --- a/src/maths/sparse/spoutput.c +++ b/src/maths/sparse/spoutput.c @@ -297,7 +297,7 @@ spPrint(MatrixPtr eMatrix, int PrintReordered, int Data, int Header) { /* Case where element exists */ if (Data) - printf(" %9.3g", (double)pElement->Real); + printf(" %9.3g", pElement->Real); else putchar('x'); @@ -328,7 +328,7 @@ spPrint(MatrixPtr eMatrix, int PrintReordered, int Data, int Header) if (pImagElements[J - StartCol] != NULL) { printf(" %8.2gj", - (double)pImagElements[J-StartCol]->Imag); + pImagElements[J-StartCol]->Imag); } else printf(" "); } @@ -520,7 +520,7 @@ spFileMatrix(MatrixPtr eMatrix, char *File, char *Label, int Reordered, } Err = fprintf ( pMatrixFile,"%d\t%d\t%-.15g\t%-.15g\n", - Row, Col, (double)pElement->Real, (double)pElement->Imag + Row, Col, pElement->Real, pElement->Imag ); if (Err < 0) return 0; pElement = pElement->NextInCol; @@ -543,7 +543,7 @@ spFileMatrix(MatrixPtr eMatrix, char *File, char *Label, int Reordered, Col = Matrix->IntToExtColMap[I]; Err = fprintf ( pMatrixFile,"%d\t%d\t%-.15g\n", - Row, Col, (double)pElement->Real + Row, Col, pElement->Real ); if (Err < 0) return 0; pElement = pElement->NextInCol; @@ -623,7 +623,7 @@ spFileVector(MatrixPtr eMatrix, char *File, RealVector RHS, RealVector iRHS) { Err = fprintf ( pMatrixFile, "%-.15g\t%-.15g\n", - (double)RHS[I], (double)iRHS[I] + RHS[I], iRHS[I] ); if (Err < 0) return 0; } @@ -632,7 +632,7 @@ spFileVector(MatrixPtr eMatrix, char *File, RealVector RHS, RealVector iRHS) { for (I = 1; I <= Size; I++) { - if (fprintf(pMatrixFile, "%-.15g\n", (double)RHS[I]) < 0) + if (fprintf(pMatrixFile, "%-.15g\n", RHS[I]) < 0) return 0; } }