diff --git a/contrib/bessel/Bessel.cpp b/contrib/bessel/Bessel.cpp
index f6c4a4a16f4c13ef20040d8bc51866b77ada36ea..0e01946ff63adfb4716204192788dd5e0a2528fd 100644
--- a/contrib/bessel/Bessel.cpp
+++ b/contrib/bessel/Bessel.cpp
@@ -11,31 +11,31 @@
 static void zbesj_(double*, double*, double*, int*, int*, double*,
     double*, int*, int*)
 {
-  std::cerr << "Bessel functions require Fortran compiler";
+  std::cerr << "Bessel functions require Fortran compiler" << std::endl;
 }
 
 static void zbesk_(double*, double*, double*, int*, int*, double*,
     double*, int*, int*)
 {
-  std::cerr << "Bessel functions require Fortran compiler";
+  std::cerr << "Bessel functions require Fortran compiler" << std::endl;
 }
 
 static void zbesy_(double*, double*, double*, int*, int*, double*,
     double*, int*, double*, double*, int*)
 {
-  std::cerr << "Bessel functions require Fortran compiler";
+  std::cerr << "Bessel functions require Fortran compiler" << std::endl;
 }
 
 static void zbesh_(double*, double*, double*, int*, int*, int*,
     double*, double*, int*, int*)
 {
-  std::cerr << "Bessel functions require Fortran compiler";
+  std::cerr << "Bessel functions require Fortran compiler" << std::endl;
 }
 
 static void zairy_(double*, double*, int*, int*, double*,
     double*, int*, int*)
 {
-  std::cerr << "Bessel functions require Fortran compiler";
+  std::cerr << "Bessel functions require Fortran compiler" << std::endl;
 }
 
 #else
@@ -59,7 +59,7 @@ extern "C" {
       double*, int*);
   void zbesh_(double*, double*, double*, int*, int*, int*, double*,
       double*, int*, int*);
-  
+
   void zairy_(double*, double*, int*, int*, double*,
       double*, int*, int*);
 }
@@ -74,24 +74,24 @@ static int BesselError(int ierr, const char *str)
     case 0 :
       return 0;
     case 1 :
-      std::cerr << "Input error in " << str;
+      std::cerr << "Input error in " << str << std::endl;
       return BESSEL_ERROR_INPUT;
     case 2 :
       return BESSEL_OVERFLOW;
     case 3 :
       if(!warn){
-        std::cout << "Half machine accuracy lost in " << str << " (large argument or order)";
+        std::cout << "Half machine accuracy lost in " << str << " (large argument or order)" << std::endl;
         warn = 1;
       }
       return BESSEL_HALF_ACCURACY;
     case 4 :
-      std::cerr << "Complete loss of significance in " << str << " (argument or order too large)";
+      std::cerr << "Complete loss of significance in " << str << " (argument or order too large)" << std::endl;
       return BESSEL_NO_ACCURACY;
     case 5 :
-      std::cerr << "Failed to converge in " << str;
+      std::cerr << "Failed to converge in " << str << std::endl;
       return BESSEL_NO_CONVERGENCE;
     default:
-      std::cout << "Unknown Bessel status in " << str << " (" << ierr << ")";
+      std::cout << "Unknown Bessel status in " << str << " (" << ierr << ")" << std::endl;
       return ierr;
   }
 }
@@ -243,7 +243,7 @@ int BesselAltSphericalHn(int type, double n, int num, double x, std::complex<dou
   return BesselError(ierr, "BesselAltSphericalHn");
 }
 
-// Airy Function and derivative (id=1) 
+// Airy Function and derivative (id=1)
 
 int AiryComplex(double xr, double xi, int id, double *valr, double *vali)
 {
@@ -251,7 +251,7 @@ int AiryComplex(double xr, double xi, int id, double *valr, double *vali)
 
   zairy_(&xr, &xi, &id, &kode, valr, vali, &nz, &ierr);
   //val = std::complex<double>(valr, vali);
-  
+
   return BesselError(ierr, "AiryComplex");
 }