Merge branch 'master' of jugit.fz-juelich.de:mlz/kww

lib/kww.c

@@ -116,31 +116,6 @@ double kwwp_lim_hig( const double b )
 }
 
 
-/*****************************************************************************/
-/*  Auxiliary functions for debugging                                        */
-/*****************************************************************************/
-
-void kww_hexprint_double( const char* name, const double x )
-{
-    union {
-        double d;
-        unsigned long u;
-    } u;
-    u.d = x;
-    printf( "  %s = %24.18f [%lx]\n", name, x, u.u );
-}
-
-void kww_hexprint_quad( const char* name, const long double x )
-{
-    union {
-        long double d;
-        unsigned long long u;
-    } u;
-    u.d = x;
-    printf( "  %s = %24.18Lf [%llx]", name, x, u.u );
-}
-
-
 /*****************************************************************************/
 /*  High-level wrapper functions                                             */
 /*****************************************************************************/
@@ -379,12 +354,6 @@ double kww_hig( const double w, const double beta,
     long double s;        // full term (with trigonometric factor)
     long double x, gl;    // local variables
 
-    if( kww_debug & 8 ) {
-        printf( "kww_hig: kappa %i mu %i deb %i\n", kappa, mu, kww_debug );
-        kww_hexprint_double( "b", beta );
-        kww_hexprint_double( "w", w );
-    }
-
     // set diagnostic variable
     kww_algorithm = 3;
 
@@ -581,11 +550,6 @@ double kww_mid( const double w, const double beta,
     else
         N = 40;
 
-    if( kww_debug & 8 ) {
-        printf( "kww_mid %i %i %i %i\n", kind, mu, N, kww_debug );
-        kww_hexprint_double( "b", beta );
-        kww_hexprint_double( "w", w );
-    }
     for ( iter=0; iter<max_iter_int; ++iter ) {
         // static initialisation of NN, ak, bk for given 'iter'
         if ( iter>iterDone[kind][j] ) {
@@ -658,11 +622,6 @@ double kww_mid( const double w, const double beta,
         kww_num_of_terms += 2*NN[j][iter]+1;
         if ( diffmode )
             S += w/sqrt(PI)/2*exp(-SQR(w)/4);
-        if( kww_debug & 8 ) {
-            printf( "iter %i N %i\n", iter, N );
-            kww_hexprint_double( "S_new", S );
-            kww_hexprint_double( "S_old", S_last );
-        }
         // termination criteria
         if      ( kww_debug & 4 )
             return -1; // we want to inspect just one sum