CUtility Class Reference

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

#include "Utility.h"

Collaboration diagram for CUtility:

Public Member Functions
	CUtility ()
	~CUtility ()

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 (CMatrixOperation::CVector vector, char *pszDebugMsg, int due=-1)
static void	ShowComplex (CComplex 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 (CLanczosMethod::LPEIGENVALUE_RESULT lpResult, int nIndex)
	Show lanczos method result. More...
static void	SetShow (bool bShow)
static void	ShowCSR (CMatrixOperation::CCSR *pCSR, char *pszDebugMsg)
	Set show message or not. More...
static void	ShowCSR (CMatrixOperation::CCSR *pCSR, int startm, int endm, int startn, int endn, char *pszDebugMsg)
	Show CSR matrix partialy. More...
static void	ShowDenseMatrix (CMatrixOperation::CDMatrix *pMatrix, int startm, int endm, int startn, int endn, char *pszDebugMsg, bool bLineFeed=false, bool bShowImaginary=false)
	Show matrix partialy. More...
static void	DumpCSR (CMatrixOperation::CCSR *pCSR, char *pszFileName, double fAtomIDStartIndex)
static void	DumpCSR (CMatrixOperation::CCSR *pCSR, char *pszFileName, double fAtomIDStartIndex, int nRowPerAtom)
	<Added by="" jhkang>=""> More...
static void	ConvertLower (char *pszBuffer)
	<Added by="" jhkang="" end>=""> More...
static void	getCurrentTime (char *pszBuffer)
	Get Current time from system. More...
static void	GetKValues (CCommandFileParser::LPINPUT_CMD_PARAM lpParam, double *pKValue[3])
	Calculate K value with K points. 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

Definition at line 16 of file Utility.h.

Constructor & Destructor Documentation

CUtility::CUtility

(

)

Definition at line 23 of file Utility.cpp.

{
}

CUtility::~CUtility

(

)

Definition at line 27 of file Utility.cpp.

{
}

Member Function Documentation

void CUtility::ConvertLower ( char *  pszBuffer )

static

<Added by="" jhkang="" end>="">

Convert string to lower case

Parameters

pszBuffer[in/out]

Target string to convert to lower case

Definition at line 415 of file Utility.cpp.

Referenced by CMaterialParam::BuildMaterialParam().

{
    int                 i, nLen;

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

Here is the caller graph for this function:

void CUtility::DumpCSR ( CMatrixOperation::CCSR *  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 296 of file Utility.cpp.

References CMPIManager::GetCurrentRank(), CMatrixOperation::CCSR::GetRowCount(), CMPIManager::GetTotalNodeCount(), CMPIManager::IsDeflationRoot(), CMPIManager::IsRootRank(), CMatrixOperation::CCSR::m_vectColumn, CMatrixOperation::CCSR::m_vectRow, CMatrixOperation::CCSR::m_vectValueImaginaryBuffer, CMatrixOperation::CCSR::m_vectValueRealBuffer, CMPIManager::ReceiveDoubleBufferSync(), and CMPIManager::SendDoubleBufferSync().

Referenced by CSPLoop::BuildHamiltonian(), and CTBMS_Solver::Launching_TBMS_Solver().

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

    if (false == CMPIManager::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( CMPIManager::IsRootRank())
        out = fopen(szFileName, "wt");
    else
    {
        CMPIManager::ReceiveDoubleBufferSync(CMPIManager::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) + 1, 
                                                            pCSR->m_vectColumn[j] + 1, 
                                                            pCSR->m_vectValueRealBuffer[j], 
                                                            pCSR->m_vectValueImaginaryBuffer[j]);
            fputs(szMsg, out);
        }
    }
    fclose(out);

#ifndef DISABLE_MPI_ROUTINE
    if( CMPIManager::GetCurrentRank() < CMPIManager::GetTotalNodeCount() - 1)
        CMPIManager::SendDoubleBufferSync(CMPIManager::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 CUtility::DumpCSR ( CMatrixOperation::CCSR *  pCSR, char *  pszFileName, double  fAtomIDStartIndex, int  nRowPerAtom  )

static

<Added by="" jhkang>="">

<Added by="" jhkang="" end>="">

Dump CSR to text file

Definition at line 355 of file Utility.cpp.

References CMPIManager::GetCurrentRank(), CMatrixOperation::CCSR::GetRowCount(), CMPIManager::GetTotalNodeCount(), CMPIManager::IsRootRank(), CMatrixOperation::CCSR::m_vectColumn, CMatrixOperation::CCSR::m_vectRow, CMatrixOperation::CCSR::m_vectValueImaginaryBuffer, CMatrixOperation::CCSR::m_vectValueRealBuffer, CMPIManager::ReceiveDoubleBufferSync(), and CMPIManager::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( CMPIManager::IsRootRank())
        out = fopen(szFileName, "wt");
    else
    {
        CMPIManager::ReceiveDoubleBufferSync(CMPIManager::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*nRowPerAtom) + 1,
                                                            pCSR->m_vectColumn[j] + 1,
                                                            pCSR->m_vectValueRealBuffer[j],
                                                            pCSR->m_vectValueImaginaryBuffer[j]);
            fputs(szMsg, out);
        }
    }
    fclose(out);

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

Here is the call graph for this function:

void CUtility::getCurrentTime ( char *  pszBuffer )

static

Get Current time from system.

Parameters

pszBuffer[out]

Charater buffer that want to save string that include time information

Definition at line 430 of file Utility.cpp.

{
    time_t          timer;
    struct tm       *t;

    timer = time(NULL);
    t = localtime(&timer);

    sprintf(pszBuffer, "%d:%d:%d", t->tm_hour, t->tm_min, t->tm_sec);
}

void CUtility::GetKValues ( CCommandFileParser::LPINPUT_CMD_PARAM  lpParam, double *  pKValue[3]  )

static

Calculate K value with K points.

Parameters

lpParam	Option parameters for program launching
pKValue[out]	Buffer for saving k values

Definition at line 445 of file Utility.cpp.

References CCommandFileParser::INPUT_CMD_PARAM::fKPoints, CCommandFileParser::INPUT_CMD_PARAM::fKValueFinal, and CCommandFileParser::INPUT_CMD_PARAM::fKValueInit.

Referenced by CTBMS_Solver::Launching_TBMS_Solver(), and CSPLoop::SolveSchroedinger().

{
    double                  fInterval[3];
    unsigned int            i, j;

    for( i = 0; i < 3; ++i)
    {
        double              *pfKValueUnit = NULL;
        pKValue[i] = (double*)malloc(sizeof(double)*(unsigned int)lpParam->fKPoints);
        pfKValueUnit = pKValue[i];
        fInterval[i] = (lpParam->fKValueFinal[i] - lpParam->fKValueInit[i]) / (lpParam->fKPoints - 1);
        pfKValueUnit[0] = lpParam->fKValueInit[i];
        if (lpParam->fKPoints > 1 )
            pfKValueUnit[(int)lpParam->fKPoints-1] = lpParam->fKValueFinal[i];
        for( j = 1; j < (unsigned int)lpParam->fKPoints-1 ; ++j )
            pfKValueUnit[j] = pfKValueUnit[j-1] + fInterval[i];
    }
}

Here is the caller graph for this function:

void CUtility::SaveResult ( CLanczosMethod::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 273 of file Utility.cpp.

References CLanczosMethod::EIGENVALUE_RESULT::nEigenValueCount, and CLanczosMethod::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 CUtility::SetShow ( bool  bShow )

inlinestatic

Definition at line 29 of file Utility.h.

References m_bShowMsg.

Referenced by InitEnvironment(), CMPIManager::InitLevel(), CTBMS_Solver::Launching_TBMS_Solver(), and CLanczosLaunching::LaunchingLanczos().

{ m_bShowMsg = bShow; }; 

Here is the caller graph for this function:

void CUtility::ShowComplex ( CComplex  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 243 of file Utility.cpp.

References CComplex::GetImaginaryNumber(), and CComplex::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 CUtility::ShowCSR ( CMatrixOperation::CCSR *  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 158 of file Utility.cpp.

References DUMY_STRING_FOR_COMPLEX_SHOWING, CMatrixOperation::CCSR::GetColumnCount(), CMatrixOperation::CCSR::GetElement(), CComplex::GetImaginaryNumber(), CComplex::GetRealNumber(), and CMatrixOperation::CCSR::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];
            CComplex    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 CUtility::ShowCSR ( CMatrixOperation::CCSR *  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 127 of file Utility.cpp.

References DUMY_STRING_FOR_COMPLEX_SHOWING, CMatrixOperation::CCSR::GetElement(), CComplex::GetImaginaryNumber(), and CComplex::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];
            CComplex    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 CUtility::ShowDenseMatrix ( CMatrixOperation::CDMatrix *  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 89 of file Utility.cpp.

References CMatrixOperation::CDMatrix::GetElement(), CComplex::GetImaginaryNumber(), and CComplex::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];
            CComplex        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 CUtility::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 258 of file Utility.cpp.

{
    char            szDebug[1024];

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

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

void CUtility::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 199 of file Utility.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 CUtility::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 59 of file Utility.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 CUtility::ShowDoubleVector ( CMatrixOperation::CVector  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 219 of file Utility.cpp.

References CMatrixOperation::CVector::GetAt(), CComplex::GetImaginaryNumber(), CComplex::GetRealNumber(), and CMatrixOperation::CVector::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 CUtility::ShowMsg ( char *  pszBuffer )

static

Show message.

Parameters

pszBuffer

Showing message

Definition at line 34 of file Utility.cpp.

Referenced by CTBMS_Solver::ApplyPhPotential(), CLanczosResultAudit::AuditResult_EV(), CLanczosResultAudit::AuditResult_WF(), CSPLoop::BuildHamiltonian(), CGeometricShape::BuildPEBiasVector(), CGeometricShape::BuildPEWaveVector(), CGeometricShape::CalculateUnitcellCount(), CGeometricShape::ConstructContactRegionOnPoissonGrid(), CSPLoop::ConstructionGeometric(), CGeometricShape::ConstructMapInfo(), CSPLoop::executeSPLoop(), InitEnvironment(), CLanczosMethod::LanczosIterationLoop(), CTBMS_Solver::Launching_TBMS_Solver(), CLanczosLaunching::LaunchingLanczos(), main(), CGeometricShape::MapElecAtomOnPoissonGrid(), CSPLoop::SetInitialPotential(), CLanczosMethod::ShowLanczosResult(), CLanczosMethod::ShowLanczosWorkingTime(), CSPLoop::SolvePoisson(), and CSPLoop::SolveSchroedinger().

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

staticprivate

Flag for display message.

Definition at line 42 of file Utility.h.

Referenced by SetShow().

---

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

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