statistics: gsl_stats_Qn0_from_sorted_data & gsl_stats_Qn_from_sorted_data

Author: mysiar

gslext.stub.php

@@ -90,4 +90,6 @@ function gsl_stats_gastwirth_from_sorted_data(array $sorted_data, int $stride):
 function gsl_stats_mad0(array $data, int $stride): float {}
 function gsl_stats_mad(array $data, int $stride): float {}
 function gsl_stats_Sn0_from_sorted_data(array $sorted_data, int $stride): float {}
-function gsl_stats_Sn_from_sorted_data(array $sorted_data, int $stride): float {}
+function gsl_stats_Sn_from_sorted_data(array $sorted_data, int $stride): float {}
+function gsl_stats_Qn0_from_sorted_data(array $sorted_data, int $stride): float {}
+function gsl_stats_Qn_from_sorted_data(array $sorted_data, int $stride): float {}

gslext_arginfo.h

@@ -1,5 +1,5 @@
 /* This is a generated file, edit the .stub.php file instead.
- * Stub hash: 30e64c54e8548bcb86f4957ae56c7255a38f13a9 */
+ * Stub hash: ff2821df6a63df9aa60e187ced5591b4928cf084 */
 
 ZEND_BEGIN_ARG_WITH_RETURN_TYPE_INFO_EX(arginfo_GSL_log1p, 0, 1, IS_DOUBLE, 0)
 	ZEND_ARG_TYPE_INFO(0, x, IS_DOUBLE, 0)
@@ -272,6 +272,10 @@ ZEND_END_ARG_INFO()
 
 #define arginfo_gsl_stats_Sn_from_sorted_data arginfo_gsl_stats_median_from_sorted_data
 
+#define arginfo_gsl_stats_Qn0_from_sorted_data arginfo_gsl_stats_median_from_sorted_data
+
+#define arginfo_gsl_stats_Qn_from_sorted_data arginfo_gsl_stats_median_from_sorted_data
+
 
 ZEND_FUNCTION(GSL_log1p);
 ZEND_FUNCTION(GSL_expm1);
@@ -351,6 +355,8 @@ ZEND_FUNCTION(gsl_stats_mad0);
 ZEND_FUNCTION(gsl_stats_mad);
 ZEND_FUNCTION(gsl_stats_Sn0_from_sorted_data);
 ZEND_FUNCTION(gsl_stats_Sn_from_sorted_data);
+ZEND_FUNCTION(gsl_stats_Qn0_from_sorted_data);
+ZEND_FUNCTION(gsl_stats_Qn_from_sorted_data);
 
 
 static const zend_function_entry ext_functions[] = {
@@ -432,5 +438,7 @@ static const zend_function_entry ext_functions[] = {
 	ZEND_FE(gsl_stats_mad, arginfo_gsl_stats_mad)
 	ZEND_FE(gsl_stats_Sn0_from_sorted_data, arginfo_gsl_stats_Sn0_from_sorted_data)
 	ZEND_FE(gsl_stats_Sn_from_sorted_data, arginfo_gsl_stats_Sn_from_sorted_data)
+	ZEND_FE(gsl_stats_Qn0_from_sorted_data, arginfo_gsl_stats_Qn0_from_sorted_data)
+	ZEND_FE(gsl_stats_Qn_from_sorted_data, arginfo_gsl_stats_Qn_from_sorted_data)
 	ZEND_FE_END
 };

inc/statistics.c

@@ -1332,3 +1332,53 @@ PHP_FUNCTION(gsl_stats_Sn_from_sorted_data)
     RETURN_DOUBLE(rv);
 }
 
+PHP_FUNCTION(gsl_stats_Qn0_from_sorted_data)
+{
+    zval *p_sorted_data, *p_stride;
+    long stride, n;
+    int *work_int;
+    double *sorted_data, *work, rv;
+
+    ZEND_PARSE_PARAMETERS_START(2, 2)
+            Z_PARAM_ZVAL(p_sorted_data)
+            Z_PARAM_ZVAL(p_stride)
+    ZEND_PARSE_PARAMETERS_END();
+
+    __zval_to_long(p_stride, &stride);
+
+    n = (long) zend_array_count(Z_ARR_P(p_sorted_data));
+
+    __alloc_double_array(&sorted_data, n);
+    __alloc_double_array(&work, 3 * n);
+    __alloc_int_array(&work_int, 5 * n);
+
+    rv = gsl_stats_Qn0_from_sorted_data(sorted_data, stride, n, work, work_int);
+
+    RETURN_DOUBLE(rv);
+}
+
+PHP_FUNCTION(gsl_stats_Qn_from_sorted_data)
+{
+    zval *p_sorted_data, *p_stride;
+    long stride, n;
+    int *work_int;
+    double *sorted_data, *work, rv;
+
+    ZEND_PARSE_PARAMETERS_START(2, 2)
+            Z_PARAM_ZVAL(p_sorted_data)
+            Z_PARAM_ZVAL(p_stride)
+    ZEND_PARSE_PARAMETERS_END();
+
+    __zval_to_long(p_stride, &stride);
+
+    n = (long) zend_array_count(Z_ARR_P(p_sorted_data));
+
+    __alloc_double_array(&sorted_data, n);
+    __alloc_double_array(&work, 3 * n);
+    __alloc_int_array(&work_int, 5 * n);
+
+    rv = gsl_stats_Qn_from_sorted_data(sorted_data, stride, n, work, work_int);
+
+    RETURN_DOUBLE(rv);
+}
+

tests/0007.statistics.phpt

@@ -260,6 +260,14 @@ var_dump("gsl_stats_Sn_from_sorted_data");
 $data = [1, 5, 11, 34, 56];
 var_dump(gsl_stats_Sn_from_sorted_data($data, 25));
 
+var_dump("gsl_stats_Qn0_from_sorted_data");
+$data = [12.6, 15.2, 16.1, 17.5, 18.3];
+var_dump(gsl_stats_Qn0_from_sorted_data($data, 1));
+
+var_dump("gsl_stats_Qn_from_sorted_data");
+$data = [12.6, 15.2, 16.1, 17.5, 18.3];
+var_dump(gsl_stats_Qn_from_sorted_data($data, 1));
+
 
 ?>
 --EXPECT--
@@ -370,4 +378,8 @@ float(1.482602218505602)
 string(30) "gsl_stats_Sn0_from_sorted_data"
 float(1.2975342124590806E-281)
 string(29) "gsl_stats_Sn_from_sorted_data"
-float(6.373052309258179E-233)
+float(6.373052309258179E-233)
+string(30) "gsl_stats_Qn0_from_sorted_data"
+float(1.0996497102983598E+248)
+string(29) "gsl_stats_Qn_from_sorted_data"
+float(2.0596179022126888E+248)