rb-gsl 1.16.0.4 → 1.16.0.5

data/ext/gsl/cqp.c

@@ -1,283 +0,0 @@
-#ifdef HAVE_GSL_GSL_CQP_H
-
-#include "include/rb_gsl.h"
-#include "gsl/gsl_cqp.h"
-
-static const gsl_cqpminimizer_type* type_by_string(VALUE t);
-static const gsl_cqpminimizer_type* get_type(VALUE t)
-{
-
-	switch (TYPE(t)) {
-	case T_STRING:
-		return type_by_string(t);
-		break;
-	default:
-		rb_raise(rb_eTypeError, "Wrong argument type %s.", rb_class2name(CLASS_OF(t)));
-	}
-}
-
-static const gsl_cqpminimizer_type* type_by_string(VALUE t)
-{
-	char *name;	
-	name = STR2CSTR(t);
-	if (strcmp(name, "mg_pdip") == 0) {
-		return gsl_cqpminimizer_mg_pdip;
-	} else {
-		rb_raise(rb_eRuntimeError, "Unknown minimizer type %s.", name);
-	}
-	return NULL; /* never reach here */
-}
-
-static VALUE rb_cqpminimizer_alloc(VALUE klass,  VALUE t, VALUE n, VALUE me, VALUE mi)
-{
-	gsl_cqpminimizer *m;
-	m = gsl_cqpminimizer_alloc(get_type(t), (size_t) FIX2INT(n), (size_t) FIX2INT(me), (size_t) FIX2INT(mi));
-	return Data_Wrap_Struct(klass, 0, gsl_cqpminimizer_free, m);
-}
-
-static VALUE rb_cqpminimizer_set(VALUE obj, VALUE data)
-{
-	gsl_cqpminimizer *m;
-	gsl_cqp_data *d;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	Data_Get_Struct(data, gsl_cqp_data, d);
-	gsl_cqpminimizer_set(m, d);
-	return Qtrue;
-}
-
-static VALUE rb_cqpminimizer_name(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return rb_str_new2(gsl_cqpminimizer_name(m));
-}
-
-static VALUE rb_cqpminimizer_iterate(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return INT2FIX(gsl_cqpminimizer_iterate(m));
-}
-
-static VALUE rb_cqpminimizer_x(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return Data_Wrap_Struct(cgsl_vector_view, 0, NULL, gsl_cqpminimizer_x(m));
-}
-
-static VALUE rb_cqpminimizer_lm_eq(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return Data_Wrap_Struct(cgsl_vector_view, 0, NULL, gsl_cqpminimizer_lm_eq(m));
-}
-static VALUE rb_cqpminimizer_lm_ineq(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return Data_Wrap_Struct(cgsl_vector_view, 0, NULL, gsl_cqpminimizer_lm_ineq(m));
-}
-static VALUE rb_cqpminimizer_f(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return rb_float_new(gsl_cqpminimizer_f(m));
-}
-static VALUE rb_cqpminimizer_gap(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return rb_float_new(gsl_cqpminimizer_gap(m));
-}
-static VALUE rb_cqpminimizer_residuals_norm(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return rb_float_new(gsl_cqpminimizer_residuals_norm(m));
-}
-/*
-static VALUE rb_cqpminimizer_minimum(VALUE obj)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return rb_float_new(gsl_cqpminimizer_minimum(m));
-}
-*/
-static VALUE rb_cqpminimizer_test_convergence(VALUE obj, VALUE g, VALUE r)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return INT2FIX(gsl_cqpminimizer_test_convergence(m, NUM2DBL(g), NUM2DBL(r)));
-}
-static VALUE rb_cqpminimizer_test_infeasibility(VALUE obj, VALUE e)
-{
-	gsl_cqpminimizer *m;
-	Data_Get_Struct(obj, gsl_cqpminimizer, m);
-	return INT2FIX(gsl_cqp_minimizer_test_infeasibility(m, NUM2DBL(e)));
-}
-
-static VALUE rb_cqp_data_alloc(VALUE klass)
-{
-	gsl_cqp_data *d;
-	d = (gsl_cqp_data*) malloc(sizeof(gsl_cqp_data));
-	return Data_Wrap_Struct(klass, 0, free, d);
-}
-
-static VALUE rb_cqp_data_Q(VALUE obj)
-{
-	gsl_cqp_data *d;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	return Data_Wrap_Struct(cgsl_matrix_view, 0, NULL, 	d->Q);
-}
-
-static VALUE rb_cqp_data_q(VALUE obj)
-{
-	gsl_cqp_data *d;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	return Data_Wrap_Struct(cgsl_vector_view, 0, NULL, 	d->q);
-}
-
-static VALUE rb_cqp_data_A(VALUE obj)
-{
-	gsl_cqp_data *d;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	return Data_Wrap_Struct(cgsl_matrix_view, 0, NULL, 	d->A);
-}
-
-static VALUE rb_cqp_data_b(VALUE obj)
-{
-	gsl_cqp_data *d;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	return Data_Wrap_Struct(cgsl_vector_view, 0, NULL, 	d->b);
-}
-
-static VALUE rb_cqp_data_C(VALUE obj)
-{
-	gsl_cqp_data *d;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	return Data_Wrap_Struct(cgsl_matrix_view, 0, NULL, 	d->C);
-}
-
-static VALUE rb_cqp_data_d(VALUE obj)
-{
-	gsl_cqp_data *d;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	return Data_Wrap_Struct(cgsl_vector_view, 0, NULL, 	d->d);
-}
-
-static VALUE rb_cqp_data_set_Q(VALUE obj, VALUE mm)
-{
-	gsl_cqp_data *d;
-	gsl_matrix *m;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	CHECK_MATRIX(mm);
-	Data_Get_Struct(mm, gsl_matrix, m);	
-	d->Q = m;
-	return Qtrue;
-}
-
-static VALUE rb_cqp_data_set_q(VALUE obj, VALUE vv)
-{
-	gsl_cqp_data *d;
-	gsl_vector *v;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	CHECK_VECTOR(vv);
-	Data_Get_Struct(vv, gsl_vector, v);
-	d->q = v;
-	return Qtrue;
-}
-
-static VALUE rb_cqp_data_set_A(VALUE obj, VALUE mm)
-{
-	gsl_cqp_data *d;
-	gsl_matrix *m;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	CHECK_MATRIX(mm);
-	Data_Get_Struct(mm, gsl_matrix, m);
-	d->A = m;
-	return Qtrue;
-}
-
-static VALUE rb_cqp_data_set_b(VALUE obj, VALUE vv)
-{
-	gsl_cqp_data *d;
-	gsl_vector *v;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	CHECK_VECTOR(vv);
-	Data_Get_Struct(vv, gsl_vector, v);
-	d->b = v;
-	return Qtrue;
-}
-
-static VALUE rb_cqp_data_set_C(VALUE obj, VALUE mm)
-{
-	gsl_cqp_data *d;
-	gsl_matrix *m;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	CHECK_MATRIX(mm);
-	Data_Get_Struct(mm, gsl_matrix, m);
-	d->C = m;
-	return Qtrue;
-}
-
-static VALUE rb_cqp_data_set_d(VALUE obj, VALUE vv)
-{
-	gsl_cqp_data *d;
-	gsl_vector *v;
-	Data_Get_Struct(obj, gsl_cqp_data, d);
-	CHECK_VECTOR(vv);
-	Data_Get_Struct(vv, gsl_vector, v);
-	d->d = v;
-	return Qtrue;
-}
-
-void Init_cqp(VALUE module)
-{
-	VALUE mCQP, cMinimizer, cData;
-	
-	mCQP = rb_define_module_under(module, "CQP");
-	cMinimizer = rb_define_class_under(mCQP, "Minimizer", cGSL_Object);
-	cData = rb_define_class_under(mCQP, "Data", cGSL_Object);	
-
-	rb_define_singleton_method(cMinimizer, "alloc", rb_cqpminimizer_alloc, 4);
-	
-	rb_define_method(cMinimizer, "set", rb_cqpminimizer_set, 1);
-	rb_define_method(cMinimizer, "name", rb_cqpminimizer_name, 0);
-	rb_define_method(cMinimizer, "iterate", rb_cqpminimizer_iterate, 0);	
-	rb_define_method(cMinimizer, "x", rb_cqpminimizer_x, 0);
-	rb_define_method(cMinimizer, "lm_eq", rb_cqpminimizer_lm_eq, 0);
-	rb_define_method(cMinimizer, "lm_ineq", rb_cqpminimizer_lm_ineq, 0);
-	rb_define_method(cMinimizer, "f", rb_cqpminimizer_f, 0);
-	rb_define_method(cMinimizer, "gap", rb_cqpminimizer_gap, 0);
-	rb_define_method(cMinimizer, "residuals_norm", rb_cqpminimizer_residuals_norm, 0);	
-/*	rb_define_method(cMinimizer, "minimum", rb_cqpminimizer_minimum, 0);	*/
-	rb_define_method(cMinimizer, "test_convergence", rb_cqpminimizer_test_convergence, 2);
-	rb_define_method(cMinimizer, "test_infeasibility", rb_cqpminimizer_test_infeasibility, 1);	
-
-	/*****/	
-	rb_define_singleton_method(cData, "alloc", rb_cqp_data_alloc, 0);	
-	
-	rb_define_method(cData, "Q", rb_cqp_data_Q, 0);
-	rb_define_method(cData, "q", rb_cqp_data_q, 0);
-	rb_define_method(cData, "A", rb_cqp_data_A, 0);
-	rb_define_method(cData, "b", rb_cqp_data_b, 0);
-	rb_define_method(cData, "C", rb_cqp_data_C, 0);
-	rb_define_method(cData, "d", rb_cqp_data_d, 0);
-	
-	rb_define_method(cData, "set_Q", rb_cqp_data_set_Q, 1);
-	rb_define_method(cData, "set_q", rb_cqp_data_set_q, 1);
-	rb_define_method(cData, "set_A", rb_cqp_data_set_A, 1);
-	rb_define_method(cData, "set_b", rb_cqp_data_set_b, 1);
-	rb_define_method(cData, "set_C", rb_cqp_data_set_C, 1);
-	rb_define_method(cData, "set_d", rb_cqp_data_set_d, 1);
-	rb_define_alias(cData, "Q=", "set_Q");
-	rb_define_alias(cData, "q=", "set_q");
-	rb_define_alias(cData, "A=", "set_A");
-	rb_define_alias(cData, "b=", "set_b");
-	rb_define_alias(cData, "C=", "set_C");
-	rb_define_alias(cData, "d=", "set_d");	
-}
-
-#endif
-

data/ext/gsl/fcmp.c

@@ -1,64 +0,0 @@
-/* sys/gsl_compare.c
- * 
- * Copyright (C) 2002 Gert Van den Eynde
- * 
- * This program is free software; you can redistribute it and/or modify
- * it under the terms of the GNU General Public License as published by
- * the Free Software Foundation; either version 2 of the License, or (at
- * your option) any later version.
- * 
- * This program is distributed in the hope that it will be useful, but
- * WITHOUT ANY WARRANTY; without even the implied warranty of
- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
- * General Public License for more details.
- * 
- * You should have received a copy of the GNU General Public License
- * along with this program; if not, write to the Free Software
- * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
- * 
- * Based on fcmp 1.2.2 Copyright (c) 1998-2000 Theodore C. Belding
- * University of Michigan Center for the Study of Complex Systems
- * Ted.Belding@umich.edu
- *
- */
-
-#ifndef GSL_1_3_LATER
-
-#include <config.h>
-#include <gsl/gsl_sys.h>
-#include <math.h>
-
-int gsl_fcmp(const double x1, const double x2, const double epsilon)
-{
-  int exponent;
-  double delta, difference;
-
-  /* Find exponent of largest absolute value */
-
-  {
-    double max = (fabs (x1) > fabs (x2)) ? x1 : x2;
-
-    frexp (max, &exponent);
-  }
-
-  /* Form a neighborhood of size  2 * delta */
-
-  delta = ldexp (epsilon, exponent);
-
-  difference = x1 - x2;
-
-  if (difference > delta)       /* x1 > x2 */
-    {
-      return 1;
-    }
-  else if (difference < -delta) /* x1 < x2 */
-    {
-      return -1;
-    }
-  else                          /* -delta <= difference <= delta */
-    {
-      return 0;                 /* x1 ~=~ x2 */
-    }
-}
-
-#endif

data/ext/gsl/fresnel.c

@@ -1,312 +0,0 @@
-/****************************************************************************
- *              fresnel.c -
- *  Calculation of Fresnel integrals by expansion to Chebyshev series
- *  Expansions are taken from the book
- *  Y.L. Luke. Mathematical functions and their approximations.
- *  Moscow, "Mir", 1980. PP. 145-149 (Russian edition)
- ****************************************************************************
- */
-/*
-	Modified for Ruby bindings
-					2006/Dec/24 Y. TSUNESADA
-*/
-#include <math.h>
-#include "include/rb_gsl.h"
-#include "include/rb_gsl_sf.h"
-
-static const double sqrt_pi_2   = 1.2533141373155002512078826424; /* sqrt(pi/2) */
-static const double sqrt_2_pi   = 0.7978845608028653558798921199; /* sqrt(2/pi) */
-static const double _1_sqrt_2pi = 0.3989422804014326779399460599; /* 1/sqrt(2*pi) */
-static const double pi_2        = 1.5707963267948966192313216916; /* pi/2 */
-
-static double f_data_a[18] =
-		{
-		  0.76435138664186000189,
-		 -0.43135547547660179313,
-		  0.43288199979726653054,
-                 -0.26973310338387111029,
-                  0.08416045320876935378,
-                 -0.01546524484461381958,
-                  0.00187855423439822018,
-                 -0.00016264977618887547,
-                  0.00001057397656383260,
-                 -0.00000053609339889243,
-                  0.00000002181658454933,
-                 -0.00000000072901621186,
-                  0.00000000002037332546,
-                 -0.00000000000048344033,
-                  0.00000000000000986533,
-                 -0.00000000000000017502,
-                  0.00000000000000000272,
-                 -0.00000000000000000004
-                };
-
-static double f_data_b[17] =
-		{
-		  0.63041404314570539241,
-		 -0.42344511405705333544,
-		  0.37617172643343656625,
-		 -0.16249489154509567415,
-		  0.03822255778633008694,
-		 -0.00564563477132190899,
-		  0.00057454951976897367,
-		 -0.00004287071532102004,
-		  0.00000245120749923299,
-		 -0.00000011098841840868,
-		  0.00000000408249731696,
-		 -0.00000000012449830219,
-		  0.00000000000320048425,
-		 -0.00000000000007032416,
-		  0.00000000000000133638,
-		 -0.00000000000000002219,
-		  0.00000000000000000032
-		};
-
-static double fresnel_cos_0_8(double x)
-{
- double x_8 = x/8.0;
- double xx = 2.0*x_8*x_8 - 1.0;
-
- double t0 = 1.0;
- double t1 = xx;
- double sumC = f_data_a[0] + f_data_a[1]*t1;
- double t2;
- int n;
- for (n=2; n < 18; n++)
- {
-  t2 = 2.0*xx*t1 - t0;
-  sumC += f_data_a[n]*t2;
-  t0 = t1; t1 = t2;
- }
- return _1_sqrt_2pi*sqrt(x)*sumC;
-}
-
-static double fresnel_sin_0_8(double x)
-{
- double x_8 = x/8.0;
- double xx = 2.0*x_8*x_8 - 1.0;
- double t0 = 1.;
- double t1 = xx;
- double ot1 = x_8;
- double ot2 = 2.0*x_8*t1 - ot1;
- double sumS = f_data_b[0]*ot1 + f_data_b[1]*ot2;
- int n;
- double t2;
- for (n=2; n < 17; n++)
- {
-  t2 = 2.0*xx*t1 - t0;
-  ot1 = ot2;
-  ot2 = 2.0*x_8*t2 - ot1;
-  sumS += f_data_b[n]*ot2;
-  t0 = t1; t1 = t2;
- }
- return _1_sqrt_2pi*sqrt(x)*sumS;
-}
-
-static double f_data_e[41] =
-		{
-		    0.97462779093296822410,
-		   -0.02424701873969321371,
-		    0.00103400906842977317,
-		   -0.00008052450246908016,
-		    0.00000905962481966582,
-		   -0.00000131016996757743,
-		    0.00000022770820391497,
-		   -0.00000004558623552026,
-		    0.00000001021567537083,
-		   -0.00000000251114508133,
-		    0.00000000066704761275,
-		   -0.00000000018931512852,
-		    0.00000000005689898935,
-		   -0.00000000001798219359,
-		    0.00000000000594162963,
-		   -0.00000000000204285065,
-		    0.00000000000072797580,
-		   -0.00000000000026797428,
-		    0.00000000000010160694,
-		   -0.00000000000003958559,
-		    0.00000000000001581262,
-		   -0.00000000000000646411,
-		    0.00000000000000269981,
-		   -0.00000000000000115038,
-		    0.00000000000000049942,
-		   -0.00000000000000022064,
-		    0.00000000000000009910,
-		   -0.00000000000000004520,
-		    0.00000000000000002092,
-		   -0.00000000000000000982,
-		    0.00000000000000000467,
-		   -0.00000000000000000225,
-		    0.00000000000000000110,
-		   -0.00000000000000000054,
-		    0.00000000000000000027,
-		   -0.00000000000000000014,
-		    0.00000000000000000007,
-		   -0.00000000000000000004,
-		    0.00000000000000000002,
-		   -0.00000000000000000001,
-		    0.00000000000000000001
-        };
-
-static double f_data_f[35] =
-		{
-		    0.99461545179407928910,
-		   -0.00524276766084297210,
-		    0.00013325864229883909,
-		   -0.00000770856452642713,
-		    0.00000070848077032045,
-		   -0.00000008812517411602,
-		    0.00000001359784717148,
-		   -0.00000000246858295747,
-		    0.00000000050925789921,
-		   -0.00000000011653400634,
-		    0.00000000002906578309,
-		   -0.00000000000779847361,
-		    0.00000000000222802542,
-		   -0.00000000000067239338,
-		    0.00000000000021296411,
-		   -0.00000000000007041482,
-		    0.00000000000002419805,
-		   -0.00000000000000861080,
-		    0.00000000000000316287,
-		   -0.00000000000000119596,
-		    0.00000000000000046444,
-		   -0.00000000000000018485,
-		    0.00000000000000007527,
-		   -0.00000000000000003131,
-		    0.00000000000000001328,
-		   -0.00000000000000000574,
-		    0.00000000000000000252,
-		   -0.00000000000000000113,
-		    0.00000000000000000051,
-		   -0.00000000000000000024,
-		    0.00000000000000000011,
-		   -0.00000000000000000005,
-		    0.00000000000000000002,
-		   -0.00000000000000000001,
-		    0.00000000000000000001
-		};
-
-static double fresnel_cos_8_inf(double x)
-{
- double xx = 128.0/(x*x) - 1.0;   /* 2.0*(8/x)^2 - 1 */
- double t0 = 1.0;
- double t1 = xx;
- double sumP = f_data_e[0] + f_data_e[1]*t1;
- double sumQ = f_data_f[0] + f_data_f[1]*t1;
- double t2;
- int n;
- for(n = 2; n < 35; n++)
- {
-   t2 = 2.0*xx*t1 - t0;
-   sumP += f_data_e[n]*t2; /*  sumP += f_data_e[n]*ChebyshevT(n,xx) */
-   sumQ += f_data_f[n]*t2; /*  sumQ += f_data_f[n]*ChebyshevT(n,xx) */
-   t0 = t1; t1 = t2;
- }
- for(n = 35; n < 41; n++)
- {
-   t2 = 2.0*xx*t1 - t0;
-   sumP += f_data_e[n]*t2; /*  sumP += f_data_e[n]*ChebyshevT(n,xx) */
-   t0 = t1; t1 = t2;
- }
- return 0.5 - _1_sqrt_2pi*(0.5*sumP*cos(x)/x - sumQ*sin(x))/sqrt(x);
-}
-
-static double fresnel_sin_8_inf(double x)
-{
- double xx = 128.0/(x*x) - 1.0;   /* 2.0*(8/x)^2 - 1 */
- double t0 = 1.0;
- double t1 = xx;
- double sumP = f_data_e[0] + f_data_e[1]*t1;
- double sumQ = f_data_f[0] + f_data_f[1]*t1;
- double t2;
- int n;
- for(n = 2; n < 35; n++)
- {
-   t2 = 2.0*xx*t1 - t0;
-   sumP += f_data_e[n]*t2; /*  sumP += f_data_e[n]*ChebyshevT(n,xx) */
-   sumQ += f_data_f[n]*t2; /*  sumQ += f_data_f[n]*ChebyshevT(n,xx) */
-   t0 = t1; t1 = t2;
- }
- for(n = 35; n < 41; n++)
- {
-   t2 = 2.0*xx*t1 - t0;
-   sumP += f_data_e[n]*t2; /*  sumQ += f_data_f[n]*ChebyshevT(n,xx) */
-   t0 = t1; t1 = t2;
- }
- return 0.5 - _1_sqrt_2pi*(0.5*sumP*sin(x)/x + sumQ*cos(x))/sqrt(x);
-}
-
-
-double fresnel_c(double x)
-{
-  double xx = x*x*pi_2;
-  double ret_val;
-  if(xx<=8.0)
-   ret_val = fresnel_cos_0_8(xx);
-  else
-   ret_val = fresnel_cos_8_inf(xx);
-  return (x<0.0) ? -ret_val : ret_val;
-}
-
-double fresnel_s(double x)
-{
-  double xx = x*x*pi_2;
-  double ret_val;
-  if(xx<=8.0)
-   ret_val = fresnel_sin_0_8(xx);
-  else
-   ret_val = fresnel_sin_8_inf(xx);
-  return (x<0.0) ? -ret_val : ret_val;
-}
-
-double fresnel_c1(double x)
-{
-  return fresnel_c(x*sqrt_2_pi);
-}
-
-double fresnel_s1(double x)
-{
-  return fresnel_s(x*sqrt_2_pi);
-}
-
-static VALUE rb_fresnel_c(VALUE obj, VALUE x)
-{
-	return rb_gsl_sf_eval1(fresnel_c, x);
-}
-static VALUE rb_fresnel_s(VALUE obj, VALUE x)
-{
-	return rb_gsl_sf_eval1(fresnel_s, x);
-}
-static VALUE rb_fresnel_c1(VALUE obj, VALUE x)
-{
-	return rb_gsl_sf_eval1(fresnel_c1, x);
-}
-static VALUE rb_fresnel_s1(VALUE obj, VALUE x)
-{
-	return rb_gsl_sf_eval1(fresnel_s1, x);
-}
-void Init_fresnel(VALUE module)
-{
-	VALUE mfresnel;
-	mfresnel = rb_define_module_under(module, "Fresnel");
-	rb_define_module_function(module, "fresnel_c", rb_fresnel_c, 1);
-	rb_define_module_function(module, "fresnel_s", rb_fresnel_s, 1);	
-	rb_define_module_function(module, "fresnel_c1", rb_fresnel_c1, 1);	
-	rb_define_module_function(module, "fresnel_s1", rb_fresnel_s1, 1);		
-	rb_define_module_function(mfresnel, "c", rb_fresnel_c, 1);
-	rb_define_module_function(mfresnel, "s", rb_fresnel_s, 1);	
-	rb_define_module_function(mfresnel, "c1", rb_fresnel_c1, 1);	
-	rb_define_module_function(mfresnel, "s1", rb_fresnel_s1, 1);	
-}
-
-
-
-
-
-
-
-
-
-
-