testOptimize_prg Program

program testOptimize_prg
	!! Test program for Optimize_mod
	!! @todo
	!! Finish tests
	use kinds_mod
	use optimize_mod
	use objective_mod
	implicit none
	
	call testObjective
	call testObjectiveN
	call testPlot
	
contains

	subroutine testObjective
		!! Verify operation of obj_t
		type(test_t)::test
		
		write(*,*) test%eval(2.0_wp)-3.0_wp
		write(*,*) test%der1(0.0_wp)
		write(*,*) test%der2(0.0_wp)
		write(*,*) test%rootNewton(4.0_wp,tol=1.0E-10_wp,maxIts=1000000)
		write(*,*) test%minNewton(4.0_wp,tol=1.0E-10_wp,maxIts=1000000)
	end subroutine testObjective

	subroutine testObjectiveN
		type(testN_t)::test
		real(wp),dimension(:),allocatable::xSD,xMN,xNM
		
		xSD = test%steepestDescent([5.0_wp,5.0_wp])
		xMN = test%minNewton([5.0_wp,5.0_wp])
		xNM = test%nelderMead([5.0_wp,5.0_wp])
		write(*,*) xSD
		write(*,*) xMN
		write(*,*) xNM
	end subroutine testObjectiveN

	subroutine testPlot
		type(testN_t)::test
		real(wp),dimension(:),allocatable::x,y
		real(wp),dimension(:,:),allocatable::f
		integer::N,i,j
		
		real(wp),dimension(:),allocatable::xm
		
		N = 100
		x = linspace(0.0_wp,7.0_wp,N)
		y = linspace(0.0_wp,7.0_wp,N)
		allocate( f(N,N) )
		
		do j=1,N
			do i=1,N
				f(i,j) = test%eval([ x(i) , y(j) ])
			end do
		end do
		
		call setup(fileName='testsOptimize-%n.svg',figSize=[400,350])
		
		kLog = 1
		xm = test%steepestDescent([5.0_wp,5.0_wp])
		call figure()
		call subplot(1,1,1,aspect=span(y)/span(x))
		call xylim(mixval(x),mixval(y))
		call contourf(x,y,f,30)
		call plot(xLog(:kLog-1,1),xLog(:kLog-1,2),lineStyle='-',markStyle='x',markColor='k')
		call colorbar2(f,30)
		call ticks()
		call labels('x','y','Steepest Descent ['//intToChar(kLog-1)//']')
		
		kLog = 1
		xm = test%minNewton([5.0_wp,5.0_wp])
		call figure()
		call subplot(1,1,1,aspect=span(y)/span(x))
		call xylim(mixval(x),mixval(y))
		call contourf(x,y,f,30)
		call plot(xLog(:kLog-1,1),xLog(:kLog-1,2),lineStyle='-',markStyle='x',markColor='k')
		call colorbar2(f,30)
		call ticks()
		call labels('x','y','Newton-Raphson ['//intToChar(kLog-1)//']')
		
		kLog = 1
		xm = test%nelderMead([5.0_wp,5.0_wp])
		call figure()
		call subplot(1,1,1,aspect=span(y)/span(x))
		call xylim(mixval(x),mixval(y))
		call contourf(x,y,f,30)
		call plot(xLog(:kLog-1,1),xLog(:kLog-1,2),lineStyle='-',markStyle='x',markColor='k')
		call colorbar2(f,30)
		call ticks()
		call labels('x','y','Nelder-Mead Simplex ['//intToChar(kLog-1)//']')
		
		call show()
		
	end subroutine testPlot

end program testOptimize_prg