#include "SteepestDescent.h"
main(){
	Vec x, y, z, gradf,p, VecListPointer[3];
	float f,x1=1.1,x2=1.2, alpha=100.0;
	int i,n=2,test;
	
	
	
	x.n=n;
	x.comp=(float *) malloc(n*sizeof(float));
	x.comp[0]=x1;
	x.comp[1]=x2;
	
	y.n=n;
	y.comp=(float *) malloc(n*sizeof(float));
	y.comp[0]=1.2;
	y.comp[1]=123.4;
	
	z.n=n;
	z.comp=(float *) malloc(n*sizeof(float));
	z.comp[0]=12.9;
	z.comp[1]=1293.4;
	
		
	VecListPointer[0]=x;
	VecListPointer[1]=y;
	VecListPointer[2]=z;
	f = F(x);
	printf("f=%f\n",f);
	gradf =GradF(x);
	p.n = n;
	p.comp=(float *) malloc(n*sizeof(float));
	for(i=0;i<n; i++){
		p.comp[i]=-gradf.comp[i];
		}
	printf("Dimension=%i and gradf={%f,%f}\n",gradf.n, gradf.comp[0],gradf.comp[1]);
	test=WolfeTest(
	x, 
	alpha, 
	p, 
	f, 
	gradf,
	F,
	GradF
	); 
	printf("test=%i\n",test);
	
	PrintScalar(VecListPointer, 3, F);
}