main.cpp

计算IsingModel的磁化强度

#include <windows.h>
#include <stdio.h>
#include <math.h>
#include "RandomGen250.h"
#include "ErrorReport.h"
#include "IsingModel.h"

CRandom CRandomGenerator;
CLog CResultLog;

void main()
{
	//unsigned long ulTotal,ulDrop,ulInterval;
	//printf( "input step: total, drop, interval:" );
	//scanf( "%u,%u,%u", &ulTotal, &ulDrop, &ulInterval );

	unsigned int uiLength;
	/*printf( "length:" );
	scanf( "%u", &uiLength );*/

	LARGE_INTEGER liCounter;
	QueryPerformanceCounter( &liCounter );
	CRandomGenerator.SetRandomSeed( (int)liCounter.LowPart );

	//CResultLog.SetOutput( "result.txt" );

	char strFilename[255];

	for( uiLength = 10; uiLength <= 60; uiLength+= 10 )
	{
		sprintf( strFilename, "PBC_L=%u.txt", uiLength );
		CResultLog.SetOutput( strFilename );
		for( double jkt = 0.01; jkt < 1.0; jkt += 0.01 )
		{
			CResultLog.WriteDoubleTab( jkt );
			IsingModel2D( uiLength, BC_TYPE_PBC, 500000, 5000, 320, jkt );
		}
		sprintf( strFilename, "FBC_L=%u.txt", uiLength );
		CResultLog.SetOutput( strFilename );
		for( jkt = 0.01; jkt < 1.0; jkt += 0.01 )
		{
			CResultLog.WriteDoubleTab( jkt );
			IsingModel2D( uiLength, BC_TYPE_FBC, 500000, 5000, 320, jkt );
		}
	}

	//IsingModel2D( uiLength, BC_TYPE_PBC, 10000, 5000, 20, 0.44 );
}