--- imach/src/imach.c	2017/07/20 13:35:01	1.279
+++ imach/src/imach.c	2018/02/21 07:58:13	1.280
@@ -1,6 +1,11 @@
-/* $Id: imach.c,v 1.279 2017/07/20 13:35:01 brouard Exp $
+/* $Id: imach.c,v 1.280 2018/02/21 07:58:13 brouard Exp $
   $State: Exp $
   $Log: imach.c,v $
+  Revision 1.280  2018/02/21 07:58:13  brouard
+  Summary: 0.99r15
+
+  New Makefile with recent VirtualBox 5.26. Bug in sqrt negatve in imach.c
+
   Revision 1.279  2017/07/20 13:35:01  brouard
   Summary: temporary working
 
@@ -1040,12 +1045,12 @@ typedef struct {
 #define ODIRSEPARATOR '\\'
 #endif
 
-/* $Id: imach.c,v 1.279 2017/07/20 13:35:01 brouard Exp $ */
+/* $Id: imach.c,v 1.280 2018/02/21 07:58:13 brouard Exp $ */
 /* $State: Exp $ */
 #include "version.h"
 char version[]=__IMACH_VERSION__;
 char copyright[]="February 2016,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121), Intel Software 2015-2018";
-char fullversion[]="$Revision: 1.279 $ $Date: 2017/07/20 13:35:01 $"; 
+char fullversion[]="$Revision: 1.280 $ $Date: 2018/02/21 07:58:13 $"; 
 char strstart[80];
 char optionfilext[10], optionfilefiname[FILENAMELENGTH];
 int erreur=0, nberr=0, nbwarn=0; /* Error number, number of errors number of warnings  */
@@ -6642,7 +6647,12 @@ To be simple, these graphs help to under
 		 }
 								
 		 /* Eigen vectors */
-		 v11=(1./sqrt(1+(v1-lc1)*(v1-lc1)/cv12/cv12));
+		 if(1+(v1-lc1)*(v1-lc1)/cv12/cv12 <1.e-5){
+		   printf(" Error sqrt of a negative number: %lf\n",1+(v1-lc1)*(v1-lc1)/cv12/cv12);
+		   fprintf(ficlog," Error sqrt of a negative number: %lf\n",1+(v1-lc1)*(v1-lc1)/cv12/cv12);
+		   v11=(1./sqrt(fabs(1+(v1-lc1)*(v1-lc1)/cv12/cv12)));
+		 }else
+		   v11=(1./sqrt(1+(v1-lc1)*(v1-lc1)/cv12/cv12));
 		 /*v21=sqrt(1.-v11*v11); *//* error */
 		 v21=(lc1-v1)/cv12*v11;
 		 v12=-v21;
@@ -6673,8 +6683,8 @@ To be simple, these graphs help to under
 		   fprintf(ficgp,"\nset label \"%d\" at %11.3e,%11.3e center",(int) age, mu1,mu2);
 		   fprintf(ficgp,"\n# Age %d, p%1d%1d - p%1d%1d",(int) age, k1,l1,k2,l2);
 		   fprintf(ficgp,"\nplot [-pi:pi] %11.3e+ %.3f*(%11.3e*%11.3e*cos(t)+%11.3e*%11.3e*sin(t)), %11.3e +%.3f*(%11.3e*%11.3e*cos(t)+%11.3e*%11.3e*sin(t)) not",	\
-			   mu1,std,v11,sqrt(lc1),v12,sqrt(fabs(lc2)),	\
-			   mu2,std,v21,sqrt(lc1),v22,sqrt(fabs(lc2))); /* For gnuplot only */
+			   mu1,std,v11,sqrt(fabs(lc1)),v12,sqrt(fabs(lc2)), \
+			   mu2,std,v21,sqrt(fabs(lc1)),v22,sqrt(fabs(lc2))); /* For gnuplot only */
 		 }else{
 		   first=0;
 		   fprintf(fichtmcov," %d (%.3f),",(int) age, c12);