fix compilation warningsHEADmaster

12 files changed, 24 insertions, 19 deletions

diff --git a/src/c/specialFunctions/besseli/dbesselia.c b/src/c/specialFunctions/besseli/dbesselia.c
index 14ac1ed3..46cedf80 100644
--- a/src/c/specialFunctions/besseli/dbesselia.c
+++ b/src/c/specialFunctions/besseli/dbesselia.c
@@ -26,8 +26,7 @@ double onebessi( double inp1, double inp2)
 
    if (inp1 < 0)
    {
-      double   dblank;
-      return( dblank );
+      return 0.0/0.0;
    }
    if (inp1 == 0)
       {
diff --git a/src/c/specialFunctions/besseli/sbesselia.c b/src/c/specialFunctions/besseli/sbesselia.c
index 6e185e0b..55ec43b6 100644
--- a/src/c/specialFunctions/besseli/sbesselia.c
+++ b/src/c/specialFunctions/besseli/sbesselia.c
@@ -26,8 +26,7 @@ float fonebessi( float inp1, float inp2)
 
    if (inp1 < 0)
    {
-      float   dblank;
-      return( dblank );
+      return 0.0/0.0;
    }
    if (inp1 == 0)
       {
diff --git a/src/c/specialFunctions/besselj/dbesselja.c b/src/c/specialFunctions/besselj/dbesselja.c
index 23a355f5..8054bd53 100644
--- a/src/c/specialFunctions/besselj/dbesselja.c
+++ b/src/c/specialFunctions/besselj/dbesselja.c
@@ -81,9 +81,7 @@ double bessj( int n, double x )
 
    if (n < 0)
    {
-      double   dblank;
-      //setdblack_c( &dblank );
-      return( dblank );
+      return 0.0/0.0;
    }
    ax=fabs(x);
    if (n == 0)
diff --git a/src/c/specialFunctions/besselj/sbesselja.c b/src/c/specialFunctions/besselj/sbesselja.c
index 68d2ea4f..57e3d63c 100644
--- a/src/c/specialFunctions/besselj/sbesselja.c
+++ b/src/c/specialFunctions/besselj/sbesselja.c
@@ -81,8 +81,7 @@ float fbessj( int n, float x )
 
    if (n < 0)
    {
-      float   dblank;
-      return( dblank );
+      return 0.0/0.0;
    }
    ax=fabs(x);
    if (n == 0)
diff --git a/src/c/specialFunctions/besselk/dbesselka.c b/src/c/specialFunctions/besselk/dbesselka.c
index df6c0708..1f01ee63 100644
--- a/src/c/specialFunctions/besselk/dbesselka.c
+++ b/src/c/specialFunctions/besselk/dbesselka.c
@@ -109,8 +109,7 @@ double bessk( int n, double x )
 
    if (n < 0 || x == 0.0)
    {
-      double   dblank;
-      return( dblank );
+      return 0.0/0.0;
    }
    if (n == 0)
       return( bessk0(x) );
diff --git a/src/c/specialFunctions/besselk/sbesselka.c b/src/c/specialFunctions/besselk/sbesselka.c
index bfb35eb1..56d3207b 100644
--- a/src/c/specialFunctions/besselk/sbesselka.c
+++ b/src/c/specialFunctions/besselk/sbesselka.c
@@ -109,8 +109,7 @@ float fbessk( int n, float x )
 
    if (n < 0 || x == 0.0)
    {
-      float   dblank;
-      return( dblank );
+      return 0.0/0.0;
    }
    if (n == 0)
       return( fbessk0(x) );
diff --git a/src/c/specialFunctions/bessely/dbesselya.c b/src/c/specialFunctions/bessely/dbesselya.c
index 76678629..14ea8583 100644
--- a/src/c/specialFunctions/bessely/dbesselya.c
+++ b/src/c/specialFunctions/bessely/dbesselya.c
@@ -135,8 +135,7 @@ double bessy( int n, double x )
 
    if (n < 0 || x == 0.0)
    {
-      double   dblank;
-      return( dblank );
+      return 0.0/0.0;
    }
    if (n == 0)
       return( bessy0(x) );
diff --git a/src/c/specialFunctions/bessely/sbesselya.c b/src/c/specialFunctions/bessely/sbesselya.c
index 1b42736c..2e8b08d8 100644
--- a/src/c/specialFunctions/bessely/sbesselya.c
+++ b/src/c/specialFunctions/bessely/sbesselya.c
@@ -136,8 +136,7 @@ float fbessy( int n, float x )
 
    if (n < 0 || x == 0.0)
    {
-      float   dblank;
-      return( dblank );
+      return 0.0/0.0;
    }
    if (n == 0)
       return( fbessy0(x) );
diff --git a/src/c/specialFunctions/calerf/dcalerfs.c b/src/c/specialFunctions/calerf/dcalerfs.c
index c3608653..3b2bfd7a 100644
--- a/src/c/specialFunctions/calerf/dcalerfs.c
+++ b/src/c/specialFunctions/calerf/dcalerfs.c
@@ -25,5 +25,6 @@ double dcalerfs (double inp1, double inp2)
 	return derfcs(inp1);
     if (inp2 == 2)
 	return derfcxs(inp1);
+    return 0.0/0.0;
 }
 
diff --git a/src/c/specialFunctions/calerf/scalerfs.c b/src/c/specialFunctions/calerf/scalerfs.c
index ffc48c37..20de0c34 100644
--- a/src/c/specialFunctions/calerf/scalerfs.c
+++ b/src/c/specialFunctions/calerf/scalerfs.c
@@ -25,5 +25,6 @@ float scalerfs (float inp1, double inp2)
 	return derfcs(inp1);
     if (inp2 == 2)
 	return derfcxs(inp1);
+    return 0.0/0.0;
 }
 
diff --git a/src/c/specialFunctions/erfinv/derfinvs.c b/src/c/specialFunctions/erfinv/derfinvs.c
index 3f7fe761..2dce2051 100644
--- a/src/c/specialFunctions/erfinv/derfinvs.c
+++ b/src/c/specialFunctions/erfinv/derfinvs.c
@@ -28,7 +28,8 @@ double derfinvs (double inp1)
     {
 	return 0.0/0.0;				// returns Nan
     }    
-    if ((inp1 >= -0.7) && (inp1 <= 0.7))
+
+    else if ((inp1 >= -0.7) && (inp1 <= 0.7))
     {
         double sq = inp1 * inp1;
         return (inp1 * (((a[3]*sq+a[2]) * sq+a[1]) * sq+a[0]) /  ((((b[3]*sq+b[2]) * sq+b[1]) * sq+b[0]) * sq+1));	//Inverse error formula
@@ -45,5 +46,10 @@ double derfinvs (double inp1)
         double z = sqrt(-log((1+inp1)/2));
         return (-(((c[3]*z+c[2]) * z+c[1]) * z+c[0]) / ((d[1]*z+d[0]) * z+1));		//Inverse error formula
     }
+
+    else
+    {
+	return 0.0/0.0;				// returns Nan
+    }
 }
 
diff --git a/src/c/specialFunctions/erfinv/serfinvs.c b/src/c/specialFunctions/erfinv/serfinvs.c
index 5436fe39..55ba532f 100644
--- a/src/c/specialFunctions/erfinv/serfinvs.c
+++ b/src/c/specialFunctions/erfinv/serfinvs.c
@@ -29,7 +29,8 @@ float serfinvs (float inp1)
 	double na = 0.0/0.0;	
 	return (float)na;
     }    
-    if ((inp1 >= -0.7) && (inp1 <= 0.7))
+
+    else if ((inp1 >= -0.7) && (inp1 <= 0.7))
     {
         float sq = inp1 * inp1;
         return (inp1 * (((a[3]*sq+a[2]) * sq+a[1]) * sq+a[0]) /  ((((b[3]*sq+b[2]) * sq+b[1]) * sq+b[0]) * sq+1));	//Inverse error formula
@@ -46,5 +47,10 @@ float serfinvs (float inp1)
         float z = sqrt(-log((1+inp1)/2));
         return (-(((c[3]*z+c[2]) * z+c[1]) * z+c[0]) / ((d[1]*z+d[0]) * z+1));		//Inverse error formula
     }
+
+    else
+    {
+	return 0.0/0.0;
+    }
 }