Ticket #50338: patch-gsl-mgsl_sf.c.diff

gsl/mlgsl_sf.c

@@ -258 +258 @@
 SF3(ellint_F, Double_val, Double_val, GSL_MODE_val)
 SF3(ellint_E, Double_val, Double_val, GSL_MODE_val)
 SF4(ellint_P, Double_val, Double_val, Double_val, GSL_MODE_val)
-SF4(ellint_D, Double_val, Double_val, Double_val, GSL_MODE_val)
+SF3(ellint_D, Double_val, Double_val, GSL_MODE_val)
 SF3(ellint_RC, Double_val, Double_val, GSL_MODE_val)
 SF4(ellint_RD, Double_val, Double_val, Double_val, GSL_MODE_val)
 SF4(ellint_RF, Double_val, Double_val, Double_val, GSL_MODE_val)
@@ -453 +453 @@
 
 /* Associated Legendre Polynomials and Spherical Harmonics  */
 SF3(legendre_Plm, Int_val, Int_val, Double_val) 
-CAMLprim value 
+/* CAMLprim value 
 ml_gsl_sf_legendre_Plm_array(value lmax, value m, value x, value result_array)
-{
-  gsl_sf_legendre_Plm_array(Int_val(lmax),
-                            Int_val(m),
-                            Double_val(x),
-                            Double_array_val(result_array));
-  return Val_unit;
-}
 
+ ******
+ function removed to make it work with gsl 2.1
+ */
 SF3(legendre_sphPlm, Int_val, Int_val, Double_val) 
-CAMLprim value 
+/* CAMLprim value 
 ml_gsl_sf_legendre_sphPlm_array(value lmax, value m, 
-                                value x, value result_array)
-{
-  gsl_sf_legendre_sphPlm_array(Int_val(lmax),
-                               Int_val(m),
-                               Double_val(x),
-                               Double_array_val(result_array));
-  return Val_unit;
-}
+                                                                value x, value result_array)
+ *****
+ function removed to make it work with gsl 2.1
 
+ */
 CAMLprim value
 ml_gsl_sf_legendre_array_size(value lmax, value m)
 {
@@ -482 +474 @@
   CAMLlocal1(ret);
   int gsl_ret;
 
-  gsl_ret = gsl_sf_legendre_array_size(Int_val(lmax), Int_val(m));
+  // hack for compatibility with V2.1 of gsl lib, since gsl_sf_legendre_array_size seem deprecated
+  // I'm not sure this is correct... 
+  gsl_ret = gsl_sf_legendre_array_index(Int_val(lmax), Int_val(m));
   
   ret = Val_int(gsl_ret);
   CAMLreturn(ret);