CKNMatrixDebug Class Reference

This class includes functions for matrix debugging. More...

#include "KNMatrixDebug.h"

Collaboration diagram for CKNMatrixDebug:

Public Member Functions
	CKNMatrixDebug ()
	~CKNMatrixDebug ()

Static Public Member Functions
static void	ShowDoubleMatrix (double *pMatrix, int m, int n, int startm, int endm, int startn, int endn, char *pszDebugMsg)
	Show matrix in output debugging windows - Visual studio and Windows DebugView only. More...
static void	ShowDoubleArray (double *vector, char *pszDebugMsg, int due)
	Show double array. More...
static void	ShowDoubleVector (CKNMatrixOperation::CKNVector vector, char *pszDebugMsg, int due=-1)
static void	ShowComplex (CKNComplex number, char *pszDebugMsg)
	Show vector class. More...
static void	ShowDouble (double number, char *pszDebugMsg)
	Show one double data. More...
static void	ShowMsg (char *pszBuffer)
	Show message. More...
static void	SaveResult (CKNLanczosMethod::LPEIGENVALUE_RESULT lpResult, int nIndex)
	Show lanczos method result. More...
static void	SetShow (bool bShow)
static void	ShowCSR (CKNMatrixOperation::CKNCSR *pCSR, char *pszDebugMsg)
	Set show message or not. More...
static void	ShowCSR (CKNMatrixOperation::CKNCSR *pCSR, int startm, int endm, int startn, int endn, char *pszDebugMsg)
	Show CSR matrix partialy. More...
static void	ShowDenseMatrix (CKNMatrixOperation::CKNDMatrix *pMatrix, int startm, int endm, int startn, int endn, char *pszDebugMsg, bool bLineFeed=false, bool bShowImaginary=false)
	Show matrix partialy. More...
static void	DumpCSR (CKNMatrixOperation::CKNCSR *pCSR, char *pszFileName, double fAtomIDStartIndex)
	Dump CSR to text file. More...
static void	DumpCSRBinary (CKNMatrixOperation::CKNCSR *pCSR, char *pszFileName, double fAtomIDStartIndex)
	Dump CSR to binary file. More...

Static Private Attributes
static bool	m_bShowMsg = false

Detailed Description

This class includes functions for matrix debugging.

Date

10/July/2014

Author

Kyu Nam Cho(mysto.nosp@m.us@k.nosp@m.orea..nosp@m.ac.k.nosp@m.r), Hoo Ryu(elec1.nosp@m.020@.nosp@m.gmail.nosp@m..com)

Definition at line 17 of file KNMatrixDebug.h.

Constructor & Destructor Documentation

CKNMatrixDebug::CKNMatrixDebug

(

)

Definition at line 24 of file KNMatrixDebug.cpp.

{
}

CKNMatrixDebug::~CKNMatrixDebug

(

)

Definition at line 28 of file KNMatrixDebug.cpp.

{
}

Member Function Documentation

void CKNMatrixDebug::DumpCSR ( CKNMatrixOperation::CKNCSR *  pCSR, char *  pszFileName, double  fAtomIDStartIndex  )

static

Dump CSR to text file.

Definition at line 340 of file KNMatrixDebug.cpp.

References CKNMPIManager::GetCurrentRank(), CKNMatrixOperation::CKNCSR::GetRowCount(), CKNMPIManager::GetTotalNodeCount(), CKNMPIManager::IsRootRank(), CKNMatrixOperation::CKNCSR::m_vectColumn, CKNMatrixOperation::CKNCSR::m_vectRow, CKNMPIManager::ReceiveDoubleBufferSync(), and CKNMPIManager::SendDoubleBufferSync().

{
    unsigned int        i, j, nSize;
    FILE                *out;
    char                szMsg[1024];
    char                szFileName[1024];

#ifdef _WIN32
    _mkdir("result");
    sprintf(szFileName, "result\\%s", pszFileName);
#else
    mkdir("result", 0777);
    sprintf(szFileName, "result/%s", pszFileName);
#endif

#ifdef DISABLE_MPI_ROUTINE
    out = fopen(szFileName, "wt");
#else //DISABLE_MPI_ROUTINE
    MPI_Request         req;
    double              fTemp = 1;

    if( CKNMPIManager::IsRootRank())
        out = fopen(szFileName, "wt");
    else
    {
        CKNMPIManager::ReceiveDoubleBufferSync(CKNMPIManager::GetCurrentRank()-1, &fTemp, 1, &req);
        out = fopen(szFileName, "at");
    }
#endif //DISABLE_MPI_ROUTINE

    if( NULL == out )
        return;

    nSize = pCSR->GetRowCount();
    for(i = 0; i < nSize ; ++ i)
    {
        unsigned int    nSubStart = pCSR->m_vectRow[i], nSubEnd = pCSR->m_vectRow[i + 1];

        for (j = nSubStart; j < nSubEnd; ++ j)
        {
            sprintf(szMsg, "%05d\t%05d\t%40.30f\t%40.30f\n", i + (int)(fAtomIDStartIndex*10) + 1, 
                                                            pCSR->m_vectColumn[j] + 1, 
                                                            pCSR->m_vectValueBuffer[j].GetRealNumber(), 
                                                            pCSR->m_vectValueBuffer[j].GetImaginaryNumber());
            fputs(szMsg, out);
        }
    }
    fclose(out);

#ifndef DISABLE_MPI_ROUTINE
    if( CKNMPIManager::GetCurrentRank() < CKNMPIManager::GetTotalNodeCount() - 1)
        CKNMPIManager::SendDoubleBufferSync(CKNMPIManager::GetCurrentRank()+1, &fTemp, 1, &req);
#endif //DISABLE_MPI_ROUTINE
}

void CKNMatrixDebug::DumpCSRBinary ( CKNMatrixOperation::CKNCSR *  pCSR, char *  pszFileName, double  fAtomIDStartIndex  )

static

Dump CSR to binary file.

Definition at line 292 of file KNMatrixDebug.cpp.

References CKNMPIManager::GetCurrentRank(), CKNMatrixOperation::CKNCSR::GetRowCount(), CKNMPIManager::GetTotalNodeCount(), CKNMPIManager::IsRootRank(), CKNMatrixOperation::CKNCSR::m_vectColumn, CKNMatrixOperation::CKNCSR::m_vectRow, CKNMPIManager::ReceiveDoubleBufferSync(), and CKNMPIManager::SendDoubleBufferSync().

{
    unsigned int        i, j, nSize;
    FILE                *out;

#ifdef DISABLE_MPI_ROUTINE
    out = fopen(pszFileName, "wb");
#else //DISABLE_MPI_ROUTINE
    MPI_Request         req;
    double              fTemp = 1;

    if (CKNMPIManager::IsRootRank())
        out = fopen(pszFileName, "wb");
    else
    {
        CKNMPIManager::ReceiveDoubleBufferSync(CKNMPIManager::GetCurrentRank() - 1, &fTemp, 1, &req);
        out = fopen(pszFileName, "ab");
    }
#endif //DISABLE_MPI_ROUTINE

    if (NULL == out)
        return;

    nSize = pCSR->GetRowCount();
    for (i = 0; i < nSize; ++i)
    {
        unsigned int    nSubStart = pCSR->m_vectRow[i], nSubEnd = pCSR->m_vectRow[i + 1];

        for (j = nSubStart; j < nSubEnd; ++j)
        {
            double  fTemp;
            int     nTemp = i + (int)(fAtomIDStartIndex * 10);
            fwrite(&nTemp, sizeof(int), 1, out);
            fwrite(&pCSR->m_vectColumn[j], sizeof(int), 1, out);
            fTemp = pCSR->m_vectValueBuffer[j].GetRealNumber();
            fwrite(&fTemp, sizeof(double), 1, out);
            fTemp = pCSR->m_vectValueBuffer[j].GetImaginaryNumber();
            fwrite(&fTemp, sizeof(double), 1, out);
        }
    }
    fclose(out);

#ifndef DISABLE_MPI_ROUTINE
    if (CKNMPIManager::GetCurrentRank() < CKNMPIManager::GetTotalNodeCount() - 1)
        CKNMPIManager::SendDoubleBufferSync(CKNMPIManager::GetCurrentRank() + 1, &fTemp, 1, &req);
#endif //DISABLE_MPI_ROUTINE
}

Here is the call graph for this function:

void CKNMatrixDebug::SaveResult ( CKNLanczosMethod::LPEIGENVALUE_RESULT  lpResult, int  nIndex  )

static

Show lanczos method result.

Parameters

lpResult	Calculated eigenvalue result
nIndex	The index that want to save result

Definition at line 274 of file KNMatrixDebug.cpp.

References CKNLanczosMethod::EIGENVALUE_RESULT::nEigenValueCount, and CKNLanczosMethod::EIGENVALUE_RESULT::pEigenValues.

Referenced by CKNLanczosTest::TestLanczos().

{
    char            szBuffer[1024];
    FILE            *fpResult;
    unsigned int    i;

    sprintf(szBuffer, "lanczosresult_%d.txt", nIndex);
    if (NULL != (fpResult = fopen(szBuffer, "wt")))
    {
        for (i = 0; i < lpResult->nEigenValueCount; i++)
        {
            sprintf(szBuffer, "%f\n", lpResult->pEigenValues[i]);
            fputs(szBuffer, fpResult);
        }
        fclose(fpResult);
    }
}

Here is the caller graph for this function:

static void CKNMatrixDebug::SetShow ( bool  bShow )

inlinestatic

Definition at line 30 of file KNMatrixDebug.h.

References m_bShowMsg.

{ m_bShowMsg = bShow; }; 

void CKNMatrixDebug::ShowComplex ( CKNComplex  number, char *  pszDebugMsg  )

static

Show vector class.

Show complex class

Parameters

number	The complex number that was wanted to show
pszDebugMsg	Text that was show first at matrix display

Definition at line 244 of file KNMatrixDebug.cpp.

References CKNComplex::GetImaginaryNumber(), and CKNComplex::GetRealNumber().

{
    char            szDebug[1024];

    sprintf(szDebug, "       %20.20f   + %20.20f i\n", number.GetRealNumber(), number.GetImaginaryNumber());

    ShowMsg(pszDebugMsg);
    ShowMsg("\n");
    ShowMsg(szDebug);
}

Here is the call graph for this function:

void CKNMatrixDebug::ShowCSR ( CKNMatrixOperation::CKNCSR *  pCSR, char *  pszDebugMsg  )

static

Set show message or not.

Show CSR matrix partialy

Parameters

pCSR	The matrix that was wanted to show
pszDebugMsg	Text that was show first at matrix display

Definition at line 159 of file KNMatrixDebug.cpp.

References DUMY_STRING_FOR_COMPLEX_SHOWING, CKNMatrixOperation::CKNCSR::GetColumnCount(), CKNMatrixOperation::CKNCSR::GetElement(), CKNComplex::GetImaginaryNumber(), CKNComplex::GetRealNumber(), and CKNMatrixOperation::CKNCSR::GetRowCount().

Referenced by CKNLanczosTest::TestCSRBuilding().

{
    bool                        bExisted;
    unsigned int                iMax, jMax;

#ifdef _WIN32
    iMax = min(10, pCSR->GetRowCount());
    jMax = min(10, pCSR->GetColumnCount());
#else
    iMax = std::min(10, (int)pCSR->GetRowCount());
    jMax = std::min(10, (int)pCSR->GetColumnCount());
#endif

    ShowMsg(pszDebugMsg);
    ShowMsg("\n");
    ShowMsg("\n");

    for (unsigned int i = 0; i < iMax; i++)
    {
        char szRow[1024] = "";
        for (unsigned int j = 0; j < jMax; j++)
        {
            char        szCol[1024];
            CKNComplex  element = pCSR->GetElement(i, j, bExisted);

            if (bExisted)
                sprintf(szCol, "%3.3f:%3.3f |\t\t", element.GetRealNumber(), element.GetImaginaryNumber());
            else
                sprintf(szCol, "%s |\t\t", DUMY_STRING_FOR_COMPLEX_SHOWING);
            strcat(szRow, szCol);
        }
        strcat(szRow, "\n");
        ShowMsg(szRow);
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

void CKNMatrixDebug::ShowCSR ( CKNMatrixOperation::CKNCSR *  pCSR, int  startm, int  endm, int  startn, int  endn, char *  pszDebugMsg  )

static

Show CSR matrix partialy.

Parameters

pCSR	The matrix that was wanted to show
pszDebugMsg	Text that was show first at matrix display
startm	Start row index that was wanted to show
endm	End row index that was wanted to show
startn	Start column index that was wanted to show
endn	End column index that was wanted to show

Definition at line 128 of file KNMatrixDebug.cpp.

References DUMY_STRING_FOR_COMPLEX_SHOWING, CKNMatrixOperation::CKNCSR::GetElement(), CKNComplex::GetImaginaryNumber(), and CKNComplex::GetRealNumber().

{
    bool                        bExisted;
    
    ShowMsg(pszDebugMsg);
    ShowMsg("\n");
    ShowMsg("\n");

    for (unsigned int i = startm; i <= (unsigned int)endm; i++)
    {
        char szRow[1024] = "";
        for (unsigned int j = startn; j <= (unsigned int)endn; j++)
        {
            char        szCol[1024];
            CKNComplex  element = pCSR->GetElement(i, j, bExisted);

            if (bExisted)
                sprintf(szCol, "%3.3f:%3.3f |\t\t", element.GetRealNumber(), element.GetImaginaryNumber());
            else
                sprintf(szCol, "%s |\t\t", DUMY_STRING_FOR_COMPLEX_SHOWING);
            strcat(szRow, szCol);
        }
        strcat(szRow, "\n");
        ShowMsg(szRow);
    }
}

Here is the call graph for this function:

void CKNMatrixDebug::ShowDenseMatrix ( CKNMatrixOperation::CKNDMatrix *  pMatrix, int  startm, int  endm, int  startn, int  endn, char *  pszDebugMsg, bool  bLineFeed = false, bool  bShowImaginary = false  )

static

Show matrix partialy.

Parameters

pMatrix	The matrix that was wanted to show
startm	Start row index that was wanted to show
endm	End row index that was wanted to show
startn	Start column index that was wanted to show
endn	End column index that was wanted to show
pszDebugMsg	Text that was show first at matrix display
bLineFeed	Has linefeed between row or not
bShowImaginary	Determine show imaginary part or not

Definition at line 90 of file KNMatrixDebug.cpp.

References CKNMatrixOperation::CKNDMatrix::GetElement(), CKNComplex::GetImaginaryNumber(), and CKNComplex::GetRealNumber().

{
    ShowMsg(pszDebugMsg);
    ShowMsg("\n");
    ShowMsg("\n");

    for (int i = startm; i <= endm; i++)
    {
        char szRow[1024] = "";
        for (int j = startn; j <= endn; j++)
        {
            char            szCol[1024];
            CKNComplex      temppNumber = pMatrix->GetElement(i, j);

            sprintf(szCol, "%10.10f ", temppNumber.GetRealNumber());
            strcat(szRow, szCol);
            if (bShowImaginary)
            {
                sprintf(szCol, "+ %10.10f i ", temppNumber.GetImaginaryNumber());
                strcat(szRow, szCol);
            }
            strcat(szRow, "\t");
        }
        if (bLineFeed)
            strcat(szRow, "\n");
        strcat(szRow, "\n");
        ShowMsg(szRow);
    }
}

Here is the call graph for this function:

void CKNMatrixDebug::ShowDouble ( double  number, char *  pszDebugMsg  )

static

Show one double data.

Parameters

number	The number that was wanted to show
pszDebugMsg	Text that was show first at matrix display

Definition at line 259 of file KNMatrixDebug.cpp.

{
    char            szDebug[1024];

    sprintf(szDebug, "       %20.20f \n", number);

    ShowMsg(pszDebugMsg);
    ShowMsg("\n");
    ShowMsg(szDebug);
}

void CKNMatrixDebug::ShowDoubleArray ( double *  vector, char *  pszDebugMsg, int  due  )

static

Show double array.

Parameters

vector	The vector that was wanted to show
pszDebugMsg	Text that was show first at matrix display
due	Debugger show vector element unti this index

Definition at line 200 of file KNMatrixDebug.cpp.

{
    char        szDebug[1024] = "";
    int         i;

    ShowMsg(pszDebugMsg);
    ShowMsg("\n");

    for (i = 0; i < due; i++)
    {
        sprintf(szDebug, "%5d: %20.20f\n", i, vector[i]);
        ShowMsg(szDebug);
    }
}

void CKNMatrixDebug::ShowDoubleMatrix ( double *  pMatrix, int  m, int  n, int  startm, int  endm, int  startn, int  endn, char *  pszDebugMsg  )

static

Show matrix in output debugging windows - Visual studio and Windows DebugView only.

Parameters

pMatrix	The matrix that was wanted to show
m	Numbers of row
n	Numbers of column
startm	Start row index that was wanted to show
endm	End row index that was wanted to show
startn	Start column index that was wanted to show
endn	End column index that was wanted to show
pszDebugMsg	Text that was show first at matrix display

Definition at line 60 of file KNMatrixDebug.cpp.

{
    ShowMsg(pszDebugMsg);
    ShowMsg("\n");
    ShowMsg("\n");

    for (int i = startm; i < endm; i++)
    {
        char szRow[1024] = "";
        for (int j = startn; j < endn; j++)
        {
            char        szCol[1024];

            sprintf(szCol, "%20.20f ", pMatrix[i*m + j]);
            strcat(szRow, szCol);
        }
        strcat(szRow, "\n");
        ShowMsg(szRow);
    }
}

void CKNMatrixDebug::ShowDoubleVector ( CKNMatrixOperation::CKNVector  vector, char *  pszDebugMsg, int  due = -1  )

static

Parameters

vector	The vector that was wanted to show
pszDebugMsg	Text that was show first at matrix display
due	Debugger show vector element unti this index

Definition at line 220 of file KNMatrixDebug.cpp.

References CKNMatrixOperation::CKNVector::GetAt(), CKNComplex::GetImaginaryNumber(), CKNComplex::GetRealNumber(), and CKNMatrixOperation::CKNVector::GetSize().

{
    char        szDebug[1024] = "";
    int         i, nSize;

    if (-1 == due)
        nSize = vector.GetSize();
    else
        nSize = due;

    ShowMsg(pszDebugMsg);
    ShowMsg("\n");

    for (i = 0; i < nSize; i++)
    {
        sprintf(szDebug, "%5d: %20.20f   + %20.20f i\n", i, vector.GetAt(i).GetRealNumber(), vector.GetAt(i).GetImaginaryNumber());
        ShowMsg(szDebug);
    }
}

Here is the call graph for this function:

void CKNMatrixDebug::ShowMsg ( char *  pszBuffer )

static

Show message.

Parameters

pszBuffer

Showing message

Definition at line 35 of file KNMatrixDebug.cpp.

Referenced by CKNLanczosTest::AuditResult(), CKNLanczosTest::AuditResult_EV(), CKNLanczosTest::AuditResult_Seorth(), CKNLanczosTest::AuditResult_WF(), CKNLanczosTest::CompareWithMatLabSeOrth(), CKNLanczosTest::COmpareWIthMatLabSeOrthMPI(), CKNLanczosTest::LanczosThread(), CKNLanczosTest::LanczosThreadForMPI(), CKNLanczosTest::LargeSizeMatrixMPI(), CKNGeometricConstructionLaunch::LaunchingGeometricConstructionMPI(), CKNLanczosTest::SolvingLargeSizeHamlitonian(), CKNLanczosTest::TestCSRBuilding(), CKNLanczosTest::TestLanczos(), and CKNLanczosTest::TestSimpleLanczos().

{
    if (false == m_bShowMsg)
        return;
#ifdef _WIN32
#if _MFC_VER == 0x0C00
    OutputDebugStringA(pszBuffer);
#else// _MFC_VER == 0x0C00
    printf("%s", pszBuffer);
#endif //_MFC_VER == 0x0C00
#else
    printf("%s", pszBuffer);
#endif
}

Here is the caller graph for this function:

Member Data Documentation

bool CKNMatrixDebug::m_bShowMsg = false

staticprivate

Definition at line 38 of file KNMatrixDebug.h.

Referenced by SetShow().

---

The documentation for this class was generated from the following files:

• D:/For Doxygen/KNMatrixDebug.h
• D:/For Doxygen/KNMatrixDebug.cpp