@math.gl/culling 3.6.2 → 4.0.0-alpha.1
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- package/dist/{esm/constants.js → constants.js} +0 -0
- package/dist/constants.js.map +1 -0
- package/dist/index.cjs +1049 -0
- package/dist/index.js +12 -0
- package/dist/index.js.map +1 -0
- package/dist/{esm/lib → lib}/algorithms/bounding-box-from-points.js +3 -3
- package/dist/lib/algorithms/bounding-box-from-points.js.map +1 -0
- package/dist/{esm/lib → lib}/algorithms/bounding-sphere-from-points.js +1 -2
- package/dist/lib/algorithms/bounding-sphere-from-points.js.map +1 -0
- package/dist/lib/algorithms/compute-eigen-decomposition.d.ts +1 -1
- package/dist/lib/algorithms/compute-eigen-decomposition.d.ts.map +1 -1
- package/dist/{esm/lib → lib}/algorithms/compute-eigen-decomposition.js +0 -0
- package/dist/lib/algorithms/compute-eigen-decomposition.js.map +1 -0
- package/dist/{esm/lib → lib}/bounding-volumes/axis-aligned-bounding-box.js +1 -1
- package/dist/lib/bounding-volumes/axis-aligned-bounding-box.js.map +1 -0
- package/dist/lib/bounding-volumes/bounding-sphere.d.ts.map +1 -1
- package/dist/{esm/lib → lib}/bounding-volumes/bounding-sphere.js +1 -1
- package/dist/lib/bounding-volumes/bounding-sphere.js.map +1 -0
- package/dist/{esm/lib → lib}/bounding-volumes/bounding-volume.js +0 -0
- package/dist/{es5/lib → lib}/bounding-volumes/bounding-volume.js.map +0 -0
- package/dist/lib/bounding-volumes/oriented-bounding-box.d.ts +2 -2
- package/dist/lib/bounding-volumes/oriented-bounding-box.d.ts.map +1 -1
- package/dist/{esm/lib → lib}/bounding-volumes/oriented-bounding-box.js +2 -2
- package/dist/lib/bounding-volumes/oriented-bounding-box.js.map +1 -0
- package/dist/{esm/lib → lib}/culling-volume.js +2 -5
- package/dist/lib/culling-volume.js.map +1 -0
- package/dist/lib/perspective-frustum.d.ts +1 -1
- package/dist/lib/perspective-frustum.d.ts.map +1 -1
- package/dist/{esm/lib → lib}/perspective-frustum.js +1 -1
- package/dist/lib/perspective-frustum.js.map +1 -0
- package/dist/lib/perspective-off-center-frustum.d.ts +1 -1
- package/dist/lib/perspective-off-center-frustum.d.ts.map +1 -1
- package/dist/{esm/lib → lib}/perspective-off-center-frustum.js +2 -2
- package/dist/lib/perspective-off-center-frustum.js.map +1 -0
- package/dist/{esm/lib → lib}/plane.js +3 -3
- package/dist/lib/plane.js.map +1 -0
- package/package.json +7 -6
- package/src/lib/algorithms/bounding-sphere-from-points.ts +1 -1
- package/src/lib/bounding-volumes/bounding-sphere.ts +2 -0
- package/src/lib/bounding-volumes/bounding-volume.ts +1 -1
- package/src/lib/bounding-volumes/oriented-bounding-box.ts +5 -2
- package/src/lib/culling-volume.ts +3 -3
- package/src/lib/perspective-frustum.ts +1 -1
- package/src/lib/plane.ts +4 -4
- package/dist/es5/constants.js +0 -13
- package/dist/es5/constants.js.map +0 -1
- package/dist/es5/index.js +0 -102
- package/dist/es5/index.js.map +0 -1
- package/dist/es5/lib/algorithms/bounding-box-from-points.js +0 -198
- package/dist/es5/lib/algorithms/bounding-box-from-points.js.map +0 -1
- package/dist/es5/lib/algorithms/bounding-sphere-from-points.js +0 -167
- package/dist/es5/lib/algorithms/bounding-sphere-from-points.js.map +0 -1
- package/dist/es5/lib/algorithms/compute-eigen-decomposition.js +0 -114
- package/dist/es5/lib/algorithms/compute-eigen-decomposition.js.map +0 -1
- package/dist/es5/lib/bounding-volumes/axis-aligned-bounding-box.js +0 -114
- package/dist/es5/lib/bounding-volumes/axis-aligned-bounding-box.js.map +0 -1
- package/dist/es5/lib/bounding-volumes/bounding-sphere.js +0 -148
- package/dist/es5/lib/bounding-volumes/bounding-sphere.js.map +0 -1
- package/dist/es5/lib/bounding-volumes/bounding-volume.js +0 -2
- package/dist/es5/lib/bounding-volumes/oriented-bounding-box.js +0 -253
- package/dist/es5/lib/bounding-volumes/oriented-bounding-box.js.map +0 -1
- package/dist/es5/lib/culling-volume.js +0 -154
- package/dist/es5/lib/culling-volume.js.map +0 -1
- package/dist/es5/lib/perspective-frustum.js +0 -161
- package/dist/es5/lib/perspective-frustum.js.map +0 -1
- package/dist/es5/lib/perspective-off-center-frustum.js +0 -183
- package/dist/es5/lib/perspective-off-center-frustum.js.map +0 -1
- package/dist/es5/lib/plane.js +0 -94
- package/dist/es5/lib/plane.js.map +0 -1
- package/dist/esm/constants.js.map +0 -1
- package/dist/esm/index.js +0 -12
- package/dist/esm/index.js.map +0 -1
- package/dist/esm/lib/algorithms/bounding-box-from-points.js.map +0 -1
- package/dist/esm/lib/algorithms/bounding-sphere-from-points.js.map +0 -1
- package/dist/esm/lib/algorithms/compute-eigen-decomposition.js.map +0 -1
- package/dist/esm/lib/bounding-volumes/axis-aligned-bounding-box.js.map +0 -1
- package/dist/esm/lib/bounding-volumes/bounding-sphere.js.map +0 -1
- package/dist/esm/lib/bounding-volumes/bounding-volume.js.map +0 -1
- package/dist/esm/lib/bounding-volumes/oriented-bounding-box.js.map +0 -1
- package/dist/esm/lib/culling-volume.js.map +0 -1
- package/dist/esm/lib/perspective-frustum.js.map +0 -1
- package/dist/esm/lib/perspective-off-center-frustum.js.map +0 -1
- package/dist/esm/lib/plane.js.map +0 -1
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"sources":["../../../../src/lib/algorithms/bounding-sphere-from-points.ts"],"names":["fromPointsXMin","Vector3","fromPointsYMin","fromPointsZMin","fromPointsXMax","fromPointsYMax","fromPointsZMax","fromPointsCurrentPos","fromPointsScratch","fromPointsRitterCenter","fromPointsMinBoxPt","fromPointsMaxBoxPt","fromPointsNaiveCenterScratch","volumeConstant","Math","PI","makeBoundingSphereFromPoints","positions","result","BoundingSphere","length","fromCenterRadius","currentPos","copy","xMin","yMin","zMin","xMax","yMax","zMax","position","x","y","z","xSpan","subtract","magnitudeSquared","ySpan","zSpan","diameter1","diameter2","maxSpan","ritterCenter","radiusSquared","ritterRadius","sqrt","minBoxPt","maxBoxPt","naiveCenter","add","multiplyByScalar","naiveRadius","r","magnitude","oldCenterToPointSquared","oldCenterToPoint","oldToNew","to","center","radius"],"mappings":";;;;;;;;;AAGA;;AACA;;;;;;;;AAGA,IAAMA,cAAc,GAAG,IAAIC,aAAJ,EAAvB;AACA,IAAMC,cAAc,GAAG,IAAID,aAAJ,EAAvB;AACA,IAAME,cAAc,GAAG,IAAIF,aAAJ,EAAvB;AACA,IAAMG,cAAc,GAAG,IAAIH,aAAJ,EAAvB;AACA,IAAMI,cAAc,GAAG,IAAIJ,aAAJ,EAAvB;AACA,IAAMK,cAAc,GAAG,IAAIL,aAAJ,EAAvB;AACA,IAAMM,oBAAoB,GAAG,IAAIN,aAAJ,EAA7B;AACA,IAAMO,iBAAiB,GAAG,IAAIP,aAAJ,EAA1B;AACA,IAAMQ,sBAAsB,GAAG,IAAIR,aAAJ,EAA/B;AACA,IAAMS,kBAAkB,GAAG,IAAIT,aAAJ,EAA3B;AACA,IAAMU,kBAAkB,GAAG,IAAIV,aAAJ,EAA3B;AACA,IAAMW,4BAA4B,GAAG,IAAIX,aAAJ,EAArC;AACA,IAAMY,cAAc,GAAI,MAAM,GAAP,GAAcC,IAAI,CAACC,EAA1C;;AAae,SAASC,4BAAT,CACbC,SADa,EAGG;AAAA,MADhBC,MACgB,uEADS,IAAIC,uBAAJ,EACT;;AAChB,MAAI,CAACF,SAAD,IAAcA,SAAS,CAACG,MAAV,KAAqB,CAAvC,EAA0C;AACxC,WAAOF,MAAM,CAACG,gBAAP,CAAwB,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAxB,EAAmC,CAAnC,CAAP;AACD;;AAED,MAAMC,UAAU,GAAGf,oBAAoB,CAACgB,IAArB,CAA0BN,SAAS,CAAC,CAAD,CAAnC,CAAnB;AAEA,MAAMO,IAAI,GAAGxB,cAAc,CAACuB,IAAf,CAAoBD,UAApB,CAAb;AACA,MAAMG,IAAI,GAAGvB,cAAc,CAACqB,IAAf,CAAoBD,UAApB,CAAb;AACA,MAAMI,IAAI,GAAGvB,cAAc,CAACoB,IAAf,CAAoBD,UAApB,CAAb;AAEA,MAAMK,IAAI,GAAGvB,cAAc,CAACmB,IAAf,CAAoBD,UAApB,CAAb;AACA,MAAMM,IAAI,GAAGvB,cAAc,CAACkB,IAAf,CAAoBD,UAApB,CAAb;AACA,MAAMO,IAAI,GAAGvB,cAAc,CAACiB,IAAf,CAAoBD,UAApB,CAAb;;AAbgB,6CAeOL,SAfP;AAAA;;AAAA;AAehB,wDAAkC;AAAA,UAAvBa,QAAuB;AAChCR,MAAAA,UAAU,CAACC,IAAX,CAAgBO,QAAhB;AAEA,UAAMC,CAAC,GAAGT,UAAU,CAACS,CAArB;AACA,UAAMC,CAAC,GAAGV,UAAU,CAACU,CAArB;AACA,UAAMC,CAAC,GAAGX,UAAU,CAACW,CAArB;;AAGA,UAAIF,CAAC,GAAGP,IAAI,CAACO,CAAb,EAAgB;AACdP,QAAAA,IAAI,CAACD,IAAL,CAAUD,UAAV;AACD;;AAED,UAAIS,CAAC,GAAGJ,IAAI,CAACI,CAAb,EAAgB;AACdJ,QAAAA,IAAI,CAACJ,IAAL,CAAUD,UAAV;AACD;;AAED,UAAIU,CAAC,GAAGP,IAAI,CAACO,CAAb,EAAgB;AACdP,QAAAA,IAAI,CAACF,IAAL,CAAUD,UAAV;AACD;;AAED,UAAIU,CAAC,GAAGJ,IAAI,CAACI,CAAb,EAAgB;AACdJ,QAAAA,IAAI,CAACL,IAAL,CAAUD,UAAV;AACD;;AAED,UAAIW,CAAC,GAAGP,IAAI,CAACO,CAAb,EAAgB;AACdP,QAAAA,IAAI,CAACH,IAAL,CAAUD,UAAV;AACD;;AAED,UAAIW,CAAC,GAAGJ,IAAI,CAACI,CAAb,EAAgB;AACdJ,QAAAA,IAAI,CAACN,IAAL,CAAUD,UAAV;AACD;AACF;AA9Ce;AAAA;AAAA;AAAA;AAAA;;AAiDhB,MAAMY,KAAK,GAAG1B,iBAAiB,CAACe,IAAlB,CAAuBI,IAAvB,EAA6BQ,QAA7B,CAAsCX,IAAtC,EAA4CY,gBAA5C,EAAd;AACA,MAAMC,KAAK,GAAG7B,iBAAiB,CAACe,IAAlB,CAAuBK,IAAvB,EAA6BO,QAA7B,CAAsCV,IAAtC,EAA4CW,gBAA5C,EAAd;AACA,MAAME,KAAK,GAAG9B,iBAAiB,CAACe,IAAlB,CAAuBM,IAAvB,EAA6BM,QAA7B,CAAsCT,IAAtC,EAA4CU,gBAA5C,EAAd;AAGA,MAAIG,SAAS,GAAGf,IAAhB;AACA,MAAIgB,SAAS,GAAGb,IAAhB;AACA,MAAIc,OAAO,GAAGP,KAAd;;AACA,MAAIG,KAAK,GAAGI,OAAZ,EAAqB;AACnBA,IAAAA,OAAO,GAAGJ,KAAV;AACAE,IAAAA,SAAS,GAAGd,IAAZ;AACAe,IAAAA,SAAS,GAAGZ,IAAZ;AACD;;AACD,MAAIU,KAAK,GAAGG,OAAZ,EAAqB;AACnBA,IAAAA,OAAO,GAAGH,KAAV;AACAC,IAAAA,SAAS,GAAGb,IAAZ;AACAc,IAAAA,SAAS,GAAGX,IAAZ;AACD;;AAGD,MAAMa,YAAY,GAAGjC,sBAArB;AACAiC,EAAAA,YAAY,CAACX,CAAb,GAAiB,CAACQ,SAAS,CAACR,CAAV,GAAcS,SAAS,CAACT,CAAzB,IAA8B,GAA/C;AACAW,EAAAA,YAAY,CAACV,CAAb,GAAiB,CAACO,SAAS,CAACP,CAAV,GAAcQ,SAAS,CAACR,CAAzB,IAA8B,GAA/C;AACAU,EAAAA,YAAY,CAACT,CAAb,GAAiB,CAACM,SAAS,CAACN,CAAV,GAAcO,SAAS,CAACP,CAAzB,IAA8B,GAA/C;AAGA,MAAIU,aAAa,GAAGnC,iBAAiB,CAACe,IAAlB,CAAuBiB,SAAvB,EAAkCL,QAAlC,CAA2CO,YAA3C,EAAyDN,gBAAzD,EAApB;AACA,MAAIQ,YAAY,GAAG9B,IAAI,CAAC+B,IAAL,CAAUF,aAAV,CAAnB;AAGA,MAAMG,QAAQ,GAAGpC,kBAAjB;AACAoC,EAAAA,QAAQ,CAACf,CAAT,GAAaP,IAAI,CAACO,CAAlB;AACAe,EAAAA,QAAQ,CAACd,CAAT,GAAaP,IAAI,CAACO,CAAlB;AACAc,EAAAA,QAAQ,CAACb,CAAT,GAAaP,IAAI,CAACO,CAAlB;AAEA,MAAMc,QAAQ,GAAGpC,kBAAjB;AACAoC,EAAAA,QAAQ,CAAChB,CAAT,GAAaJ,IAAI,CAACI,CAAlB;AACAgB,EAAAA,QAAQ,CAACf,CAAT,GAAaJ,IAAI,CAACI,CAAlB;AACAe,EAAAA,QAAQ,CAACd,CAAT,GAAaJ,IAAI,CAACI,CAAlB;AAEA,MAAMe,WAAW,GAAGpC,4BAA4B,CAC7CW,IADiB,CACZuB,QADY,EAEjBG,GAFiB,CAEbF,QAFa,EAGjBG,gBAHiB,CAGA,GAHA,CAApB;AAMA,MAAIC,WAAW,GAAG,CAAlB;;AA/FgB,8CAgGOlC,SAhGP;AAAA;;AAAA;AAgGhB,2DAAkC;AAAA,UAAvBa,SAAuB;AAChCR,MAAAA,UAAU,CAACC,IAAX,CAAgBO,SAAhB;AAGA,UAAMsB,CAAC,GAAG5C,iBAAiB,CAACe,IAAlB,CAAuBD,UAAvB,EAAmCa,QAAnC,CAA4Ca,WAA5C,EAAyDK,SAAzD,EAAV;;AACA,UAAID,CAAC,GAAGD,WAAR,EAAqB;AACnBA,QAAAA,WAAW,GAAGC,CAAd;AACD;;AAGD,UAAME,uBAAuB,GAAG9C,iBAAiB,CAC9Ce,IAD6B,CACxBD,UADwB,EAE7Ba,QAF6B,CAEpBO,YAFoB,EAG7BN,gBAH6B,EAAhC;;AAKA,UAAIkB,uBAAuB,GAAGX,aAA9B,EAA6C;AAC3C,YAAMY,gBAAgB,GAAGzC,IAAI,CAAC+B,IAAL,CAAUS,uBAAV,CAAzB;AAEAV,QAAAA,YAAY,GAAG,CAACA,YAAY,GAAGW,gBAAhB,IAAoC,GAAnD;AACAZ,QAAAA,aAAa,GAAGC,YAAY,GAAGA,YAA/B;AAEA,YAAMY,QAAQ,GAAGD,gBAAgB,GAAGX,YAApC;AACAF,QAAAA,YAAY,CAACX,CAAb,GAAiB,CAACa,YAAY,GAAGF,YAAY,CAACX,CAA5B,GAAgCyB,QAAQ,GAAGlC,UAAU,CAACS,CAAvD,IAA4DwB,gBAA7E;AACAb,QAAAA,YAAY,CAACV,CAAb,GAAiB,CAACY,YAAY,GAAGF,YAAY,CAACV,CAA5B,GAAgCwB,QAAQ,GAAGlC,UAAU,CAACU,CAAvD,IAA4DuB,gBAA7E;AACAb,QAAAA,YAAY,CAACT,CAAb,GAAiB,CAACW,YAAY,GAAGF,YAAY,CAACT,CAA5B,GAAgCuB,QAAQ,GAAGlC,UAAU,CAACW,CAAvD,IAA4DsB,gBAA7E;AACD;AACF;AA1He;AAAA;AAAA;AAAA;AAAA;;AA4HhB,MAAIX,YAAY,GAAGO,WAAnB,EAAgC;AAC9BT,IAAAA,YAAY,CAACe,EAAb,CAAgBvC,MAAM,CAACwC,MAAvB;AACAxC,IAAAA,MAAM,CAACyC,MAAP,GAAgBf,YAAhB;AACD,GAHD,MAGO;AACLI,IAAAA,WAAW,CAACS,EAAZ,CAAevC,MAAM,CAACwC,MAAtB;AACAxC,IAAAA,MAAM,CAACyC,MAAP,GAAgBR,WAAhB;AACD;;AAED,SAAOjC,MAAP;AACD","sourcesContent":["// This file is derived from the Cesium math library under Apache 2 license\n// See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md\n\nimport {Vector3} from '@math.gl/core';\nimport BoundingSphere from '../bounding-volumes/bounding-sphere';\n\n/* eslint-disable */\nconst fromPointsXMin = new Vector3();\nconst fromPointsYMin = new Vector3();\nconst fromPointsZMin = new Vector3();\nconst fromPointsXMax = new Vector3();\nconst fromPointsYMax = new Vector3();\nconst fromPointsZMax = new Vector3();\nconst fromPointsCurrentPos = new Vector3();\nconst fromPointsScratch = new Vector3();\nconst fromPointsRitterCenter = new Vector3();\nconst fromPointsMinBoxPt = new Vector3();\nconst fromPointsMaxBoxPt = new Vector3();\nconst fromPointsNaiveCenterScratch = new Vector3();\nconst volumeConstant = (4.0 / 3.0) * Math.PI;\n\n/**\n * Computes a tight-fitting bounding sphere enclosing a list of 3D Cartesian points.\n *\n * The bounding sphere is computed by running two algorithms, a naive algorithm and\n * Ritter's algorithm. The smaller of the two spheres is used to ensure a tight fit.\n * Bounding sphere computation article http://blogs.agi.com/insight3d/index.php/2008/02/04/a-bounding\n *\n * @param positions An array of points that the bounding sphere will enclose.\n * @param result Optional object onto which to store the result.\n * @returns The modified result parameter or a new `BoundingSphere` instance if not provided.\n */\nexport default function makeBoundingSphereFromPoints(\n positions: number[][],\n result: BoundingSphere = new BoundingSphere()\n): BoundingSphere {\n if (!positions || positions.length === 0) {\n return result.fromCenterRadius([0, 0, 0], 0);\n }\n\n const currentPos = fromPointsCurrentPos.copy(positions[0]);\n\n const xMin = fromPointsXMin.copy(currentPos);\n const yMin = fromPointsYMin.copy(currentPos);\n const zMin = fromPointsZMin.copy(currentPos);\n\n const xMax = fromPointsXMax.copy(currentPos);\n const yMax = fromPointsYMax.copy(currentPos);\n const zMax = fromPointsZMax.copy(currentPos);\n\n for (const position of positions) {\n currentPos.copy(position);\n\n const x = currentPos.x;\n const y = currentPos.y;\n const z = currentPos.z;\n\n // Store points containing the the smallest and largest components\n if (x < xMin.x) {\n xMin.copy(currentPos);\n }\n\n if (x > xMax.x) {\n xMax.copy(currentPos);\n }\n\n if (y < yMin.y) {\n yMin.copy(currentPos);\n }\n\n if (y > yMax.y) {\n yMax.copy(currentPos);\n }\n\n if (z < zMin.z) {\n zMin.copy(currentPos);\n }\n\n if (z > zMax.z) {\n zMax.copy(currentPos);\n }\n }\n\n // Compute x-, y-, and z-spans (Squared distances b/n each component's min. and max.).\n const xSpan = fromPointsScratch.copy(xMax).subtract(xMin).magnitudeSquared();\n const ySpan = fromPointsScratch.copy(yMax).subtract(yMin).magnitudeSquared();\n const zSpan = fromPointsScratch.copy(zMax).subtract(zMin).magnitudeSquared();\n\n // Set the diameter endpoints to the largest span.\n let diameter1 = xMin;\n let diameter2 = xMax;\n let maxSpan = xSpan;\n if (ySpan > maxSpan) {\n maxSpan = ySpan;\n diameter1 = yMin;\n diameter2 = yMax;\n }\n if (zSpan > maxSpan) {\n maxSpan = zSpan;\n diameter1 = zMin;\n diameter2 = zMax;\n }\n\n // Calculate the center of the initial sphere found by Ritter's algorithm\n const ritterCenter = fromPointsRitterCenter;\n ritterCenter.x = (diameter1.x + diameter2.x) * 0.5;\n ritterCenter.y = (diameter1.y + diameter2.y) * 0.5;\n ritterCenter.z = (diameter1.z + diameter2.z) * 0.5;\n\n // Calculate the radius of the initial sphere found by Ritter's algorithm\n let radiusSquared = fromPointsScratch.copy(diameter2).subtract(ritterCenter).magnitudeSquared();\n let ritterRadius = Math.sqrt(radiusSquared);\n\n // Find the center of the sphere found using the Naive method.\n const minBoxPt = fromPointsMinBoxPt;\n minBoxPt.x = xMin.x;\n minBoxPt.y = yMin.y;\n minBoxPt.z = zMin.z;\n\n const maxBoxPt = fromPointsMaxBoxPt;\n maxBoxPt.x = xMax.x;\n maxBoxPt.y = yMax.y;\n maxBoxPt.z = zMax.z;\n\n const naiveCenter = fromPointsNaiveCenterScratch\n .copy(minBoxPt)\n .add(maxBoxPt)\n .multiplyByScalar(0.5);\n\n // Begin 2nd pass to find naive radius and modify the ritter sphere.\n let naiveRadius = 0;\n for (const position of positions) {\n currentPos.copy(position);\n\n // Find the furthest point from the naive center to calculate the naive radius.\n const r = fromPointsScratch.copy(currentPos).subtract(naiveCenter).magnitude();\n if (r > naiveRadius) {\n naiveRadius = r;\n }\n\n // Make adjustments to the Ritter Sphere to include all points.\n const oldCenterToPointSquared = fromPointsScratch\n .copy(currentPos)\n .subtract(ritterCenter)\n .magnitudeSquared();\n\n if (oldCenterToPointSquared > radiusSquared) {\n const oldCenterToPoint = Math.sqrt(oldCenterToPointSquared);\n // Calculate new radius to include the point that lies outside\n ritterRadius = (ritterRadius + oldCenterToPoint) * 0.5;\n radiusSquared = ritterRadius * ritterRadius;\n // Calculate center of new Ritter sphere\n const oldToNew = oldCenterToPoint - ritterRadius;\n ritterCenter.x = (ritterRadius * ritterCenter.x + oldToNew * currentPos.x) / oldCenterToPoint;\n ritterCenter.y = (ritterRadius * ritterCenter.y + oldToNew * currentPos.y) / oldCenterToPoint;\n ritterCenter.z = (ritterRadius * ritterCenter.z + oldToNew * currentPos.z) / oldCenterToPoint;\n }\n }\n\n if (ritterRadius < naiveRadius) {\n ritterCenter.to(result.center);\n result.radius = ritterRadius;\n } else {\n naiveCenter.to(result.center);\n result.radius = naiveRadius;\n }\n\n return result;\n}\n"],"file":"bounding-sphere-from-points.js"}
|
|
@@ -1,114 +0,0 @@
|
|
|
1
|
-
"use strict";
|
|
2
|
-
|
|
3
|
-
Object.defineProperty(exports, "__esModule", {
|
|
4
|
-
value: true
|
|
5
|
-
});
|
|
6
|
-
exports.default = computeEigenDecomposition;
|
|
7
|
-
|
|
8
|
-
var _core = require("@math.gl/core");
|
|
9
|
-
|
|
10
|
-
var scratchMatrix = new _core.Matrix3();
|
|
11
|
-
var scratchUnitary = new _core.Matrix3();
|
|
12
|
-
var scratchDiagonal = new _core.Matrix3();
|
|
13
|
-
var jMatrix = new _core.Matrix3();
|
|
14
|
-
var jMatrixTranspose = new _core.Matrix3();
|
|
15
|
-
|
|
16
|
-
function computeEigenDecomposition(matrix) {
|
|
17
|
-
var result = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : {};
|
|
18
|
-
var EIGEN_TOLERANCE = _core._MathUtils.EPSILON20;
|
|
19
|
-
var EIGEN_MAX_SWEEPS = 10;
|
|
20
|
-
var count = 0;
|
|
21
|
-
var sweep = 0;
|
|
22
|
-
var unitaryMatrix = scratchUnitary;
|
|
23
|
-
var diagonalMatrix = scratchDiagonal;
|
|
24
|
-
unitaryMatrix.identity();
|
|
25
|
-
diagonalMatrix.copy(matrix);
|
|
26
|
-
var epsilon = EIGEN_TOLERANCE * computeFrobeniusNorm(diagonalMatrix);
|
|
27
|
-
|
|
28
|
-
while (sweep < EIGEN_MAX_SWEEPS && offDiagonalFrobeniusNorm(diagonalMatrix) > epsilon) {
|
|
29
|
-
shurDecomposition(diagonalMatrix, jMatrix);
|
|
30
|
-
jMatrixTranspose.copy(jMatrix).transpose();
|
|
31
|
-
diagonalMatrix.multiplyRight(jMatrix);
|
|
32
|
-
diagonalMatrix.multiplyLeft(jMatrixTranspose);
|
|
33
|
-
unitaryMatrix.multiplyRight(jMatrix);
|
|
34
|
-
|
|
35
|
-
if (++count > 2) {
|
|
36
|
-
++sweep;
|
|
37
|
-
count = 0;
|
|
38
|
-
}
|
|
39
|
-
}
|
|
40
|
-
|
|
41
|
-
result.unitary = unitaryMatrix.toTarget(result.unitary);
|
|
42
|
-
result.diagonal = diagonalMatrix.toTarget(result.diagonal);
|
|
43
|
-
return result;
|
|
44
|
-
}
|
|
45
|
-
|
|
46
|
-
function computeFrobeniusNorm(matrix) {
|
|
47
|
-
var norm = 0.0;
|
|
48
|
-
|
|
49
|
-
for (var i = 0; i < 9; ++i) {
|
|
50
|
-
var temp = matrix[i];
|
|
51
|
-
norm += temp * temp;
|
|
52
|
-
}
|
|
53
|
-
|
|
54
|
-
return Math.sqrt(norm);
|
|
55
|
-
}
|
|
56
|
-
|
|
57
|
-
var rowVal = [1, 0, 0];
|
|
58
|
-
var colVal = [2, 2, 1];
|
|
59
|
-
|
|
60
|
-
function offDiagonalFrobeniusNorm(matrix) {
|
|
61
|
-
var norm = 0.0;
|
|
62
|
-
|
|
63
|
-
for (var i = 0; i < 3; ++i) {
|
|
64
|
-
var temp = matrix[scratchMatrix.getElementIndex(colVal[i], rowVal[i])];
|
|
65
|
-
norm += 2.0 * temp * temp;
|
|
66
|
-
}
|
|
67
|
-
|
|
68
|
-
return Math.sqrt(norm);
|
|
69
|
-
}
|
|
70
|
-
|
|
71
|
-
function shurDecomposition(matrix, result) {
|
|
72
|
-
var tolerance = _core._MathUtils.EPSILON15;
|
|
73
|
-
var maxDiagonal = 0.0;
|
|
74
|
-
var rotAxis = 1;
|
|
75
|
-
|
|
76
|
-
for (var i = 0; i < 3; ++i) {
|
|
77
|
-
var temp = Math.abs(matrix[scratchMatrix.getElementIndex(colVal[i], rowVal[i])]);
|
|
78
|
-
|
|
79
|
-
if (temp > maxDiagonal) {
|
|
80
|
-
rotAxis = i;
|
|
81
|
-
maxDiagonal = temp;
|
|
82
|
-
}
|
|
83
|
-
}
|
|
84
|
-
|
|
85
|
-
var p = rowVal[rotAxis];
|
|
86
|
-
var q = colVal[rotAxis];
|
|
87
|
-
var c = 1.0;
|
|
88
|
-
var s = 0.0;
|
|
89
|
-
|
|
90
|
-
if (Math.abs(matrix[scratchMatrix.getElementIndex(q, p)]) > tolerance) {
|
|
91
|
-
var qq = matrix[scratchMatrix.getElementIndex(q, q)];
|
|
92
|
-
var pp = matrix[scratchMatrix.getElementIndex(p, p)];
|
|
93
|
-
var qp = matrix[scratchMatrix.getElementIndex(q, p)];
|
|
94
|
-
var tau = (qq - pp) / 2.0 / qp;
|
|
95
|
-
var t;
|
|
96
|
-
|
|
97
|
-
if (tau < 0.0) {
|
|
98
|
-
t = -1.0 / (-tau + Math.sqrt(1.0 + tau * tau));
|
|
99
|
-
} else {
|
|
100
|
-
t = 1.0 / (tau + Math.sqrt(1.0 + tau * tau));
|
|
101
|
-
}
|
|
102
|
-
|
|
103
|
-
c = 1.0 / Math.sqrt(1.0 + t * t);
|
|
104
|
-
s = t * c;
|
|
105
|
-
}
|
|
106
|
-
|
|
107
|
-
_core.Matrix3.IDENTITY.to(result);
|
|
108
|
-
|
|
109
|
-
result[scratchMatrix.getElementIndex(p, p)] = result[scratchMatrix.getElementIndex(q, q)] = c;
|
|
110
|
-
result[scratchMatrix.getElementIndex(q, p)] = s;
|
|
111
|
-
result[scratchMatrix.getElementIndex(p, q)] = -s;
|
|
112
|
-
return result;
|
|
113
|
-
}
|
|
114
|
-
//# sourceMappingURL=compute-eigen-decomposition.js.map
|
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"sources":["../../../../src/lib/algorithms/compute-eigen-decomposition.ts"],"names":["scratchMatrix","Matrix3","scratchUnitary","scratchDiagonal","jMatrix","jMatrixTranspose","computeEigenDecomposition","matrix","result","EIGEN_TOLERANCE","_MathUtils","EPSILON20","EIGEN_MAX_SWEEPS","count","sweep","unitaryMatrix","diagonalMatrix","identity","copy","epsilon","computeFrobeniusNorm","offDiagonalFrobeniusNorm","shurDecomposition","transpose","multiplyRight","multiplyLeft","unitary","toTarget","diagonal","norm","i","temp","Math","sqrt","rowVal","colVal","getElementIndex","tolerance","EPSILON15","maxDiagonal","rotAxis","abs","p","q","c","s","qq","pp","qp","tau","t","IDENTITY","to"],"mappings":";;;;;;;AAGA;;AAEA,IAAMA,aAAa,GAAG,IAAIC,aAAJ,EAAtB;AACA,IAAMC,cAAc,GAAG,IAAID,aAAJ,EAAvB;AACA,IAAME,eAAe,GAAG,IAAIF,aAAJ,EAAxB;AAEA,IAAMG,OAAO,GAAG,IAAIH,aAAJ,EAAhB;AACA,IAAMI,gBAAgB,GAAG,IAAIJ,aAAJ,EAAzB;;AAqCe,SAASK,yBAAT,CACbC,MADa,EAIO;AAAA,MADpBC,MACoB,uEADS,EACT;AACpB,MAAMC,eAAe,GAAGC,iBAAWC,SAAnC;AACA,MAAMC,gBAAgB,GAAG,EAAzB;AAEA,MAAIC,KAAK,GAAG,CAAZ;AACA,MAAIC,KAAK,GAAG,CAAZ;AAEA,MAAMC,aAAa,GAAGb,cAAtB;AACA,MAAMc,cAAc,GAAGb,eAAvB;AAEAY,EAAAA,aAAa,CAACE,QAAd;AACAD,EAAAA,cAAc,CAACE,IAAf,CAAoBX,MAApB;AAEA,MAAMY,OAAO,GAAGV,eAAe,GAAGW,oBAAoB,CAACJ,cAAD,CAAtD;;AAEA,SAAOF,KAAK,GAAGF,gBAAR,IAA4BS,wBAAwB,CAACL,cAAD,CAAxB,GAA2CG,OAA9E,EAAuF;AACrFG,IAAAA,iBAAiB,CAACN,cAAD,EAAiBZ,OAAjB,CAAjB;AAEAC,IAAAA,gBAAgB,CAACa,IAAjB,CAAsBd,OAAtB,EAA+BmB,SAA/B;AAEAP,IAAAA,cAAc,CAACQ,aAAf,CAA6BpB,OAA7B;AACAY,IAAAA,cAAc,CAACS,YAAf,CAA4BpB,gBAA5B;AACAU,IAAAA,aAAa,CAACS,aAAd,CAA4BpB,OAA5B;;AAEA,QAAI,EAAES,KAAF,GAAU,CAAd,EAAiB;AACf,QAAEC,KAAF;AACAD,MAAAA,KAAK,GAAG,CAAR;AACD;AACF;;AAEDL,EAAAA,MAAM,CAACkB,OAAP,GAAiBX,aAAa,CAACY,QAAd,CAAuBnB,MAAM,CAACkB,OAA9B,CAAjB;AACAlB,EAAAA,MAAM,CAACoB,QAAP,GAAkBZ,cAAc,CAACW,QAAf,CAAwBnB,MAAM,CAACoB,QAA/B,CAAlB;AAEA,SAAOpB,MAAP;AACD;;AAED,SAASY,oBAAT,CAA8Bb,MAA9B,EAAuD;AACrD,MAAIsB,IAAI,GAAG,GAAX;;AACA,OAAK,IAAIC,CAAC,GAAG,CAAb,EAAgBA,CAAC,GAAG,CAApB,EAAuB,EAAEA,CAAzB,EAA4B;AAC1B,QAAMC,IAAI,GAAGxB,MAAM,CAACuB,CAAD,CAAnB;AACAD,IAAAA,IAAI,IAAIE,IAAI,GAAGA,IAAf;AACD;;AACD,SAAOC,IAAI,CAACC,IAAL,CAAUJ,IAAV,CAAP;AACD;;AAED,IAAMK,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAf;AACA,IAAMC,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAf;;AAIA,SAASd,wBAAT,CAAkCd,MAAlC,EAA2D;AACzD,MAAIsB,IAAI,GAAG,GAAX;;AACA,OAAK,IAAIC,CAAC,GAAG,CAAb,EAAgBA,CAAC,GAAG,CAApB,EAAuB,EAAEA,CAAzB,EAA4B;AAC1B,QAAMC,IAAI,GAAGxB,MAAM,CAACP,aAAa,CAACoC,eAAd,CAA8BD,MAAM,CAACL,CAAD,CAApC,EAAyCI,MAAM,CAACJ,CAAD,CAA/C,CAAD,CAAnB;AACAD,IAAAA,IAAI,IAAI,MAAME,IAAN,GAAaA,IAArB;AACD;;AACD,SAAOC,IAAI,CAACC,IAAL,CAAUJ,IAAV,CAAP;AACD;;AAUD,SAASP,iBAAT,CAA2Bf,MAA3B,EAA4CC,MAA5C,EAAsE;AACpE,MAAM6B,SAAS,GAAG3B,iBAAW4B,SAA7B;AAEA,MAAIC,WAAW,GAAG,GAAlB;AACA,MAAIC,OAAO,GAAG,CAAd;;AAGA,OAAK,IAAIV,CAAC,GAAG,CAAb,EAAgBA,CAAC,GAAG,CAApB,EAAuB,EAAEA,CAAzB,EAA4B;AAC1B,QAAMC,IAAI,GAAGC,IAAI,CAACS,GAAL,CAASlC,MAAM,CAACP,aAAa,CAACoC,eAAd,CAA8BD,MAAM,CAACL,CAAD,CAApC,EAAyCI,MAAM,CAACJ,CAAD,CAA/C,CAAD,CAAf,CAAb;;AACA,QAAIC,IAAI,GAAGQ,WAAX,EAAwB;AACtBC,MAAAA,OAAO,GAAGV,CAAV;AACAS,MAAAA,WAAW,GAAGR,IAAd;AACD;AACF;;AAED,MAAMW,CAAC,GAAGR,MAAM,CAACM,OAAD,CAAhB;AACA,MAAMG,CAAC,GAAGR,MAAM,CAACK,OAAD,CAAhB;AAEA,MAAII,CAAC,GAAG,GAAR;AACA,MAAIC,CAAC,GAAG,GAAR;;AAEA,MAAIb,IAAI,CAACS,GAAL,CAASlC,MAAM,CAACP,aAAa,CAACoC,eAAd,CAA8BO,CAA9B,EAAiCD,CAAjC,CAAD,CAAf,IAAwDL,SAA5D,EAAuE;AACrE,QAAMS,EAAE,GAAGvC,MAAM,CAACP,aAAa,CAACoC,eAAd,CAA8BO,CAA9B,EAAiCA,CAAjC,CAAD,CAAjB;AACA,QAAMI,EAAE,GAAGxC,MAAM,CAACP,aAAa,CAACoC,eAAd,CAA8BM,CAA9B,EAAiCA,CAAjC,CAAD,CAAjB;AACA,QAAMM,EAAE,GAAGzC,MAAM,CAACP,aAAa,CAACoC,eAAd,CAA8BO,CAA9B,EAAiCD,CAAjC,CAAD,CAAjB;AAEA,QAAMO,GAAG,GAAG,CAACH,EAAE,GAAGC,EAAN,IAAY,GAAZ,GAAkBC,EAA9B;AACA,QAAIE,CAAJ;;AAEA,QAAID,GAAG,GAAG,GAAV,EAAe;AACbC,MAAAA,CAAC,GAAG,CAAC,GAAD,IAAQ,CAACD,GAAD,GAAOjB,IAAI,CAACC,IAAL,CAAU,MAAMgB,GAAG,GAAGA,GAAtB,CAAf,CAAJ;AACD,KAFD,MAEO;AACLC,MAAAA,CAAC,GAAG,OAAOD,GAAG,GAAGjB,IAAI,CAACC,IAAL,CAAU,MAAMgB,GAAG,GAAGA,GAAtB,CAAb,CAAJ;AACD;;AAEDL,IAAAA,CAAC,GAAG,MAAMZ,IAAI,CAACC,IAAL,CAAU,MAAMiB,CAAC,GAAGA,CAApB,CAAV;AACAL,IAAAA,CAAC,GAAGK,CAAC,GAAGN,CAAR;AACD;;AAGD3C,gBAAQkD,QAAR,CAAiBC,EAAjB,CAAoB5C,MAApB;;AACAA,EAAAA,MAAM,CAACR,aAAa,CAACoC,eAAd,CAA8BM,CAA9B,EAAiCA,CAAjC,CAAD,CAAN,GAA8ClC,MAAM,CAACR,aAAa,CAACoC,eAAd,CAA8BO,CAA9B,EAAiCA,CAAjC,CAAD,CAAN,GAA8CC,CAA5F;AACApC,EAAAA,MAAM,CAACR,aAAa,CAACoC,eAAd,CAA8BO,CAA9B,EAAiCD,CAAjC,CAAD,CAAN,GAA8CG,CAA9C;AACArC,EAAAA,MAAM,CAACR,aAAa,CAACoC,eAAd,CAA8BM,CAA9B,EAAiCC,CAAjC,CAAD,CAAN,GAA8C,CAACE,CAA/C;AAEA,SAAOrC,MAAP;AACD","sourcesContent":["// This file is derived from the Cesium math library under Apache 2 license\n// See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md\n\nimport {Matrix3, _MathUtils} from '@math.gl/core';\n\nconst scratchMatrix = new Matrix3();\nconst scratchUnitary = new Matrix3();\nconst scratchDiagonal = new Matrix3();\n\nconst jMatrix = new Matrix3();\nconst jMatrixTranspose = new Matrix3();\n\nexport type EigenDecomposition = {\n unitary: Matrix3;\n diagonal: Matrix3;\n};\n\n/**\n * Computes the eigenvectors and eigenvalues of a symmetric matrix.\n *\n * - Returns a diagonal matrix and unitary matrix such that:\n * `matrix = unitary matrix * diagonal matrix * transpose(unitary matrix)`\n * - The values along the diagonal of the diagonal matrix are the eigenvalues. The columns\n * of the unitary matrix are the corresponding eigenvectors.\n * - This routine was created based upon Matrix Computations, 3rd ed., by Golub and Van Loan,\n * section 8.4.3 The Classical Jacobi Algorithm\n *\n * @param matrix The 3x3 matrix to decompose into diagonal and unitary matrix. Expected to be symmetric.\n * @param result Optional object with unitary and diagonal properties which are matrices onto which to store the result.\n * @returns An object with unitary and diagonal properties which are the unitary and diagonal matrices, respectively.\n *\n * @example\n * const a = //... symmetric matrix\n * const result = {\n * unitary : new Matrix3(),\n * diagonal : new Matrix3()\n * };\n * computeEigenDecomposition(a, result);\n *\n * const unitaryTranspose = Matrix3.transpose(result.unitary, new Matrix3());\n * const b = Matrix3.multiply(result.unitary, result.diagonal, new Matrix3());\n * Matrix3.multiply(b, unitaryTranspose, b); // b is now equal to a\n *\n * const lambda = result.diagonal.getColumn(0, new Vector3()).x; // first eigenvalue\n * const v = result.unitary.getColumn(0, new Vector3()); // first eigenvector\n * const c = v.multiplyByScalar(lambda); // equal to v.transformByMatrix3(a)\n */\nexport default function computeEigenDecomposition(\n matrix: number[],\n // @ts-expect-error accept empty object type\n result: EigenDecomposition = {}\n): EigenDecomposition {\n const EIGEN_TOLERANCE = _MathUtils.EPSILON20;\n const EIGEN_MAX_SWEEPS = 10;\n\n let count = 0;\n let sweep = 0;\n\n const unitaryMatrix = scratchUnitary;\n const diagonalMatrix = scratchDiagonal;\n\n unitaryMatrix.identity();\n diagonalMatrix.copy(matrix);\n\n const epsilon = EIGEN_TOLERANCE * computeFrobeniusNorm(diagonalMatrix);\n\n while (sweep < EIGEN_MAX_SWEEPS && offDiagonalFrobeniusNorm(diagonalMatrix) > epsilon) {\n shurDecomposition(diagonalMatrix, jMatrix);\n\n jMatrixTranspose.copy(jMatrix).transpose();\n\n diagonalMatrix.multiplyRight(jMatrix);\n diagonalMatrix.multiplyLeft(jMatrixTranspose);\n unitaryMatrix.multiplyRight(jMatrix);\n\n if (++count > 2) {\n ++sweep;\n count = 0;\n }\n }\n\n result.unitary = unitaryMatrix.toTarget(result.unitary);\n result.diagonal = diagonalMatrix.toTarget(result.diagonal);\n\n return result;\n}\n\nfunction computeFrobeniusNorm(matrix: Matrix3): number {\n let norm = 0.0;\n for (let i = 0; i < 9; ++i) {\n const temp = matrix[i];\n norm += temp * temp;\n }\n return Math.sqrt(norm);\n}\n\nconst rowVal = [1, 0, 0];\nconst colVal = [2, 2, 1];\n\n// Computes the \"off-diagonal\" Frobenius norm.\n// Assumes matrix is symmetric.\nfunction offDiagonalFrobeniusNorm(matrix: Matrix3): number {\n let norm = 0.0;\n for (let i = 0; i < 3; ++i) {\n const temp = matrix[scratchMatrix.getElementIndex(colVal[i], rowVal[i])];\n norm += 2.0 * temp * temp;\n }\n return Math.sqrt(norm);\n}\n\n// The routine takes a matrix, which is assumed to be symmetric, and\n// finds the largest off-diagonal term, and then creates\n// a matrix (result) which can be used to help reduce it\n//\n// This routine was created based upon Matrix Computations, 3rd ed., by Golub and Van Loan,\n// section 8.4.2 The 2by2 Symmetric Schur Decomposition.\n//\n// eslint-disable-next-line max-statements\nfunction shurDecomposition(matrix: Matrix3, result: Matrix3): Matrix3 {\n const tolerance = _MathUtils.EPSILON15;\n\n let maxDiagonal = 0.0;\n let rotAxis = 1;\n\n // find pivot (rotAxis) based on max diagonal of matrix\n for (let i = 0; i < 3; ++i) {\n const temp = Math.abs(matrix[scratchMatrix.getElementIndex(colVal[i], rowVal[i])]);\n if (temp > maxDiagonal) {\n rotAxis = i;\n maxDiagonal = temp;\n }\n }\n\n const p = rowVal[rotAxis];\n const q = colVal[rotAxis];\n\n let c = 1.0;\n let s = 0.0;\n\n if (Math.abs(matrix[scratchMatrix.getElementIndex(q, p)]) > tolerance) {\n const qq = matrix[scratchMatrix.getElementIndex(q, q)];\n const pp = matrix[scratchMatrix.getElementIndex(p, p)];\n const qp = matrix[scratchMatrix.getElementIndex(q, p)];\n\n const tau = (qq - pp) / 2.0 / qp;\n let t;\n\n if (tau < 0.0) {\n t = -1.0 / (-tau + Math.sqrt(1.0 + tau * tau));\n } else {\n t = 1.0 / (tau + Math.sqrt(1.0 + tau * tau));\n }\n\n c = 1.0 / Math.sqrt(1.0 + t * t);\n s = t * c;\n }\n\n // Copy into result\n Matrix3.IDENTITY.to(result);\n result[scratchMatrix.getElementIndex(p, p)] = result[scratchMatrix.getElementIndex(q, q)] = c;\n result[scratchMatrix.getElementIndex(q, p)] = s;\n result[scratchMatrix.getElementIndex(p, q)] = -s;\n\n return result;\n}\n"],"file":"compute-eigen-decomposition.js"}
|
|
@@ -1,114 +0,0 @@
|
|
|
1
|
-
"use strict";
|
|
2
|
-
|
|
3
|
-
var _interopRequireDefault = require("@babel/runtime/helpers/interopRequireDefault");
|
|
4
|
-
|
|
5
|
-
Object.defineProperty(exports, "__esModule", {
|
|
6
|
-
value: true
|
|
7
|
-
});
|
|
8
|
-
exports.default = void 0;
|
|
9
|
-
|
|
10
|
-
var _classCallCheck2 = _interopRequireDefault(require("@babel/runtime/helpers/classCallCheck"));
|
|
11
|
-
|
|
12
|
-
var _createClass2 = _interopRequireDefault(require("@babel/runtime/helpers/createClass"));
|
|
13
|
-
|
|
14
|
-
var _defineProperty2 = _interopRequireDefault(require("@babel/runtime/helpers/defineProperty"));
|
|
15
|
-
|
|
16
|
-
var _core = require("@math.gl/core");
|
|
17
|
-
|
|
18
|
-
var _constants = require("../../constants");
|
|
19
|
-
|
|
20
|
-
var scratchVector = new _core.Vector3();
|
|
21
|
-
var scratchNormal = new _core.Vector3();
|
|
22
|
-
|
|
23
|
-
var AxisAlignedBoundingBox = function () {
|
|
24
|
-
function AxisAlignedBoundingBox() {
|
|
25
|
-
var minimum = arguments.length > 0 && arguments[0] !== undefined ? arguments[0] : [0, 0, 0];
|
|
26
|
-
var maximum = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : [0, 0, 0];
|
|
27
|
-
var center = arguments.length > 2 ? arguments[2] : undefined;
|
|
28
|
-
(0, _classCallCheck2.default)(this, AxisAlignedBoundingBox);
|
|
29
|
-
(0, _defineProperty2.default)(this, "center", void 0);
|
|
30
|
-
(0, _defineProperty2.default)(this, "halfDiagonal", void 0);
|
|
31
|
-
(0, _defineProperty2.default)(this, "minimum", void 0);
|
|
32
|
-
(0, _defineProperty2.default)(this, "maximum", void 0);
|
|
33
|
-
center = center || scratchVector.copy(minimum).add(maximum).scale(0.5);
|
|
34
|
-
this.center = new _core.Vector3(center);
|
|
35
|
-
this.halfDiagonal = new _core.Vector3(maximum).subtract(this.center);
|
|
36
|
-
this.minimum = new _core.Vector3(minimum);
|
|
37
|
-
this.maximum = new _core.Vector3(maximum);
|
|
38
|
-
}
|
|
39
|
-
|
|
40
|
-
(0, _createClass2.default)(AxisAlignedBoundingBox, [{
|
|
41
|
-
key: "clone",
|
|
42
|
-
value: function clone() {
|
|
43
|
-
return new AxisAlignedBoundingBox(this.minimum, this.maximum, this.center);
|
|
44
|
-
}
|
|
45
|
-
}, {
|
|
46
|
-
key: "equals",
|
|
47
|
-
value: function equals(right) {
|
|
48
|
-
return this === right || Boolean(right) && this.minimum.equals(right.minimum) && this.maximum.equals(right.maximum);
|
|
49
|
-
}
|
|
50
|
-
}, {
|
|
51
|
-
key: "transform",
|
|
52
|
-
value: function transform(_transform) {
|
|
53
|
-
this.center.transformAsPoint(_transform);
|
|
54
|
-
this.halfDiagonal.transform(_transform);
|
|
55
|
-
this.minimum.transform(_transform);
|
|
56
|
-
this.maximum.transform(_transform);
|
|
57
|
-
return this;
|
|
58
|
-
}
|
|
59
|
-
}, {
|
|
60
|
-
key: "intersectPlane",
|
|
61
|
-
value: function intersectPlane(plane) {
|
|
62
|
-
var halfDiagonal = this.halfDiagonal;
|
|
63
|
-
var normal = scratchNormal.from(plane.normal);
|
|
64
|
-
var e = halfDiagonal.x * Math.abs(normal.x) + halfDiagonal.y * Math.abs(normal.y) + halfDiagonal.z * Math.abs(normal.z);
|
|
65
|
-
var s = this.center.dot(normal) + plane.distance;
|
|
66
|
-
|
|
67
|
-
if (s - e > 0) {
|
|
68
|
-
return _constants.INTERSECTION.INSIDE;
|
|
69
|
-
}
|
|
70
|
-
|
|
71
|
-
if (s + e < 0) {
|
|
72
|
-
return _constants.INTERSECTION.OUTSIDE;
|
|
73
|
-
}
|
|
74
|
-
|
|
75
|
-
return _constants.INTERSECTION.INTERSECTING;
|
|
76
|
-
}
|
|
77
|
-
}, {
|
|
78
|
-
key: "distanceTo",
|
|
79
|
-
value: function distanceTo(point) {
|
|
80
|
-
return Math.sqrt(this.distanceSquaredTo(point));
|
|
81
|
-
}
|
|
82
|
-
}, {
|
|
83
|
-
key: "distanceSquaredTo",
|
|
84
|
-
value: function distanceSquaredTo(point) {
|
|
85
|
-
var offset = scratchVector.from(point).subtract(this.center);
|
|
86
|
-
var halfDiagonal = this.halfDiagonal;
|
|
87
|
-
var distanceSquared = 0.0;
|
|
88
|
-
var d;
|
|
89
|
-
d = Math.abs(offset.x) - halfDiagonal.x;
|
|
90
|
-
|
|
91
|
-
if (d > 0) {
|
|
92
|
-
distanceSquared += d * d;
|
|
93
|
-
}
|
|
94
|
-
|
|
95
|
-
d = Math.abs(offset.y) - halfDiagonal.y;
|
|
96
|
-
|
|
97
|
-
if (d > 0) {
|
|
98
|
-
distanceSquared += d * d;
|
|
99
|
-
}
|
|
100
|
-
|
|
101
|
-
d = Math.abs(offset.z) - halfDiagonal.z;
|
|
102
|
-
|
|
103
|
-
if (d > 0) {
|
|
104
|
-
distanceSquared += d * d;
|
|
105
|
-
}
|
|
106
|
-
|
|
107
|
-
return distanceSquared;
|
|
108
|
-
}
|
|
109
|
-
}]);
|
|
110
|
-
return AxisAlignedBoundingBox;
|
|
111
|
-
}();
|
|
112
|
-
|
|
113
|
-
exports.default = AxisAlignedBoundingBox;
|
|
114
|
-
//# sourceMappingURL=axis-aligned-bounding-box.js.map
|
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"sources":["../../../../src/lib/bounding-volumes/axis-aligned-bounding-box.ts"],"names":["scratchVector","Vector3","scratchNormal","AxisAlignedBoundingBox","minimum","maximum","center","copy","add","scale","halfDiagonal","subtract","right","Boolean","equals","transform","transformAsPoint","plane","normal","from","e","x","Math","abs","y","z","s","dot","distance","INTERSECTION","INSIDE","OUTSIDE","INTERSECTING","point","sqrt","distanceSquaredTo","offset","distanceSquared","d"],"mappings":";;;;;;;;;;;;;;;AACA;;AAEA;;AAEA,IAAMA,aAAa,GAAG,IAAIC,aAAJ,EAAtB;AACA,IAAMC,aAAa,GAAG,IAAID,aAAJ,EAAtB;;IAQqBE,sB;AAgBnB,oCAIE;AAAA,QAHAC,OAGA,uEAH6B,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAG7B;AAAA,QAFAC,OAEA,uEAF6B,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAE7B;AAAA,QADAC,MACA;AAAA;AAAA;AAAA;AAAA;AAAA;AAEAA,IAAAA,MAAM,GAAGA,MAAM,IAAIN,aAAa,CAACO,IAAd,CAAmBH,OAAnB,EAA4BI,GAA5B,CAAgCH,OAAhC,EAAyCI,KAAzC,CAA+C,GAA/C,CAAnB;AACA,SAAKH,MAAL,GAAc,IAAIL,aAAJ,CAAYK,MAAZ,CAAd;AACA,SAAKI,YAAL,GAAoB,IAAIT,aAAJ,CAAYI,OAAZ,EAAqBM,QAArB,CAA8B,KAAKL,MAAnC,CAApB;AAOA,SAAKF,OAAL,GAAe,IAAIH,aAAJ,CAAYG,OAAZ,CAAf;AAOA,SAAKC,OAAL,GAAe,IAAIJ,aAAJ,CAAYI,OAAZ,CAAf;AACD;;;;WAOD,iBAAgC;AAC9B,aAAO,IAAIF,sBAAJ,CAA2B,KAAKC,OAAhC,EAAyC,KAAKC,OAA9C,EAAuD,KAAKC,MAA5D,CAAP;AACD;;;WASD,gBAAOM,KAAP,EAA+C;AAC7C,aACE,SAASA,KAAT,IACCC,OAAO,CAACD,KAAD,CAAP,IAAkB,KAAKR,OAAL,CAAaU,MAAb,CAAoBF,KAAK,CAACR,OAA1B,CAAlB,IAAwD,KAAKC,OAAL,CAAaS,MAAb,CAAoBF,KAAK,CAACP,OAA1B,CAF3D;AAID;;;WAOD,mBAAUU,UAAV,EAA8C;AAC5C,WAAKT,MAAL,CAAYU,gBAAZ,CAA6BD,UAA7B;AAEA,WAAKL,YAAL,CAAkBK,SAAlB,CAA4BA,UAA5B;AACA,WAAKX,OAAL,CAAaW,SAAb,CAAuBA,UAAvB;AACA,WAAKV,OAAL,CAAaU,SAAb,CAAuBA,UAAvB;AACA,aAAO,IAAP;AACD;;;WAKD,wBAAeE,KAAf,EAAqC;AACnC,UAAOP,YAAP,GAAuB,IAAvB,CAAOA,YAAP;AACA,UAAMQ,MAAM,GAAGhB,aAAa,CAACiB,IAAd,CAAmBF,KAAK,CAACC,MAAzB,CAAf;AACA,UAAME,CAAC,GACLV,YAAY,CAACW,CAAb,GAAiBC,IAAI,CAACC,GAAL,CAASL,MAAM,CAACG,CAAhB,CAAjB,GACAX,YAAY,CAACc,CAAb,GAAiBF,IAAI,CAACC,GAAL,CAASL,MAAM,CAACM,CAAhB,CADjB,GAEAd,YAAY,CAACe,CAAb,GAAiBH,IAAI,CAACC,GAAL,CAASL,MAAM,CAACO,CAAhB,CAHnB;AAIA,UAAMC,CAAC,GAAG,KAAKpB,MAAL,CAAYqB,GAAZ,CAAgBT,MAAhB,IAA0BD,KAAK,CAACW,QAA1C;;AAEA,UAAIF,CAAC,GAAGN,CAAJ,GAAQ,CAAZ,EAAe;AACb,eAAOS,wBAAaC,MAApB;AACD;;AAED,UAAIJ,CAAC,GAAGN,CAAJ,GAAQ,CAAZ,EAAe;AAEb,eAAOS,wBAAaE,OAApB;AACD;;AAED,aAAOF,wBAAaG,YAApB;AACD;;;WAGD,oBAAWC,KAAX,EAA6C;AAC3C,aAAOX,IAAI,CAACY,IAAL,CAAU,KAAKC,iBAAL,CAAuBF,KAAvB,CAAV,CAAP;AACD;;;WAGD,2BAAkBA,KAAlB,EAAoD;AAClD,UAAMG,MAAM,GAAGpC,aAAa,CAACmB,IAAd,CAAmBc,KAAnB,EAA0BtB,QAA1B,CAAmC,KAAKL,MAAxC,CAAf;AACA,UAAOI,YAAP,GAAuB,IAAvB,CAAOA,YAAP;AAEA,UAAI2B,eAAe,GAAG,GAAtB;AACA,UAAIC,CAAJ;AAEAA,MAAAA,CAAC,GAAGhB,IAAI,CAACC,GAAL,CAASa,MAAM,CAACf,CAAhB,IAAqBX,YAAY,CAACW,CAAtC;;AACA,UAAIiB,CAAC,GAAG,CAAR,EAAW;AACTD,QAAAA,eAAe,IAAIC,CAAC,GAAGA,CAAvB;AACD;;AAEDA,MAAAA,CAAC,GAAGhB,IAAI,CAACC,GAAL,CAASa,MAAM,CAACZ,CAAhB,IAAqBd,YAAY,CAACc,CAAtC;;AACA,UAAIc,CAAC,GAAG,CAAR,EAAW;AACTD,QAAAA,eAAe,IAAIC,CAAC,GAAGA,CAAvB;AACD;;AAEDA,MAAAA,CAAC,GAAGhB,IAAI,CAACC,GAAL,CAASa,MAAM,CAACX,CAAhB,IAAqBf,YAAY,CAACe,CAAtC;;AACA,UAAIa,CAAC,GAAG,CAAR,EAAW;AACTD,QAAAA,eAAe,IAAIC,CAAC,GAAGA,CAAvB;AACD;;AAED,aAAOD,eAAP;AACD","sourcesContent":["import {BoundingVolume} from './bounding-volume';\nimport {Vector3} from '@math.gl/core';\nimport Plane from '../plane';\nimport {INTERSECTION} from '../../constants';\n\nconst scratchVector = new Vector3();\nconst scratchNormal = new Vector3();\n\n/**\n * An axis aligned bounding box - aligned with coordinate axes\n * @see BoundingVolume\n * @see BoundingRectangle\n * @see OrientedBoundingBox\n */\nexport default class AxisAlignedBoundingBox implements BoundingVolume {\n /** The center point of the bounding box. */\n readonly center: Vector3;\n /** The positive half diagonal of the bounding box. */\n readonly halfDiagonal: Vector3;\n /** The minimum point defining the bounding box. [0, 0, 0] for empty box */\n readonly minimum: Vector3;\n /** The maximum point defining the bounding box. [0, 0, 0] for empty box */\n readonly maximum: Vector3;\n\n /**\n * Creates an instance of an AxisAlignedBoundingBox from the minimum and maximum points along the x, y, and z axes.\n * @param minimum=[0, 0, 0] The minimum point along the x, y, and z axes.\n * @param maximum=[0, 0, 0] The maximum point along the x, y, and z axes.\n * @param center The center of the box; automatically computed if not supplied.\n */\n constructor(\n minimum: readonly number[] = [0, 0, 0],\n maximum: readonly number[] = [0, 0, 0],\n center?: readonly number[]\n ) {\n // If center was not defined, compute it.\n center = center || scratchVector.copy(minimum).add(maximum).scale(0.5);\n this.center = new Vector3(center);\n this.halfDiagonal = new Vector3(maximum).subtract(this.center);\n\n /**\n * The minimum point defining the bounding box.\n * @type {Vector3}\n * @default {@link 0, 0, 0}\n */\n this.minimum = new Vector3(minimum);\n\n /**\n * The maximum point defining the bounding box.\n * @type {Vector3}\n * @default {@link 0, 0, 0}\n */\n this.maximum = new Vector3(maximum);\n }\n\n /**\n * Duplicates a AxisAlignedBoundingBox instance.\n *\n * @returns {AxisAlignedBoundingBox} A new AxisAlignedBoundingBox instance.\n */\n clone(): AxisAlignedBoundingBox {\n return new AxisAlignedBoundingBox(this.minimum, this.maximum, this.center);\n }\n\n /**\n * Compares the provided AxisAlignedBoundingBox componentwise and returns\n * <code>true</code> if they are equal, <code>false</code> otherwise.\n *\n * @param {AxisAlignedBoundingBox} [right] The second AxisAlignedBoundingBox to compare with.\n * @returns {Boolean} <code>true</code> if left and right are equal, <code>false</code> otherwise.\n */\n equals(right: AxisAlignedBoundingBox): boolean {\n return (\n this === right ||\n (Boolean(right) && this.minimum.equals(right.minimum) && this.maximum.equals(right.maximum))\n );\n }\n\n /**\n * Applies a 4x4 affine transformation matrix to a bounding sphere.\n * @param transform The transformation matrix to apply to the bounding sphere.\n * @returns itself, i.e. the modified BoundingVolume.\n */\n transform(transform: readonly number[]): this {\n this.center.transformAsPoint(transform);\n // TODO - this.halfDiagonal.transformAsVector(transform);\n this.halfDiagonal.transform(transform);\n this.minimum.transform(transform);\n this.maximum.transform(transform);\n return this;\n }\n\n /**\n * Determines which side of a plane a box is located.\n */\n intersectPlane(plane: Plane): number {\n const {halfDiagonal} = this;\n const normal = scratchNormal.from(plane.normal);\n const e =\n halfDiagonal.x * Math.abs(normal.x) +\n halfDiagonal.y * Math.abs(normal.y) +\n halfDiagonal.z * Math.abs(normal.z);\n const s = this.center.dot(normal) + plane.distance; // signed distance from center\n\n if (s - e > 0) {\n return INTERSECTION.INSIDE;\n }\n\n if (s + e < 0) {\n // Not in front because normals point inward\n return INTERSECTION.OUTSIDE;\n }\n\n return INTERSECTION.INTERSECTING;\n }\n\n /** Computes the estimated distance from the closest point on a bounding box to a point. */\n distanceTo(point: readonly number[]): number {\n return Math.sqrt(this.distanceSquaredTo(point));\n }\n\n /** Computes the estimated distance squared from the closest point on a bounding box to a point. */\n distanceSquaredTo(point: readonly number[]): number {\n const offset = scratchVector.from(point).subtract(this.center);\n const {halfDiagonal} = this;\n\n let distanceSquared = 0.0;\n let d;\n\n d = Math.abs(offset.x) - halfDiagonal.x;\n if (d > 0) {\n distanceSquared += d * d;\n }\n\n d = Math.abs(offset.y) - halfDiagonal.y;\n if (d > 0) {\n distanceSquared += d * d;\n }\n\n d = Math.abs(offset.z) - halfDiagonal.z;\n if (d > 0) {\n distanceSquared += d * d;\n }\n\n return distanceSquared;\n }\n}\n"],"file":"axis-aligned-bounding-box.js"}
|
|
@@ -1,148 +0,0 @@
|
|
|
1
|
-
"use strict";
|
|
2
|
-
|
|
3
|
-
var _interopRequireDefault = require("@babel/runtime/helpers/interopRequireDefault");
|
|
4
|
-
|
|
5
|
-
var _typeof = require("@babel/runtime/helpers/typeof");
|
|
6
|
-
|
|
7
|
-
Object.defineProperty(exports, "__esModule", {
|
|
8
|
-
value: true
|
|
9
|
-
});
|
|
10
|
-
exports.default = void 0;
|
|
11
|
-
|
|
12
|
-
var _classCallCheck2 = _interopRequireDefault(require("@babel/runtime/helpers/classCallCheck"));
|
|
13
|
-
|
|
14
|
-
var _createClass2 = _interopRequireDefault(require("@babel/runtime/helpers/createClass"));
|
|
15
|
-
|
|
16
|
-
var _defineProperty2 = _interopRequireDefault(require("@babel/runtime/helpers/defineProperty"));
|
|
17
|
-
|
|
18
|
-
var _core = require("@math.gl/core");
|
|
19
|
-
|
|
20
|
-
var mat4 = _interopRequireWildcard(require("gl-matrix/mat4"));
|
|
21
|
-
|
|
22
|
-
var _constants = require("../../constants");
|
|
23
|
-
|
|
24
|
-
function _getRequireWildcardCache(nodeInterop) { if (typeof WeakMap !== "function") return null; var cacheBabelInterop = new WeakMap(); var cacheNodeInterop = new WeakMap(); return (_getRequireWildcardCache = function _getRequireWildcardCache(nodeInterop) { return nodeInterop ? cacheNodeInterop : cacheBabelInterop; })(nodeInterop); }
|
|
25
|
-
|
|
26
|
-
function _interopRequireWildcard(obj, nodeInterop) { if (!nodeInterop && obj && obj.__esModule) { return obj; } if (obj === null || _typeof(obj) !== "object" && typeof obj !== "function") { return { default: obj }; } var cache = _getRequireWildcardCache(nodeInterop); if (cache && cache.has(obj)) { return cache.get(obj); } var newObj = {}; var hasPropertyDescriptor = Object.defineProperty && Object.getOwnPropertyDescriptor; for (var key in obj) { if (key !== "default" && Object.prototype.hasOwnProperty.call(obj, key)) { var desc = hasPropertyDescriptor ? Object.getOwnPropertyDescriptor(obj, key) : null; if (desc && (desc.get || desc.set)) { Object.defineProperty(newObj, key, desc); } else { newObj[key] = obj[key]; } } } newObj.default = obj; if (cache) { cache.set(obj, newObj); } return newObj; }
|
|
27
|
-
|
|
28
|
-
var scratchVector = new _core.Vector3();
|
|
29
|
-
var scratchVector2 = new _core.Vector3();
|
|
30
|
-
|
|
31
|
-
var BoundingSphere = function () {
|
|
32
|
-
function BoundingSphere() {
|
|
33
|
-
var center = arguments.length > 0 && arguments[0] !== undefined ? arguments[0] : [0, 0, 0];
|
|
34
|
-
var radius = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : 0.0;
|
|
35
|
-
(0, _classCallCheck2.default)(this, BoundingSphere);
|
|
36
|
-
(0, _defineProperty2.default)(this, "center", void 0);
|
|
37
|
-
(0, _defineProperty2.default)(this, "radius", void 0);
|
|
38
|
-
this.radius = -0;
|
|
39
|
-
this.center = new _core.Vector3();
|
|
40
|
-
this.fromCenterRadius(center, radius);
|
|
41
|
-
}
|
|
42
|
-
|
|
43
|
-
(0, _createClass2.default)(BoundingSphere, [{
|
|
44
|
-
key: "fromCenterRadius",
|
|
45
|
-
value: function fromCenterRadius(center, radius) {
|
|
46
|
-
this.center.from(center);
|
|
47
|
-
this.radius = radius;
|
|
48
|
-
return this;
|
|
49
|
-
}
|
|
50
|
-
}, {
|
|
51
|
-
key: "fromCornerPoints",
|
|
52
|
-
value: function fromCornerPoints(corner, oppositeCorner) {
|
|
53
|
-
oppositeCorner = scratchVector.from(oppositeCorner);
|
|
54
|
-
this.center = new _core.Vector3().from(corner).add(oppositeCorner).scale(0.5);
|
|
55
|
-
this.radius = this.center.distance(oppositeCorner);
|
|
56
|
-
return this;
|
|
57
|
-
}
|
|
58
|
-
}, {
|
|
59
|
-
key: "equals",
|
|
60
|
-
value: function equals(right) {
|
|
61
|
-
return this === right || Boolean(right) && this.center.equals(right.center) && this.radius === right.radius;
|
|
62
|
-
}
|
|
63
|
-
}, {
|
|
64
|
-
key: "clone",
|
|
65
|
-
value: function clone() {
|
|
66
|
-
return new BoundingSphere(this.center, this.radius);
|
|
67
|
-
}
|
|
68
|
-
}, {
|
|
69
|
-
key: "union",
|
|
70
|
-
value: function union(boundingSphere) {
|
|
71
|
-
var leftCenter = this.center;
|
|
72
|
-
var leftRadius = this.radius;
|
|
73
|
-
var rightCenter = boundingSphere.center;
|
|
74
|
-
var rightRadius = boundingSphere.radius;
|
|
75
|
-
var toRightCenter = scratchVector.copy(rightCenter).subtract(leftCenter);
|
|
76
|
-
var centerSeparation = toRightCenter.magnitude();
|
|
77
|
-
|
|
78
|
-
if (leftRadius >= centerSeparation + rightRadius) {
|
|
79
|
-
return this.clone();
|
|
80
|
-
}
|
|
81
|
-
|
|
82
|
-
if (rightRadius >= centerSeparation + leftRadius) {
|
|
83
|
-
return boundingSphere.clone();
|
|
84
|
-
}
|
|
85
|
-
|
|
86
|
-
var halfDistanceBetweenTangentPoints = (leftRadius + centerSeparation + rightRadius) * 0.5;
|
|
87
|
-
scratchVector2.copy(toRightCenter).scale((-leftRadius + halfDistanceBetweenTangentPoints) / centerSeparation).add(leftCenter);
|
|
88
|
-
this.center.copy(scratchVector2);
|
|
89
|
-
this.radius = halfDistanceBetweenTangentPoints;
|
|
90
|
-
return this;
|
|
91
|
-
}
|
|
92
|
-
}, {
|
|
93
|
-
key: "expand",
|
|
94
|
-
value: function expand(point) {
|
|
95
|
-
var scratchPoint = scratchVector.from(point);
|
|
96
|
-
var radius = scratchPoint.subtract(this.center).magnitude();
|
|
97
|
-
|
|
98
|
-
if (radius > this.radius) {
|
|
99
|
-
this.radius = radius;
|
|
100
|
-
}
|
|
101
|
-
|
|
102
|
-
return this;
|
|
103
|
-
}
|
|
104
|
-
}, {
|
|
105
|
-
key: "transform",
|
|
106
|
-
value: function transform(_transform) {
|
|
107
|
-
this.center.transform(_transform);
|
|
108
|
-
var scale = mat4.getScaling(scratchVector, _transform);
|
|
109
|
-
this.radius = Math.max(scale[0], Math.max(scale[1], scale[2])) * this.radius;
|
|
110
|
-
return this;
|
|
111
|
-
}
|
|
112
|
-
}, {
|
|
113
|
-
key: "distanceSquaredTo",
|
|
114
|
-
value: function distanceSquaredTo(point) {
|
|
115
|
-
var d = this.distanceTo(point);
|
|
116
|
-
return d * d;
|
|
117
|
-
}
|
|
118
|
-
}, {
|
|
119
|
-
key: "distanceTo",
|
|
120
|
-
value: function distanceTo(point) {
|
|
121
|
-
var scratchPoint = scratchVector.from(point);
|
|
122
|
-
var delta = scratchPoint.subtract(this.center);
|
|
123
|
-
return Math.max(0, delta.len() - this.radius);
|
|
124
|
-
}
|
|
125
|
-
}, {
|
|
126
|
-
key: "intersectPlane",
|
|
127
|
-
value: function intersectPlane(plane) {
|
|
128
|
-
var center = this.center;
|
|
129
|
-
var radius = this.radius;
|
|
130
|
-
var normal = plane.normal;
|
|
131
|
-
var distanceToPlane = normal.dot(center) + plane.distance;
|
|
132
|
-
|
|
133
|
-
if (distanceToPlane < -radius) {
|
|
134
|
-
return _constants.INTERSECTION.OUTSIDE;
|
|
135
|
-
}
|
|
136
|
-
|
|
137
|
-
if (distanceToPlane < radius) {
|
|
138
|
-
return _constants.INTERSECTION.INTERSECTING;
|
|
139
|
-
}
|
|
140
|
-
|
|
141
|
-
return _constants.INTERSECTION.INSIDE;
|
|
142
|
-
}
|
|
143
|
-
}]);
|
|
144
|
-
return BoundingSphere;
|
|
145
|
-
}();
|
|
146
|
-
|
|
147
|
-
exports.default = BoundingSphere;
|
|
148
|
-
//# sourceMappingURL=bounding-sphere.js.map
|
|
@@ -1 +0,0 @@
|
|
|
1
|
-
{"version":3,"sources":["../../../../src/lib/bounding-volumes/bounding-sphere.ts"],"names":["scratchVector","Vector3","scratchVector2","BoundingSphere","center","radius","fromCenterRadius","from","corner","oppositeCorner","add","scale","distance","right","Boolean","equals","boundingSphere","leftCenter","leftRadius","rightCenter","rightRadius","toRightCenter","copy","subtract","centerSeparation","magnitude","clone","halfDistanceBetweenTangentPoints","point","scratchPoint","transform","mat4","getScaling","Math","max","d","distanceTo","delta","len","plane","normal","distanceToPlane","dot","INTERSECTION","OUTSIDE","INTERSECTING","INSIDE"],"mappings":";;;;;;;;;;;;;;;;;AAGA;;AACA;;AACA;;;;;;AAIA,IAAMA,aAAa,GAAG,IAAIC,aAAJ,EAAtB;AACA,IAAMC,cAAc,GAAG,IAAID,aAAJ,EAAvB;;IAGqBE,c;AAKnB,4BAAyE;AAAA,QAA7DC,MAA6D,uEAAjC,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAiC;AAAA,QAAtBC,MAAsB,uEAAL,GAAK;AAAA;AAAA;AAAA;AACvE,SAAKA,MAAL,GAAc,CAAC,CAAf;AACA,SAAKD,MAAL,GAAc,IAAIH,aAAJ,EAAd;AACA,SAAKK,gBAAL,CAAsBF,MAAtB,EAA8BC,MAA9B;AACD;;;;WAGD,0BAAiBD,MAAjB,EAA4CC,MAA5C,EAAkE;AAChE,WAAKD,MAAL,CAAYG,IAAZ,CAAiBH,MAAjB;AACA,WAAKC,MAAL,GAAcA,MAAd;AACA,aAAO,IAAP;AACD;;;WAMD,0BAAiBG,MAAjB,EAA4CC,cAA5C,EAAqF;AACnFA,MAAAA,cAAc,GAAGT,aAAa,CAACO,IAAd,CAAmBE,cAAnB,CAAjB;AACA,WAAKL,MAAL,GAAc,IAAIH,aAAJ,GAAcM,IAAd,CAAmBC,MAAnB,EAA2BE,GAA3B,CAA+BD,cAA/B,EAA+CE,KAA/C,CAAqD,GAArD,CAAd;AACA,WAAKN,MAAL,GAAc,KAAKD,MAAL,CAAYQ,QAAZ,CAAqBH,cAArB,CAAd;AACA,aAAO,IAAP;AACD;;;WAGD,gBAAOI,KAAP,EAAuC;AACrC,aACE,SAASA,KAAT,IACCC,OAAO,CAACD,KAAD,CAAP,IAAkB,KAAKT,MAAL,CAAYW,MAAZ,CAAmBF,KAAK,CAACT,MAAzB,CAAlB,IAAsD,KAAKC,MAAL,KAAgBQ,KAAK,CAACR,MAF/E;AAID;;;WAGD,iBAAwB;AACtB,aAAO,IAAIF,cAAJ,CAAmB,KAAKC,MAAxB,EAAgC,KAAKC,MAArC,CAAP;AACD;;;WAGD,eAAMW,cAAN,EAAsD;AACpD,UAAMC,UAAU,GAAG,KAAKb,MAAxB;AACA,UAAMc,UAAU,GAAG,KAAKb,MAAxB;AACA,UAAMc,WAAW,GAAGH,cAAc,CAACZ,MAAnC;AACA,UAAMgB,WAAW,GAAGJ,cAAc,CAACX,MAAnC;AAEA,UAAMgB,aAAa,GAAGrB,aAAa,CAACsB,IAAd,CAAmBH,WAAnB,EAAgCI,QAAhC,CAAyCN,UAAzC,CAAtB;AACA,UAAMO,gBAAgB,GAAGH,aAAa,CAACI,SAAd,EAAzB;;AAEA,UAAIP,UAAU,IAAIM,gBAAgB,GAAGJ,WAArC,EAAkD;AAEhD,eAAO,KAAKM,KAAL,EAAP;AACD;;AAED,UAAIN,WAAW,IAAII,gBAAgB,GAAGN,UAAtC,EAAkD;AAEhD,eAAOF,cAAc,CAACU,KAAf,EAAP;AACD;;AAGD,UAAMC,gCAAgC,GAAG,CAACT,UAAU,GAAGM,gBAAb,GAAgCJ,WAAjC,IAAgD,GAAzF;AAGAlB,MAAAA,cAAc,CACXoB,IADH,CACQD,aADR,EAEGV,KAFH,CAES,CAAC,CAACO,UAAD,GAAcS,gCAAf,IAAmDH,gBAF5D,EAGGd,GAHH,CAGOO,UAHP;AAKA,WAAKb,MAAL,CAAYkB,IAAZ,CAAiBpB,cAAjB;AACA,WAAKG,MAAL,GAAcsB,gCAAd;AAEA,aAAO,IAAP;AACD;;;WAGD,gBAAOC,KAAP,EAAuC;AACrC,UAAMC,YAAY,GAAG7B,aAAa,CAACO,IAAd,CAAmBqB,KAAnB,CAArB;AACA,UAAMvB,MAAM,GAAGwB,YAAY,CAACN,QAAb,CAAsB,KAAKnB,MAA3B,EAAmCqB,SAAnC,EAAf;;AACA,UAAIpB,MAAM,GAAG,KAAKA,MAAlB,EAA0B;AACxB,aAAKA,MAAL,GAAcA,MAAd;AACD;;AACD,aAAO,IAAP;AACD;;;WAUD,mBAAUyB,UAAV,EAA8C;AAC5C,WAAK1B,MAAL,CAAY0B,SAAZ,CAAsBA,UAAtB;AACA,UAAMnB,KAAK,GAAGoB,IAAI,CAACC,UAAL,CAAgBhC,aAAhB,EAA+B8B,UAA/B,CAAd;AACA,WAAKzB,MAAL,GAAc4B,IAAI,CAACC,GAAL,CAASvB,KAAK,CAAC,CAAD,CAAd,EAAmBsB,IAAI,CAACC,GAAL,CAASvB,KAAK,CAAC,CAAD,CAAd,EAAmBA,KAAK,CAAC,CAAD,CAAxB,CAAnB,IAAmD,KAAKN,MAAtE;AACA,aAAO,IAAP;AACD;;;WAGD,2BAAkBuB,KAAlB,EAAyD;AACvD,UAAMO,CAAC,GAAG,KAAKC,UAAL,CAAgBR,KAAhB,CAAV;AACA,aAAOO,CAAC,GAAGA,CAAX;AACD;;;WAGD,oBAAWP,KAAX,EAAkD;AAChD,UAAMC,YAAY,GAAG7B,aAAa,CAACO,IAAd,CAAmBqB,KAAnB,CAArB;AACA,UAAMS,KAAK,GAAGR,YAAY,CAACN,QAAb,CAAsB,KAAKnB,MAA3B,CAAd;AACA,aAAO6B,IAAI,CAACC,GAAL,CAAS,CAAT,EAAYG,KAAK,CAACC,GAAN,KAAc,KAAKjC,MAA/B,CAAP;AACD;;;WAGD,wBAAekC,KAAf,EAAqC;AACnC,UAAMnC,MAAM,GAAG,KAAKA,MAApB;AACA,UAAMC,MAAM,GAAG,KAAKA,MAApB;AACA,UAAMmC,MAAM,GAAGD,KAAK,CAACC,MAArB;AACA,UAAMC,eAAe,GAAGD,MAAM,CAACE,GAAP,CAAWtC,MAAX,IAAqBmC,KAAK,CAAC3B,QAAnD;;AAGA,UAAI6B,eAAe,GAAG,CAACpC,MAAvB,EAA+B;AAC7B,eAAOsC,wBAAaC,OAApB;AACD;;AAED,UAAIH,eAAe,GAAGpC,MAAtB,EAA8B;AAC5B,eAAOsC,wBAAaE,YAApB;AACD;;AAED,aAAOF,wBAAaG,MAApB;AACD","sourcesContent":["// This file is derived from the Cesium math library under Apache 2 license\n// See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md\n\nimport {NumericArray, Vector3} from '@math.gl/core';\nimport * as mat4 from 'gl-matrix/mat4';\nimport {INTERSECTION} from '../../constants';\nimport {BoundingVolume} from './bounding-volume';\nimport Plane from '../plane';\n\nconst scratchVector = new Vector3();\nconst scratchVector2 = new Vector3();\n\n/** A BoundingSphere */\nexport default class BoundingSphere implements BoundingVolume {\n center: Vector3;\n radius: number;\n\n /** Creates a bounding sphere */\n constructor(center: readonly number[] = [0, 0, 0], radius: number = 0.0) {\n this.radius = -0;\n this.center = new Vector3();\n this.fromCenterRadius(center, radius);\n }\n\n /** Sets the bounding sphere from `center` and `radius`. */\n fromCenterRadius(center: readonly number[], radius: number): this {\n this.center.from(center);\n this.radius = radius;\n return this;\n }\n\n /**\n * Computes a bounding sphere from the corner points of an axis-aligned bounding box. The sphere\n * tightly and fully encompasses the box.\n */\n fromCornerPoints(corner: readonly number[], oppositeCorner: readonly number[]): this {\n oppositeCorner = scratchVector.from(oppositeCorner);\n this.center = new Vector3().from(corner).add(oppositeCorner).scale(0.5);\n this.radius = this.center.distance(oppositeCorner);\n return this;\n }\n\n /** Compares the provided BoundingSphere component wise */\n equals(right: BoundingSphere): boolean {\n return (\n this === right ||\n (Boolean(right) && this.center.equals(right.center) && this.radius === right.radius)\n );\n }\n\n /** Duplicates a BoundingSphere instance. */\n clone(): BoundingSphere {\n return new BoundingSphere(this.center, this.radius);\n }\n\n /** Computes a bounding sphere that contains both the left and right bounding spheres. */\n union(boundingSphere: BoundingSphere): BoundingSphere {\n const leftCenter = this.center;\n const leftRadius = this.radius;\n const rightCenter = boundingSphere.center;\n const rightRadius = boundingSphere.radius;\n\n const toRightCenter = scratchVector.copy(rightCenter).subtract(leftCenter);\n const centerSeparation = toRightCenter.magnitude();\n\n if (leftRadius >= centerSeparation + rightRadius) {\n // Left sphere wins.\n return this.clone();\n }\n\n if (rightRadius >= centerSeparation + leftRadius) {\n // Right sphere wins.\n return boundingSphere.clone();\n }\n\n // There are two tangent points, one on far side of each sphere.\n const halfDistanceBetweenTangentPoints = (leftRadius + centerSeparation + rightRadius) * 0.5;\n\n // Compute the center point halfway between the two tangent points.\n scratchVector2\n .copy(toRightCenter)\n .scale((-leftRadius + halfDistanceBetweenTangentPoints) / centerSeparation)\n .add(leftCenter);\n\n this.center.copy(scratchVector2);\n this.radius = halfDistanceBetweenTangentPoints;\n\n return this;\n }\n\n /** Computes a bounding sphere by enlarging the provided sphere to contain the provided point. */\n expand(point: readonly number[]): this {\n const scratchPoint = scratchVector.from(point);\n const radius = scratchPoint.subtract(this.center).magnitude();\n if (radius > this.radius) {\n this.radius = radius;\n }\n return this;\n }\n\n // BoundingVolume interface\n\n /**\n * Applies a 4x4 affine transformation matrix to a bounding sphere.\n * @param sphere The bounding sphere to apply the transformation to.\n * @param transform The transformation matrix to apply to the bounding sphere.\n * @returns self.\n */\n transform(transform: readonly number[]): this {\n this.center.transform(transform);\n const scale = mat4.getScaling(scratchVector, transform);\n this.radius = Math.max(scale[0], Math.max(scale[1], scale[2])) * this.radius;\n return this;\n }\n\n /** Computes the estimated distance squared from the closest point on a bounding sphere to a point. */\n distanceSquaredTo(point: Readonly<NumericArray>): number {\n const d = this.distanceTo(point);\n return d * d;\n }\n\n /** Computes the estimated distance from the closest point on a bounding sphere to a point. */\n distanceTo(point: Readonly<NumericArray>): number {\n const scratchPoint = scratchVector.from(point);\n const delta = scratchPoint.subtract(this.center);\n return Math.max(0, delta.len() - this.radius);\n }\n\n /** Determines which side of a plane a sphere is located. */\n intersectPlane(plane: Plane): number {\n const center = this.center;\n const radius = this.radius;\n const normal = plane.normal;\n const distanceToPlane = normal.dot(center) + plane.distance;\n\n // The center point is negative side of the plane normal\n if (distanceToPlane < -radius) {\n return INTERSECTION.OUTSIDE;\n }\n // The center point is positive side of the plane, but radius extends beyond it; partial overlap\n if (distanceToPlane < radius) {\n return INTERSECTION.INTERSECTING;\n }\n // The center point and radius is positive side of the plane\n return INTERSECTION.INSIDE;\n }\n}\n"],"file":"bounding-sphere.js"}
|