robust-predicates 2.0.2 → 3.0.1

package/esm/util.js

@@ -2,7 +2,8 @@ export const epsilon = 1.1102230246251565e-16;
 export const splitter = 134217729;
 export const resulterrbound = (3 + 8 * epsilon) * epsilon;
 
-export function fast_expansion_sum_zeroelim(elen, e, flen, f, h) {
+// fast_expansion_sum_zeroelim routine from oritinal code
+export function sum(elen, e, flen, f, h) {
     let Q, Qnew, hh, bvirt;
     let enow = e[0];
     let fnow = f[0];
@@ -16,7 +17,7 @@ export function fast_expansion_sum_zeroelim(elen, e, flen, f, h) {
         fnow = f[++findex];
     }
     let hindex = 0;
-    if ((eindex < elen) && (findex < flen)) {
+    if (eindex < elen && findex < flen) {
         if ((fnow > enow) === (fnow > -enow)) {
             Qnew = enow + Q;
             hh = Q - (Qnew - enow);
@@ -30,7 +31,7 @@ export function fast_expansion_sum_zeroelim(elen, e, flen, f, h) {
         if (hh !== 0) {
             h[hindex++] = hh;
         }
-        while ((eindex < elen) && (findex < flen)) {
+        while (eindex < elen && findex < flen) {
             if ((fnow > enow) === (fnow > -enow)) {
                 Qnew = Q + enow;
                 bvirt = Qnew - Q;
@@ -74,24 +75,30 @@ export function fast_expansion_sum_zeroelim(elen, e, flen, f, h) {
     return hindex;
 }
 
-export function scale_expansion_zeroelim(elen, e, b, h) {
+export function sum_three(alen, a, blen, b, clen, c, tmp, out) {
+    return sum(sum(alen, a, blen, b, tmp), tmp, clen, c, out);
+}
+
+// scale_expansion_zeroelim routine from oritinal code
+export function scale(elen, e, b, h) {
     let Q, sum, hh, product1, product0;
     let bvirt, c, ahi, alo, bhi, blo;
 
     c = splitter * b;
     bhi = c - (c - b);
     blo = b - bhi;
-    Q = e[0] * b;
-    c = splitter * e[0];
-    ahi = c - (c - e[0]);
-    alo = e[0] - ahi;
+    let enow = e[0];
+    Q = enow * b;
+    c = splitter * enow;
+    ahi = c - (c - enow);
+    alo = enow - ahi;
     hh = alo * blo - (Q - ahi * bhi - alo * bhi - ahi * blo);
     let hindex = 0;
     if (hh !== 0) {
         h[hindex++] = hh;
     }
-    for (let eindex = 1; eindex < elen; eindex++) {
-        const enow = e[eindex];
+    for (let i = 1; i < elen; i++) {
+        enow = e[i];
         product1 = enow * b;
         c = splitter * enow;
         ahi = c - (c - enow);
@@ -109,12 +116,17 @@ export function scale_expansion_zeroelim(elen, e, b, h) {
             h[hindex++] = hh;
         }
     }
-    if (Q !== 0.0 || hindex === 0) {
+    if (Q !== 0 || hindex === 0) {
         h[hindex++] = Q;
     }
     return hindex;
 }
 
+export function negate(elen, e) {
+    for (let i = 0; i < elen; i++) e[i] = -e[i];
+    return elen;
+}
+
 export function estimate(elen, e) {
     let Q = e[0];
     for (let i = 1; i < elen; i++) Q += e[i];

package/index.d.ts

@@ -0,0 +1,49 @@
+export as namespace predicates;
+
+export function orient2d(
+    ax: number, ay: number,
+    bx: number, by: number,
+    cx: number, cy: number): number;
+
+export function orient2dfast(
+    ax: number, ay: number,
+    bx: number, by: number,
+    cx: number, cy: number): number;
+
+export function incircle(
+    ax: number, ay: number,
+    bx: number, by: number,
+    cx: number, cy: number,
+    dx: number, dy: number): number;
+
+export function incirclefast(
+    ax: number, ay: number,
+    bx: number, by: number,
+    cx: number, cy: number,
+    dx: number, dy: number): number;
+
+export function orient3d(
+    ax: number, ay: number, az: number,
+    bx: number, by: number, bz: number,
+    cx: number, cy: number, cz: number,
+    dx: number, dy: number, dz: number): number;
+
+export function orient3dfast(
+    ax: number, ay: number, az: number,
+    bx: number, by: number, bz: number,
+    cx: number, cy: number, cz: number,
+    dx: number, dy: number, dz: number): number;
+
+export function insphere(
+    ax: number, ay: number, az: number,
+    bx: number, by: number, bz: number,
+    cx: number, cy: number, cz: number,
+    dx: number, dy: number, dz: number,
+    ex: number, ey: number, ez: number): number;
+
+export function inspherefast(
+    ax: number, ay: number, az: number,
+    bx: number, by: number, bz: number,
+    cx: number, cy: number, cz: number,
+    dx: number, dy: number, dz: number,
+    ex: number, ey: number, ez: number): number;

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "robust-predicates",
-  "version": "2.0.2",
+  "version": "3.0.1",
   "description": "Fast robust predicates for computational geometry",
   "keywords": [
     "computational geometry",
@@ -8,46 +8,55 @@
   ],
   "author": "Vladimir Agafonkin",
   "license": "Unlicense",
-  "main": "umd/predicates.js",
+  "type": "module",
+  "main": "index.js",
   "unpkg": "umd/predicates.min.js",
   "module": "index.js",
+  "exports": "./index.js",
   "types": "index.d.ts",
   "scripts": {
     "build": "mkdirp esm && node compile.js",
     "lint": "eslint *.js src test/test.js",
-    "test": "npm run lint && npm run build && node -r esm test/test.js",
-    "cov": "rm -rf node_modules/.cache/esm && nyc --require esm -r lcov -r text node test/test.js",
-    "bench": "node -r esm bench.js",
+    "test": "npm run lint && npm run build && node test/test.js",
+    "cov": "c8 node test/test.js",
+    "bench": "node bench.js",
     "prepublishOnly": "npm run test && rollup -c"
   },
   "devDependencies": {
-    "eslint": "^6.3.0",
+    "c8": "^7.7.0",
+    "eslint": "^7.23.0",
     "eslint-config-mourner": "^3.0.0",
-    "esm": "^3.2.25",
-    "mkdirp": "^0.5.1",
+    "mkdirp": "^1.0.4",
     "nextafter": "^1.0.0",
-    "nyc": "^14.1.1",
     "robust-in-sphere": "^1.1.3",
     "robust-orientation": "^1.1.3",
-    "rollup": "^1.20.3",
-    "rollup-plugin-terser": "^5.1.1",
-    "tape": "^4.11.0",
-    "terser": "^4.2.1"
+    "rollup": "^2.44.0",
+    "rollup-plugin-terser": "^7.0.2",
+    "tape": "^5.2.2",
+    "terser": "^5.6.1"
   },
   "files": [
     "index.js",
+    "index.d.ts",
     "esm",
     "umd"
   ],
+  "repository": {
+    "type": "git",
+    "url": "https://github.com/mourner/robust-predicates.git"
+  },
   "eslintConfig": {
     "extends": "mourner",
+    "parserOptions": {
+      "ecmaVersion": 2020
+    },
     "rules": {
       "camelcase": 0,
       "new-cap": 0,
       "no-unused-vars": [
         2,
         {
-          "varsIgnorePattern": "splitter|bvirt|c|[ab]hi|[ab]lo|_[ijk0]"
+          "varsIgnorePattern": "splitter|bvirt|c|[ab]hi|[ab]lo|_[ijk0]|u3|[st][01]"
         }
       ],
       "no-lonely-if": 0
@@ -56,16 +65,13 @@
       "$Fast_Two_Sum": false,
       "$Two_Sum": false,
       "$Two_Diff_Tail": false,
-      "$Two_Diff": false,
       "$Split": false,
       "$Two_Product": false,
       "$Two_Product_Presplit": false,
-      "$Square": false,
-      "$Two_One_Sum": false,
-      "$Two_One_Diff": false,
-      "$Two_Two_Sum": false,
-      "$Two_Two_Diff": false,
-      "$Two_One_Product": false
+      "$Two_One_Product": false,
+      "$Cross_Product": false,
+      "$Square_Sum": false,
+      "$Two_Product_Sum": false
     }
   }
 }