--- imach/src/imach.c 2015/04/30 08:27:53 1.188 +++ imach/src/imach.c 2015/04/30 14:45:16 1.189 @@ -1,6 +1,9 @@ -/* $Id: imach.c,v 1.188 2015/04/30 08:27:53 brouard Exp $ +/* $Id: imach.c,v 1.189 2015/04/30 14:45:16 brouard Exp $ $State: Exp $ $Log: imach.c,v $ + Revision 1.189 2015/04/30 14:45:16 brouard + Summary: 0.98q2 + Revision 1.188 2015/04/30 08:27:53 brouard *** empty log message *** @@ -678,11 +681,11 @@ typedef struct { #define ODIRSEPARATOR '\\' #endif -/* $Id: imach.c,v 1.188 2015/04/30 08:27:53 brouard Exp $ */ +/* $Id: imach.c,v 1.189 2015/04/30 14:45:16 brouard Exp $ */ /* $State: Exp $ */ -char version[]="Imach version 0.98q1, April 2015,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121), Intel Software 2015"; -char fullversion[]="$Revision: 1.188 $ $Date: 2015/04/30 08:27:53 $"; +char version[]="Imach version 0.98q2, April 2015,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121), Intel Software 2015"; +char fullversion[]="$Revision: 1.189 $ $Date: 2015/04/30 14:45:16 $"; char strstart[80]; char optionfilext[10], optionfilefiname[FILENAMELENGTH]; int erreur=0, nberr=0, nbwarn=0; /* Error number, number of errors number of warnings */ @@ -1566,16 +1569,16 @@ void linmin(double p[], double xi[], int printf("retour brent fret=%.12e xmin=%.12e\n",*fret,xmin); fprintf(ficlog,"retour brent fret=%.12e xmin=%.12e\n",*fret,xmin); #endif - printf("linmin end "); + /* printf("linmin end "); */ for (j=1;j<=n;j++) { /* printf(" before xi[%d]=%12.8f", j,xi[j]); */ xi[j] *= xmin; /* xi rescaled by xmin: if xmin=-1.237 and xi=(1,0,...,0) xi=(-1.237,0,...,0) */ - if(xxs <1.0) - printf(" after xi[%d]=%12.8f, xmin=%12.8f, ax=%12.8f, xx=%12.8f, bx=%12.8f, xxs=%12.8f", j,xi[j], xmin, ax, xx, bx,xxs ); + /* if(xxs <1.0) */ + /* printf(" after xi[%d]=%12.8f, xmin=%12.8f, ax=%12.8f, xx=%12.8f, bx=%12.8f, xxs=%12.8f", j,xi[j], xmin, ax, xx, bx,xxs ); */ p[j] += xi[j]; /* Parameters values are updated accordingly */ } - printf("\n"); - printf("Comparing last *frec(xmin)=%12.8f from Brent and frec(0.)=%12.8f \n", (*func)(p)); + /* printf("\n"); */ + /* printf("Comparing last *frec(xmin)=%12.8f from Brent and frec(0.)=%12.8f \n", *fret, (*func)(p)); */ free_vector(xicom,1,n); free_vector(pcom,1,n); }