//?Contour.cpp:?implementation?of?the?CContour?class.
//
//////////////////////////////////////////////////////////////////////

#include?“stdafx.h“
#include?“stdlib.h“
#include?“math.h“
#include?“Contour.h“

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

double?TestFunction（double?xdouble?y）
{??
	return?0.5*（cos（x+3.14/4）+sin（y+3.14/4））;?
};

//////////////////////////////////////////////////////////////////////
//?Construction/Destruction
//////////////////////////////////////////////////////////////////////

CContour::CContour（）
{
	m_iColFir=m_iRowFir=32;
	m_iColSec=m_iRowSec=256;
	m_dDx=m_dDy=0;
	m_pFieldFcn=NULL;
	m_pLimits[0]=m_pLimits[2]=0;
	m_pLimits[1]=m_pLimits[3]=5.;
	m_ppFnData=NULL;

	//?temporary?stuff
	m_pFieldFcn=TestFunction;
	m_vPlanes.resize（20）;
	for?（UINT?i=0;i	{
		m_vPlanes[i]=（i-m_vPlanes.size（）/2.0）*0.1;
	}
}

CContour::~CContour（）
{
	CleanMemory（）;
}

void?CContour::InitMemory（）
{
	if?（!m_ppFnData）
	{
		m_ppFnData=new?CFnStr*[m_iColSec+1];
		for?（int?i=0;i		{
			m_ppFnData[i]=NULL;
		}
	}
}

void?CContour::CleanMemory（）
{
	if?（m_ppFnData）
	{
		int?i;
		for?（i=0;i		{
			if?（m_ppFnData[i]）
				delete[]?（m_ppFnData[i]）;
		}
		delete[]?m_ppFnData;
		m_ppFnData=NULL;
	}
}

void?CContour::Generate（）
{

	int?i?j;
	int?x3?x4?y3?y4?x?y?oldx3?xlow;
	const?int?cols=m_iColSec+1;
	const?int?rows=m_iRowSec+1;
	double?xoffyoff;
	
	//?Initialize?memroy?if?needed
	InitMemory（）;

	m_dDx?=?（m_pLimits[1]-m_pLimits[0]）/（double）（m_iColSec）;
	xoff?=?m_pLimits[0];
	m_dDy?=?（m_pLimits[3]-m_pLimits[2]）/（double）（m_iRowSec）;
	yoff?=?m_pLimits[2];

	xlow?=?0;
	oldx3?=?0;
	x3?=?（cols-1）/m_iRowFir;
	x4?=?（?2*（cols-1）?）/m_iRowFir;
	for?（x?=?oldx3;?x?<=?x4;?x++）?
	{	??/*?allocate?new?columns?needed
		*/
		if?（x?>=?cols）
			break;
		if?（m_ppFnData[x]==NULL）
			m_ppFnData[x]?=?new?CFnStr[rows];

		for?（y?=?0;?y?			FnctData（xy）->m_sTopLen?=?-1;
	}

	y4?=?0;
	for?（j?=?0;?j?	{
		y3?=?y4;
		y4?=?（（j+1）*（rows-1））/m_iColFir;
		Cntr1（oldx3?x3?y3?y4）;
	}

	for?（i?=?1;?i?	{
		y4?=?0;
		for?（j?=?0;?j?		{
			y3?=?y4;
			y4?=?（（j+1）*（rows-1））/m_iColFir;
			Cntr1（x3?x4?y3?y4）;
		}

		y4?=?0;
		for?（j?=?0;?j?		{
			y3?=?y4;
			y4?=?（（j+1）*（rows-1））/m_iColFir;
			Pass2（oldx3x3y3y4）;
		}

		if?（i?		{	?/*?re-use?columns?no?longer?needed?*/
			oldx3?=?x3;
			x3?=?x4;
			x4?=?（（i+2）*（cols-1））/m_iRowFir;
			for?（x?=?x3+1;?x?<=?x4;?x++）?
			{
				if?（xlow?				{
					if?（m_ppFnData[x]）
						delete[]?m_ppFnData[x];
					m_ppFnData[x]?=?m_ppFnData[xlow];
					m_ppFnData[?xlow++?]?=?NULL;
				}?
				else
					if?（m_ppFnData[x]==NULL）
						m_ppFnData[x]?=?new?CFnStr[rows];

				for?（y?=?0;?