/***********************************************
	Test Routine (Rosenblock function)
	f(x0,x1) = 100*(x1-x0^2)^2 + (1-x0)^2
	xinit = (-1.2, 1.0)
	xopt = (1, 1)   fopt = 0
************************************************/

#include <stdio.h>
#include <math.h>
#include "opt.h"

static double Rosenblock(int n, double x[])
{
	double	x0, x1;
	
	x0 = x[0];
	x1 = x[1];
	return 100*(x1-x0*x0)*(x1-x0*x0) + (1-x0)*(1-x0);
}

static double *grad_Rosenblock(int n, double x[], double grad[])
{
	double	x0, x1;
	
	x0 = x[0];
	x1 = x[1];
	grad[0] = 200*(x1-x0*x0)*(-2*x0) + 2*(1-x0)*(-1);
	grad[1] = 200*(x1-x0*x0);
	return grad;
}

main()
{
	int	n = 2;
	double	x[2];
	double	fopt;
	int	timeout = 1000;
	double	eps = 1.0e-20;
	
	x[0] = -1.20;
	x[1] = 1.00;
	
	if (QuasiNewtonMethod(n, Rosenblock, grad_Rosenblock, x, &fopt, timeout, eps) == success) {
		printf("converge ");
	} else {
		printf("timeout  ");
	}
	printf(" [ %lf %lf ]  %lf\n", x[0], x[1], fopt);
}