quote:

---

So what else files such as internal.c I should add to that System folder?

---

quote:

---

Now it work. only 8.42s to do so

---

 
bool MakeMatrixCircle(Matrix &mMat, double dOutsideVal = 0)
{
	if(mMat.IsValid())
	{
		int nCols=mMat.GetNumCols();
		int nRows=mMat.GetNumRows();
		int i, nCells;
		if(nRows>1 && nCols>1)
		{
			vector vctCol;
			nCells=nCols;
			if(nRows!=nCols)
			{
				if(nRows<nCols)
				{
					for(i=nRows;i<nCols;i++)
					{
						if(mMat.GetColumn(vctCol,i)==false) return (false);
						vctCol*=0;
						vctCol+=dOutsideVal;
						if(mMat.SetColumn(vctCol,i)!=0) return (false);
					}
					nCells=nRows;
				}
				else
				{
					for(i=nCols;i<nRows;i++)
					{
						if(mMat.GetRow(vctCol,i)==false) return (false);
						vctCol*=0;
						vctCol+=dOutsideVal;
						if(mMat.SetRow(vctCol,i)!=0) return (false);
					}
					nCells=nCols;
				}
			}
			double Radius=nCells-1;
			Radius/=2;
			double Rsq=Radius^2;
			double *RSQ;
			RSQ=&Rsq;
			int nHalf;
			vector vctMaskMult, vctMaskAdd;
			if(is_odd(nCells)==true)
			{
				nHalf=1+((nCells-1)/2);
			}
			else
			{
				nHalf=nCells/2;
			}
			int k=nCols-1;
			double *repVal;
			double repValue=dOutsideVal;
			repVal=&repValue;
			progressBox progBox("Modifying Matrix ... ",PBOX_TOPMOST);
			progBox.SetRange(0,nHalf);
			vctMaskMult.SetSize(nCells);
			vctMaskAdd.SetSize(nCells);
			double *vctAdd, *vctAdd1, *vctMult, *vctMult1;
			for(i=0;i<nHalf;i++)
			{
				progBox.Set(i);
				if(mMat.GetColumn(vctCol,i)==false) return (false);
				vctAdd=&vctMaskAdd[0];
				vctMult=&vctMaskMult[0];
				vctAdd1=&vctMaskAdd[nCells-1];
				vctMult1=&vctMaskMult[nCells-1];				
				double cVal=i+1;
				double *RCval, rcVal=(Radius-cVal)^2;
				RCval=&rcVal;
				for(int j=0;j<nHalf;j++)
				{
					double rVal=j+1;
					if((*RSQ)<=(((Radius-rVal)^2)+(*RCval)))
					{
						*vctMult=0;
						*vctMult1=0;
						*vctAdd=*repVal;
						*vctAdd1=*repVal;
					}
					else
					{
						*vctMult=1;
						*vctAdd=0;
						*vctMult1=1;
						*vctAdd1=0;
					}
					vctMult++;
					vctAdd++;
					vctMult1--;
					vctAdd1--;					
				}
				vctCol*=vctMaskMult;
				vctCol+=vctMaskAdd;
				if(mMat.SetColumn(vctCol,i)!=0) return (false);
				if(mMat.GetColumn(vctCol,k)==false) return (false);
				vctCol*=vctMaskMult;
				vctCol+=vctMaskAdd;
				if(mMat.SetColumn(vctCol,k)!=0) return (false);
				k--;
			}
			return (true);
		}
	}
	return (false);
}

void mkcircle(double dOutsideVal = 0)
{
	MatrixLayer ml = Project.ActiveLayer();
	
	Matrix mm(ml);
	int nRows = mm.GetNumRows();
	int nCols = mm.GetNumCols();
	int nrc = 0.5 + nRows/2.0; // nearest int
	int ncc = 0.5 + nCols/2.0;
	double rr = min(nrc, ncc);
	
	for(int nr = 0; nr < nRows; nr++)
	{
		for(int nc = 0; nc < mm.GetNumCols(); nc++)
		{
			double xsqr = (nr - nrc)^2;
			double ysqr = (nc - ncc)^2;
			if(sqrt(xsqr + ysqr) > rr)
			{
				mm[nr][nc] = dOutsideVal;
			}
		}
	}
}