/* note #undef's at end of file */
#define NRANSI
#include <iostream>
#include <iostream.h>
#include "nrutil.h"
#define TOL 2.0e-4

int ncom;
float *pcom,*xicom,(*nrfunc)(float []);

void linmin(float p[], float xi[], int n, float *fret, float (*func)(float []))
{
	float brent (float ax, float bx, float cx,
		     float (*f)(float), float tol, float *xmin);
	void  mnbrak(float *ax, float *bx, float *cx, float *fa, float *fb,
		     float *fc, float (*func)(float));
	float f1dim (float x);	
	int j;
	float xx,xmin,fx,fb,fa,bx,ax;

	ncom=n;
	pcom=nrvector(0,n-1);
	xicom=nrvector(0,n-1);
	nrfunc=func;
	for (j=0;j<n;j++) {
	  pcom[j]=p[j];
	  xicom[j]=xi[j];
	}
	ax=0.0;
	xx=1.0;
	mnbrak(&ax,&xx,&bx,&fa,&fx,&fb,f1dim);
	*fret=brent(ax,xx,bx,f1dim,TOL,&xmin);
	for (j=0;j<n;j++) {
	  xi[j] *= xmin;
	  p[j] += xi[j];
	}
	//cout << "xmin="<< xmin<< "\n";
	//cout << "linmin: theta=";
	//cout << "linmin 0 : p=";
	//for(j=0;j<n;j++){cout << p[j] << " ";}cout << "\n";	

	
	free_nrvector(xicom,0,n-1);
	free_nrvector(pcom,0,n-1);
}
#undef TOL
#undef NRANSI