CKNIPCCUtility Class Reference

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

#include "KNIPCCUtility.h"

Collaboration diagram for CKNIPCCUtility:

Public Member Functions
	CKNIPCCUtility ()
	~CKNIPCCUtility ()

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 void	ConvertLower (char *pszBuffer)
	Convert string to lower case. More...

Static Private Attributes
static bool	m_bShowMsg = false
	Flag for display message. More...

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), Hoon Ryu(elec1.nosp@m.020@.nosp@m.gmail.nosp@m..com)

Definition at line 17 of file KNIPCCUtility.h.

Constructor & Destructor Documentation

CKNIPCCUtility::CKNIPCCUtility

(

)

Definition at line 24 of file KNIPCCUtility.cpp.

{
}

CKNIPCCUtility::~CKNIPCCUtility

(

)

Definition at line 28 of file KNIPCCUtility.cpp.

{
}

Member Function Documentation

void CKNIPCCUtility::ConvertLower ( char *  pszBuffer )

static

Convert string to lower case.

Parameters

pszBuffer[in/out]

Target string to convert to lower case

Definition at line 411 of file KNIPCCUtility.cpp.

Referenced by CKNMaterialParam::BuildMaterialParam().

{
    int                 i, nLen;

    nLen = strlen(pszBuffer);
    for (i = 0; i < nLen; ++i)
    {
        if (pszBuffer[i] < 'Z')
            pszBuffer[i] += 32;
    }
}

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

static

Dump CSR to text file.

Parameters

pCSR	CSR matrix
pszFileName	Save file name
fAtomIDStartIndex	Start index in shape

Definition at line 350 of file KNIPCCUtility.cpp.

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

Referenced by CKNTBMS_Solver::Launching_TBMS_Solver().

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

    if( false == CKNMPIManager::IsDeflationRoot() )
        return;

#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_vectValueRealBuffer[j], 
                                                            pCSR->m_vectValueImaginaryBuffer[j]);
            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
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

static

Dump CSR to binary file.

Parameters

pCSR	CSR matrix
pszFileName	Save file name
fAtomIDStartIndex	Start index in shape

Definition at line 297 of file KNIPCCUtility.cpp.

References CKNMPIManager::GetCurrentRank(), CKNMatrixOperation::CKNCSR::GetRowCount(), CKNMPIManager::GetTotalNodeCount(), CKNMPIManager::IsRootRank(), CKNMatrixOperation::CKNCSR::m_vectColumn, CKNMatrixOperation::CKNCSR::m_vectRow, CKNMatrixOperation::CKNCSR::m_vectValueImaginaryBuffer, CKNMatrixOperation::CKNCSR::m_vectValueRealBuffer, 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_vectValueRealBuffer[j];
            fwrite(&fTemp, sizeof(double), 1, out);
            fTemp = pCSR->m_vectValueImaginaryBuffer[j];
            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 CKNIPCCUtility::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 KNIPCCUtility.cpp.

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

{
    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);
    }
}

static void CKNIPCCUtility::SetShow ( bool  bShow )

inlinestatic

Definition at line 30 of file KNIPCCUtility.h.

References m_bShowMsg.

Referenced by CKNMPIManager::InitLevel(), CKNTBMS_Solver::Launching_TBMS_Solver(), and CKNLanczosLaunching::LaunchingLanczos().

{ m_bShowMsg = bShow; }; 

Here is the caller graph for this function:

void CKNIPCCUtility::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 KNIPCCUtility.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 CKNIPCCUtility::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 KNIPCCUtility.cpp.

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

{
    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:

void CKNIPCCUtility::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 KNIPCCUtility.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 CKNIPCCUtility::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 KNIPCCUtility.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 CKNIPCCUtility::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 KNIPCCUtility.cpp.

{
    char            szDebug[1024];

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

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

void CKNIPCCUtility::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 KNIPCCUtility.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 CKNIPCCUtility::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 KNIPCCUtility.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 CKNIPCCUtility::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 KNIPCCUtility.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 CKNIPCCUtility::ShowMsg ( char *  pszBuffer )

static

Show message.

Parameters

pszBuffer

Showing message

Definition at line 35 of file KNIPCCUtility.cpp.

Referenced by CKNTBMS_Solver::ApplyPhPotential(), CKNLanczosResultAudit::AuditResult_EV(), CKNLanczosResultAudit::AuditResult_WF(), CKNGeometricShape::CalculateUnitcellCount(), CKNGeometricShape::ConstructMapInfo(), CKNLanczosMethod::LanczosIterationLoop(), CKNTBMS_Solver::Launching_TBMS_Solver(), CKNLanczosLaunching::LaunchingLanczos(), CKNMPIManager::LoadBlancingForLanczos(), and CKNLanczosMethod::ShowLanczosResult().

{
    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 CKNIPCCUtility::m_bShowMsg = false

staticprivate

Flag for display message.

Definition at line 39 of file KNIPCCUtility.h.

Referenced by SetShow().

---

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

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