/*
 *  MATRIX OUTPUT MODULE
 *
 *  Author:                     Advisor:
 *      Kenneth S. Kundert          Alberto Sangiovanni-Vincentelli
 *      UC Berkeley
 *
 *  This file contains the output-to-file and output-to-screen routines for
 *  the matrix package.
 *
 *  >>> User accessible functions contained in this file:
 *  spPrint
 *  spFileMatrix
 *  spFileVector
 *  spFileStats
 *
 *  >>> Other functions contained in this file:
 */


/*
 *  Revision and copyright information.
 *
 *  Copyright (c) 1985,86,87,88
 *  by Kenneth S. Kundert and the University of California.
 *
 *  Permission to use, copy, modify, and distribute this software and
 *  its documentation for any purpose and without fee is hereby granted,
 *  provided that the copyright notices appear in all copies and
 *  supporting documentation and that the authors and the University of
 *  California are properly credited.  The authors and the University of
 *  California make no representations as to the suitability of this
 *  software for any purpose.  It is provided `as is', without express
 *  or implied warranty.
 */

#include <stdio.h>
#include "spmalloc.h"
#include "localization.h"
#include "charEncoding.h"
/*
 *  IMPORTS
 *
 *  >>> Import descriptions:
 *  spConfig.h
 *     Macros that customize the sparse matrix routines.
 *  spmatrix.h
 *     Macros and declarations to be imported by the user.
 *  spDefs.h
 *     Matrix type and macro definitions for the sparse matrix routines.
 */

#define spINSIDE_SPARSE
#include "spConfig.h"
#include "spmatrix.h"
#include "spDefs.h"





#if DOCUMENTATION

/*
 *  PRINT MATRIX
 *
 *  Formats and send the matrix to standard output.  Some elementary
 *  statistics are also output.  The matrix is output in a format that is
 *  readable by people.
 *
 *  >>> Arguments:
 *  Matrix  <input>  (char *)
 *      Pointer to matrix.
 *  PrintReordered  <input>  (int)
 *      Indicates whether the matrix should be printed out in its original
 *      form, as input by the user, or whether it should be printed in its
 *      reordered form, as used by the matrix routines.  A zero indicates that
 *      the matrix should be printed as inputed, a one indicates that it
 *      should be printed reordered.
 *  Data  <input>  (int)
 *      Boolean flag that when false indicates that output should be
 *      compressed such that only the existence of an element should be
 *      indicated rather than giving the actual value.  Thus 11 times as
 *      many can be printed on a row.  A zero signifies that the matrix
 *      should be printed compressed. A one indicates that the matrix
 *      should be printed in all its glory.
 *  Header  <input>  (int)
 *      Flag indicating that extra information should be given, such as row
 *      and column numbers.
 *
 *  >>> Local variables:
 *  Col  (int)
 *      Column being printed.
 *  ElementCount  (int)
 *      Variable used to count the number of nonzero elements in the matrix.
 *  LargestElement  (RealNumber)
 *      The magnitude of the largest element in the matrix.
 *  LargestDiag  (RealNumber)
 *      The magnitude of the largest diagonal in the matrix.
 *  Magnitude  (RealNumber)
 *      The absolute value of the matrix element being printed.
 *  PrintOrdToIntColMap  (int [])
 *      A translation array that maps the order that columns will be
 *      printed in (if not PrintReordered) to the internal column numbers.
 *  PrintOrdToIntRowMap  (int [])
 *      A translation array that maps the order that rows will be
 *      printed in (if not PrintReordered) to the internal row numbers.
 *  pElement  (ElementPtr)
 *      Pointer to the element in the matrix that is to be printed.
 *  pImagElements  (ElementPtr [ ])
 *      Array of pointers to elements in the matrix.  These pointers point
 *      to the elements whose real values have just been printed.  They are
 *      used to quickly access those same elements so their imaginary values
 *      can be printed.
 *  Row  (int)
 *      Row being printed.
 *  Size  (int)
 *      The size of the matrix.
 *  SmallestDiag  (RealNumber)
 *      The magnitude of the smallest diagonal in the matrix.
 *  SmallestElement  (RealNumber)
 *      The magnitude of the smallest element in the matrix excluding zero
 *      elements.
 *  StartCol  (int)
 *      The column number of the first column to be printed in the group of
 *      columns currently being printed.
 *  StopCol  (int)
 *      The column number of the last column to be printed in the group of
 *      columns currently being printed.
 *  Top  (int)
 *      The largest expected external row or column number.
 */

void spPrint(char *eMatrix, int PrintReordered, int Data, int Header )
{
    register  int  J = 0;
    int I = 0, Row = 0, Col = 0, Size = 0, Top_ = 0, StartCol = 1, StopCol = 0, Columns = 0, ElementCount = 0;
    double  Magnitude = 0., SmallestDiag = 0., SmallestElement = 0.;
    double  LargestElement = 0.0, LargestDiag = 0.0;
    MatrixPtr  Matrix = (MatrixPtr)eMatrix;
    ElementPtr  pElement = NULL, pImagElements[PRINTER_WIDTH / 10 + 1];
    int  *PrintOrdToIntRowMap = NULL, *PrintOrdToIntColMap = NULL;

    /* Begin `spPrint'. */
    ASSERT( IS_SPARSE( Matrix ) );
    Size = Matrix->Size;

    /* Create a packed external to internal row and column translation array. */
# if TRANSLATE
    Top_ = Matrix->AllocatedExtSize;
#else
    Top_ = Matrix->AllocatedSize;
#endif
    SPCALLOC( PrintOrdToIntRowMap, int, Top_ + 1 );
    SPCALLOC( PrintOrdToIntColMap, int, Top_ + 1 );
    if ( PrintOrdToIntRowMap == NULL OR PrintOrdToIntColMap == NULL)
    {
        Matrix->Error = spNO_MEMORY;
        return;
    }
    for (I = 1; I <= Size; I++)
    {
        PrintOrdToIntRowMap[ Matrix->IntToExtRowMap[I] ] = I;
        PrintOrdToIntColMap[ Matrix->IntToExtColMap[I] ] = I;
    }

    /* Pack the arrays. */
    for (J = 1, I = 1; I <= Top_; I++)
    {
        if (PrintOrdToIntRowMap[I] != 0)
        {
            PrintOrdToIntRowMap[ J++ ] = PrintOrdToIntRowMap[ I ];
        }
    }
    for (J = 1, I = 1; I <= Top_; I++)
    {
        if (PrintOrdToIntColMap[I] != 0)
        {
            PrintOrdToIntColMap[ J++ ] = PrintOrdToIntColMap[ I ];
        }
    }

    /* Print header. */
    if (Header)
    {
        printf(_("MATRIX SUMMARY\n\n"));
        printf(_("Size of matrix = %1u x %1u.\n"), Size, Size);
        if ( Matrix->Reordered AND PrintReordered )
        {
            printf(_("Matrix has been reordered.\n"));
        }
        putchar('\n');

        if ( Matrix->Factored )
        {
            printf(_("Matrix after factorization:\n"));
        }
        else
        {
            printf(_("Matrix before factorization:\n"));
        }

        SmallestElement = LARGEST_REAL;
        SmallestDiag = SmallestElement;
    }

    /* Determine how many columns to use. */
    Columns = PRINTER_WIDTH;
    if (Header)
    {
        Columns -= 5;
    }
    if (Data)
    {
        Columns = (Columns + 1) / 10;
    }

    /*
     * Print matrix by printing groups of complete columns until all the columns
     * are printed.
     */
    J = 0;
    while ( J <= Size )

        /* Calculate index of last column to printed in this group. */
    {
        StopCol = StartCol + Columns - 1;
        if (StopCol > Size)
        {
            StopCol = Size;
        }

        /* Label the columns. */
        if (Header)
        {
            if (Data)
            {
                printf("    ");
                for (I = StartCol; I <= StopCol; I++)
                {
                    if (PrintReordered)
                    {
                        Col = I;
                    }
                    else
                    {
                        Col = PrintOrdToIntColMap[I];
                    }
                    printf(" %9d", Matrix->IntToExtColMap[ Col ]);
                }
                printf("\n\n");
            }
            else
            {
                if (PrintReordered)
                {
                    printf(_("Columns %1d to %1d.\n"), StartCol, StopCol);
                }
                else
                {
                    printf(_("Columns %1d to %1d.\n"),
                           Matrix->IntToExtColMap[ PrintOrdToIntColMap[StartCol] ],
                           Matrix->IntToExtColMap[ PrintOrdToIntColMap[StopCol] ]);
                }
            }
        }

        /* Print every row ...  */
        for (I = 1; I <= Size; I++)
        {
            if (PrintReordered)
            {
                Row = I;
            }
            else
            {
                Row = PrintOrdToIntRowMap[I];
            }

            if (Header)
            {
                if (PrintReordered AND NOT Data)
                {
                    printf("%4d", I);
                }
                else
                {
                    printf("%4d", Matrix->IntToExtRowMap[ Row ]);
                }
                if (NOT Data)
                {
                    putchar(' ');
                }
            }

            /* ... in each column of the group. */
            for (J = StartCol; J <= StopCol; J++)
            {
                if (PrintReordered)
                {
                    Col = J;
                }
                else
                {
                    Col = PrintOrdToIntColMap[J];
                }

                pElement = Matrix->FirstInCol[Col];
                while (pElement != NULL AND pElement->Row != Row)
                {
                    pElement = pElement->NextInCol;
                }

                if (Data)
                {
                    pImagElements[J - StartCol] = pElement;
                }

                if (pElement != NULL)

                    /* Case where element exists */
                {
                    if (Data)
                    {
                        printf(" %9.3lg", (double)pElement->Real);
                    }
                    else
                    {
                        putchar('x');
                    }

                    /* Update status variables */
                    if ( (Magnitude = ELEMENT_MAG(pElement)) > LargestElement )
                    {
                        LargestElement = Magnitude;
                    }
                    if ((Magnitude < SmallestElement) AND (Magnitude != 0.0))
                    {
                        SmallestElement = Magnitude;
                    }
                    ElementCount++;
                }

                /* Case where element is structurally zero */
                else
                {
                    if (Data)
                    {
                        printf("       ...");
                    }
                    else
                    {
                        putchar('.');
                    }
                }
            }
            putchar('\n');

#if spCOMPLEX
            if (Matrix->Complex AND Data)
            {
                printf("    ");
                for (J = StartCol; J <= StopCol; J++)
                {
                    if (pImagElements[J - StartCol] != NULL)
                    {
                        printf(" %8.2lgj",
                               (double)pImagElements[J - StartCol]->Imag);
                    }
                    else
                    {
                        printf("          ");
                    }
                }
                putchar('\n');
            }
#endif /* spCOMPLEX */
        }

        /* Calculate index of first column in next group. */
        StartCol = StopCol;
        StartCol++;
        putchar('\n');
    }
    if (Header)
    {
        printf(_("\nLargest element in matrix = %-1.4lg.\n"), LargestElement);
        printf(_("Smallest element in matrix = %-1.4lg.\n"), SmallestElement);

        /* Search for largest and smallest diagonal values */
        for (I = 1; I <= Size; I++)
        {
            if (Matrix->Diag[I] != NULL)
            {
                Magnitude = ELEMENT_MAG( Matrix->Diag[I] );
                if ( Magnitude > LargestDiag )
                {
                    LargestDiag = Magnitude;
                }
                if ( Magnitude < SmallestDiag )
                {
                    SmallestDiag = Magnitude;
                }
            }
        }

        /* Print the largest and smallest diagonal values */
        if ( Matrix->Factored )
        {
            printf(_("\nLargest diagonal element = %-1.4lg.\n"), LargestDiag);
            printf(_("Smallest diagonal element = %-1.4lg.\n"), SmallestDiag);
        }
        else
        {
            printf(_("\nLargest pivot element = %-1.4lg.\n"), LargestDiag);
            printf(_("Smallest pivot element = %-1.4lg.\n"), SmallestDiag);
        }

        /* Calculate and print sparsity and number of fill-ins created. */
        printf(_("\nDensity = %2.2lf%%.\n"), ((double)(ElementCount * 100)) /
               ((double)(Size * Size)));
        if (NOT Matrix->NeedsOrdering)
        {
            printf(_("Number of fill-ins = %1d.\n"), Matrix->Fillins);
        }
    }
    putchar('\n');
    (void)fflush(stdout);

    SPFREE(PrintOrdToIntColMap);
    SPFREE(PrintOrdToIntRowMap);
    return;
}











/*
 *  OUTPUT MATRIX TO FILE
 *
 *  Writes matrix to file in format suitable to be read back in by the
 *  matrix test program.
 *
 *  >>> Returns:
 *  One is returned if routine was successful, otherwise zero is returned.
 *  The calling function can query errno (the system global error variable)
 *  as to the reason why this routine failed.
 *
 *  >>> Arguments:
 *  Matrix  <input>  (char *)
 *      Pointer to matrix.
 *  File  <input>  (char *)
 *      Name of file into which matrix is to be written.
 *  Label  <input>  (char *)
 *      String that is transferred to file and is used as a label.
 *  Reordered  <input> (SPBOOLEAN)
 *      Specifies whether matrix should be output in reordered form,
 *      or in original order.
 *  Data  <input> (SPBOOLEAN)
 *      Indicates that the element values should be output along with
 *      the indices for each element.  This parameter must be true if
 *      matrix is to be read by the sparse test program.
 *  Header  <input> (SPBOOLEAN)
 *      Indicates that header is desired.  This parameter must be true if
 *      matrix is to be read by the sparse test program.
 *
 *  >>> Local variables:
 *  Col  (int)
 *      The original column number of the element being output.
 *  pElement  (ElementPtr)
 *      Pointer to an element in the matrix.
 *  pMatrixFile  (FILE *)
 *      File pointer to the matrix file.
 *  Row  (int)
 *      The original row number of the element being output.
 *  Size  (int)
 *      The size of the matrix.
 */

int spFileMatrix(char * eMatrix, char * File, char * Label, int  Reordered, int Data, int Header )
{
    MatrixPtr  Matrix = (MatrixPtr)eMatrix;
    register  int  I = 0, Size = 0;
    register  ElementPtr  pElement = NULL;
    int  Row = 0, Col = 0, Err = 0;
    FILE  *pMatrixFile = NULL;

    /* Begin `spFileMatrix'. */
    ASSERT( IS_SPARSE( Matrix ) );

    /* Open file matrix file in write mode. */
    wcfopen(pMatrixFile, File, "w");
    if (pMatrixFile == NULL)
    {
        return 0;
    }

    /* Output header. */
    Size = Matrix->Size;
    if (Header)
    {
        if (Matrix->Factored AND Data)
        {
            Err = fprintf
                  (   pMatrixFile,
                      _("Warning : The following matrix is factored in to LU form.\n")
                  );
        }
        if (Err < 0)
        {
            return 0;
        }
        if (fprintf(pMatrixFile, "%s\n", Label) < 0)
        {
            return 0;
        }
        Err = fprintf( pMatrixFile, "%d\t%s\n", Size,
                       (Matrix->Complex ? "complex" : "real"));
        if (Err < 0)
        {
            return 0;
        }
    }

    /* Output matrix. */
    if (NOT Data)
    {
        for (I = 1; I <= Size; I++)
        {
            pElement = Matrix->FirstInCol[I];
            while (pElement != NULL)
            {
                if (Reordered)
                {
                    Row = pElement->Row;
                    Col = I;
                }
                else
                {
                    Row = Matrix->IntToExtRowMap[pElement->Row];
                    Col = Matrix->IntToExtColMap[I];
                }
                pElement = pElement->NextInCol;
                if (fprintf(pMatrixFile, "%d\t%d\n", Row, Col) < 0)
                {
                    return 0;
                }
            }
        }
        /* Output terminator, a line of zeros. */
        if (Header)
            if (fprintf(pMatrixFile, "0\t0\n") < 0)
            {
                return 0;
            }
    }

#if spCOMPLEX
    if (Data AND Matrix->Complex)
    {
        for (I = 1; I <= Size; I++)
        {
            pElement = Matrix->FirstInCol[I];
            while (pElement != NULL)
            {
                if (Reordered)
                {
                    Row = pElement->Row;
                    Col = I;
                }
                else
                {
                    Row = Matrix->IntToExtRowMap[pElement->Row];
                    Col = Matrix->IntToExtColMap[I];
                }
                Err = fprintf
                      (   pMatrixFile, "%d\t%d\t%-.15lg\t%-.15lg\n",
                          Row, Col, (double)pElement->Real, (double)pElement->Imag
                      );
                if (Err < 0)
                {
                    return 0;
                }
                pElement = pElement->NextInCol;
            }
        }
        /* Output terminator, a line of zeros. */
        if (Header)
            if (fprintf(pMatrixFile, "0\t0\t0.0\t0.0\n") < 0)
            {
                return 0;
            }

    }
#endif /* spCOMPLEX */

#if REAL
    if (Data AND NOT Matrix->Complex)
    {
        for (I = 1; I <= Size; I++)
        {
            pElement = Matrix->FirstInCol[I];
            while (pElement != NULL)
            {
                Row = Matrix->IntToExtRowMap[pElement->Row];
                Col = Matrix->IntToExtColMap[I];
                Err = fprintf
                      (   pMatrixFile, "%d\t%d\t%-.15lg\n",
                          Row, Col, (double)pElement->Real
                      );
                if (Err < 0)
                {
                    return 0;
                }
                pElement = pElement->NextInCol;
            }
        }
        /* Output terminator, a line of zeros. */
        if (Header)
            if (fprintf(pMatrixFile, "0\t0\t0.0\n") < 0)
            {
                return 0;
            }

    }
#endif /* REAL */

    /* Close file. */
    if (fclose(pMatrixFile) < 0)
    {
        return 0;
    }
    return 1;
}







/*
 *  OUTPUT SOURCE VECTOR TO FILE
 *
 *  Writes vector to file in format suitable to be read back in by the
 *  matrix test program.  This routine should be executed after the function
 *  spFileMatrix.
 *
 *  >>> Returns:
 *  One is returned if routine was successful, otherwise zero is returned.
 *  The calling function can query errno (the system global error variable)
 *  as to the reason why this routine failed.
 *
 *  >>> Arguments:
 *  Matrix  <input>  (char *)
 *      Pointer to matrix.
 *  File  <input>  (char *)
 *      Name of file into which matrix is to be written.
 *  RHS  <input>  (RealNumber [])
 *      Right-hand side vector. This is only the real portion if
 *      spSEPARATED_COMPLEX_VECTORS is true.
 *  iRHS  <input>  (RealNumber [])
 *      Right-hand side vector, imaginary portion.  Not necessary if matrix
 *      is real or if spSEPARATED_COMPLEX_VECTORS is set false.
 *
 *  >>> Local variables:
 *  pMatrixFile  (FILE *)
 *      File pointer to the matrix file.
 *  Size  (int)
 *      The size of the matrix.
 *
 *  >>> Obscure Macros
 *  IMAG_RHS
 *      Replaces itself with `, iRHS' if the options spCOMPLEX and
 *      spSEPARATED_COMPLEX_VECTORS are set, otherwise it disappears
 *      without a trace.
 */

int spFileVector(char *eMatrix, char *File, RealVector RHS IMAG_RHS )
{
    MatrixPtr  Matrix = (MatrixPtr)eMatrix;
    register  int  I, Size, Err;
    FILE  *pMatrixFile;

    /* Begin `spFileVector'. */
    ASSERT( IS_SPARSE( Matrix ) AND RHS != NULL)

    /* Open File in append mode. */
    wcfopen(pMatrixFile, File, "a");
    if (pMatrixFile == NULL)
    {
        return 0;
    }

    /* Correct array pointers for ARRAY_OFFSET. */
#if NOT ARRAY_OFFSET
#if spCOMPLEX
    if (Matrix->Complex)
    {
#if spSEPARATED_COMPLEX_VECTORS
        ASSERT(iRHS != NULL)
        --RHS;
        --iRHS;
#else
        RHS -= 2;
#endif
    }
    else
#endif /* spCOMPLEX */
        --RHS;
#endif /* NOT ARRAY_OFFSET */


    /* Output vector. */
    Size = Matrix->Size;
#if spCOMPLEX
    if (Matrix->Complex)
    {
#if spSEPARATED_COMPLEX_VECTORS
        for (I = 1; I <= Size; I++)
        {
            Err = fprintf
                  (   pMatrixFile, "%-.15lg\t%-.15lg\n",
                      (double)RHS[I], (double)iRHS[I]
                  );
            if (Err < 0)
            {
                return 0;
            }
        }
#else
        for (I = 1; I <= Size; I++)
        {
            Err = fprintf
                  (   pMatrixFile, "%-.15lg\t%-.15lg\n",
                      (double)RHS[2 * I], (double)RHS[2 * I + 1]
                  );
            if (Err < 0)
            {
                return 0;
            }
        }
#endif
    }
#endif /* spCOMPLEX */
#if REAL AND spCOMPLEX
    else
#endif
#if REAL
    {
        for (I = 1; I <= Size; I++)
        {
            if (fprintf(pMatrixFile, "%-.15lg\n", (double)RHS[I]) < 0)
            {
                return 0;
            }
        }
    }
#endif /* REAL */

    /* Close file. */
    if (fclose(pMatrixFile) < 0)
    {
        return 0;
    }
    return 1;
}









/*
 *  OUTPUT STATISTICS TO FILE
 *
 *  Writes useful information concerning the matrix to a file.  Should be
 *  executed after the matrix is factored.
 *
 *  >>> Returns:
 *  One is returned if routine was successful, otherwise zero is returned.
 *  The calling function can query errno (the system global error variable)
 *  as to the reason why this routine failed.
 *
 *  >>> Arguments:
 *  Matrix  <input>  (char *)
 *      Pointer to matrix.
 *  File  <input>  (char *)
 *      Name of file into which matrix is to be written.
 *  Label  <input>  (char *)
 *      String that is transferred to file and is used as a label.
 *
 *  >>> Local variables:
 *  Data  (RealNumber)
 *      The value of the matrix element being output.
 *  LargestElement  (RealNumber)
 *      The largest element in the matrix.
 *  NumberOfElements  (int)
 *      Number of nonzero elements in the matrix.
 *  pElement  (ElementPtr)
 *      Pointer to an element in the matrix.
 *  pStatsFile  (FILE *)
 *      File pointer to the statistics file.
 *  Size  (int)
 *      The size of the matrix.
 *  SmallestElement  (RealNumber)
 *      The smallest element in the matrix excluding zero elements.
 */

int spFileStats(char *eMatrix, char *File, char *Label )
{
    MatrixPtr  Matrix = (MatrixPtr)eMatrix;
    register  int  Size, I;
    register  ElementPtr  pElement;
    int NumberOfElements;
    RealNumber  Data, LargestElement, SmallestElement;
    FILE  *pStatsFile;

    /* Begin `spFileStats'. */
    ASSERT( IS_SPARSE( Matrix ) );

    /* Open File in append mode. */
    wcfopen(pStatsFile, File, "a");
    if (pStatsFile == NULL)
    {
        return 0;
    }

    /* Output statistics. */
    Size = Matrix->Size;
    if (NOT Matrix->Factored)
    {
        fprintf(pStatsFile, _("Matrix has not been factored.\n"));
    }
    fprintf(pStatsFile, _("|||  Starting new matrix  |||\n"));
    fprintf(pStatsFile, "%s\n", Label);
    if (Matrix->Complex)
    {
        fprintf(pStatsFile, _("Matrix is complex.\n"));
    }
    else
    {
        fprintf(pStatsFile, _("Matrix is real.\n"));
    }
    fprintf(pStatsFile, "     Size = %d\n", Size);

    /* Search matrix. */
    NumberOfElements = 0;
    LargestElement = 0.0;
    SmallestElement = LARGEST_REAL;

    for (I = 1; I <= Size; I++)
    {
        pElement = Matrix->FirstInCol[I];
        while (pElement != NULL)
        {
            NumberOfElements++;
            Data = ELEMENT_MAG(pElement);
            if (Data > LargestElement)
            {
                LargestElement = Data;
            }
            if (Data < SmallestElement AND Data != 0.0)
            {
                SmallestElement = Data;
            }
            pElement = pElement->NextInCol;
        }
    }

    SmallestElement = Min( SmallestElement, LargestElement );

    /* Output remaining statistics. */
    fprintf(pStatsFile, _("     Initial number of elements = %d\n"),
            NumberOfElements - Matrix->Fillins);
    fprintf(pStatsFile,
            _("     Initial average number of elements per row = %lf\n"),
            (double)(NumberOfElements - Matrix->Fillins) / (double)Size);
    fprintf(pStatsFile, _("     Fill-ins = %d\n"), Matrix->Fillins);
    fprintf(pStatsFile, _("     Average number of fill-ins per row = %lf%%\n"),
            (double)Matrix->Fillins / (double)Size);
    fprintf(pStatsFile, _("     Total number of elements = %d\n"),
            NumberOfElements);
    fprintf(pStatsFile, _("     Average number of elements per row = %lf\n"),
            (double)NumberOfElements / (double)Size);
    fprintf(pStatsFile, _("     Density = %lf%%\n"),
            (double)(100.0 * NumberOfElements) / (double)(Size * Size));
    fprintf(pStatsFile, _("     Relative Threshold = %e\n"), Matrix->RelThreshold);
    fprintf(pStatsFile, _("     Absolute Threshold = %e\n"), Matrix->AbsThreshold);
    fprintf(pStatsFile, _("     Largest Element = %e\n"), LargestElement);
    fprintf(pStatsFile, _("     Smallest Element = %e\n\n\n"), SmallestElement);

    /* Close file. */
    (void)fclose(pStatsFile);
    return 1;
}
#endif /* DOCUMENTATION */