hypergeometricdlg.cpp

A program for calculation of hypergeometric distribution

// HypergeometricDlg.cpp : implementation file
//

#include "stdafx.h"
#include "Hypergeometric.h"
#include "HypergeometricDlg.h"
#include "HyperGeometricProbability.h"

#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

int nFactorialFunctionCall;
int nDPFactorialFunctionCall;
int nDPFactorialFunctionCall2;
int nDPFactorialFunctionCall3;

/////////////////////////////////////////////////////////////////////////////
// CAboutDlg dialog used for App About

class CAboutDlg : public CDialog
{
public:
	CAboutDlg();

// Dialog Data
	//{{AFX_DATA(CAboutDlg)
	enum { IDD = IDD_ABOUTBOX };
	//}}AFX_DATA

	// ClassWizard generated virtual function overrides
	//{{AFX_VIRTUAL(CAboutDlg)
	protected:
	virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV support
	//}}AFX_VIRTUAL

// Implementation
protected:
	//{{AFX_MSG(CAboutDlg)
	//}}AFX_MSG
	DECLARE_MESSAGE_MAP()
};

CAboutDlg::CAboutDlg() : CDialog(CAboutDlg::IDD)
{
	//{{AFX_DATA_INIT(CAboutDlg)
	//}}AFX_DATA_INIT
}

void CAboutDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CAboutDlg)
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
	//{{AFX_MSG_MAP(CAboutDlg)
		// No message handlers
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CHypergeometricDlg dialog

CHypergeometricDlg::CHypergeometricDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CHypergeometricDlg::IDD, pParent)
{
	//{{AFX_DATA_INIT(CHypergeometricDlg)
	m_doubeResult = 0.0;
	m_nC = 0;
	m_nG = 0;
	m_nK = 0;
	m_nN = 0;
	m_dwCalcTime = 0;
	//}}AFX_DATA_INIT
	// Note that LoadIcon does not require a subsequent DestroyIcon in Win32
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
	m_vecMemFactorial.push_back(new CLargeNumber("1"));
	m_vecMemFactorialTmp.push_back(new string("1"));
	m_vecMemFactorialTmp2.push_back( "1" );
	nFactorialFunctionCall = 0;
	nDPFactorialFunctionCall = 0;
	nDPFactorialFunctionCall2 = 0;
}

void CHypergeometricDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CHypergeometricDlg)
	DDX_Control(pDX, IDC_PROGRESSBAR, m_ctrlProgress);
	DDX_Text(pDX, IDC_EDIT_P, m_doubeResult);
	DDX_Text(pDX, IDC_EDIT_C, m_nC);
	DDX_Text(pDX, IDC_EDIT_G, m_nG);
	DDX_Text(pDX, IDC_EDIT_K, m_nK);
	DDX_Text(pDX, IDC_EDIT_N, m_nN);
	DDX_Text(pDX, IDC_EDIT_CALCTIME, m_dwCalcTime);
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CHypergeometricDlg, CDialog)
	//{{AFX_MSG_MAP(CHypergeometricDlg)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	ON_BN_CLICKED(IDC_BTN_CALCULATE, OnBtnCalculate)
	ON_BN_CLICKED(IDC_BTN_DPCALCULATE, OnBtnDpcalculate)
	ON_WM_DESTROY()
	ON_BN_CLICKED(IDC_BTN_DPCALCULATE2, OnBtnDpcalculate2)
	ON_BN_CLICKED(IDC_BTN_DPCALCULATE3, OnBtnDpcalculate3)
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CHypergeometricDlg message handlers

BOOL CHypergeometricDlg::OnInitDialog()
{
	CDialog::OnInitDialog();

	// Add "About..." menu item to system menu.

	// IDM_ABOUTBOX must be in the system command range.
	ASSERT((IDM_ABOUTBOX & 0xFFF0) == IDM_ABOUTBOX);
	ASSERT(IDM_ABOUTBOX < 0xF000);

	CMenu* pSysMenu = GetSystemMenu(FALSE);
	if (pSysMenu != NULL)
	{
		CString strAboutMenu;
		strAboutMenu.LoadString(IDS_ABOUTBOX);
		if (!strAboutMenu.IsEmpty())
		{
			pSysMenu->AppendMenu(MF_SEPARATOR);
			pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);
		}
	}

	// Set the icon for this dialog.  The framework does this automatically
	//  when the application's main window is not a dialog
	SetIcon(m_hIcon, TRUE);			// Set big icon
	SetIcon(m_hIcon, FALSE);		// Set small icon
	
	// TODO: Add extra initialization here
	
	return TRUE;  // return TRUE  unless you set the focus to a control
}

void CHypergeometricDlg::OnSysCommand(UINT nID, LPARAM lParam)
{
	if ((nID & 0xFFF0) == IDM_ABOUTBOX)
	{
		CAboutDlg dlgAbout;
		dlgAbout.DoModal();
	}
	else
	{
		CDialog::OnSysCommand(nID, lParam);
	}
}

// If you add a minimize button to your dialog, you will need the code below
//  to draw the icon.  For MFC applications using the document/view model,
//  this is automatically done for you by the framework.

void CHypergeometricDlg::OnPaint() 
{
	if (IsIconic())
	{
		CPaintDC dc(this); // device context for painting

		SendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);

		// Center icon in client rectangle
		int cxIcon = GetSystemMetrics(SM_CXICON);
		int cyIcon = GetSystemMetrics(SM_CYICON);
		CRect rect;
		GetClientRect(&rect);
		int x = (rect.Width() - cxIcon + 1) / 2;
		int y = (rect.Height() - cyIcon + 1) / 2;

		// Draw the icon
		dc.DrawIcon(x, y, m_hIcon);
	}
	else
	{
		CDialog::OnPaint();
	}
}

// The system calls this to obtain the cursor to display while the user drags
//  the minimized window.
HCURSOR CHypergeometricDlg::OnQueryDragIcon()
{
	return (HCURSOR) m_hIcon;
}

CLargeNumber CHypergeometricDlg::Factorial(int iNumber)
{
	CLargeNumber oLN("1");
	if(iNumber > 1)
	{
		for(int i=2; i<=iNumber; i++)
		{
			oLN *= i;
			TRACE1("%d\n", nFactorialFunctionCall++);
		}
	}
	return oLN;
}

CLargeNumber CHypergeometricDlg::DPFactorial2(int iNumber)
{
	int nSize = m_vecMemFactorialTmp2.size();
	if ( iNumber < nSize )
		return CLargeNumber( (m_vecMemFactorialTmp2[iNumber]) );
	CLargeNumber oLN = CLargeNumber( m_vecMemFactorialTmp2[nSize-1] );

	for(int i = nSize; i<=iNumber; i++)
	{
		oLN *= i;
		m_vecMemFactorialTmp2.push_back( oLN.ToString() );
		TRACE1("%d\n", nDPFactorialFunctionCall2++);
	}

	return oLN;
}


CLargeNumber CHypergeometricDlg::DPFactorial(int iNumber)
{
	int nSize = m_vecMemFactorialTmp.size();
	if ( iNumber < nSize )
		return CLargeNumber( *(m_vecMemFactorialTmp[iNumber]) );
	CLargeNumber oLN = CLargeNumber( *m_vecMemFactorialTmp[nSize-1] );

	for(int i = nSize; i<=iNumber; i++)
	{
		oLN *= i;
		m_vecMemFactorialTmp.push_back( new string(oLN.ToString()) );
		TRACE1("%d\n", nDPFactorialFunctionCall++);
	}

	return oLN;
}

double CHypergeometricDlg::PValue(int nG, int nC, int nN, int nK)
{
	double Result = 1;
	CLargeNumber lnTemp1 = Combination( nG, nN );
	CLargeNumber lnTemp2;
	double doubleTemp1;
	int nTemp2 = 10000;
	CLargeNumber lnTemp3("10000");
	for( int i = 0; i < nK - 1; i++ )
	{
		m_ctrlProgress.StepIt();
		lnTemp2 = Combination( nC, i ) * Combination( nG - nC, nN - i );
		lnTemp2 *= lnTemp3;
		lnTemp2 /= lnTemp1 ;

		doubleTemp1 = atoi(lnTemp2.ToString().data());
		doubleTemp1 /= nTemp2;
		
		Result -= doubleTemp1;
	}
	return Result;
}

double CHypergeometricDlg::DPPValue(int nG, int nC, int nN, int nK)
{
	double Result = 1;
	CLargeNumber lnTemp1 = DPCombination( nG, nN );
	CLargeNumber lnTemp2;
	double doubleTemp1;
	int nTemp2 = 10000;
	CLargeNumber lnTemp3("10000");
	for( int i = 0; i < nK - 1; i++ )
	{
		m_ctrlProgress.StepIt();
		lnTemp2 = DPCombination( nC, i ) * DPCombination( nG - nC, nN - i );
		lnTemp2 *= lnTemp3;
		lnTemp2 /= lnTemp1 ;

		doubleTemp1 = atoi(lnTemp2.ToString().data());
		doubleTemp1 /= nTemp2;
		
		Result -= doubleTemp1;
	}
	return Result;
}

double CHypergeometricDlg::DPPValue2(int nG, int nC, int nN, int nK)
{
	double Result = 1;
	CLargeNumber lnTemp1 = DPCombination2( nG, nN );
	CLargeNumber lnTemp2;
	double doubleTemp1;
	int nTemp2 = 10000;
	CLargeNumber lnTemp3("10000");
	for( int i = 0; i < nK - 1; i++ )
	{
		m_ctrlProgress.StepIt();
		lnTemp2 = DPCombination2( nC, i ) * DPCombination2( nG - nC, nN - i );
		lnTemp2 *= lnTemp3;
		lnTemp2 /= lnTemp1 ;

		doubleTemp1 = atoi(lnTemp2.ToString().data());
		doubleTemp1 /= nTemp2;
		
		Result -= doubleTemp1;
	}
	return Result;
}

void CHypergeometricDlg::OnBtnCalculate() 
{
	// TODO: Add your control notification handler code here
	UpdateData();
	
	m_ctrlProgress.SetRange( 0, m_nK - 1 );
	m_ctrlProgress.SetPos(0);
	m_ctrlProgress.SetStep(1);
	DWORD dwStart = GetTickCount();
	m_doubeResult = PValue( m_nG, m_nC, m_nN, m_nK );
	DWORD dwEnd = GetTickCount();

	m_dwCalcTime = dwEnd - dwStart;

	UpdateData( FALSE );

	nFactorialFunctionCall = 0;
	nDPFactorialFunctionCall = 0;
	nDPFactorialFunctionCall2 = 0;
}

void CHypergeometricDlg::OnBtnDpcalculate() 
{
	UpdateData();
	
	m_ctrlProgress.SetRange( 0, m_nK - 1 );
	m_ctrlProgress.SetPos(0);
	m_ctrlProgress.SetStep(1);

	DWORD dwStart = GetTickCount();

	m_doubeResult = DPPValue( m_nG, m_nC, m_nN, m_nK );
	
	DWORD dwEnd = GetTickCount();
	
	m_dwCalcTime = dwEnd - dwStart;
	UpdateData( FALSE );	

	nFactorialFunctionCall = 0;
	nDPFactorialFunctionCall = 0;
	nDPFactorialFunctionCall2 = 0;
}

void CHypergeometricDlg::OnBtnDpcalculate2() 
{
	UpdateData();
	
	m_ctrlProgress.SetRange( 0, m_nK - 1 );
	m_ctrlProgress.SetPos(0);
	m_ctrlProgress.SetStep(1);
	
	DWORD dwStart = GetTickCount();
	
	m_doubeResult = DPPValue2( m_nG, m_nC, m_nN, m_nK );
	
	DWORD dwEnd = GetTickCount();
	
	m_dwCalcTime = dwEnd - dwStart;
	UpdateData( FALSE );

	nFactorialFunctionCall = 0;
	nDPFactorialFunctionCall = 0;
	nDPFactorialFunctionCall2 = 0;
}

void CHypergeometricDlg::OnBtnDpcalculate3() 
{
	UpdateData();
	m_ctrlProgress.SetRange( 0, m_nK - 1 );
	m_ctrlProgress.SetPos(0);
	m_ctrlProgress.SetStep(1);

	DWORD dwStart = GetTickCount();
	
	double doubleResult = 1.0;
	for( int i = 0 ; i < m_nK - 1; i++ )
	{
		HyperGeometricProbability hgp( m_nG, m_nN, i, m_nC );
		doubleResult -= hgp.GetResult();
	}
	m_doubeResult = doubleResult;

	DWORD dwEnd = GetTickCount();
	
	m_dwCalcTime = dwEnd - dwStart;
	UpdateData( FALSE );

	nFactorialFunctionCall = 0;
	nDPFactorialFunctionCall = 0;
	nDPFactorialFunctionCall2 = 0;
}

void CHypergeometricDlg::OnDestroy() 
{
	CDialog::OnDestroy();
	
	// TODO: Add your message handler code here
	for( int i = 0; i < m_vecMemFactorial.size(); i++ )
	{
		delete m_vecMemFactorial[i];
	}
	for( i = 0; i < m_vecMemFactorialTmp.size(); i++ )
	{
		delete m_vecMemFactorialTmp[i];
	}
}