w1.cpp

分子动力学的VC源码,可以方面的用语化学与材料计算方面的开发,希望有用.

// W1.cpp : implementation file
//

#include "stdafx.h"
#include "mds.h"
#include "W1.h"
#include "dlg.h"
//#include "MainFrm.h"
#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

/////////////////////////////////////////////////////////////////////////////
// W1

IMPLEMENT_DYNCREATE(W1, CMDIChildWnd)

W1::W1()
{
}

W1::~W1()
{
}
BOOL W1::Create(LPCTSTR szTitle, LONG style /* = 0 */,
	const RECT& rect /* = rectDefault */,
	CMDIFrameWnd* parent /* = NULL */)
{

	LPCTSTR lpszBounceClass =
		AfxRegisterWndClass(CS_HREDRAW | CS_VREDRAW,
			LoadCursor(NULL, IDC_UPARROW),
			(HBRUSH)(COLOR_WINDOW+1),
			NULL);

	return CMDIChildWnd::Create(lpszBounceClass, szTitle, style, rect, parent);



}

BEGIN_MESSAGE_MAP(W1, CMDIChildWnd)
	//{{AFX_MSG_MAP(W1)
	ON_WM_PAINT()
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// W1 message handlers

void W1::OnPaint() 
{
		dlg DLG;


	CPaintDC dc(this); 



    int i,j;

    double R=0.361*40;
    double r=3.5;
    double X[20000],Y[20000];
    int m_nPenWidth1=1;
    int m_nPenWidth=1;

	int n0,n1,n2,n3;

    FILE *fp;

	   fp=fopen("k1.txt","r");
       fscanf(fp,"%d",&n0);
       fclose(fp);

       fp=fopen("lx.txt","r");
       fscanf(fp,"%d",&n1);
       fclose(fp);

       fp=fopen("wy.txt","r");
       fscanf(fp,"%d",&n2);
       fclose(fp);

       fp=fopen("hz.txt","r");
       fscanf(fp,"%d",&n3);
       fclose(fp);

/*
	n0=MF.k1;
	n1=DLG.lx;
	n2=DLG.wy;
	n3=DLG.hz;
*/


        CPen  penBlue(PS_SOLID,m_nPenWidth1,RGB(80,0,255));
		CPen*   pOldPen=dc.SelectObject(&penBlue);        


			
        dc.MoveTo(210,60);
		dc.LineTo(210+R*n2,60);

        dc.MoveTo(210+R*n2,60);
		dc.LineTo(210+R*n2,60+R*n3);

        dc.MoveTo(210+R*n2,60+R*n3);
		dc.LineTo(210,60+R*n3);

        dc.MoveTo(210,60+R*n3);
		dc.LineTo(210,60);



       fp=fopen("output.txt","r");
       for(i=1;i<=n0;i++)
	   {
        	fscanf(fp,"%le%le",&X[i],&Y[i]);

			CRect	rcEllipse1(X[i]*40+210-r,Y[i]*40+60-r, X[i]*40+210+r, Y[i]*40+60+r);
			dc.Ellipse(rcEllipse1);

	   }
       fclose(fp);

/*
	    CPen  penRed(PS_SOLID,m_nPenWidth,RGB(255,0,0));
		CPen*  pOldPen1=dc.SelectObject(&penRed);  

	   fp=fopen("proutput.txt","r");
       for(i=1;i<=n0;i++)
	   {
        	fscanf(fp,"%le%le",&X[i],&Y[i]);

			CRect	rcEllipse1(X[i]*40+210-r,Y[i]*40+60-r, X[i]*40+210+r, Y[i]*40+60+r);
			dc.Ellipse(rcEllipse1);

	   }
       fclose(fp);

*/
}