@stdlib/stats-base-dists-f-kurtosis 0.2.2 → 0.3.1

package/NOTICE

@@ -1 +1 @@
-Copyright (c) 2016-2024 The Stdlib Authors.
+Copyright (c) 2016-2026 The Stdlib Authors.

package/README.md

@@ -78,7 +78,7 @@ var kurtosis = require( '@stdlib/stats-base-dists-f-kurtosis' );
 
 #### kurtosis( d1, d2 )
 
-Returns the [excess kurtosis][kurtosis] of a [F][f-distribution] distribution with parameters `d1` (numerator degrees of freedom) and `d2` (denominator degrees of freedom).
+Returns the [excess kurtosis][kurtosis] of an [F][f-distribution] distribution with parameters `d1` (numerator degrees of freedom) and `d2` (denominator degrees of freedom).
 
 ```javascript
 var v = kurtosis( 4.0, 9.0 );
@@ -145,20 +145,109 @@ v = kurtosis( 3.0, -1.0 );
 <!-- eslint no-undef: "error" -->
 
 ```javascript
-var randu = require( '@stdlib/random-base-randu' );
+var uniform = require( '@stdlib/random-array-uniform' );
+var logEachMap = require( '@stdlib/console-log-each-map' );
 var EPS = require( '@stdlib/constants-float64-eps' );
 var kurtosis = require( '@stdlib/stats-base-dists-f-kurtosis' );
 
-var d1;
-var d2;
-var v;
-var i;
+var opts = {
+    'dtype': 'float64'
+};
+var d1 = uniform( 10, EPS, 10.0, opts );
+var d2 = uniform( 10, EPS, 20.0, opts );
 
-for ( i = 0; i < 10; i++ ) {
-    d1 = ( randu()*10.0 ) + EPS;
-    d2 = ( randu()*20.0 ) + EPS;
-    v = kurtosis( d1, d2 );
-    console.log( 'd1: %d, d2: %d, skew(X;d1,d2): %d', d1.toFixed( 4 ), d2.toFixed( 4 ), v.toFixed( 4 ) );
+logEachMap( 'd1: %0.4f, d2: %0.4f, Kurt(X;d1,d2): %0.4f', d1, d2, kurtosis );
+```
+
+</section>
+
+<!-- /.examples -->
+
+<!-- C interface documentation. -->
+
+* * *
+
+<section class="c">
+
+## C APIs
+
+<!-- Section to include introductory text. Make sure to keep an empty line after the intro `section` element and another before the `/section` close. -->
+
+<section class="intro">
+
+</section>
+
+<!-- /.intro -->
+
+<!-- C usage documentation. -->
+
+<section class="usage">
+
+### Usage
+
+```c
+#include "stdlib/stats/base/dists/f/kurtosis.h"
+```
+
+#### stdlib_base_dists_f_kurtosis( d1, d2 )
+
+Evaluates the [excess kurtosis][kurtosis] of an [F][f-distribution] distribution with parameters `d1` (numerator degrees of freedom) and `d2` (denominator degrees of freedom).
+
+```c
+double out = stdlib_base_dists_f_kurtosis( 3.0, 9.0 );
+// returns ~124.667
+```
+
+The function accepts the following arguments:
+
+-   **d1**: `[in] double` numerator degrees of freedom.
+-   **d2**: `[in] double` denominator degrees of freedom.
+
+```c
+double stdlib_base_dists_f_kurtosis( const double d1, const double d2 );
+```
+
+</section>
+
+<!-- /.usage -->
+
+<!-- C API usage notes. Make sure to keep an empty line after the `section` element and another before the `/section` close. -->
+
+<section class="notes">
+
+</section>
+
+<!-- /.notes -->
+
+<!-- C API usage examples. -->
+
+<section class="examples">
+
+### Examples
+
+```c
+#include "stdlib/stats/base/dists/f/kurtosis.h"
+#include "stdlib/constants/float64/eps.h"
+#include <stdlib.h>
+#include <stdio.h>
+
+static double random_uniform( const double min, const double max ) {
+    double v = (double)rand() / ( (double)RAND_MAX + 1.0 );
+    return min + ( v*(max-min) );
+}
+
+int main( void ) {
+    double d1;
+    double d2;
+    double y;
+    int i;
+
+    for ( i = 0; i < 25; i++ ) {
+        d1 = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+        d1 = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 10.0 );
+        y = stdlib_base_dists_f_kurtosis( d1, d2 );
+        printf( "d1: %lf, d2: %lf, Kurt(X;d1,d2): %lf\n", d1, d2, y );
+    }
 }
 ```
 
@@ -166,6 +255,10 @@ for ( i = 0; i < 10; i++ ) {
 
 <!-- /.examples -->
 
+</section>
+
+<!-- /.c -->
+
 <!-- Section to include cited references. If references are included, add a horizontal rule *before* the section. Make sure to keep an empty line after the `section` element and another before the `/section` close. -->
 
 <section class="references">
@@ -208,7 +301,7 @@ See [LICENSE][stdlib-license].
 
 ## Copyright
 
-Copyright &copy; 2016-2024. The Stdlib [Authors][stdlib-authors].
+Copyright &copy; 2016-2026. The Stdlib [Authors][stdlib-authors].
 
 </section>
 
@@ -221,8 +314,8 @@ Copyright &copy; 2016-2024. The Stdlib [Authors][stdlib-authors].
 [npm-image]: http://img.shields.io/npm/v/@stdlib/stats-base-dists-f-kurtosis.svg
 [npm-url]: https://npmjs.org/package/@stdlib/stats-base-dists-f-kurtosis
 
-[test-image]: https://github.com/stdlib-js/stats-base-dists-f-kurtosis/actions/workflows/test.yml/badge.svg?branch=v0.2.2
-[test-url]: https://github.com/stdlib-js/stats-base-dists-f-kurtosis/actions/workflows/test.yml?query=branch:v0.2.2
+[test-image]: https://github.com/stdlib-js/stats-base-dists-f-kurtosis/actions/workflows/test.yml/badge.svg?branch=v0.3.1
+[test-url]: https://github.com/stdlib-js/stats-base-dists-f-kurtosis/actions/workflows/test.yml?query=branch:v0.3.1
 
 [coverage-image]: https://img.shields.io/codecov/c/github/stdlib-js/stats-base-dists-f-kurtosis/main.svg
 [coverage-url]: https://codecov.io/github/stdlib-js/stats-base-dists-f-kurtosis?branch=main
@@ -234,8 +327,8 @@ Copyright &copy; 2016-2024. The Stdlib [Authors][stdlib-authors].
 
 -->
 
-[chat-image]: https://img.shields.io/gitter/room/stdlib-js/stdlib.svg
-[chat-url]: https://app.gitter.im/#/room/#stdlib-js_stdlib:gitter.im
+[chat-image]: https://img.shields.io/badge/zulip-join_chat-brightgreen.svg
+[chat-url]: https://stdlib.zulipchat.com
 
 [stdlib]: https://github.com/stdlib-js/stdlib
 

package/dist/index.js.map

@@ -1,7 +1,7 @@
 {
   "version": 3,
   "sources": ["../lib/main.js", "../lib/index.js"],
-  "sourcesContent": ["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n// MODULES //\n\nvar pow = require( '@stdlib/math-base-special-pow' );\n\n\n// MAIN //\n\n/**\n* Returns the excess kurtosis of an F distribution.\n*\n* @param {PositiveNumber} d1 - numerator degrees of freedom\n* @param {PositiveNumber} d2 - denominator degrees of freedom\n* @returns {PositiveNumber} kurtosis\n*\n* @example\n* var v = kurtosis( 3.0, 9.0 );\n* // returns ~124.667\n*\n* @example\n* var v = kurtosis( 4.0, 12.0 );\n* // returns ~26.143\n*\n* @example\n* var v = kurtosis( 8.0, 9.0 );\n* // returns ~100.167\n*\n* @example\n* var v = kurtosis( 1.0, 8.0 );\n* // returns NaN\n*\n* @example\n* var v = kurtosis( 1.0, -0.1 );\n* // returns NaN\n*\n* @example\n* var v = kurtosis( -0.1, 1.0 );\n* // returns NaN\n*\n* @example\n* var v = kurtosis( 2.0, NaN );\n* // returns NaN\n*\n* @example\n* var v = kurtosis( NaN, 2.0 );\n* // returns NaN\n*/\nfunction kurtosis( d1, d2 ) {\n\tvar denom;\n\tvar num;\n\n\tif ( d1 <= 0.0 || d2 <= 8.0 ) {\n\t\treturn NaN;\n\t}\n\tnum = ( d1 * ( ( 5.0*d2 ) - 22.0 ) * ( d1+d2-2.0 ) ) +\n\t\t( ( d2-4.0 ) * pow( d2-2.0, 2.0 ) );\n\tdenom = d1 * ( d2-6.0 ) * ( d2-8.0 ) * ( d1+d2-2.0 );\n\treturn 12.0 * num / denom;\n}\n\n\n// EXPORTS //\n\nmodule.exports = kurtosis;\n", "/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n/**\n* F distribution excess kurtosis.\n*\n* @module @stdlib/stats-base-dists-f-kurtosis\n*\n* @example\n* var kurtosis = require( '@stdlib/stats-base-dists-f-kurtosis' );\n*\n* var v = kurtosis( 3.0, 9.0 );\n* // returns ~124.667\n*\n* v = kurtosis( 4.0, 12.0 );\n* // returns ~26.143\n*\n* v = kurtosis( 8.0, 9.0 );\n* // returns ~100.167\n*/\n\n// MODULES //\n\nvar main = require( './main.js' );\n\n\n// EXPORTS //\n\nmodule.exports = main;\n"],
+  "sourcesContent": ["/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n// MODULES //\n\nvar pow = require( '@stdlib/math-base-special-pow' );\n\n\n// MAIN //\n\n/**\n* Returns the excess kurtosis of an F distribution.\n*\n* @param {PositiveNumber} d1 - numerator degrees of freedom\n* @param {PositiveNumber} d2 - denominator degrees of freedom\n* @returns {PositiveNumber} excess kurtosis\n*\n* @example\n* var v = kurtosis( 3.0, 9.0 );\n* // returns ~124.667\n*\n* @example\n* var v = kurtosis( 4.0, 12.0 );\n* // returns ~26.143\n*\n* @example\n* var v = kurtosis( 8.0, 9.0 );\n* // returns ~100.167\n*\n* @example\n* var v = kurtosis( 1.0, 8.0 );\n* // returns NaN\n*\n* @example\n* var v = kurtosis( 1.0, -0.1 );\n* // returns NaN\n*\n* @example\n* var v = kurtosis( -0.1, 1.0 );\n* // returns NaN\n*\n* @example\n* var v = kurtosis( 2.0, NaN );\n* // returns NaN\n*\n* @example\n* var v = kurtosis( NaN, 2.0 );\n* // returns NaN\n*/\nfunction kurtosis( d1, d2 ) {\n\tvar denom;\n\tvar num;\n\n\tif ( d1 <= 0.0 || d2 <= 8.0 ) {\n\t\treturn NaN;\n\t}\n\tnum = ( d1 * ( ( 5.0*d2 ) - 22.0 ) * ( d1+d2-2.0 ) ) +\n\t\t( ( d2-4.0 ) * pow( d2-2.0, 2.0 ) );\n\tdenom = d1 * ( d2-6.0 ) * ( d2-8.0 ) * ( d1+d2-2.0 );\n\treturn 12.0 * num / denom;\n}\n\n\n// EXPORTS //\n\nmodule.exports = kurtosis;\n", "/**\n* @license Apache-2.0\n*\n* Copyright (c) 2018 The Stdlib Authors.\n*\n* Licensed under the Apache License, Version 2.0 (the \"License\");\n* you may not use this file except in compliance with the License.\n* You may obtain a copy of the License at\n*\n*    http://www.apache.org/licenses/LICENSE-2.0\n*\n* Unless required by applicable law or agreed to in writing, software\n* distributed under the License is distributed on an \"AS IS\" BASIS,\n* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n* See the License for the specific language governing permissions and\n* limitations under the License.\n*/\n\n'use strict';\n\n/**\n* F distribution excess kurtosis.\n*\n* @module @stdlib/stats-base-dists-f-kurtosis\n*\n* @example\n* var kurtosis = require( '@stdlib/stats-base-dists-f-kurtosis' );\n*\n* var v = kurtosis( 3.0, 9.0 );\n* // returns ~124.667\n*\n* v = kurtosis( 4.0, 12.0 );\n* // returns ~26.143\n*\n* v = kurtosis( 8.0, 9.0 );\n* // returns ~100.167\n*/\n\n// MODULES //\n\nvar main = require( './main.js' );\n\n\n// EXPORTS //\n\nmodule.exports = main;\n"],
   "mappings": "uGAAA,IAAAA,EAAAC,EAAA,SAAAC,EAAAC,EAAA,cAsBA,IAAIC,EAAM,QAAS,+BAAgC,EA4CnD,SAASC,EAAUC,EAAIC,EAAK,CAC3B,IAAIC,EACAC,EAEJ,OAAKH,GAAM,GAAOC,GAAM,EAChB,KAERE,EAAQH,GAAS,EAAIC,EAAO,KAAWD,EAAGC,EAAG,IACxCA,EAAG,GAAQH,EAAKG,EAAG,EAAK,CAAI,EACjCC,EAAQF,GAAOC,EAAG,IAAUA,EAAG,IAAUD,EAAGC,EAAG,GACxC,GAAOE,EAAMD,EACrB,CAKAL,EAAO,QAAUE,IC1CjB,IAAIK,EAAO,IAKX,OAAO,QAAUA",
   "names": ["require_main", "__commonJSMin", "exports", "module", "pow", "kurtosis", "d1", "d2", "denom", "num", "main"]
 }

package/include/stdlib/stats/base/dists/f/kurtosis.h

@@ -0,0 +1,38 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+#ifndef STDLIB_STATS_BASE_DISTS_F_KURTOSIS_H
+#define STDLIB_STATS_BASE_DISTS_F_KURTOSIS_H
+
+/*
+* If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler.
+*/
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+/**
+* Returns the excess kurtosis of an F distribution.
+*/
+double stdlib_base_dists_f_kurtosis( const double d1, const double d2 );
+
+#ifdef __cplusplus
+}
+#endif
+
+#endif // !STDLIB_STATS_BASE_DISTS_F_KURTOSIS_H

package/lib/main.js

@@ -30,7 +30,7 @@ var pow = require( '@stdlib/math-base-special-pow' );
 *
 * @param {PositiveNumber} d1 - numerator degrees of freedom
 * @param {PositiveNumber} d2 - denominator degrees of freedom
-* @returns {PositiveNumber} kurtosis
+* @returns {PositiveNumber} excess kurtosis
 *
 * @example
 * var v = kurtosis( 3.0, 9.0 );

package/lib/native.js

@@ -0,0 +1,75 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+'use strict';
+
+// MODULES //
+
+var addon = require( './../src/addon.node' );
+
+
+// MAIN //
+
+/**
+* Returns the excess kurtosis of an F distribution.
+*
+* @private
+* @param {PositiveNumber} d1 - numerator degrees of freedom
+* @param {PositiveNumber} d2 - denominator degrees of freedom
+* @returns {PositiveNumber} excess kurtosis
+*
+* @example
+* var v = kurtosis( 3.0, 9.0 );
+* // returns ~124.667
+*
+* @example
+* var v = kurtosis( 4.0, 12.0 );
+* // returns ~26.143
+*
+* @example
+* var v = kurtosis( 8.0, 9.0 );
+* // returns ~100.167
+*
+* @example
+* var v = kurtosis( 1.0, 8.0 );
+* // returns NaN
+*
+* @example
+* var v = kurtosis( 1.0, -0.1 );
+* // returns NaN
+*
+* @example
+* var v = kurtosis( -0.1, 1.0 );
+* // returns NaN
+*
+* @example
+* var v = kurtosis( 2.0, NaN );
+* // returns NaN
+*
+* @example
+* var v = kurtosis( NaN, 2.0 );
+* // returns NaN
+*/
+function kurtosis( d1, d2 ) {
+	return addon( d1, d2 );
+}
+
+
+// EXPORTS //
+
+module.exports = kurtosis;

package/manifest.json

@@ -0,0 +1,78 @@
+{
+  "options": {
+    "task": "build",
+    "wasm": false
+  },
+  "fields": [
+    {
+      "field": "src",
+      "resolve": true,
+      "relative": true
+    },
+    {
+      "field": "include",
+      "resolve": true,
+      "relative": true
+    },
+    {
+      "field": "libraries",
+      "resolve": false,
+      "relative": false
+    },
+    {
+      "field": "libpath",
+      "resolve": true,
+      "relative": false
+    }
+  ],
+  "confs": [
+    {
+      "task": "build",
+      "wasm": false,
+      "src": [
+        "./src/main.c"
+      ],
+      "include": [
+        "./include"
+      ],
+      "libraries": [],
+      "libpath": [],
+      "dependencies": [
+        "@stdlib/math-base-napi-binary",
+        "@stdlib/math-base-special-pow"
+      ]
+    },
+    {
+      "task": "benchmark",
+      "wasm": false,
+      "src": [
+        "./src/main.c"
+      ],
+      "include": [
+        "./include"
+      ],
+      "libraries": [],
+      "libpath": [],
+      "dependencies": [
+        "@stdlib/math-base-special-pow",
+        "@stdlib/constants-float64-eps"
+      ]
+    },
+    {
+      "task": "examples",
+      "wasm": false,
+      "src": [
+        "./src/main.c"
+      ],
+      "include": [
+        "./include"
+      ],
+      "libraries": [],
+      "libpath": [],
+      "dependencies": [
+        "@stdlib/math-base-special-pow",
+        "@stdlib/constants-float64-eps"
+      ]
+    }
+  ]
+}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@stdlib/stats-base-dists-f-kurtosis",
-  "version": "0.2.2",
+  "version": "0.3.1",
   "description": "F distribution excess kurtosis.",
   "license": "Apache-2.0",
   "author": {
@@ -14,9 +14,12 @@
     }
   ],
   "main": "./lib",
+  "gypfile": false,
   "directories": {
     "doc": "./docs",
+    "include": "./include",
     "lib": "./lib",
+    "src": "./src",
     "dist": "./dist"
   },
   "types": "./docs/types",
@@ -30,7 +33,9 @@
     "url": "https://github.com/stdlib-js/stdlib/issues"
   },
   "dependencies": {
-    "@stdlib/math-base-special-pow": "^0.3.0"
+    "@stdlib/math-base-napi-binary": "^0.3.3",
+    "@stdlib/math-base-special-pow": "^0.3.0",
+    "@stdlib/utils-library-manifest": "^0.2.4"
   },
   "devDependencies": {},
   "engines": {

package/src/addon.c

@@ -0,0 +1,22 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+#include "stdlib/stats/base/dists/f/kurtosis.h"
+#include "stdlib/math/base/napi/binary.h"
+
+STDLIB_MATH_BASE_NAPI_MODULE_DD_D( stdlib_base_dists_f_kurtosis )

package/src/main.c

@@ -0,0 +1,43 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+#include "stdlib/stats/base/dists/f/kurtosis.h"
+#include "stdlib/math/base/special/pow.h"
+
+/**
+* Returns the excess kurtosis of an F distribution.
+*
+* @param d1    numerator degrees of freedom
+* @param d2    denominator degrees of freedom
+* @return      excess kurtosis
+*
+* @example
+* double y = stdlib_base_dists_f_kurtosis( 3.0, 9.0 );
+* // returns ~124.667
+*/
+double stdlib_base_dists_f_kurtosis( const double d1, const double d2 ) {
+	double denom;
+	double num;
+
+	if ( d1 <= 0.0 || d2 <= 8.0 ) {
+		return 0.0/0.0; // NaN
+	}
+	num = ( d1 * ( ( 5.0*d2 ) - 22.0 ) * ( d1+d2-2.0 ) ) + ( ( d2-4.0 ) * stdlib_base_pow( d2-2.0, 2.0 ) );
+	denom = d1 * ( d2-6.0 ) * ( d2-8.0 ) * ( d1+d2-2.0 );
+	return 12.0 * num / denom;
+}