#pragma rtGlobals=1		// Use modern global access method.
Menu "simulation"
	"Calc SAW or Gauss Compiled"
	"-"
	submenu "Parts"
		"calculate gaussian"
		"Calculate SAW"
	end
end

Macro CalculateGaussian(name,npt)
	string/g gname
	variable/g gnpt
	string name=gname
	variable npt=gnpt
	prompt name,"Enter the base name for x, y z waves"
	prompt npt,"Enter number of steps in chain (10,000)"
	
	silent 1
	gname=name;gnpt=npt
	variable counter=1//npt=10000,
	string xpos="xpos_"+name,ypos="ypos_"+name,zpos="zpos_"+name
	make/n=(npt)/o $xpos=0, $ypos=0, $zpos=0
	variable xtst=enoise(100),ytst=enoise(100),ztst=enoise(100),addint=enoise(100)
	do
		xtst=enoise(100);ytst=enoise(100);ztst=enoise(100);addint=enoise(100)
		addint=(addint/abs(addint))
		do
			if((xtst>ytst)&&(xtst>ztst))
				$xpos[counter]=addint+$xpos[counter-1]
				$ypos[counter]=$ypos[counter-1]
				$zpos[counter]=$zpos[counter-1]
				break
			endif
			if((ytst>xtst)&&(ytst>ztst))
				$ypos[counter]=addint+$ypos[counter-1]
				$xpos[counter]=$xpos[counter-1]
				$zpos[counter]=$zpos[counter-1]
				break
			endif
			if((ztst>ytst)&&(ztst>xtst))
				$zpos[counter]=addint+$zpos[counter-1]
				$ypos[counter]=$ypos[counter-1]
				$xpos[counter]=$xpos[counter-1]
				break
			endif
			break
		while(0)
		counter+=1
		print counter
	while(counter<npt)
ENdMacro

macro GetName_xpos(xpos)
	string/g gxpos
	string xpos=gxpos
	prompt xpos,"Enter the xwave to get name from:",popup,wavelist("xpos_*", ";", "")
	
	silent 1
	gxpos=xpos
	string/g gname=xpos[5,50]
	variable/g gnpt=numpnts($xpos)
end

Macro CalcSAWorGaussCompiled(name,npt,SAWorGauss)
	string/g gname
	variable/g gnpt,gSAWorGauss
	string name=gname
	variable npt=gnpt,SAWorGauss=gSAWorGauss
	prompt name,"Enter the base name for x, y z waves (not for just plot)"
	prompt npt,"Enter number of steps in chain (10,000)"
	prompt SAWorGauss,"SAW or Gaussian walk or Just plot results?",popup,"SAW;Gaussian Walk;Just Plot Results"
	
	silent 1
	if(SAWorGauss==3)
		GetName_xpos()//returns gname and gnpt
		name=gname
		npt=gnpt
	endif
	gname=name;gnpt=npt;gSAWorGauss=SAWorGauss
	variable counter=1 // npt=10000,
	string xpos="xpos_"+name,ypos="ypos_"+name,zpos="zpos_"+name,rpos="rpos_"+name
	make/n=(npt)/o $xpos, $ypos, $zpos,$rpos
	if(SAWorGauss==1)
		CalcSAW($xpos,$ypos,$zpos,npt)
	else
		if(SAWorGauss==2)
			CalcGauss($xpos,$ypos,$zpos,npt)
		endif
	endif
	
	//2d plot of results
	if (exists("CreateSurfer") !=4)
		DoAlert 0, "Surface Plotter XOP must be installed"
		return
	endif
	CreateSurfer
	MoveWindow 37.5,57.5,428.25,318.5 
	ModifySurfer  FactoryDefaults, Update=0
	ModifySurfer/N=Surface1
	ModifySurfer srcWave=($xpos,$ypos,$zpos)
	ModifySurfer  srcType=4,plotType=5
	ModifySurfer  setControlView=3
	ModifySurfer  theta=169.333,  phi=298.833,  zScale=1,  xStep=2,  yStep=2
	ModifySurfer  frame=895,  drawFrame=1
	ModifySurfer  drawBox=1
	ModifySurfer palette=Rainbow
	ModifySurfer markerRGB={0,39168,0}
	ModifySurfer  marker=8,  markerSize=2
	ModifySurfer  Update=1
	
	$rpos=sqrt($xpos^2+$ypos^2+$zpos^2)
	print "End-to-end distance = "+num2str($rpos[npt-1])
	
	//plot eted
	Display /W=(5.25,41.75,474.75,396.5) $rpos
	ModifyGraph grid=2
	ModifyGraph log=1
	ModifyGraph tick=2
	ModifyGraph mirror=1
	ModifyGraph standoff=0
	Label left "End-toEnd-Distance"
	Label bottom "N"
	SetAxis left 1,100
	SetAxis bottom 1,10000
	ShowTools
	SetDrawLayer UserFront
	SetDrawEnv xcoord= bottom,ycoord= left
	DrawLine 1.01859138805411,1.01124030028642,10185.9138805411,102.260694492338
	SetDrawEnv xcoord= bottom,ycoord= left
	DrawLine 1,1,100000,1000
EndMacro

FUnction CalcGauss(xpos,ypos,zpos,npt)
	wave xpos,ypos,zpos
	variable npt
	
	variable countprint=0
	xpos=0;ypos=0;zpos=0
	variable xtst=enoise(100),ytst=enoise(100),ztst=enoise(100),addint=enoise(100)
	variable countB=0,counter=1,CountC=0,val1=0,val2=0
	do
		xtst=enoise(100);ytst=enoise(100);ztst=enoise(100);addint=enoise(100)
		addint=(addint/abs(addint))
		do
			//val1=counter
			if((xtst>ytst)&&(xtst>ztst))
				xpos[counter]=addint+xpos[counter-1]
				ypos[counter]=ypos[counter-1]
				zpos[counter]=zpos[counter-1]
				break
			endif
			if((ytst>xtst)&&(ytst>ztst))
				ypos[counter]=addint+ypos[counter-1]
				xpos[counter]=xpos[counter-1]
				zpos[counter]=zpos[counter-1]
				break
			endif
			if((ztst>ytst)&&(ztst>xtst))
				zpos[counter]=addint+zpos[counter-1]
				ypos[counter]=ypos[counter-1]
				xpos[counter]=xpos[counter-1]
				break
			endif
			break
		while(0)
		counter+=1
		countprint+=1
		if(countprint>300)
			print counter
			countprint=0
		endif
	while(counter<npt)
End

Function CalcSAW(xpos,ypos,zpos,npt)
	wave xpos,ypos,zpos
	variable npt
	
	variable countprint=0
	xpos=0;ypos=0;zpos=0
	variable xtst=enoise(100),ytst=enoise(100),ztst=enoise(100),addint=enoise(100)
	variable countB=0,counter=1,CountC=0,val1=0,val2=0
	do
		xtst=enoise(100);ytst=enoise(100);ztst=enoise(100);addint=enoise(100)
		addint=(addint/abs(addint))
		do
			val1=counter
			if((xtst>ytst)&&(xtst>ztst))
				xpos[counter]=addint+xpos[counter-1]
				ypos[counter]=ypos[counter-1]
				zpos[counter]=zpos[counter-1]
				do
					if(xpos[countB]==xpos[counter]&&ypos[countB]==ypos[counter]&&zpos[countB]==zpos[counter])
						xpos[counter]=xpos[counter-1]
						counter-=1
					endif
					countB+=1
				while(countB<counter)
				countB=0
				break
			endif
			if((ytst>xtst)&&(ytst>ztst))
				ypos[counter]=addint+ypos[counter-1]
				xpos[counter]=xpos[counter-1]
				zpos[counter]=zpos[counter-1]
				do
					if(xpos[countB]==xpos[counter]&&ypos[countB]==ypos[counter]&&zpos[countB]==zpos[counter])
						ypos[counter]=ypos[counter-1]
						counter-=1
					endif
					countB+=1
				while(countB<counter)
				countB=0
				break
			endif
			if((ztst>ytst)&&(ztst>xtst))
				zpos[counter]=addint+zpos[counter-1]
				ypos[counter]=ypos[counter-1]
				xpos[counter]=xpos[counter-1]
				do
					if(xpos[countB]==xpos[counter]&&ypos[countB]==ypos[counter]&&zpos[countB]==zpos[counter])
						zpos[counter]=zpos[counter-1]
						counter-=1
					endif
					countB+=1
				while(countB<counter)
				countB=0
				break
			endif
			break
		while(0)
		counter+=1
		val2=counter
		if(val1==val2)
			countC+=1
			if(countC>216)//then you have reached a deadend and have to back track to continue, say 100 back
					//216 is 6^3
				counter-=100
				countC=0
			endif
		else
			countC=0
		endif
		countprint+=1
		if(countprint>300)
			print counter
			countprint=0
		endif
	while(counter<npt)
End

Macro CalculateSAW(name,npt)
	string/g gname
	variable/g gnpt
	string name=gname
	variable npt=gnpt
	prompt name,"Enter the base name for x, y z waves"
	prompt npt,"Enter number of steps in chain (10,000)"
	
	silent 1
	gname=name;gnpt=npt
	variable counter=1 // npt=10000,
	string xpos="xpos_"+name,ypos="ypos_"+name,zpos="zpos_"+name
	make/n=(npt)/o $xpos=0, $ypos=0, $zpos=0
	variable xtst=enoise(100),ytst=enoise(100),ztst=enoise(100),addint=enoise(100)
	variable countB=0
	do
		xtst=enoise(100);ytst=enoise(100);ztst=enoise(100);addint=enoise(100)
		addint=(addint/abs(addint))
		do
			if((xtst>ytst)&&(xtst>ztst))
				$xpos[counter]=addint+$xpos[counter-1]
				$ypos[counter]=$ypos[counter-1]
				$zpos[counter]=$zpos[counter-1]
				do
					if($xpos[countB]==$xpos[counter]&&$ypos[countB]==$ypos[counter]&&$zpos[countB]==$zpos[counter])
						$xpos[counter]=$xpos[counter-1]
						counter-=1
					endif
					countB+=1
				while(countB<counter)
				countB=0
				break
			endif
			if((ytst>xtst)&&(ytst>ztst))
				$ypos[counter]=addint+$ypos[counter-1]
				$xpos[counter]=$xpos[counter-1]
				$zpos[counter]=$zpos[counter-1]
				do
					if($xpos[countB]==$xpos[counter]&&$ypos[countB]==$ypos[counter]&&$zpos[countB]==$zpos[counter])
						$ypos[counter]=$ypos[counter-1]
						counter-=1
					endif
					countB+=1
				while(countB<counter)
				countB=0
				break
			endif
			if((ztst>ytst)&&(ztst>xtst))
				$zpos[counter]=addint+$zpos[counter-1]
				$ypos[counter]=$ypos[counter-1]
				$xpos[counter]=$xpos[counter-1]
				do
					if($xpos[countB]==$xpos[counter]&&$ypos[countB]==$ypos[counter]&&$zpos[countB]==$zpos[counter])
						$zpos[counter]=$zpos[counter-1]
						counter-=1
					endif
					countB+=1
				while(countB<counter)
				countB=0
				break
			endif
			break
		while(0)
		counter+=1
		print counter
	while(counter<npt)
ENdMacro