#pragma rtGlobals=1		// Use modern global access method.

Menu "Dynamics"
	"Calc Correlation Function Panel"
	"Calculate Position Panel"
End

Function NoTImeSteps(ctrlName,varNum,varStr,varName) : SetVariableControl
	String ctrlName
	Variable varNum
	String varStr
	String varName
	
	nvar gnosteps,gtimemax
	Make/n=(gnosteps)/o qxpos_time=p*gtimemax/(gnosteps-1),rxpos_time
End

Function MaxTime(ctrlName,varNum,varStr,varName) : SetVariableControl
	String ctrlName
	Variable varNum
	String varStr
	String varName

	nvar gnosteps,gtimemax
	Make/n=(gnosteps)/o qxpos_time=p*gtimemax/(gnosteps-1),rxpos_time
End

Function CalculateXpos(ctrlName) : ButtonControl
	String ctrlName

	nvar gnosteps,gtimemax,gneda,gks,gDiffCoeff
	Make/n=(gnosteps)/o qxpos_time=p*gtimemax/(gnosteps-1),rxpos_time
	Rxpos_time[0]=0
	variable counter=0,tau=gneda/gks
	//First Calculate the flutuations for all time steps
	Make/n=(gnosteps)/o g_of_t
	do
		g_of_t[counter]=gDiffcoeff*enoise(1)
		counter+=1
	while(counter<gnosteps)
	counter=1
	variable counter2=0
	//then calculate the xpos by integral for each point
	do
		Rxpos_time[counter]=0
		do
			//use memory function and add up the history of g(t) pulses
			Rxpos_time[counter]+=g_of_t[counter2]*exp(-(counter-counter2)/tau)
			counter2+=1
		while(counter2<counter)
		counter2=0
		counter+=1
	while(counter<gnosteps)
End

Function CalculateCorrelationFUn(ctrlName) : ButtonControl
	String ctrlName

	nvar gnosteps
	wave rxpos_time,qxpos_time
	make/n=(gnosteps)/o R_TimeCorrFun
	
	variable counter=1
	variable counter2=0,counter3=0
	//then calculate the xpos by integral for each point
	R_TimeCorrFun[0]=1
	do
		R_TimeCorrFun[counter]=0
		do
			R_TimeCorrFun[counter]+=rxpos_time[counter2]*rxpos_time[counter2+counter]
			counter2+=1
		while((counter2+counter)<gnosteps)
		//R_TimeCorrFun[counter]/=counter2
		R_TimeCorrFun[counter]/=(gnosteps-counter)//from http://www.nyu.edu/classes/tuckerman/stat.mech/lectures/lecture_21/node4.html
		counter2=0
		counter+=1
	while(counter<gnosteps)
	//do
	//	R_TimeCorrFun[counter]=0
	//	do
	//		//use memory function and add up the history of g(t) pulses
	//		do
	//			R_TimeCorrFun[counter]+=rxpos_time[counter3]*rxpos_time[counter2]
	//			counter3+=1
	//		while(counter3<gnosteps)
	//		counter3=0
	//		counter2+=1
	//	while(counter2<counter)
	//	counter2=0
	//	counter+=1
	//while(counter<gnosteps)
End

Window CalculatePositionPanel() : Graph
	PauseUpdate; Silent 1		// building window...
		
	variable/g gks,gneda,gdiffcoeff,gTimeMax,gnosteps
	make/n=(gnosteps)/o R_TimeCorrFun
	make/n=(gnosteps)/o Rxpos_time, qxpos_time

	Display /W=(327.75,48.5,756,398) rxpos_time vs qxpos_time
	ModifyGraph grid=2
	ModifyGraph tick=2
	ModifyGraph zero(left)=1
	ModifyGraph mirror=1
	ModifyGraph font="Times New Roman"
	ModifyGraph minor=1
	ModifyGraph fSize=15
	ModifyGraph fStyle=1
	ModifyGraph standoff=0
	Label left "\\F'Times'\\Z18\\f01x Position"
	Label bottom "\\F'Times'\\Z18\\f01Time"
	ShowInfo
	Legend/N=text0/J "\\s(Rxpos_time) Rxpos_time"
	ShowTools
	SetVariable SpringConst,pos={8,307},size={150,16},title="Spring Constant"
	SetVariable SpringConst,value= gks
	SetVariable FricFactor,pos={8,328},size={150,16},title="Friction Factor"
	SetVariable FricFactor,value= gneda
	SetVariable DiffusionCoef,pos={8,349},size={100,16},title="D",value= gdiffcoeff
	SetVariable TimeMax,pos={8,370},size={150,16},proc=MaxTime,title="Maximum Time"
	SetVariable TimeMax,value= gtimemax
	SetVariable NumberofTimeSteps,pos={8,392},size={180,16},proc=NoTImeSteps,title="Number of Time Steps"
	SetVariable NumberofTimeSteps,value= gnosteps
	Button CalculateXposvsTime,pos={10,415},size={130,20},proc=CalculateXpos,title="Calculate X Position"
EndMacro

Window CalcCorrelationFunctionPanel() : Graph
	PauseUpdate; Silent 1		// building window...
		
	variable/g gks,gneda,gdiffcoeff,gTimeMax,gnosteps
	make/n=(gnosteps)/o R_TimeCorrFun
	make/n=(gnosteps)/o Rxpos_time, qxpos_time

	Display /W=(5.25,41.75,338.25,388.25) R_TimeCorrFun vs qxpos_time
	ModifyGraph mode=4
	ModifyGraph gaps=0
	ModifyGraph grid=2
	ModifyGraph tick=2
	ModifyGraph zero(left)=1
	ModifyGraph mirror=1
	ModifyGraph font="Times New Roman"
	ModifyGraph minor=1
	ModifyGraph fSize=15
	ModifyGraph fStyle=1
	ModifyGraph standoff=0
	Label left "\\F'Times'\\Z18\\f01Time Correlation Function"
	Label bottom "\\F'Times'\\Z18\\f01Time"
	SetAxis left 0,2e-05
	SetAxis bottom 0,1000000000
	ShowInfo
	Legend/N=text0/J "\\s(R_TimeCorrFun) R_TimeCorrFun"
	ShowTools
	Button CalculateCorrelationFunc,pos={12,398},size={190,20},proc=CalculateCorrelationFUn,title="Calculate Correlation Function"
EndMacro