gsl 1.15.3 → 1.16.0.6

data/{include → ext/gsl_native/include}/rb_gsl_histogram3d.h

@@ -3,11 +3,11 @@
 
 typedef struct {
   size_t nx, ny, nz;
-  double * xrange ;
-  double * yrange ;
-  double * zrange ;
-  double * bin ;
-} mygsl_histogram3d ;
+  double * xrange;
+  double * yrange;
+  double * zrange;
+  double * bin;
+} mygsl_histogram3d;
 
 typedef struct {
   gsl_histogram2d h;
@@ -15,56 +15,56 @@ typedef struct {
 
 
 mygsl_histogram3d* mygsl_histogram3d_alloc(const size_t nx, const size_t ny,
-					   const size_t nz);
+                                           const size_t nz);
 void mygsl_histogram3d_free (mygsl_histogram3d * h);
-mygsl_histogram3d* mygsl_histogram3d_calloc_uniform(const size_t nx, 
-						    const size_t ny,
-						    const size_t nz,
-						    const double xmin,
-						    const double xmax,
-						    const double ymin,
-						    const double ymax,
-						    const double zmin,
-						    const double zmax);
-mygsl_histogram3d* mygsl_histogram3d_calloc(const size_t nx, 
-					    const size_t ny,
-					    const size_t nz);
-int mygsl_histogram3d_set_ranges_uniform (mygsl_histogram3d * h, 
-					  double xmin, double xmax,
-					  double ymin, double ymax,
-					  double zmin, double zmax);
-int mygsl_histogram3d_set_ranges (mygsl_histogram3d * h, 
-				  const double xrange[], size_t xsize,
-				  const double yrange[], size_t ysize,
-				  const double zrange[], size_t zsize);
+mygsl_histogram3d* mygsl_histogram3d_calloc_uniform(const size_t nx,
+                                                    const size_t ny,
+                                                    const size_t nz,
+                                                    const double xmin,
+                                                    const double xmax,
+                                                    const double ymin,
+                                                    const double ymax,
+                                                    const double zmin,
+                                                    const double zmax);
+mygsl_histogram3d* mygsl_histogram3d_calloc(const size_t nx,
+                                            const size_t ny,
+                                            const size_t nz);
+int mygsl_histogram3d_set_ranges_uniform (mygsl_histogram3d * h,
+                                          double xmin, double xmax,
+                                          double ymin, double ymax,
+                                          double zmin, double zmax);
+int mygsl_histogram3d_set_ranges (mygsl_histogram3d * h,
+                                  const double xrange[], size_t xsize,
+                                  const double yrange[], size_t ysize,
+                                  const double zrange[], size_t zsize);
 int mygsl_histogram3d_memcpy(mygsl_histogram3d * dest, const mygsl_histogram3d * src);
 mygsl_histogram3d* mygsl_histogram3d_clone(const mygsl_histogram3d * src);
 int mygsl_histogram3d_fread(FILE * stream, mygsl_histogram3d * h);
 int mygsl_histogram3d_fwrite(FILE * stream, const mygsl_histogram3d * h);
 int mygsl_histogram3d_increment(mygsl_histogram3d * h, double x, double y, double z);
-int mygsl_histogram3d_accumulate (mygsl_histogram3d * h, 
-				  double x, double y, double z, double weight);
-int mygsl_histogram3d_increment2(mygsl_histogram3d * h, 
-				 double x, double y, double z);
-int mygsl_histogram3d_accumulate2(mygsl_histogram3d * h, 
-				  double x, double y, double z, double weight);
-double mygsl_histogram3d_get (const mygsl_histogram3d * h, const size_t i, 
-			      const size_t j, const size_t k);
-int mygsl_histogram3d_get_xrange(const mygsl_histogram3d * h, const size_t i,  
-				 double *xlower, double *xupper);
-int mygsl_histogram3d_get_yrange(const mygsl_histogram3d * h, const size_t j,  
-				 double *ylower, double *yupper);
-int mygsl_histogram3d_get_zrange(const mygsl_histogram3d * h, const size_t k,  
-				 double *zlower, double *zupper);
+int mygsl_histogram3d_accumulate (mygsl_histogram3d * h,
+                                  double x, double y, double z, double weight);
+int mygsl_histogram3d_increment2(mygsl_histogram3d * h,
+                                 double x, double y, double z);
+int mygsl_histogram3d_accumulate2(mygsl_histogram3d * h,
+                                  double x, double y, double z, double weight);
+double mygsl_histogram3d_get (const mygsl_histogram3d * h, const size_t i,
+                              const size_t j, const size_t k);
+int mygsl_histogram3d_get_xrange(const mygsl_histogram3d * h, const size_t i,
+                                 double *xlower, double *xupper);
+int mygsl_histogram3d_get_yrange(const mygsl_histogram3d * h, const size_t j,
+                                 double *ylower, double *yupper);
+int mygsl_histogram3d_get_zrange(const mygsl_histogram3d * h, const size_t k,
+                                 double *zlower, double *zupper);
 int mygsl_histogram3d_find (const mygsl_histogram3d * h,
-			    const double x, const double y, const double z,
-			    size_t * i, size_t * j, size_t *k);
+                            const double x, const double y, const double z,
+                            size_t * i, size_t * j, size_t *k);
 gsl_histogram2d* mygsl_histogram3d_xyproject(const mygsl_histogram3d * h3,
-					     size_t kstart, size_t kend);
+                                             size_t kstart, size_t kend);
 gsl_histogram2d* mygsl_histogram3d_xzproject(const mygsl_histogram3d * h3,
-					     size_t jstart, size_t jend);
+                                             size_t jstart, size_t jend);
 gsl_histogram2d* mygsl_histogram3d_yzproject(const mygsl_histogram3d * h3,
-					     size_t istart, size_t iend);
+                                             size_t istart, size_t iend);
 int mygsl_histogram3d_scale(mygsl_histogram3d * h, double scale);
 int mygsl_histogram3d_shift(mygsl_histogram3d * h, double shift);
 double mygsl_histogram3d_xmax(const mygsl_histogram3d * h);
@@ -74,11 +74,11 @@ double mygsl_histogram3d_ymin(const mygsl_histogram3d * h);
 double mygsl_histogram3d_zmax(const mygsl_histogram3d * h);
 double mygsl_histogram3d_zmin(const mygsl_histogram3d * h);
 double mygsl_histogram3d_max_val(const mygsl_histogram3d * h);
-void mygsl_histogram3d_max_bin(const mygsl_histogram3d * h, 
-			       size_t *imax_out, size_t *jmax_out, size_t *kmax_out);
+void mygsl_histogram3d_max_bin(const mygsl_histogram3d * h,
+                               size_t *imax_out, size_t *jmax_out, size_t *kmax_out);
 double mygsl_histogram3d_min_val(const mygsl_histogram3d * h);
-void mygsl_histogram3d_min_bin(const mygsl_histogram3d * h, 
-			       size_t *imin_out, size_t *jmin_out, size_t *kmin_out);
+void mygsl_histogram3d_min_bin(const mygsl_histogram3d * h,
+                               size_t *imin_out, size_t *jmin_out, size_t *kmin_out);
 double mygsl_histogram3d_sum (const mygsl_histogram3d * h);
 double mygsl_histogram3d_xmean (const mygsl_histogram3d * h);
 double mygsl_histogram3d_ymean (const mygsl_histogram3d * h);
@@ -88,7 +88,7 @@ double mygsl_histogram3d_ysigma(const mygsl_histogram3d * h);
 double mygsl_histogram3d_zsigma(const mygsl_histogram3d * h);
 void mygsl_histogram3d_reset(mygsl_histogram3d * h);
 int mygsl_histogram3d_equal_bins_p(const mygsl_histogram3d * h1,
-				   const mygsl_histogram3d * h2);
+                                   const mygsl_histogram3d * h2);
 int mygsl_histogram3d_add(mygsl_histogram3d * h1, const mygsl_histogram3d * h2);
 int mygsl_histogram3d_sub(mygsl_histogram3d * h1, const mygsl_histogram3d * h2);
 int mygsl_histogram3d_mul(mygsl_histogram3d * h1, const mygsl_histogram3d * h2);

data/{include → ext/gsl_native/include}/rb_gsl_integration.h

@@ -12,6 +12,6 @@
 #ifndef ___RB_GSL_INTEGRATION_H___
 #define ___RB_GSL_INTEGRATION_H___
 
-#include "gsl/gsl_integration.h"
+#include <gsl/gsl_integration.h>
 
 #endif

data/{include → ext/gsl_native/include}/rb_gsl_interp.h

@@ -13,14 +13,9 @@
 #ifndef ___RB_GSL_INTERP_H___
 #define ___RB_GSL_INTERP_H___
 
-#include "ruby.h"
 #include "rb_gsl.h"
 #include <gsl/gsl_interp.h>
 #include <gsl/gsl_spline.h>
-#ifdef HAVE_NARRAY_H
-#include "narray.h"
-#endif
-
 
 typedef struct {
   gsl_interp *p;

data/{include → ext/gsl_native/include}/rb_gsl_linalg.h

@@ -12,12 +12,8 @@
 #ifndef ___RB_GSL_LINALG_H___
 #define ___RB_GSL_LINALG_H___
 
-#include "gsl/gsl_linalg.h"
-#include "gsl/gsl_math.h"
-
-#ifdef HAVE_NARRAY_H
-#include "rb_gsl_with_narray.h"
-#endif
+#include <gsl/gsl_linalg.h>
+#include <gsl/gsl_math.h>
 
 VALUE rb_gsl_linalg_complex_LU_decomp(int argc, VALUE *argv, VALUE obj);
 VALUE rb_gsl_linalg_complex_LU_decomp2(int argc, VALUE *argv, VALUE obj);

data/{include → ext/gsl_native/include}/rb_gsl_math.h

@@ -13,14 +13,8 @@
 #ifndef ___RB_GSL_MATH_H___
 #define ___RB_GSL_MATH_H___
 
-#include "ruby.h"
-#include <gsl/gsl_math.h>
 #include "rb_gsl.h"
 
-#ifndef GSL_1_3_LATER
-int gsl_fcmp (const double x1, const double x2, const double epsilon);
-#endif
-
 VALUE rb_gsl_math_complex_eval(gsl_complex (*func)(gsl_complex), VALUE obj);
 
 #endif

data/{include → ext/gsl_native/include}/rb_gsl_odeiv.h

@@ -13,9 +13,6 @@
 #ifndef ___RB_GSL_ODEIV_H___
 #define ___RB_GSL_ODEIV_H___
 
-#include "ruby.h"
-#include <gsl/gsl_odeiv.h>
 #include "rb_gsl.h"
-#include "rb_gsl_array.h"
 
 #endif

data/{include → ext/gsl_native/include}/rb_gsl_poly.h

@@ -13,11 +13,7 @@
 #ifndef ___RB_GSL_POLY_H___
 #define ___RB_GSL_POLY_H___
 
-#include "ruby.h"
-#include <gsl/gsl_vector.h>
 #include <gsl/gsl_poly.h>
-#include <gsl/gsl_complex.h>
-#include <gsl/gsl_complex_math.h>
 #include "rb_gsl_complex.h"
 #include "rb_gsl_array.h"
 
@@ -31,7 +27,7 @@ EXTERN VALUE cgsl_rational;
 typedef gsl_vector gsl_poly;
 typedef gsl_vector_int gsl_poly_int;
 /*
-typedef struct ___gsl_rational 
+typedef struct ___gsl_rational
 {
   VALUE num, den;
   gsl_poly *pnum;
@@ -40,7 +36,7 @@ typedef struct ___gsl_rational
 */
 
 int gsl_poly_conv(const double *a, size_t na, const double *b, size_t nb,
-		  double *c, size_t *nc);
+                  double *c, size_t *nc);
 
 gsl_vector* gsl_poly_deconv_vector(const gsl_vector *c, const gsl_vector *a, gsl_vector **r);
 gsl_vector* gsl_poly_deriv(const gsl_vector *v);
@@ -55,7 +51,7 @@ VALUE rb_gsl_poly_complex_solve2(int argc, VALUE *argv, VALUE obj);
 
 
 int gsl_poly_int_conv(const int *a, size_t na, const int *b, size_t nb,
-		  int *c, size_t *nc);
+                      int *c, size_t *nc);
 
 gsl_vector_int* gsl_poly_int_deconv_vector(const gsl_vector_int *c, const gsl_vector_int *a, gsl_vector_int **r);
 gsl_vector_int* gsl_poly_int_deriv(const gsl_vector_int *v);

data/{include → ext/gsl_native/include}/rb_gsl_rational.h

@@ -13,16 +13,9 @@
 #ifndef ___RB_GSL_RATIONAL_H___
 #define ___RB_GSL_RATIONAL_H___
 
-#include "ruby.h"
-#include <gsl/gsl_vector.h>
-#include <gsl/gsl_poly.h>
-#include <gsl/gsl_complex.h>
-#include <gsl/gsl_complex_math.h>
-#include "rb_gsl_complex.h"
-#include "rb_gsl_array.h"
 #include "rb_gsl_poly.h"
 
-typedef struct ___gsl_rational 
+typedef struct ___gsl_rational
 {
   VALUE num, den;
   gsl_poly *pnum;

data/{include → ext/gsl_native/include}/rb_gsl_rng.h

@@ -12,7 +12,6 @@
 #ifndef ___RB_GSL_RNG_H___
 #define ___RB_GSL_RNG_H___
 
-#include "ruby.h"
 #include <gsl/gsl_rng.h>
 #include "rb_gsl.h"
 

data/{include → ext/gsl_native/include}/rb_gsl_root.h

@@ -13,9 +13,9 @@
 #define ___RB_GSL_ROOT_H___
 
 #include <gsl/gsl_errno.h>
-#include <gsl/gsl_math.h>
 #include <gsl/gsl_roots.h>
 #include "rb_gsl.h"
+
 EXTERN VALUE cgsl_fsolver;
 EXTERN VALUE cgsl_fdfsolver;
 

data/{include → ext/gsl_native/include}/rb_gsl_sf.h

@@ -11,13 +11,9 @@
 #ifndef ___RB_GSL_SF_H___
 #define ___RB_GSL_SF_H___
 
-#include "ruby.h"
 #include <gsl/gsl_sf.h>
-#include <gsl/gsl_errno.h>
-#include <gsl/gsl_vector.h>
-#include <gsl/gsl_math.h>
 #include <gsl/gsl_complex.h>
-#include "rb_gsl_config.h"
+#include <gsl/gsl_sf_mathieu.h>
 #include "rb_gsl.h"
 
 EXTERN VALUE cgsl_sf_result, cgsl_sf_result_e10;
@@ -28,56 +24,56 @@ VALUE rb_gsl_sf_eval1(double (*func)(double), VALUE argv);
 VALUE rb_gsl_sf_eval_int_double(double (*func)(int, double), VALUE jj, VALUE argv);
 VALUE rb_gsl_sf_eval_double_double(double (*func)(double, double), VALUE ff, VALUE argv);
 VALUE rb_gsl_sf_eval1_uint(double (*func)(unsigned int), VALUE argv);
-VALUE rb_gsl_sf_eval_double_m(double (*func)(double, gsl_mode_t), VALUE argv, VALUE m);;
-VALUE rb_gsl_sf_eval_int_double_double(double (*func)(int, double, double), VALUE jj, 
-				       VALUE ff, VALUE argv);
+VALUE rb_gsl_sf_eval_double_m(double (*func)(double, gsl_mode_t), VALUE argv, VALUE m);
+VALUE rb_gsl_sf_eval_int_double_double(double (*func)(int, double, double), VALUE jj,
+                                       VALUE ff, VALUE argv);
 VALUE rb_gsl_sf_eval_int_int_double(double (*func)(int, int, double), VALUE jj,
-				    VALUE jj2, VALUE argv);
-VALUE rb_gsl_sf_eval_double3(double (*func)(double, double, double), 
-					   VALUE ff, VALUE ff2, VALUE argv);
-VALUE rb_gsl_sf_eval_double4(double (*func)(double, double, double, double), 
-			     VALUE ff, VALUE ff2, VALUE ff3, VALUE argv);
+                                    VALUE jj2, VALUE argv);
+VALUE rb_gsl_sf_eval_double3(double (*func)(double, double, double),
+                             VALUE ff, VALUE ff2, VALUE argv);
+VALUE rb_gsl_sf_eval_double4(double (*func)(double, double, double, double),
+                             VALUE ff, VALUE ff2, VALUE ff3, VALUE argv);
 VALUE rb_gsl_sf_eval_double_int(double (*func)(double, int), VALUE argv, VALUE jj);
 VALUE rb_gsl_sf_eval1_int(double (*func)(int), VALUE argv);
 VALUE rb_gsl_sf_eval_e(int (*func)(double, gsl_sf_result*), VALUE x);
-VALUE rb_gsl_sf_eval_e_int_double(int (*func)(int, double, gsl_sf_result*), 
-				  VALUE n, VALUE x);
-VALUE rb_gsl_sf_eval_e_int_int_double(int (*func)(int, int, double, gsl_sf_result*), 
-				      VALUE n1, VALUE n2, VALUE x);
-VALUE rb_gsl_sf_eval_e_double2(int (*func)(double, double, gsl_sf_result*), 
-			       VALUE x1, VALUE x2);
+VALUE rb_gsl_sf_eval_e_int_double(int (*func)(int, double, gsl_sf_result*),
+                                  VALUE n, VALUE x);
+VALUE rb_gsl_sf_eval_e_int_int_double(int (*func)(int, int, double, gsl_sf_result*),
+                                      VALUE n1, VALUE n2, VALUE x);
+VALUE rb_gsl_sf_eval_e_double2(int (*func)(double, double, gsl_sf_result*),
+                               VALUE x1, VALUE x2);
 VALUE rb_gsl_sf_eval_e_uint(int (*func)(unsigned int, gsl_sf_result*), VALUE x);
-VALUE rb_gsl_sf_eval_e_int_uint(int (*func)(int, unsigned int, gsl_sf_result*), 
-				VALUE n, VALUE x);
-VALUE rb_gsl_sf_eval_e_double_uint(int (*func)(double, unsigned int, gsl_sf_result*), 
-				   VALUE y, VALUE x);
-VALUE rb_gsl_sf_eval_e_m(int (*func)(double, gsl_mode_t, gsl_sf_result*), 
-			 VALUE x, VALUE m);
-VALUE rb_gsl_sf_eval_e_double2_m(int (*func)(double, double, gsl_mode_t, gsl_sf_result*), 
-				 VALUE x1, VALUE x2, VALUE m);
-VALUE rb_gsl_sf_eval_e_double3_m(int (*func)(double, double, double, gsl_mode_t, gsl_sf_result*), 
-				 VALUE x1, VALUE x2, VALUE x3, VALUE m);
-VALUE rb_gsl_sf_eval_double2_m(double (*func)(double, double, gsl_mode_t), 
-			       VALUE argv, VALUE x2, VALUE m);
-VALUE rb_gsl_sf_eval_double3_m(double (*func)(double, double, double, gsl_mode_t), 
-			       VALUE argv, VALUE x2, VALUE x3, VALUE m);
+VALUE rb_gsl_sf_eval_e_int_uint(int (*func)(int, unsigned int, gsl_sf_result*),
+                                VALUE n, VALUE x);
+VALUE rb_gsl_sf_eval_e_double_uint(int (*func)(double, unsigned int, gsl_sf_result*),
+                                   VALUE y, VALUE x);
+VALUE rb_gsl_sf_eval_e_m(int (*func)(double, gsl_mode_t, gsl_sf_result*),
+                         VALUE x, VALUE m);
+VALUE rb_gsl_sf_eval_e_double2_m(int (*func)(double, double, gsl_mode_t, gsl_sf_result*),
+                                 VALUE x1, VALUE x2, VALUE m);
+VALUE rb_gsl_sf_eval_e_double3_m(int (*func)(double, double, double, gsl_mode_t, gsl_sf_result*),
+                                 VALUE x1, VALUE x2, VALUE x3, VALUE m);
+VALUE rb_gsl_sf_eval_double2_m(double (*func)(double, double, gsl_mode_t),
+                               VALUE argv, VALUE x2, VALUE m);
+VALUE rb_gsl_sf_eval_double3_m(double (*func)(double, double, double, gsl_mode_t),
+                               VALUE argv, VALUE x2, VALUE x3, VALUE m);
 
-VALUE rb_gsl_sf_eval_e_double4_m(int (*func)(double, double, double, double, gsl_mode_t, gsl_sf_result*), 
-				 VALUE x1, VALUE x2, VALUE x3, VALUE x4, VALUE m);
+VALUE rb_gsl_sf_eval_e_double4_m(int (*func)(double, double, double, double, gsl_mode_t, gsl_sf_result*),
+                                 VALUE x1, VALUE x2, VALUE x3, VALUE x4, VALUE m);
 
 VALUE rb_gsl_sf_eval_e_int(int (*func)(int, gsl_sf_result*), VALUE x);
-VALUE rb_gsl_sf_eval_e_double3(int (*func)(double, double, double, gsl_sf_result*), 
-			      VALUE x1, VALUE x2, VALUE x3);
-VALUE rb_gsl_sf_eval_e_int_double2(int (*func)(int, double, double, gsl_sf_result*), 
-				   VALUE n, VALUE x1, VALUE x2);
-VALUE rb_gsl_sf_eval_e_double2(int (*func)(double, double, gsl_sf_result*), 
-			       VALUE x1, VALUE x2);
+VALUE rb_gsl_sf_eval_e_double3(int (*func)(double, double, double, gsl_sf_result*),
+                               VALUE x1, VALUE x2, VALUE x3);
+VALUE rb_gsl_sf_eval_e_int_double2(int (*func)(int, double, double, gsl_sf_result*),
+                                   VALUE n, VALUE x1, VALUE x2);
+VALUE rb_gsl_sf_eval_e_double2(int (*func)(double, double, gsl_sf_result*),
+                               VALUE x1, VALUE x2);
 
 VALUE eval_sf(double (*func)(double, gsl_mode_t), VALUE argv);
 
 VALUE rb_gsl_sf_eval_double4_m(double (*func)(double, double, double, double,
-					      gsl_mode_t), 
-			       VALUE argv, VALUE x2, VALUE x3, VALUE x4, VALUE m);
+                                              gsl_mode_t),
+                               VALUE argv, VALUE x2, VALUE x3, VALUE x4, VALUE m);
 
 VALUE rb_gsl_sf_eval_complex(double (*f)(double), VALUE obj);
 
@@ -109,11 +105,6 @@ void Init_gsl_sf_synchrotron(VALUE module);
 void Init_gsl_sf_transport(VALUE module);
 void Init_gsl_sf_trigonometric(VALUE module);
 void Init_gsl_sf_zeta(VALUE module);
-
-
-#ifdef GSL_1_9_LATER
-#include <gsl/gsl_sf_mathieu.h>
 void Init_sf_mathieu(VALUE module);
-#endif
 
 #endif

data/{include → ext/gsl_native/include}/rb_gsl_statistics.h

@@ -12,6 +12,6 @@
 #ifndef ___RB_GSL_STATISTICS_H___
 #define ___RB_GSL_STATISTICS_H___
 
-#include "gsl/gsl_statistics.h"
+#include <gsl/gsl_statistics.h>
 
 #endif

data/{include → ext/gsl_native/include}/rb_gsl_tensor.h

@@ -3,8 +3,6 @@
 
 #ifdef HAVE_TENSOR_TENSOR_H
 #include "rb_gsl.h"
-#include "rb_gsl_common.h"
-//#include <tensor/tensor.h>
 #include <tensor/tensor.h>
 
 EXTERN VALUE cgsl_tensor, cgsl_tensor_int;

data/{include → ext/gsl_native/include}/rb_gsl_with_narray.h

@@ -1,10 +1,11 @@
 #ifndef RB_GSL_WITH_NARRAY_H
 #define RB_GSL_WITH_NARRAY_H
 
+#include <gsl/gsl_vector.h>
+
 #ifdef HAVE_NARRAY_H
 
 #include "narray.h"
-#include "gsl/gsl_vector.h"
 
 gsl_vector* make_cvector_from_narray(VALUE);
 void cvector_set_from_narray(gsl_vector*, VALUE);
@@ -25,5 +26,6 @@ extern VALUE cNVector, cNMatrix;
 
 gsl_vector_complex* na_to_gv_complex(VALUE na);
 gsl_vector_complex_view* na_to_gv_complex_view(VALUE na);
+
 #endif // HAVE_NARRAY_H
 #endif // RB_GSL_WITH_NARRAY_H

data/{include → ext/gsl_native/include}/templates_on.h

@@ -173,7 +173,7 @@
 #error unknown BASE_ directive in source.h
 #endif
 
-#define CONCAT2x(a,b) a ## _ ## b 
+#define CONCAT2x(a,b) a ## _ ## b
 #define CONCAT2(a,b) CONCAT2x(a,b)
 #define CONCAT3x(a,b,c) a ## _ ## b ## _ ## c
 #define CONCAT3(a,b,c) CONCAT3x(a,b,c)