@math.gl/geospatial 3.5.6 → 3.6.0-alpha.2

package/dist/es5/ellipsoid/helpers/ellipsoid-transform.js.map

@@ -1 +1 @@
-{"version":3,"sources":["../../../../src/ellipsoid/helpers/ellipsoid-transform.js"],"names":["EPSILON14","scratchOrigin","Vector3","VECTOR_PRODUCT_LOCAL_FRAME","up","south","north","west","east","down","degeneratePositionLocalFrame","scratchAxisVectors","scratchVector1","scratchVector2","scratchVector3","localFrameToFixedFrame","ellipsoid","firstAxis","secondAxis","thirdAxis","cartesianOrigin","result","thirdAxisInferred","firstAxisVector","secondAxisVector","thirdAxisVector","origin","copy","atPole","x","y","sign","Math","z","fromArray","scale","set","normalize","geodeticSurfaceNormal","cross"],"mappings":";;;;;;;AAAA;;AAEA,IAAMA,SAAS,GAAG,KAAlB;AAEA,IAAMC,aAAa,GAAG,IAAIC,aAAJ,EAAtB;AAGA,IAAMC,0BAA0B,GAAG;AACjCC,EAAAA,EAAE,EAAE;AACFC,IAAAA,KAAK,EAAE,MADL;AAEFC,IAAAA,KAAK,EAAE,MAFL;AAGFC,IAAAA,IAAI,EAAE,OAHJ;AAIFC,IAAAA,IAAI,EAAE;AAJJ,GAD6B;AAOjCC,EAAAA,IAAI,EAAE;AACJJ,IAAAA,KAAK,EAAE,MADH;AAEJC,IAAAA,KAAK,EAAE,MAFH;AAGJC,IAAAA,IAAI,EAAE,OAHF;AAIJC,IAAAA,IAAI,EAAE;AAJF,GAP2B;AAajCH,EAAAA,KAAK,EAAE;AACLD,IAAAA,EAAE,EAAE,MADC;AAELK,IAAAA,IAAI,EAAE,MAFD;AAGLF,IAAAA,IAAI,EAAE,MAHD;AAILC,IAAAA,IAAI,EAAE;AAJD,GAb0B;AAmBjCF,EAAAA,KAAK,EAAE;AACLF,IAAAA,EAAE,EAAE,MADC;AAELK,IAAAA,IAAI,EAAE,MAFD;AAGLF,IAAAA,IAAI,EAAE,IAHD;AAILC,IAAAA,IAAI,EAAE;AAJD,GAnB0B;AAyBjCD,EAAAA,IAAI,EAAE;AACJH,IAAAA,EAAE,EAAE,OADA;AAEJK,IAAAA,IAAI,EAAE,OAFF;AAGJH,IAAAA,KAAK,EAAE,MAHH;AAIJD,IAAAA,KAAK,EAAE;AAJH,GAzB2B;AA+BjCG,EAAAA,IAAI,EAAE;AACJJ,IAAAA,EAAE,EAAE,OADA;AAEJK,IAAAA,IAAI,EAAE,OAFF;AAGJH,IAAAA,KAAK,EAAE,IAHH;AAIJD,IAAAA,KAAK,EAAE;AAJH;AA/B2B,CAAnC;AAuCA,IAAMK,4BAA4B,GAAG;AACnCJ,EAAAA,KAAK,EAAE,CAAC,CAAC,CAAF,EAAK,CAAL,EAAQ,CAAR,CAD4B;AAEnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAF6B;AAGnCJ,EAAAA,EAAE,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAH+B;AAInCC,EAAAA,KAAK,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAJ4B;AAKnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAC,CAAL,EAAQ,CAAR,CAL6B;AAMnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAC,CAAR;AAN6B,CAArC;AASA,IAAME,kBAAkB,GAAG;AACzBH,EAAAA,IAAI,EAAE,IAAIN,aAAJ,EADmB;AAEzBI,EAAAA,KAAK,EAAE,IAAIJ,aAAJ,EAFkB;AAGzBE,EAAAA,EAAE,EAAE,IAAIF,aAAJ,EAHqB;AAIzBK,EAAAA,IAAI,EAAE,IAAIL,aAAJ,EAJmB;AAKzBG,EAAAA,KAAK,EAAE,IAAIH,aAAJ,EALkB;AAMzBO,EAAAA,IAAI,EAAE,IAAIP,aAAJ;AANmB,CAA3B;AASA,IAAMU,cAAc,GAAG,IAAIV,aAAJ,EAAvB;AACA,IAAMW,cAAc,GAAG,IAAIX,aAAJ,EAAvB;AACA,IAAMY,cAAc,GAAG,IAAIZ,aAAJ,EAAvB;;AAKe,SAASa,sBAAT,CACbC,SADa,EAEbC,SAFa,EAGbC,UAHa,EAIbC,SAJa,EAKbC,eALa,EAMbC,MANa,EAOb;AACA,MAAMC,iBAAiB,GACrBnB,0BAA0B,CAACc,SAAD,CAA1B,IAAyCd,0BAA0B,CAACc,SAAD,CAA1B,CAAsCC,UAAtC,CAD3C;AAGA,oBAAOI,iBAAiB,KAAK,CAACH,SAAD,IAAcA,SAAS,KAAKG,iBAAjC,CAAxB;AAEA,MAAIC,eAAJ;AACA,MAAIC,gBAAJ;AACA,MAAIC,eAAJ;AAEA,MAAMC,MAAM,GAAGzB,aAAa,CAAC0B,IAAd,CAAmBP,eAAnB,CAAf;AAGA,MAAMQ,MAAM,GAAG,kBAAcF,MAAM,CAACG,CAArB,EAAwB,GAAxB,EAA6B7B,SAA7B,KAA2C,kBAAc0B,MAAM,CAACI,CAArB,EAAwB,GAAxB,EAA6B9B,SAA7B,CAA1D;;AAEA,MAAI4B,MAAJ,EAAY;AAEV,QAAMG,IAAI,GAAGC,IAAI,CAACD,IAAL,CAAUL,MAAM,CAACO,CAAjB,CAAb;AAEAV,IAAAA,eAAe,GAAGX,cAAc,CAACsB,SAAf,CAAyBxB,4BAA4B,CAACO,SAAD,CAArD,CAAlB;;AACA,QAAIA,SAAS,KAAK,MAAd,IAAwBA,SAAS,KAAK,MAA1C,EAAkD;AAChDM,MAAAA,eAAe,CAACY,KAAhB,CAAsBJ,IAAtB;AACD;;AAEDP,IAAAA,gBAAgB,GAAGX,cAAc,CAACqB,SAAf,CAAyBxB,4BAA4B,CAACQ,UAAD,CAArD,CAAnB;;AACA,QAAIA,UAAU,KAAK,MAAf,IAAyBA,UAAU,KAAK,MAA5C,EAAoD;AAClDM,MAAAA,gBAAgB,CAACW,KAAjB,CAAuBJ,IAAvB;AACD;;AAEDN,IAAAA,eAAe,GAAGX,cAAc,CAACoB,SAAf,CAAyBxB,4BAA4B,CAACS,SAAD,CAArD,CAAlB;;AACA,QAAIA,SAAS,KAAK,MAAd,IAAwBA,SAAS,KAAK,MAA1C,EAAkD;AAChDM,MAAAA,eAAe,CAACU,KAAhB,CAAsBJ,IAAtB;AACD;AACF,GAlBD,MAkBO;AAEL,QAAO3B,EAAP,GAA0BO,kBAA1B,CAAOP,EAAP;AAAA,QAAWI,IAAX,GAA0BG,kBAA1B,CAAWH,IAAX;AAAA,QAAiBF,KAAjB,GAA0BK,kBAA1B,CAAiBL,KAAjB;AAEAE,IAAAA,IAAI,CAAC4B,GAAL,CAAS,CAACV,MAAM,CAACI,CAAjB,EAAoBJ,MAAM,CAACG,CAA3B,EAA8B,GAA9B,EAAmCQ,SAAnC;AACArB,IAAAA,SAAS,CAACsB,qBAAV,CAAgCZ,MAAhC,EAAwCtB,EAAxC;AACAE,IAAAA,KAAK,CAACqB,IAAN,CAAWvB,EAAX,EAAemC,KAAf,CAAqB/B,IAArB;AAEA,QAAOC,IAAP,GAA4BE,kBAA5B,CAAOF,IAAP;AAAA,QAAaF,IAAb,GAA4BI,kBAA5B,CAAaJ,IAAb;AAAA,QAAmBF,KAAnB,GAA4BM,kBAA5B,CAAmBN,KAAnB;AAEAI,IAAAA,IAAI,CAACkB,IAAL,CAAUvB,EAAV,EAAc+B,KAAd,CAAoB,CAAC,CAArB;AACA5B,IAAAA,IAAI,CAACoB,IAAL,CAAUnB,IAAV,EAAgB2B,KAAhB,CAAsB,CAAC,CAAvB;AACA9B,IAAAA,KAAK,CAACsB,IAAN,CAAWrB,KAAX,EAAkB6B,KAAlB,CAAwB,CAAC,CAAzB;AAGAZ,IAAAA,eAAe,GAAGZ,kBAAkB,CAACM,SAAD,CAApC;AACAO,IAAAA,gBAAgB,GAAGb,kBAAkB,CAACO,UAAD,CAArC;AACAO,IAAAA,eAAe,GAAGd,kBAAkB,CAACQ,SAAD,CAApC;AACD;;AAGDE,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACM,CAA5B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACO,CAA5B;AACAT,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACU,CAA5B;AACAZ,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAY,GAAZ;AACAA,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACK,CAA7B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACM,CAA7B;AACAT,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACS,CAA7B;AACAZ,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAY,GAAZ;AACAA,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYI,eAAe,CAACI,CAA5B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYI,eAAe,CAACK,CAA5B;AACAT,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaI,eAAe,CAACQ,CAA7B;AACAZ,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAa,GAAb;AACAA,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACG,CAApB;AACAR,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACI,CAApB;AACAT,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACO,CAApB;AACAZ,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAa,GAAb;AACA,SAAOA,MAAP;AACD","sourcesContent":["import {Vector3, assert, equals as equalsEpsilon} from '@math.gl/core';\n\nconst EPSILON14 = 1e-14;\n\nconst scratchOrigin = new Vector3();\n\n// Caclulate third axis from given two axii\nconst VECTOR_PRODUCT_LOCAL_FRAME = {\n  up: {\n    south: 'east',\n    north: 'west',\n    west: 'south',\n    east: 'north'\n  },\n  down: {\n    south: 'west',\n    north: 'east',\n    west: 'north',\n    east: 'south'\n  },\n  south: {\n    up: 'west',\n    down: 'east',\n    west: 'down',\n    east: 'up'\n  },\n  north: {\n    up: 'east',\n    down: 'west',\n    west: 'up',\n    east: 'down'\n  },\n  west: {\n    up: 'north',\n    down: 'south',\n    north: 'down',\n    south: 'up'\n  },\n  east: {\n    up: 'south',\n    down: 'north',\n    north: 'up',\n    south: 'down'\n  }\n};\n\nconst degeneratePositionLocalFrame = {\n  north: [-1, 0, 0],\n  east: [0, 1, 0],\n  up: [0, 0, 1],\n  south: [1, 0, 0],\n  west: [0, -1, 0],\n  down: [0, 0, -1]\n};\n\nconst scratchAxisVectors = {\n  east: new Vector3(),\n  north: new Vector3(),\n  up: new Vector3(),\n  west: new Vector3(),\n  south: new Vector3(),\n  down: new Vector3()\n};\n\nconst scratchVector1 = new Vector3();\nconst scratchVector2 = new Vector3();\nconst scratchVector3 = new Vector3();\n\n// Computes a 4x4 transformation matrix from a reference frame\n// centered at the provided origin to the provided ellipsoid's fixed reference frame.\n// eslint-disable-next-line max-statements, max-params, complexity\nexport default function localFrameToFixedFrame(\n  ellipsoid,\n  firstAxis,\n  secondAxis,\n  thirdAxis,\n  cartesianOrigin,\n  result\n) {\n  const thirdAxisInferred =\n    VECTOR_PRODUCT_LOCAL_FRAME[firstAxis] && VECTOR_PRODUCT_LOCAL_FRAME[firstAxis][secondAxis];\n  // firstAxis and secondAxis must be east, north, up, west, south or down.');\n  assert(thirdAxisInferred && (!thirdAxis || thirdAxis === thirdAxisInferred));\n\n  let firstAxisVector;\n  let secondAxisVector;\n  let thirdAxisVector;\n\n  const origin = scratchOrigin.copy(cartesianOrigin);\n\n  // If x and y are zero, assume origin is at a pole, which is a special case.\n  const atPole = equalsEpsilon(origin.x, 0.0, EPSILON14) && equalsEpsilon(origin.y, 0.0, EPSILON14);\n\n  if (atPole) {\n    // Look up axis value and adjust\n    const sign = Math.sign(origin.z);\n\n    firstAxisVector = scratchVector1.fromArray(degeneratePositionLocalFrame[firstAxis]);\n    if (firstAxis !== 'east' && firstAxis !== 'west') {\n      firstAxisVector.scale(sign);\n    }\n\n    secondAxisVector = scratchVector2.fromArray(degeneratePositionLocalFrame[secondAxis]);\n    if (secondAxis !== 'east' && secondAxis !== 'west') {\n      secondAxisVector.scale(sign);\n    }\n\n    thirdAxisVector = scratchVector3.fromArray(degeneratePositionLocalFrame[thirdAxis]);\n    if (thirdAxis !== 'east' && thirdAxis !== 'west') {\n      thirdAxisVector.scale(sign);\n    }\n  } else {\n    // Calculate all axis\n    const {up, east, north} = scratchAxisVectors;\n\n    east.set(-origin.y, origin.x, 0.0).normalize();\n    ellipsoid.geodeticSurfaceNormal(origin, up);\n    north.copy(up).cross(east);\n\n    const {down, west, south} = scratchAxisVectors;\n\n    down.copy(up).scale(-1);\n    west.copy(east).scale(-1);\n    south.copy(north).scale(-1);\n\n    // Pick three axis based on desired orientation\n    firstAxisVector = scratchAxisVectors[firstAxis];\n    secondAxisVector = scratchAxisVectors[secondAxis];\n    thirdAxisVector = scratchAxisVectors[thirdAxis];\n  }\n\n  // TODO - assuming the result is column-major\n  result[0] = firstAxisVector.x;\n  result[1] = firstAxisVector.y;\n  result[2] = firstAxisVector.z;\n  result[3] = 0.0;\n  result[4] = secondAxisVector.x;\n  result[5] = secondAxisVector.y;\n  result[6] = secondAxisVector.z;\n  result[7] = 0.0;\n  result[8] = thirdAxisVector.x;\n  result[9] = thirdAxisVector.y;\n  result[10] = thirdAxisVector.z;\n  result[11] = 0.0;\n  result[12] = origin.x;\n  result[13] = origin.y;\n  result[14] = origin.z;\n  result[15] = 1.0;\n  return result;\n}\n"],"file":"ellipsoid-transform.js"}
+{"version":3,"sources":["../../../../src/ellipsoid/helpers/ellipsoid-transform.ts"],"names":["EPSILON14","scratchOrigin","Vector3","VECTOR_PRODUCT_LOCAL_FRAME","up","south","north","west","east","down","degeneratePositionLocalFrame","scratchAxisVectors","scratchVector1","scratchVector2","scratchVector3","localFrameToFixedFrame","ellipsoid","firstAxis","secondAxis","thirdAxis","cartesianOrigin","result","thirdAxisInferred","firstAxisVector","secondAxisVector","thirdAxisVector","origin","copy","atPole","x","y","sign","Math","z","fromArray","scale","set","normalize","geodeticSurfaceNormal","cross"],"mappings":";;;;;;;AAAA;;AAIA,IAAMA,SAAS,GAAG,KAAlB;AAEA,IAAMC,aAAa,GAAG,IAAIC,aAAJ,EAAtB;AAGA,IAAMC,0BAA0B,GAAG;AACjCC,EAAAA,EAAE,EAAE;AACFC,IAAAA,KAAK,EAAE,MADL;AAEFC,IAAAA,KAAK,EAAE,MAFL;AAGFC,IAAAA,IAAI,EAAE,OAHJ;AAIFC,IAAAA,IAAI,EAAE;AAJJ,GAD6B;AAOjCC,EAAAA,IAAI,EAAE;AACJJ,IAAAA,KAAK,EAAE,MADH;AAEJC,IAAAA,KAAK,EAAE,MAFH;AAGJC,IAAAA,IAAI,EAAE,OAHF;AAIJC,IAAAA,IAAI,EAAE;AAJF,GAP2B;AAajCH,EAAAA,KAAK,EAAE;AACLD,IAAAA,EAAE,EAAE,MADC;AAELK,IAAAA,IAAI,EAAE,MAFD;AAGLF,IAAAA,IAAI,EAAE,MAHD;AAILC,IAAAA,IAAI,EAAE;AAJD,GAb0B;AAmBjCF,EAAAA,KAAK,EAAE;AACLF,IAAAA,EAAE,EAAE,MADC;AAELK,IAAAA,IAAI,EAAE,MAFD;AAGLF,IAAAA,IAAI,EAAE,IAHD;AAILC,IAAAA,IAAI,EAAE;AAJD,GAnB0B;AAyBjCD,EAAAA,IAAI,EAAE;AACJH,IAAAA,EAAE,EAAE,OADA;AAEJK,IAAAA,IAAI,EAAE,OAFF;AAGJH,IAAAA,KAAK,EAAE,MAHH;AAIJD,IAAAA,KAAK,EAAE;AAJH,GAzB2B;AA+BjCG,EAAAA,IAAI,EAAE;AACJJ,IAAAA,EAAE,EAAE,OADA;AAEJK,IAAAA,IAAI,EAAE,OAFF;AAGJH,IAAAA,KAAK,EAAE,IAHH;AAIJD,IAAAA,KAAK,EAAE;AAJH;AA/B2B,CAAnC;AAuCA,IAAMK,4BAA4B,GAAG;AACnCJ,EAAAA,KAAK,EAAE,CAAC,CAAC,CAAF,EAAK,CAAL,EAAQ,CAAR,CAD4B;AAEnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAF6B;AAGnCJ,EAAAA,EAAE,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAH+B;AAInCC,EAAAA,KAAK,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAJ4B;AAKnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAC,CAAL,EAAQ,CAAR,CAL6B;AAMnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAC,CAAR;AAN6B,CAArC;AASA,IAAME,kBAAkB,GAAG;AACzBH,EAAAA,IAAI,EAAE,IAAIN,aAAJ,EADmB;AAEzBI,EAAAA,KAAK,EAAE,IAAIJ,aAAJ,EAFkB;AAGzBE,EAAAA,EAAE,EAAE,IAAIF,aAAJ,EAHqB;AAIzBK,EAAAA,IAAI,EAAE,IAAIL,aAAJ,EAJmB;AAKzBG,EAAAA,KAAK,EAAE,IAAIH,aAAJ,EALkB;AAMzBO,EAAAA,IAAI,EAAE,IAAIP,aAAJ;AANmB,CAA3B;AASA,IAAMU,cAAc,GAAG,IAAIV,aAAJ,EAAvB;AACA,IAAMW,cAAc,GAAG,IAAIX,aAAJ,EAAvB;AACA,IAAMY,cAAc,GAAG,IAAIZ,aAAJ,EAAvB;;AAOe,SAASa,sBAAT,CACbC,SADa,EAEbC,SAFa,EAGbC,UAHa,EAIbC,SAJa,EAKbC,eALa,EAMbC,MANa,EAOH;AACV,MAAMC,iBAAiB,GACrBnB,0BAA0B,CAACc,SAAD,CAA1B,IACCd,0BAA0B,CAACc,SAAD,CAA1B,CAAsCC,UAAtC,CAFH;AAIA,oBAAOI,iBAAiB,KAAK,CAACH,SAAD,IAAcA,SAAS,KAAKG,iBAAjC,CAAxB;AAEA,MAAIC,eAAJ;AACA,MAAIC,gBAAJ;AACA,MAAIC,eAAJ;AAEA,MAAMC,MAAM,GAAGzB,aAAa,CAAC0B,IAAd,CAAmBP,eAAnB,CAAf;AAGA,MAAMQ,MAAM,GAAG,kBAAcF,MAAM,CAACG,CAArB,EAAwB,GAAxB,EAA6B7B,SAA7B,KAA2C,kBAAc0B,MAAM,CAACI,CAArB,EAAwB,GAAxB,EAA6B9B,SAA7B,CAA1D;;AAEA,MAAI4B,MAAJ,EAAY;AAEV,QAAMG,IAAI,GAAGC,IAAI,CAACD,IAAL,CAAUL,MAAM,CAACO,CAAjB,CAAb;AAEAV,IAAAA,eAAe,GAAGX,cAAc,CAACsB,SAAf,CAAyBxB,4BAA4B,CAACO,SAAD,CAArD,CAAlB;;AACA,QAAIA,SAAS,KAAK,MAAd,IAAwBA,SAAS,KAAK,MAA1C,EAAkD;AAChDM,MAAAA,eAAe,CAACY,KAAhB,CAAsBJ,IAAtB;AACD;;AAEDP,IAAAA,gBAAgB,GAAGX,cAAc,CAACqB,SAAf,CAAyBxB,4BAA4B,CAACQ,UAAD,CAArD,CAAnB;;AACA,QAAIA,UAAU,KAAK,MAAf,IAAyBA,UAAU,KAAK,MAA5C,EAAoD;AAClDM,MAAAA,gBAAgB,CAACW,KAAjB,CAAuBJ,IAAvB;AACD;;AAEDN,IAAAA,eAAe,GAAGX,cAAc,CAACoB,SAAf,CAAyBxB,4BAA4B,CAACS,SAAD,CAArD,CAAlB;;AACA,QAAIA,SAAS,KAAK,MAAd,IAAwBA,SAAS,KAAK,MAA1C,EAAkD;AAChDM,MAAAA,eAAe,CAACU,KAAhB,CAAsBJ,IAAtB;AACD;AACF,GAlBD,MAkBO;AAEL,QAAO3B,EAAP,GAA0BO,kBAA1B,CAAOP,EAAP;AAAA,QAAWI,IAAX,GAA0BG,kBAA1B,CAAWH,IAAX;AAAA,QAAiBF,KAAjB,GAA0BK,kBAA1B,CAAiBL,KAAjB;AAEAE,IAAAA,IAAI,CAAC4B,GAAL,CAAS,CAACV,MAAM,CAACI,CAAjB,EAAoBJ,MAAM,CAACG,CAA3B,EAA8B,GAA9B,EAAmCQ,SAAnC;AACArB,IAAAA,SAAS,CAACsB,qBAAV,CAAgCZ,MAAhC,EAAwCtB,EAAxC;AACAE,IAAAA,KAAK,CAACqB,IAAN,CAAWvB,EAAX,EAAemC,KAAf,CAAqB/B,IAArB;AAEA,QAAOC,IAAP,GAA4BE,kBAA5B,CAAOF,IAAP;AAAA,QAAaF,IAAb,GAA4BI,kBAA5B,CAAaJ,IAAb;AAAA,QAAmBF,KAAnB,GAA4BM,kBAA5B,CAAmBN,KAAnB;AAEAI,IAAAA,IAAI,CAACkB,IAAL,CAAUvB,EAAV,EAAc+B,KAAd,CAAoB,CAAC,CAArB;AACA5B,IAAAA,IAAI,CAACoB,IAAL,CAAUnB,IAAV,EAAgB2B,KAAhB,CAAsB,CAAC,CAAvB;AACA9B,IAAAA,KAAK,CAACsB,IAAN,CAAWrB,KAAX,EAAkB6B,KAAlB,CAAwB,CAAC,CAAzB;AAGAZ,IAAAA,eAAe,GAAGZ,kBAAkB,CAACM,SAAD,CAApC;AACAO,IAAAA,gBAAgB,GAAGb,kBAAkB,CAACO,UAAD,CAArC;AACAO,IAAAA,eAAe,GAAGd,kBAAkB,CAACQ,SAAD,CAApC;AACD;;AAGDE,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACM,CAA5B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACO,CAA5B;AACAT,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACU,CAA5B;AACAZ,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAY,GAAZ;AACAA,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACK,CAA7B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACM,CAA7B;AACAT,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACS,CAA7B;AACAZ,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAY,GAAZ;AACAA,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYI,eAAe,CAACI,CAA5B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYI,eAAe,CAACK,CAA5B;AACAT,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaI,eAAe,CAACQ,CAA7B;AACAZ,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAa,GAAb;AACAA,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACG,CAApB;AACAR,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACI,CAApB;AACAT,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACO,CAApB;AACAZ,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAa,GAAb;AACA,SAAOA,MAAP;AACD","sourcesContent":["import {Vector3, assert, equals as equalsEpsilon} from '@math.gl/core';\n\nimport type Ellipsoid from '../ellipsoid';\n\nconst EPSILON14 = 1e-14;\n\nconst scratchOrigin = new Vector3();\n\n// Caclulate third axis from given two axii\nconst VECTOR_PRODUCT_LOCAL_FRAME = {\n  up: {\n    south: 'east',\n    north: 'west',\n    west: 'south',\n    east: 'north'\n  },\n  down: {\n    south: 'west',\n    north: 'east',\n    west: 'north',\n    east: 'south'\n  },\n  south: {\n    up: 'west',\n    down: 'east',\n    west: 'down',\n    east: 'up'\n  },\n  north: {\n    up: 'east',\n    down: 'west',\n    west: 'up',\n    east: 'down'\n  },\n  west: {\n    up: 'north',\n    down: 'south',\n    north: 'down',\n    south: 'up'\n  },\n  east: {\n    up: 'south',\n    down: 'north',\n    north: 'up',\n    south: 'down'\n  }\n} as const;\n\nconst degeneratePositionLocalFrame = {\n  north: [-1, 0, 0],\n  east: [0, 1, 0],\n  up: [0, 0, 1],\n  south: [1, 0, 0],\n  west: [0, -1, 0],\n  down: [0, 0, -1]\n} as const;\n\nconst scratchAxisVectors = {\n  east: new Vector3(),\n  north: new Vector3(),\n  up: new Vector3(),\n  west: new Vector3(),\n  south: new Vector3(),\n  down: new Vector3()\n};\n\nconst scratchVector1 = new Vector3();\nconst scratchVector2 = new Vector3();\nconst scratchVector3 = new Vector3();\n\ntype Axis = 'up' | 'down' | 'north' | 'east' | 'south' | 'west';\n\n// Computes a 4x4 transformation matrix from a reference frame\n// centered at the provided origin to the provided ellipsoid's fixed reference frame.\n// eslint-disable-next-line max-statements, max-params, complexity\nexport default function localFrameToFixedFrame(\n  ellipsoid: Ellipsoid,\n  firstAxis: Axis,\n  secondAxis: Axis,\n  thirdAxis: Axis,\n  cartesianOrigin: number[],\n  result: number[]\n): number[] {\n  const thirdAxisInferred =\n    VECTOR_PRODUCT_LOCAL_FRAME[firstAxis] &&\n    (VECTOR_PRODUCT_LOCAL_FRAME[firstAxis][secondAxis] as Axis);\n  // firstAxis and secondAxis must be east, north, up, west, south or down.');\n  assert(thirdAxisInferred && (!thirdAxis || thirdAxis === thirdAxisInferred));\n\n  let firstAxisVector: Vector3;\n  let secondAxisVector: Vector3;\n  let thirdAxisVector: Vector3;\n\n  const origin = scratchOrigin.copy(cartesianOrigin);\n\n  // If x and y are zero, assume origin is at a pole, which is a special case.\n  const atPole = equalsEpsilon(origin.x, 0.0, EPSILON14) && equalsEpsilon(origin.y, 0.0, EPSILON14);\n\n  if (atPole) {\n    // Look up axis value and adjust\n    const sign = Math.sign(origin.z);\n\n    firstAxisVector = scratchVector1.fromArray(degeneratePositionLocalFrame[firstAxis]);\n    if (firstAxis !== 'east' && firstAxis !== 'west') {\n      firstAxisVector.scale(sign);\n    }\n\n    secondAxisVector = scratchVector2.fromArray(degeneratePositionLocalFrame[secondAxis]);\n    if (secondAxis !== 'east' && secondAxis !== 'west') {\n      secondAxisVector.scale(sign);\n    }\n\n    thirdAxisVector = scratchVector3.fromArray(degeneratePositionLocalFrame[thirdAxis]);\n    if (thirdAxis !== 'east' && thirdAxis !== 'west') {\n      thirdAxisVector.scale(sign);\n    }\n  } else {\n    // Calculate all axis\n    const {up, east, north} = scratchAxisVectors;\n\n    east.set(-origin.y, origin.x, 0.0).normalize();\n    ellipsoid.geodeticSurfaceNormal(origin, up);\n    north.copy(up).cross(east);\n\n    const {down, west, south} = scratchAxisVectors;\n\n    down.copy(up).scale(-1);\n    west.copy(east).scale(-1);\n    south.copy(north).scale(-1);\n\n    // Pick three axis based on desired orientation\n    firstAxisVector = scratchAxisVectors[firstAxis];\n    secondAxisVector = scratchAxisVectors[secondAxis];\n    thirdAxisVector = scratchAxisVectors[thirdAxis];\n  }\n\n  // TODO - assuming the result is column-major\n  result[0] = firstAxisVector.x;\n  result[1] = firstAxisVector.y;\n  result[2] = firstAxisVector.z;\n  result[3] = 0.0;\n  result[4] = secondAxisVector.x;\n  result[5] = secondAxisVector.y;\n  result[6] = secondAxisVector.z;\n  result[7] = 0.0;\n  result[8] = thirdAxisVector.x;\n  result[9] = thirdAxisVector.y;\n  result[10] = thirdAxisVector.z;\n  result[11] = 0.0;\n  result[12] = origin.x;\n  result[13] = origin.y;\n  result[14] = origin.z;\n  result[15] = 1.0;\n  return result;\n}\n"],"file":"ellipsoid-transform.js"}

package/dist/es5/ellipsoid/helpers/scale-to-geodetic-surface.js

@@ -1,7 +1,5 @@
 "use strict";
 
-var _typeof = require("@babel/runtime/helpers/typeof");
-
 Object.defineProperty(exports, "__esModule", {
   value: true
 });
@@ -9,25 +7,19 @@ exports.default = scaleToGeodeticSurface;
 
 var _core = require("@math.gl/core");
 
-var vec3 = _interopRequireWildcard(require("gl-matrix/vec3"));
-
-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); }
-
-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; }
-
 var scratchVector = new _core.Vector3();
 var scaleToGeodeticSurfaceIntersection = new _core.Vector3();
 var scaleToGeodeticSurfaceGradient = new _core.Vector3();
 
 function scaleToGeodeticSurface(cartesian, ellipsoid) {
-  var result = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : new _core.Vector3();
+  var result = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : [];
   var oneOverRadii = ellipsoid.oneOverRadii,
       oneOverRadiiSquared = ellipsoid.oneOverRadiiSquared,
       centerToleranceSquared = ellipsoid.centerToleranceSquared;
   scratchVector.from(cartesian);
-  var positionX = cartesian.x;
-  var positionY = cartesian.y;
-  var positionZ = cartesian.z;
+  var positionX = scratchVector.x;
+  var positionY = scratchVector.y;
+  var positionZ = scratchVector.z;
   var oneOverRadiiX = oneOverRadii.x;
   var oneOverRadiiY = oneOverRadii.y;
   var oneOverRadiiZ = oneOverRadii.z;
@@ -53,7 +45,7 @@ function scaleToGeodeticSurface(cartesian, ellipsoid) {
   var oneOverRadiiSquaredZ = oneOverRadiiSquared.z;
   var gradient = scaleToGeodeticSurfaceGradient;
   gradient.set(intersection.x * oneOverRadiiSquaredX * 2.0, intersection.y * oneOverRadiiSquaredY * 2.0, intersection.z * oneOverRadiiSquaredZ * 2.0);
-  var lambda = (1.0 - ratio) * cartesian.len() / (0.5 * gradient.len());
+  var lambda = (1.0 - ratio) * scratchVector.len() / (0.5 * gradient.len());
   var correction = 0.0;
   var xMultiplier;
   var yMultiplier;

package/dist/es5/ellipsoid/helpers/scale-to-geodetic-surface.js.map

@@ -1 +1 @@
-{"version":3,"sources":["../../../../src/ellipsoid/helpers/scale-to-geodetic-surface.js"],"names":["scratchVector","Vector3","scaleToGeodeticSurfaceIntersection","scaleToGeodeticSurfaceGradient","scaleToGeodeticSurface","cartesian","ellipsoid","result","oneOverRadii","oneOverRadiiSquared","centerToleranceSquared","from","positionX","x","positionY","y","positionZ","z","oneOverRadiiX","oneOverRadiiY","oneOverRadiiZ","x2","y2","z2","squaredNorm","ratio","Math","sqrt","Number","isFinite","undefined","intersection","copy","scale","to","oneOverRadiiSquaredX","oneOverRadiiSquaredY","oneOverRadiiSquaredZ","gradient","set","lambda","len","correction","xMultiplier","yMultiplier","zMultiplier","func","xMultiplier2","yMultiplier2","zMultiplier2","xMultiplier3","yMultiplier3","zMultiplier3","denominator","derivative","abs","_MathUtils","EPSILON12"],"mappings":";;;;;;;;;AACA;;AACA;;;;;;AAEA,IAAMA,aAAa,GAAG,IAAIC,aAAJ,EAAtB;AACA,IAAMC,kCAAkC,GAAG,IAAID,aAAJ,EAA3C;AACA,IAAME,8BAA8B,GAAG,IAAIF,aAAJ,EAAvC;;AAKe,SAASG,sBAAT,CAAgCC,SAAhC,EAA2CC,SAA3C,EAA8E;AAAA,MAAxBC,MAAwB,uEAAf,IAAIN,aAAJ,EAAe;AAC3F,MAAOO,YAAP,GAAoEF,SAApE,CAAOE,YAAP;AAAA,MAAqBC,mBAArB,GAAoEH,SAApE,CAAqBG,mBAArB;AAAA,MAA0CC,sBAA1C,GAAoEJ,SAApE,CAA0CI,sBAA1C;AAEAV,EAAAA,aAAa,CAACW,IAAd,CAAmBN,SAAnB;AAEA,MAAMO,SAAS,GAAGP,SAAS,CAACQ,CAA5B;AACA,MAAMC,SAAS,GAAGT,SAAS,CAACU,CAA5B;AACA,MAAMC,SAAS,GAAGX,SAAS,CAACY,CAA5B;AAEA,MAAMC,aAAa,GAAGV,YAAY,CAACK,CAAnC;AACA,MAAMM,aAAa,GAAGX,YAAY,CAACO,CAAnC;AACA,MAAMK,aAAa,GAAGZ,YAAY,CAACS,CAAnC;AAEA,MAAMI,EAAE,GAAGT,SAAS,GAAGA,SAAZ,GAAwBM,aAAxB,GAAwCA,aAAnD;AACA,MAAMI,EAAE,GAAGR,SAAS,GAAGA,SAAZ,GAAwBK,aAAxB,GAAwCA,aAAnD;AACA,MAAMI,EAAE,GAAGP,SAAS,GAAGA,SAAZ,GAAwBI,aAAxB,GAAwCA,aAAnD;AAGA,MAAMI,WAAW,GAAGH,EAAE,GAAGC,EAAL,GAAUC,EAA9B;AACA,MAAME,KAAK,GAAGC,IAAI,CAACC,IAAL,CAAU,MAAMH,WAAhB,CAAd;;AAGA,MAAI,CAACI,MAAM,CAACC,QAAP,CAAgBJ,KAAhB,CAAL,EAA6B;AAC3B,WAAOK,SAAP;AACD;;AAGD,MAAMC,YAAY,GAAG7B,kCAArB;AACA6B,EAAAA,YAAY,CAACC,IAAb,CAAkB3B,SAAlB,EAA6B4B,KAA7B,CAAmCR,KAAnC;;AAGA,MAAID,WAAW,GAAGd,sBAAlB,EAA0C;AACxC,WAAOqB,YAAY,CAACG,EAAb,CAAgB3B,MAAhB,CAAP;AACD;;AAED,MAAM4B,oBAAoB,GAAG1B,mBAAmB,CAACI,CAAjD;AACA,MAAMuB,oBAAoB,GAAG3B,mBAAmB,CAACM,CAAjD;AACA,MAAMsB,oBAAoB,GAAG5B,mBAAmB,CAACQ,CAAjD;AAIA,MAAMqB,QAAQ,GAAGnC,8BAAjB;AACAmC,EAAAA,QAAQ,CAACC,GAAT,CACER,YAAY,CAAClB,CAAb,GAAiBsB,oBAAjB,GAAwC,GAD1C,EAEEJ,YAAY,CAAChB,CAAb,GAAiBqB,oBAAjB,GAAwC,GAF1C,EAGEL,YAAY,CAACd,CAAb,GAAiBoB,oBAAjB,GAAwC,GAH1C;AAOA,MAAIG,MAAM,GAAI,CAAC,MAAMf,KAAP,IAAgBpB,SAAS,CAACoC,GAAV,EAAjB,IAAqC,MAAMH,QAAQ,CAACG,GAAT,EAA3C,CAAb;AACA,MAAIC,UAAU,GAAG,GAAjB;AAEA,MAAIC,WAAJ;AACA,MAAIC,WAAJ;AACA,MAAIC,WAAJ;AACA,MAAIC,IAAJ;;AAEA,KAAG;AACDN,IAAAA,MAAM,IAAIE,UAAV;AAEAC,IAAAA,WAAW,GAAG,OAAO,MAAMH,MAAM,GAAGL,oBAAtB,CAAd;AACAS,IAAAA,WAAW,GAAG,OAAO,MAAMJ,MAAM,GAAGJ,oBAAtB,CAAd;AACAS,IAAAA,WAAW,GAAG,OAAO,MAAML,MAAM,GAAGH,oBAAtB,CAAd;AAEA,QAAMU,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AACA,QAAMK,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AACA,QAAMK,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AAEA,QAAMK,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AACA,QAAMQ,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AACA,QAAMQ,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AAEAC,IAAAA,IAAI,GAAGzB,EAAE,GAAG0B,YAAL,GAAoBzB,EAAE,GAAG0B,YAAzB,GAAwCzB,EAAE,GAAG0B,YAA7C,GAA4D,GAAnE;AAIA,QAAMI,WAAW,GACfhC,EAAE,GAAG6B,YAAL,GAAoBf,oBAApB,GACAb,EAAE,GAAG6B,YAAL,GAAoBf,oBADpB,GAEAb,EAAE,GAAG6B,YAAL,GAAoBf,oBAHtB;AAKA,QAAMiB,UAAU,GAAG,CAAC,GAAD,GAAOD,WAA1B;AAEAX,IAAAA,UAAU,GAAGI,IAAI,GAAGQ,UAApB;AACD,GA3BD,QA2BS5B,IAAI,CAAC6B,GAAL,CAAST,IAAT,IAAiBU,iBAAWC,SA3BrC;;AA6BA,SAAOzD,aAAa,CAACiC,KAAd,CAAoB,CAACU,WAAD,EAAcC,WAAd,EAA2BC,WAA3B,CAApB,EAA6DX,EAA7D,CAAgE3B,MAAhE,CAAP;AACD","sourcesContent":["/* eslint-disable */\nimport {Vector3, assert, _MathUtils} from '@math.gl/core';\nimport * as vec3 from 'gl-matrix/vec3';\n\nconst scratchVector = new Vector3();\nconst scaleToGeodeticSurfaceIntersection = new Vector3();\nconst scaleToGeodeticSurfaceGradient = new Vector3();\n\n// Scales the provided Cartesian position along the geodetic surface normal\n// so that it is on the surface of this ellipsoid.  If the position is\n// at the center of the ellipsoid, this function returns undefined.\nexport default function scaleToGeodeticSurface(cartesian, ellipsoid, result = new Vector3()) {\n  const {oneOverRadii, oneOverRadiiSquared, centerToleranceSquared} = ellipsoid;\n\n  scratchVector.from(cartesian);\n\n  const positionX = cartesian.x;\n  const positionY = cartesian.y;\n  const positionZ = cartesian.z;\n\n  const oneOverRadiiX = oneOverRadii.x;\n  const oneOverRadiiY = oneOverRadii.y;\n  const oneOverRadiiZ = oneOverRadii.z;\n\n  const x2 = positionX * positionX * oneOverRadiiX * oneOverRadiiX;\n  const y2 = positionY * positionY * oneOverRadiiY * oneOverRadiiY;\n  const z2 = positionZ * positionZ * oneOverRadiiZ * oneOverRadiiZ;\n\n  // Compute the squared ellipsoid norm.\n  const squaredNorm = x2 + y2 + z2;\n  const ratio = Math.sqrt(1.0 / squaredNorm);\n\n  // When very close to center or at center\n  if (!Number.isFinite(ratio)) {\n    return undefined;\n  }\n\n  // As an initial approximation, assume that the radial intersection is the projection point.\n  const intersection = scaleToGeodeticSurfaceIntersection;\n  intersection.copy(cartesian).scale(ratio);\n\n  // If the position is near the center, the iteration will not converge.\n  if (squaredNorm < centerToleranceSquared) {\n    return intersection.to(result);\n  }\n\n  const oneOverRadiiSquaredX = oneOverRadiiSquared.x;\n  const oneOverRadiiSquaredY = oneOverRadiiSquared.y;\n  const oneOverRadiiSquaredZ = oneOverRadiiSquared.z;\n\n  // Use the gradient at the intersection point in place of the true unit normal.\n  // The difference in magnitude will be absorbed in the multiplier.\n  const gradient = scaleToGeodeticSurfaceGradient;\n  gradient.set(\n    intersection.x * oneOverRadiiSquaredX * 2.0,\n    intersection.y * oneOverRadiiSquaredY * 2.0,\n    intersection.z * oneOverRadiiSquaredZ * 2.0\n  );\n\n  // Compute the initial guess at the normal vector multiplier, lambda.\n  let lambda = ((1.0 - ratio) * cartesian.len()) / (0.5 * gradient.len());\n  let correction = 0.0;\n\n  let xMultiplier;\n  let yMultiplier;\n  let zMultiplier;\n  let func;\n\n  do {\n    lambda -= correction;\n\n    xMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredX);\n    yMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredY);\n    zMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredZ);\n\n    const xMultiplier2 = xMultiplier * xMultiplier;\n    const yMultiplier2 = yMultiplier * yMultiplier;\n    const zMultiplier2 = zMultiplier * zMultiplier;\n\n    const xMultiplier3 = xMultiplier2 * xMultiplier;\n    const yMultiplier3 = yMultiplier2 * yMultiplier;\n    const zMultiplier3 = zMultiplier2 * zMultiplier;\n\n    func = x2 * xMultiplier2 + y2 * yMultiplier2 + z2 * zMultiplier2 - 1.0;\n\n    // \"denominator\" here refers to the use of this expression in the velocity and acceleration\n    // computations in the sections to follow.\n    const denominator =\n      x2 * xMultiplier3 * oneOverRadiiSquaredX +\n      y2 * yMultiplier3 * oneOverRadiiSquaredY +\n      z2 * zMultiplier3 * oneOverRadiiSquaredZ;\n\n    const derivative = -2.0 * denominator;\n\n    correction = func / derivative;\n  } while (Math.abs(func) > _MathUtils.EPSILON12);\n\n  return scratchVector.scale([xMultiplier, yMultiplier, zMultiplier]).to(result);\n}\n"],"file":"scale-to-geodetic-surface.js"}
+{"version":3,"sources":["../../../../src/ellipsoid/helpers/scale-to-geodetic-surface.ts"],"names":["scratchVector","Vector3","scaleToGeodeticSurfaceIntersection","scaleToGeodeticSurfaceGradient","scaleToGeodeticSurface","cartesian","ellipsoid","result","oneOverRadii","oneOverRadiiSquared","centerToleranceSquared","from","positionX","x","positionY","y","positionZ","z","oneOverRadiiX","oneOverRadiiY","oneOverRadiiZ","x2","y2","z2","squaredNorm","ratio","Math","sqrt","Number","isFinite","undefined","intersection","copy","scale","to","oneOverRadiiSquaredX","oneOverRadiiSquaredY","oneOverRadiiSquaredZ","gradient","set","lambda","len","correction","xMultiplier","yMultiplier","zMultiplier","func","xMultiplier2","yMultiplier2","zMultiplier2","xMultiplier3","yMultiplier3","zMultiplier3","denominator","derivative","abs","_MathUtils","EPSILON12"],"mappings":";;;;;;;AACA;;AAGA,IAAMA,aAAa,GAAG,IAAIC,aAAJ,EAAtB;AACA,IAAMC,kCAAkC,GAAG,IAAID,aAAJ,EAA3C;AACA,IAAME,8BAA8B,GAAG,IAAIF,aAAJ,EAAvC;;AAKe,SAASG,sBAAT,CACbC,SADa,EAEbC,SAFa,EAIH;AAAA,MADVC,MACU,uEADS,EACT;AACV,MAAOC,YAAP,GAAoEF,SAApE,CAAOE,YAAP;AAAA,MAAqBC,mBAArB,GAAoEH,SAApE,CAAqBG,mBAArB;AAAA,MAA0CC,sBAA1C,GAAoEJ,SAApE,CAA0CI,sBAA1C;AAEAV,EAAAA,aAAa,CAACW,IAAd,CAAmBN,SAAnB;AAEA,MAAMO,SAAS,GAAGZ,aAAa,CAACa,CAAhC;AACA,MAAMC,SAAS,GAAGd,aAAa,CAACe,CAAhC;AACA,MAAMC,SAAS,GAAGhB,aAAa,CAACiB,CAAhC;AAEA,MAAMC,aAAa,GAAGV,YAAY,CAACK,CAAnC;AACA,MAAMM,aAAa,GAAGX,YAAY,CAACO,CAAnC;AACA,MAAMK,aAAa,GAAGZ,YAAY,CAACS,CAAnC;AAEA,MAAMI,EAAE,GAAGT,SAAS,GAAGA,SAAZ,GAAwBM,aAAxB,GAAwCA,aAAnD;AACA,MAAMI,EAAE,GAAGR,SAAS,GAAGA,SAAZ,GAAwBK,aAAxB,GAAwCA,aAAnD;AACA,MAAMI,EAAE,GAAGP,SAAS,GAAGA,SAAZ,GAAwBI,aAAxB,GAAwCA,aAAnD;AAGA,MAAMI,WAAW,GAAGH,EAAE,GAAGC,EAAL,GAAUC,EAA9B;AACA,MAAME,KAAK,GAAGC,IAAI,CAACC,IAAL,CAAU,MAAMH,WAAhB,CAAd;;AAGA,MAAI,CAACI,MAAM,CAACC,QAAP,CAAgBJ,KAAhB,CAAL,EAA6B;AAC3B,WAAOK,SAAP;AACD;;AAGD,MAAMC,YAAY,GAAG7B,kCAArB;AACA6B,EAAAA,YAAY,CAACC,IAAb,CAAkB3B,SAAlB,EAA6B4B,KAA7B,CAAmCR,KAAnC;;AAGA,MAAID,WAAW,GAAGd,sBAAlB,EAA0C;AACxC,WAAOqB,YAAY,CAACG,EAAb,CAAgB3B,MAAhB,CAAP;AACD;;AAED,MAAM4B,oBAAoB,GAAG1B,mBAAmB,CAACI,CAAjD;AACA,MAAMuB,oBAAoB,GAAG3B,mBAAmB,CAACM,CAAjD;AACA,MAAMsB,oBAAoB,GAAG5B,mBAAmB,CAACQ,CAAjD;AAIA,MAAMqB,QAAQ,GAAGnC,8BAAjB;AACAmC,EAAAA,QAAQ,CAACC,GAAT,CACER,YAAY,CAAClB,CAAb,GAAiBsB,oBAAjB,GAAwC,GAD1C,EAEEJ,YAAY,CAAChB,CAAb,GAAiBqB,oBAAjB,GAAwC,GAF1C,EAGEL,YAAY,CAACd,CAAb,GAAiBoB,oBAAjB,GAAwC,GAH1C;AAOA,MAAIG,MAAM,GAAI,CAAC,MAAMf,KAAP,IAAgBzB,aAAa,CAACyC,GAAd,EAAjB,IAAyC,MAAMH,QAAQ,CAACG,GAAT,EAA/C,CAAb;AACA,MAAIC,UAAU,GAAG,GAAjB;AAEA,MAAIC,WAAJ;AACA,MAAIC,WAAJ;AACA,MAAIC,WAAJ;AACA,MAAIC,IAAJ;;AAEA,KAAG;AACDN,IAAAA,MAAM,IAAIE,UAAV;AAEAC,IAAAA,WAAW,GAAG,OAAO,MAAMH,MAAM,GAAGL,oBAAtB,CAAd;AACAS,IAAAA,WAAW,GAAG,OAAO,MAAMJ,MAAM,GAAGJ,oBAAtB,CAAd;AACAS,IAAAA,WAAW,GAAG,OAAO,MAAML,MAAM,GAAGH,oBAAtB,CAAd;AAEA,QAAMU,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AACA,QAAMK,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AACA,QAAMK,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AAEA,QAAMK,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AACA,QAAMQ,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AACA,QAAMQ,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AAEAC,IAAAA,IAAI,GAAGzB,EAAE,GAAG0B,YAAL,GAAoBzB,EAAE,GAAG0B,YAAzB,GAAwCzB,EAAE,GAAG0B,YAA7C,GAA4D,GAAnE;AAIA,QAAMI,WAAW,GACfhC,EAAE,GAAG6B,YAAL,GAAoBf,oBAApB,GACAb,EAAE,GAAG6B,YAAL,GAAoBf,oBADpB,GAEAb,EAAE,GAAG6B,YAAL,GAAoBf,oBAHtB;AAKA,QAAMiB,UAAU,GAAG,CAAC,GAAD,GAAOD,WAA1B;AAEAX,IAAAA,UAAU,GAAGI,IAAI,GAAGQ,UAApB;AACD,GA3BD,QA2BS5B,IAAI,CAAC6B,GAAL,CAAST,IAAT,IAAiBU,iBAAWC,SA3BrC;;AA6BA,SAAOzD,aAAa,CAACiC,KAAd,CAAoB,CAACU,WAAD,EAAcC,WAAd,EAA2BC,WAA3B,CAApB,EAA6DX,EAA7D,CAAgE3B,MAAhE,CAAP;AACD","sourcesContent":["/* eslint-disable */\nimport {Vector3, _MathUtils} from '@math.gl/core';\nimport type Ellipsoid from '../ellipsoid';\n\nconst scratchVector = new Vector3();\nconst scaleToGeodeticSurfaceIntersection = new Vector3();\nconst scaleToGeodeticSurfaceGradient = new Vector3();\n\n// Scales the provided Cartesian position along the geodetic surface normal\n// so that it is on the surface of this ellipsoid.  If the position is\n// at the center of the ellipsoid, this function returns undefined.\nexport default function scaleToGeodeticSurface(\n  cartesian: number[],\n  ellipsoid: Ellipsoid,\n  result: number[] = []\n): number[] {\n  const {oneOverRadii, oneOverRadiiSquared, centerToleranceSquared} = ellipsoid;\n\n  scratchVector.from(cartesian);\n\n  const positionX = scratchVector.x;\n  const positionY = scratchVector.y;\n  const positionZ = scratchVector.z;\n\n  const oneOverRadiiX = oneOverRadii.x;\n  const oneOverRadiiY = oneOverRadii.y;\n  const oneOverRadiiZ = oneOverRadii.z;\n\n  const x2 = positionX * positionX * oneOverRadiiX * oneOverRadiiX;\n  const y2 = positionY * positionY * oneOverRadiiY * oneOverRadiiY;\n  const z2 = positionZ * positionZ * oneOverRadiiZ * oneOverRadiiZ;\n\n  // Compute the squared ellipsoid norm.\n  const squaredNorm = x2 + y2 + z2;\n  const ratio = Math.sqrt(1.0 / squaredNorm);\n\n  // When very close to center or at center\n  if (!Number.isFinite(ratio)) {\n    return undefined;\n  }\n\n  // As an initial approximation, assume that the radial intersection is the projection point.\n  const intersection = scaleToGeodeticSurfaceIntersection;\n  intersection.copy(cartesian).scale(ratio);\n\n  // If the position is near the center, the iteration will not converge.\n  if (squaredNorm < centerToleranceSquared) {\n    return intersection.to(result);\n  }\n\n  const oneOverRadiiSquaredX = oneOverRadiiSquared.x;\n  const oneOverRadiiSquaredY = oneOverRadiiSquared.y;\n  const oneOverRadiiSquaredZ = oneOverRadiiSquared.z;\n\n  // Use the gradient at the intersection point in place of the true unit normal.\n  // The difference in magnitude will be absorbed in the multiplier.\n  const gradient = scaleToGeodeticSurfaceGradient;\n  gradient.set(\n    intersection.x * oneOverRadiiSquaredX * 2.0,\n    intersection.y * oneOverRadiiSquaredY * 2.0,\n    intersection.z * oneOverRadiiSquaredZ * 2.0\n  );\n\n  // Compute the initial guess at the normal vector multiplier, lambda.\n  let lambda = ((1.0 - ratio) * scratchVector.len()) / (0.5 * gradient.len());\n  let correction = 0.0;\n\n  let xMultiplier;\n  let yMultiplier;\n  let zMultiplier;\n  let func;\n\n  do {\n    lambda -= correction;\n\n    xMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredX);\n    yMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredY);\n    zMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredZ);\n\n    const xMultiplier2 = xMultiplier * xMultiplier;\n    const yMultiplier2 = yMultiplier * yMultiplier;\n    const zMultiplier2 = zMultiplier * zMultiplier;\n\n    const xMultiplier3 = xMultiplier2 * xMultiplier;\n    const yMultiplier3 = yMultiplier2 * yMultiplier;\n    const zMultiplier3 = zMultiplier2 * zMultiplier;\n\n    func = x2 * xMultiplier2 + y2 * yMultiplier2 + z2 * zMultiplier2 - 1.0;\n\n    // \"denominator\" here refers to the use of this expression in the velocity and acceleration\n    // computations in the sections to follow.\n    const denominator =\n      x2 * xMultiplier3 * oneOverRadiiSquaredX +\n      y2 * yMultiplier3 * oneOverRadiiSquaredY +\n      z2 * zMultiplier3 * oneOverRadiiSquaredZ;\n\n    const derivative = -2.0 * denominator;\n\n    correction = func / derivative;\n  } while (Math.abs(func) > _MathUtils.EPSILON12);\n\n  return scratchVector.scale([xMultiplier, yMultiplier, zMultiplier]).to(result);\n}\n"],"file":"scale-to-geodetic-surface.js"}

package/dist/es5/index.js.map

@@ -1 +1 @@
-{"version":3,"sources":["../../src/index.js"],"names":[],"mappings":";;;;;;;;;;;;;;;;;;;;AAAA;;AACA","sourcesContent":["export {default as Ellipsoid} from './ellipsoid/ellipsoid';\nexport {isWGS84} from './type-utils';\n"],"file":"index.js"}
+{"version":3,"sources":["../../src/index.ts"],"names":[],"mappings":";;;;;;;;;;;;;;;;;;;;AAAA;;AACA","sourcesContent":["export {default as Ellipsoid} from './ellipsoid/ellipsoid';\nexport {isWGS84} from './type-utils';\n"],"file":"index.js"}

package/dist/es5/type-utils.js

@@ -15,68 +15,69 @@ var _core = require("@math.gl/core");
 
 var _constants = require("./constants");
 
-var noop = function noop(x) {
+function identity(x) {
   return x;
-};
+}
 
 var scratchVector = new _core.Vector3();
 
-function fromCartographic(cartographic, result) {
-  var map = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : noop;
+function fromCartographic(cartographic) {
+  var result = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : [];
+  var map = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : identity;
 
-  if ((0, _core.isArray)(cartographic)) {
-    result[0] = map(cartographic[0]);
-    result[1] = map(cartographic[1]);
-    result[2] = cartographic[2];
-  } else if ('longitude' in cartographic) {
+  if ('longitude' in cartographic) {
     result[0] = map(cartographic.longitude);
     result[1] = map(cartographic.latitude);
     result[2] = cartographic.height;
-  } else {
+  } else if ('x' in cartographic) {
     result[0] = map(cartographic.x);
     result[1] = map(cartographic.y);
     result[2] = cartographic.z;
+  } else {
+    result[0] = map(cartographic[0]);
+    result[1] = map(cartographic[1]);
+    result[2] = cartographic[2];
   }
 
   return result;
 }
 
 function fromCartographicToRadians(cartographic) {
-  var vector = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : scratchVector;
-  return fromCartographic(cartographic, vector, _core.config._cartographicRadians ? noop : _core.toRadians);
+  var vector = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : [];
+  return fromCartographic(cartographic, vector, _core.config._cartographicRadians ? identity : _core.toRadians);
 }
 
 function fromCartographicToDegrees(cartographic) {
-  var vector = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : scratchVector;
-  return fromCartographic(cartographic, vector, _core.config._cartographicRadians ? _core.toDegrees : noop);
+  var vector = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : [];
+  return fromCartographic(cartographic, vector, _core.config._cartographicRadians ? _core.toDegrees : identity);
 }
 
 function toCartographic(vector, cartographic) {
-  var map = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : noop;
+  var map = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : identity;
 
-  if ((0, _core.isArray)(cartographic)) {
-    cartographic[0] = map(vector[0]);
-    cartographic[1] = map(vector[1]);
-    cartographic[2] = vector[2];
-  } else if ('longitude' in cartographic) {
+  if ('longitude' in cartographic) {
     cartographic.longitude = map(vector[0]);
     cartographic.latitude = map(vector[1]);
     cartographic.height = vector[2];
-  } else {
+  } else if ('x' in cartographic) {
     cartographic.x = map(vector[0]);
     cartographic.y = map(vector[1]);
     cartographic.z = vector[2];
+  } else {
+    cartographic[0] = map(vector[0]);
+    cartographic[1] = map(vector[1]);
+    cartographic[2] = vector[2];
   }
 
   return cartographic;
 }
 
 function toCartographicFromRadians(vector, cartographic) {
-  return toCartographic(vector, cartographic, _core.config._cartographicRadians ? noop : _core.toDegrees);
+  return toCartographic(vector, cartographic, _core.config._cartographicRadians ? identity : _core.toDegrees);
 }
 
 function toCartographicFromDegrees(vector, cartographic) {
-  return toCartographic(vector, cartographic, _core.config._cartographicRadians ? _core.toRadians : noop);
+  return toCartographic(vector, cartographic, _core.config._cartographicRadians ? _core.toRadians : identity);
 }
 
 function isWGS84(vector) {

package/dist/es5/type-utils.js.map

@@ -1 +1 @@
-{"version":3,"sources":["../../src/type-utils.js"],"names":["noop","x","scratchVector","Vector3","fromCartographic","cartographic","result","map","longitude","latitude","height","y","z","fromCartographicToRadians","vector","config","_cartographicRadians","toRadians","fromCartographicToDegrees","toDegrees","toCartographic","toCartographicFromRadians","toCartographicFromDegrees","isWGS84","from","oneOverRadiiSquared","WGS84_CONSTANTS","centerToleranceSquared","x2","y2","z2","Math","abs"],"mappings":";;;;;;;;;;;;;AAGA;;AACA;;AAEA,IAAMA,IAAI,GAAG,SAAPA,IAAO,CAACC,CAAD;AAAA,SAAOA,CAAP;AAAA,CAAb;;AAEA,IAAMC,aAAa,GAAG,IAAIC,aAAJ,EAAtB;;AAEO,SAASC,gBAAT,CAA0BC,YAA1B,EAAwCC,MAAxC,EAA4D;AAAA,MAAZC,GAAY,uEAANP,IAAM;;AACjE,MAAI,mBAAQK,YAAR,CAAJ,EAA2B;AACzBC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAAC,CAAD,CAAb,CAAf;AACAC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAAC,CAAD,CAAb,CAAf;AACAC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYD,YAAY,CAAC,CAAD,CAAxB;AACD,GAJD,MAIO,IAAI,eAAeA,YAAnB,EAAiC;AACtCC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACG,SAAd,CAAf;AACAF,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACI,QAAd,CAAf;AACAH,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYD,YAAY,CAACK,MAAzB;AACD,GAJM,MAIA;AACLJ,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACJ,CAAd,CAAf;AACAK,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACM,CAAd,CAAf;AACAL,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYD,YAAY,CAACO,CAAzB;AACD;;AACD,SAAON,MAAP;AACD;;AAEM,SAASO,yBAAT,CAAmCR,YAAnC,EAAyE;AAAA,MAAxBS,MAAwB,uEAAfZ,aAAe;AAC9E,SAAOE,gBAAgB,CAACC,YAAD,EAAeS,MAAf,EAAuBC,aAAOC,oBAAP,GAA8BhB,IAA9B,GAAqCiB,eAA5D,CAAvB;AACD;;AAEM,SAASC,yBAAT,CAAmCb,YAAnC,EAAyE;AAAA,MAAxBS,MAAwB,uEAAfZ,aAAe;AAC9E,SAAOE,gBAAgB,CAACC,YAAD,EAAeS,MAAf,EAAuBC,aAAOC,oBAAP,GAA8BG,eAA9B,GAA0CnB,IAAjE,CAAvB;AACD;;AAEM,SAASoB,cAAT,CAAwBN,MAAxB,EAAgCT,YAAhC,EAA0D;AAAA,MAAZE,GAAY,uEAANP,IAAM;;AAC/D,MAAI,mBAAQK,YAAR,CAAJ,EAA2B;AACzBA,IAAAA,YAAY,CAAC,CAAD,CAAZ,GAAkBE,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAArB;AACAT,IAAAA,YAAY,CAAC,CAAD,CAAZ,GAAkBE,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAArB;AACAT,IAAAA,YAAY,CAAC,CAAD,CAAZ,GAAkBS,MAAM,CAAC,CAAD,CAAxB;AACD,GAJD,MAIO,IAAI,eAAeT,YAAnB,EAAiC;AACtCA,IAAAA,YAAY,CAACG,SAAb,GAAyBD,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAA5B;AACAT,IAAAA,YAAY,CAACI,QAAb,GAAwBF,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAA3B;AACAT,IAAAA,YAAY,CAACK,MAAb,GAAsBI,MAAM,CAAC,CAAD,CAA5B;AACD,GAJM,MAIA;AACLT,IAAAA,YAAY,CAACJ,CAAb,GAAiBM,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAApB;AACAT,IAAAA,YAAY,CAACM,CAAb,GAAiBJ,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAApB;AACAT,IAAAA,YAAY,CAACO,CAAb,GAAiBE,MAAM,CAAC,CAAD,CAAvB;AACD;;AACD,SAAOT,YAAP;AACD;;AAEM,SAASgB,yBAAT,CAAmCP,MAAnC,EAA2CT,YAA3C,EAAyD;AAC9D,SAAOe,cAAc,CAACN,MAAD,EAAST,YAAT,EAAuBU,aAAOC,oBAAP,GAA8BhB,IAA9B,GAAqCmB,eAA5D,CAArB;AACD;;AAEM,SAASG,yBAAT,CAAmCR,MAAnC,EAA2CT,YAA3C,EAAyD;AAC9D,SAAOe,cAAc,CAACN,MAAD,EAAST,YAAT,EAAuBU,aAAOC,oBAAP,GAA8BC,eAA9B,GAA0CjB,IAAjE,CAArB;AACD;;AAEM,SAASuB,OAAT,CAAiBT,MAAjB,EAAyB;AAC9B,MAAI,CAACA,MAAL,EAAa;AACX,WAAO,KAAP;AACD;;AACDZ,EAAAA,aAAa,CAACsB,IAAd,CAAmBV,MAAnB;AACA,MAAOW,mBAAP,GAAsDC,0BAAtD,CAAOD,mBAAP;AAAA,MAA4BE,sBAA5B,GAAsDD,0BAAtD,CAA4BC,sBAA5B;AACA,MAAMC,EAAE,GAAGd,MAAM,CAAC,CAAD,CAAN,GAAYA,MAAM,CAAC,CAAD,CAAlB,GAAwBW,mBAAmB,CAAC,CAAD,CAAtD;AACA,MAAMI,EAAE,GAAGf,MAAM,CAAC,CAAD,CAAN,GAAYA,MAAM,CAAC,CAAD,CAAlB,GAAwBW,mBAAmB,CAAC,CAAD,CAAtD;AACA,MAAMK,EAAE,GAAGhB,MAAM,CAAC,CAAD,CAAN,GAAYA,MAAM,CAAC,CAAD,CAAlB,GAAwBW,mBAAmB,CAAC,CAAD,CAAtD;AACA,SAAOM,IAAI,CAACC,GAAL,CAASJ,EAAE,GAAGC,EAAL,GAAUC,EAAV,GAAe,CAAxB,IAA6BH,sBAApC;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, isArray, toRadians, toDegrees, config} from '@math.gl/core';\nimport {WGS84_CONSTANTS} from './constants';\n\nconst noop = (x) => x;\n\nconst scratchVector = new Vector3();\n\nexport function fromCartographic(cartographic, result, map = noop) {\n  if (isArray(cartographic)) {\n    result[0] = map(cartographic[0]);\n    result[1] = map(cartographic[1]);\n    result[2] = cartographic[2];\n  } else if ('longitude' in cartographic) {\n    result[0] = map(cartographic.longitude);\n    result[1] = map(cartographic.latitude);\n    result[2] = cartographic.height;\n  } else {\n    result[0] = map(cartographic.x);\n    result[1] = map(cartographic.y);\n    result[2] = cartographic.z;\n  }\n  return result;\n}\n\nexport function fromCartographicToRadians(cartographic, vector = scratchVector) {\n  return fromCartographic(cartographic, vector, config._cartographicRadians ? noop : toRadians);\n}\n\nexport function fromCartographicToDegrees(cartographic, vector = scratchVector) {\n  return fromCartographic(cartographic, vector, config._cartographicRadians ? toDegrees : noop);\n}\n\nexport function toCartographic(vector, cartographic, map = noop) {\n  if (isArray(cartographic)) {\n    cartographic[0] = map(vector[0]);\n    cartographic[1] = map(vector[1]);\n    cartographic[2] = vector[2];\n  } else if ('longitude' in cartographic) {\n    cartographic.longitude = map(vector[0]);\n    cartographic.latitude = map(vector[1]);\n    cartographic.height = vector[2];\n  } else {\n    cartographic.x = map(vector[0]);\n    cartographic.y = map(vector[1]);\n    cartographic.z = vector[2];\n  }\n  return cartographic;\n}\n\nexport function toCartographicFromRadians(vector, cartographic) {\n  return toCartographic(vector, cartographic, config._cartographicRadians ? noop : toDegrees);\n}\n\nexport function toCartographicFromDegrees(vector, cartographic) {\n  return toCartographic(vector, cartographic, config._cartographicRadians ? toRadians : noop);\n}\n\nexport function isWGS84(vector) {\n  if (!vector) {\n    return false;\n  }\n  scratchVector.from(vector);\n  const {oneOverRadiiSquared, centerToleranceSquared} = WGS84_CONSTANTS;\n  const x2 = vector[0] * vector[0] * oneOverRadiiSquared[0];\n  const y2 = vector[1] * vector[1] * oneOverRadiiSquared[1];\n  const z2 = vector[2] * vector[2] * oneOverRadiiSquared[2];\n  return Math.abs(x2 + y2 + z2 - 1) < centerToleranceSquared;\n}\n"],"file":"type-utils.js"}
+{"version":3,"sources":["../../src/type-utils.ts"],"names":["identity","x","scratchVector","Vector3","fromCartographic","cartographic","result","map","longitude","latitude","height","y","z","fromCartographicToRadians","vector","config","_cartographicRadians","toRadians","fromCartographicToDegrees","toDegrees","toCartographic","toCartographicFromRadians","toCartographicFromDegrees","isWGS84","from","oneOverRadiiSquared","WGS84_CONSTANTS","centerToleranceSquared","x2","y2","z2","Math","abs"],"mappings":";;;;;;;;;;;;;AAIA;;AACA;;AAgBA,SAASA,QAAT,CAAkBC,CAAlB,EAAqC;AACnC,SAAOA,CAAP;AACD;;AAED,IAAMC,aAAa,GAAG,IAAIC,aAAJ,EAAtB;;AAQO,SAASC,gBAAT,CACLC,YADK,EAIK;AAAA,MAFVC,MAEU,uEAFD,EAEC;AAAA,MADVC,GACU,uEADJP,QACI;;AACV,MAAI,eAAeK,YAAnB,EAAiC;AAC/BC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACG,SAAd,CAAf;AACAF,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACI,QAAd,CAAf;AACAH,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYD,YAAY,CAACK,MAAzB;AACD,GAJD,MAIO,IAAI,OAAOL,YAAX,EAAyB;AAC9BC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACJ,CAAd,CAAf;AACAK,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACM,CAAd,CAAf;AACAL,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYD,YAAY,CAACO,CAAzB;AACD,GAJM,MAIA;AACLN,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAAC,CAAD,CAAb,CAAf;AACAC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAAC,CAAD,CAAb,CAAf;AACAC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYD,YAAY,CAAC,CAAD,CAAxB;AACD;;AACD,SAAOC,MAAP;AACD;;AAOM,SAASO,yBAAT,CACLR,YADK,EAGK;AAAA,MADVS,MACU,uEADD,EACC;AACV,SAAOV,gBAAgB,CAACC,YAAD,EAAeS,MAAf,EAAuBC,aAAOC,oBAAP,GAA8BhB,QAA9B,GAAyCiB,eAAhE,CAAvB;AACD;;AAOM,SAASC,yBAAT,CACLb,YADK,EAGK;AAAA,MADVS,MACU,uEADD,EACC;AACV,SAAOV,gBAAgB,CAACC,YAAD,EAAeS,MAAf,EAAuBC,aAAOC,oBAAP,GAA8BG,eAA9B,GAA0CnB,QAAjE,CAAvB;AACD;;AAEM,SAASoB,cAAT,CACLN,MADK,EAELT,YAFK,EAIF;AAAA,MADHE,GACG,uEAD0BP,QAC1B;;AACH,MAAI,eAAeK,YAAnB,EAAiC;AAC/BA,IAAAA,YAAY,CAACG,SAAb,GAAyBD,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAA5B;AACAT,IAAAA,YAAY,CAACI,QAAb,GAAwBF,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAA3B;AACAT,IAAAA,YAAY,CAACK,MAAb,GAAsBI,MAAM,CAAC,CAAD,CAA5B;AACD,GAJD,MAIO,IAAI,OAAOT,YAAX,EAAyB;AAC9BA,IAAAA,YAAY,CAACJ,CAAb,GAAiBM,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAApB;AACAT,IAAAA,YAAY,CAACM,CAAb,GAAiBJ,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAApB;AACAT,IAAAA,YAAY,CAACO,CAAb,GAAiBE,MAAM,CAAC,CAAD,CAAvB;AACD,GAJM,MAIA;AACLT,IAAAA,YAAY,CAAC,CAAD,CAAZ,GAAkBE,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAArB;AACAT,IAAAA,YAAY,CAAC,CAAD,CAAZ,GAAkBE,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAArB;AACAT,IAAAA,YAAY,CAAC,CAAD,CAAZ,GAAkBS,MAAM,CAAC,CAAD,CAAxB;AACD;;AACD,SAAOT,YAAP;AACD;;AAEM,SAASgB,yBAAT,CACLP,MADK,EAELT,YAFK,EAGF;AACH,SAAOe,cAAc,CAACN,MAAD,EAAST,YAAT,EAAuBU,aAAOC,oBAAP,GAA8BhB,QAA9B,GAAyCmB,eAAhE,CAArB;AACD;;AAEM,SAASG,yBAAT,CACLR,MADK,EAELT,YAFK,EAGF;AACH,SAAOe,cAAc,CAACN,MAAD,EAAST,YAAT,EAAuBU,aAAOC,oBAAP,GAA8BC,eAA9B,GAA0CjB,QAAjE,CAArB;AACD;;AAGM,SAASuB,OAAT,CAAiBT,MAAjB,EAA0D;AAC/D,MAAI,CAACA,MAAL,EAAa;AACX,WAAO,KAAP;AACD;;AACDZ,EAAAA,aAAa,CAACsB,IAAd,CAAmBV,MAAnB;AACA,MAAOW,mBAAP,GAAsDC,0BAAtD,CAAOD,mBAAP;AAAA,MAA4BE,sBAA5B,GAAsDD,0BAAtD,CAA4BC,sBAA5B;AACA,MAAMC,EAAE,GAAGd,MAAM,CAAC,CAAD,CAAN,GAAYA,MAAM,CAAC,CAAD,CAAlB,GAAwBW,mBAAmB,CAAC,CAAD,CAAtD;AACA,MAAMI,EAAE,GAAGf,MAAM,CAAC,CAAD,CAAN,GAAYA,MAAM,CAAC,CAAD,CAAlB,GAAwBW,mBAAmB,CAAC,CAAD,CAAtD;AACA,MAAMK,EAAE,GAAGhB,MAAM,CAAC,CAAD,CAAN,GAAYA,MAAM,CAAC,CAAD,CAAlB,GAAwBW,mBAAmB,CAAC,CAAD,CAAtD;AACA,SAAOM,IAAI,CAACC,GAAL,CAASJ,EAAE,GAAGC,EAAL,GAAUC,EAAV,GAAe,CAAxB,IAA6BH,sBAApC;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 type {NumericArray} from '@math.gl/core';\nimport {Vector3, toRadians, toDegrees, config} from '@math.gl/core';\nimport {WGS84_CONSTANTS} from './constants';\n\ntype LngLatHeightObject = {\n  longitude: number;\n  latitude: number;\n  height: number;\n};\n\ntype XYZObject = {\n  x: number;\n  y: number;\n  z: number;\n};\n\ntype Cartographic = LngLatHeightObject | XYZObject | NumericArray;\n\nfunction identity(x: number): number {\n  return x;\n}\n\nconst scratchVector = new Vector3();\n\nexport function fromCartographic(cartographic: Cartographic): number[];\nexport function fromCartographic<NumArrayT>(\n  cartographic: Cartographic,\n  result: NumArrayT,\n  map?: (x: number) => number\n): NumArrayT;\nexport function fromCartographic(\n  cartographic: Cartographic,\n  result = [] as number[],\n  map = identity\n): number[] {\n  if ('longitude' in cartographic) {\n    result[0] = map(cartographic.longitude);\n    result[1] = map(cartographic.latitude);\n    result[2] = cartographic.height;\n  } else if ('x' in cartographic) {\n    result[0] = map(cartographic.x);\n    result[1] = map(cartographic.y);\n    result[2] = cartographic.z;\n  } else {\n    result[0] = map(cartographic[0]);\n    result[1] = map(cartographic[1]);\n    result[2] = cartographic[2];\n  }\n  return result;\n}\n\nexport function fromCartographicToRadians(cartographic: Cartographic, result?: number[]): number[];\nexport function fromCartographicToRadians<TArray>(\n  cartographic: Cartographic,\n  result: TArray\n): TArray;\nexport function fromCartographicToRadians(\n  cartographic: Cartographic,\n  vector = [] as number[]\n): number[] {\n  return fromCartographic(cartographic, vector, config._cartographicRadians ? identity : toRadians);\n}\n\nexport function fromCartographicToDegrees(cartographic: Cartographic, result?: number[]): number[];\nexport function fromCartographicToDegrees<TArray>(\n  cartographic: Cartographic,\n  result: TArray\n): TArray;\nexport function fromCartographicToDegrees(\n  cartographic: Cartographic,\n  vector = [] as number[]\n): number[] {\n  return fromCartographic(cartographic, vector, config._cartographicRadians ? toDegrees : identity);\n}\n\nexport function toCartographic<T extends Cartographic>(\n  vector: Readonly<NumericArray>,\n  cartographic: T,\n  map: (x: number) => number = identity\n): T {\n  if ('longitude' in cartographic) {\n    cartographic.longitude = map(vector[0]);\n    cartographic.latitude = map(vector[1]);\n    cartographic.height = vector[2];\n  } else if ('x' in cartographic) {\n    cartographic.x = map(vector[0]);\n    cartographic.y = map(vector[1]);\n    cartographic.z = vector[2];\n  } else {\n    cartographic[0] = map(vector[0]);\n    cartographic[1] = map(vector[1]);\n    cartographic[2] = vector[2];\n  }\n  return cartographic;\n}\n\nexport function toCartographicFromRadians<T extends Cartographic>(\n  vector: Readonly<NumericArray>,\n  cartographic: T\n): T {\n  return toCartographic(vector, cartographic, config._cartographicRadians ? identity : toDegrees);\n}\n\nexport function toCartographicFromDegrees<T extends Cartographic>(\n  vector: Readonly<NumericArray>,\n  cartographic: T\n): T {\n  return toCartographic(vector, cartographic, config._cartographicRadians ? toRadians : identity);\n}\n\n// Estimates if a vector is close to the surface of the WGS84 Ellipsoid\nexport function isWGS84(vector: Readonly<NumericArray>): boolean {\n  if (!vector) {\n    return false;\n  }\n  scratchVector.from(vector);\n  const {oneOverRadiiSquared, centerToleranceSquared} = WGS84_CONSTANTS;\n  const x2 = vector[0] * vector[0] * oneOverRadiiSquared[0];\n  const y2 = vector[1] * vector[1] * oneOverRadiiSquared[1];\n  const z2 = vector[2] * vector[2] * oneOverRadiiSquared[2];\n  return Math.abs(x2 + y2 + z2 - 1) < centerToleranceSquared;\n}\n\n/*\n\nexport function fromCartographic(cartographic: Cartographic, result?: number[]): number[];\nexport function fromCartographic(cartographic: Cartographic, result: TypedArray): TypedArray;\nexport function fromCartographicToRadians(cartographic: Cartographic, result?: number[]): number[];\nexport function fromCartographicToRadians(\n  cartographic: Cartographic,\n  result: TypedArray\n): TypedArray;\nexport function fromCartographicToDegrees(cartographic: Cartographic, result?: number[]): number[];\nexport function fromCartographicToDegrees(\n  cartographic: Cartographic,\n  result: TypedArray\n): TypedArray;\n\nexport function toCartographic(vector: number[] | TypedArray, result: Cartographic): number[];\nexport function toCartographicFromRadians(\n  vector: number[] | TypedArray,\n  result: Cartographic\n): number[];\nexport function toCartographicFromDegrees(\n  vector: number[] | TypedArray,\n  result: Cartographic\n): number[];\n\n// Estimates if a vector is close to the surface of the WGS84 Ellipsoid\nexport function isWGS84(vector: number[] | TypedArray): boolean;\n*/\n"],"file":"type-utils.js"}

package/dist/esm/constants.js.map

@@ -1 +1 @@
-{"version":3,"sources":["../../src/constants.js"],"names":["WGS84_RADIUS_X","WGS84_RADIUS_Y","WGS84_RADIUS_Z","WGS84_CONSTANTS","radii","radiiSquared","oneOverRadii","oneOverRadiiSquared","maximumRadius","Math","max","centerToleranceSquared"],"mappings":"AAGA,OAAO,MAAMA,cAAc,GAAG,SAAvB;AACP,OAAO,MAAMC,cAAc,GAAG,SAAvB;AACP,OAAO,MAAMC,cAAc,GAAG,kBAAvB;AAIP,OAAO,MAAMC,eAAe,GAAG;AAC7BC,EAAAA,KAAK,EAAE,CAACJ,cAAD,EAAiBC,cAAjB,EAAiCC,cAAjC,CADsB;AAE7BG,EAAAA,YAAY,EAAE,CACZL,cAAc,GAAGA,cADL,EAEZC,cAAc,GAAGA,cAFL,EAGZC,cAAc,GAAGA,cAHL,CAFe;AAO7BI,EAAAA,YAAY,EAAE,CAAC,MAAMN,cAAP,EAAuB,MAAMC,cAA7B,EAA6C,MAAMC,cAAnD,CAPe;AAQ7BK,EAAAA,mBAAmB,EAAE,CACnB,OAAOP,cAAc,GAAGA,cAAxB,CADmB,EAEnB,OAAOC,cAAc,GAAGA,cAAxB,CAFmB,EAGnB,OAAOC,cAAc,GAAGA,cAAxB,CAHmB,CARQ;AAa7BM,EAAAA,aAAa,EAAEC,IAAI,CAACC,GAAL,CAASV,cAAT,EAAyBC,cAAzB,EAAyCC,cAAzC,CAbc;AAc7BS,EAAAA,sBAAsB,EAAE;AAdK,CAAxB","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\nexport const WGS84_RADIUS_X = 6378137.0;\nexport const WGS84_RADIUS_Y = 6378137.0;\nexport const WGS84_RADIUS_Z = 6356752.3142451793;\n\n// Pre-calculated ellipsoid defaults to avoid utils depending on `ellipsoid.js`\n\nexport const WGS84_CONSTANTS = {\n  radii: [WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z],\n  radiiSquared: [\n    WGS84_RADIUS_X * WGS84_RADIUS_X,\n    WGS84_RADIUS_Y * WGS84_RADIUS_Y,\n    WGS84_RADIUS_Z * WGS84_RADIUS_Z\n  ],\n  oneOverRadii: [1.0 / WGS84_RADIUS_X, 1.0 / WGS84_RADIUS_Y, 1.0 / WGS84_RADIUS_Z],\n  oneOverRadiiSquared: [\n    1.0 / (WGS84_RADIUS_X * WGS84_RADIUS_X),\n    1.0 / (WGS84_RADIUS_Y * WGS84_RADIUS_Y),\n    1.0 / (WGS84_RADIUS_Z * WGS84_RADIUS_Z)\n  ],\n  maximumRadius: Math.max(WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z),\n  centerToleranceSquared: 1e-1 // EPSILON1;\n};\n"],"file":"constants.js"}
+{"version":3,"sources":["../../src/constants.ts"],"names":["WGS84_RADIUS_X","WGS84_RADIUS_Y","WGS84_RADIUS_Z","WGS84_CONSTANTS","radii","radiiSquared","oneOverRadii","oneOverRadiiSquared","maximumRadius","Math","max","centerToleranceSquared"],"mappings":"AAGA,OAAO,MAAMA,cAAc,GAAG,SAAvB;AACP,OAAO,MAAMC,cAAc,GAAG,SAAvB;AACP,OAAO,MAAMC,cAAc,GAAG,kBAAvB;AAIP,OAAO,MAAMC,eAAe,GAAG;AAC7BC,EAAAA,KAAK,EAAE,CAACJ,cAAD,EAAiBC,cAAjB,EAAiCC,cAAjC,CADsB;AAE7BG,EAAAA,YAAY,EAAE,CACZL,cAAc,GAAGA,cADL,EAEZC,cAAc,GAAGA,cAFL,EAGZC,cAAc,GAAGA,cAHL,CAFe;AAO7BI,EAAAA,YAAY,EAAE,CAAC,MAAMN,cAAP,EAAuB,MAAMC,cAA7B,EAA6C,MAAMC,cAAnD,CAPe;AAQ7BK,EAAAA,mBAAmB,EAAE,CACnB,OAAOP,cAAc,GAAGA,cAAxB,CADmB,EAEnB,OAAOC,cAAc,GAAGA,cAAxB,CAFmB,EAGnB,OAAOC,cAAc,GAAGA,cAAxB,CAHmB,CARQ;AAa7BM,EAAAA,aAAa,EAAEC,IAAI,CAACC,GAAL,CAASV,cAAT,EAAyBC,cAAzB,EAAyCC,cAAzC,CAbc;AAc7BS,EAAAA,sBAAsB,EAAE;AAdK,CAAxB","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\nexport const WGS84_RADIUS_X = 6378137.0;\nexport const WGS84_RADIUS_Y = 6378137.0;\nexport const WGS84_RADIUS_Z = 6356752.3142451793;\n\n// Pre-calculated ellipsoid defaults to avoid utils depending on `ellipsoid.js`\n\nexport const WGS84_CONSTANTS = {\n  radii: [WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z],\n  radiiSquared: [\n    WGS84_RADIUS_X * WGS84_RADIUS_X,\n    WGS84_RADIUS_Y * WGS84_RADIUS_Y,\n    WGS84_RADIUS_Z * WGS84_RADIUS_Z\n  ],\n  oneOverRadii: [1.0 / WGS84_RADIUS_X, 1.0 / WGS84_RADIUS_Y, 1.0 / WGS84_RADIUS_Z],\n  oneOverRadiiSquared: [\n    1.0 / (WGS84_RADIUS_X * WGS84_RADIUS_X),\n    1.0 / (WGS84_RADIUS_Y * WGS84_RADIUS_Y),\n    1.0 / (WGS84_RADIUS_Z * WGS84_RADIUS_Z)\n  ],\n  maximumRadius: Math.max(WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z),\n  centerToleranceSquared: 1e-1 // EPSILON1;\n};\n"],"file":"constants.js"}

package/dist/esm/ellipsoid/ellipsoid.js

@@ -1,4 +1,5 @@
-import { Vector3, Matrix4, toRadians, toDegrees, assert, equals, _MathUtils } from '@math.gl/core';
+import _defineProperty from "@babel/runtime/helpers/esm/defineProperty";
+import { Vector3, Matrix4, assert, equals, _MathUtils } from '@math.gl/core';
 import * as vec3 from 'gl-matrix/vec3';
 import { WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z } from '../constants';
 import { fromCartographicToRadians, toCartographicFromRadians } from '../type-utils';
@@ -12,12 +13,25 @@ const scratchHeight = new Vector3();
 const scratchCartesian = new Vector3();
 let wgs84;
 export default class Ellipsoid {
-  static get WGS84() {
-    wgs84 = wgs84 || new Ellipsoid(WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z);
-    return wgs84;
-  }
-
   constructor(x = 0.0, y = 0.0, z = 0.0) {
+    _defineProperty(this, "radii", void 0);
+
+    _defineProperty(this, "radiiSquared", void 0);
+
+    _defineProperty(this, "radiiToTheFourth", void 0);
+
+    _defineProperty(this, "oneOverRadii", void 0);
+
+    _defineProperty(this, "oneOverRadiiSquared", void 0);
+
+    _defineProperty(this, "minimumRadius", void 0);
+
+    _defineProperty(this, "maximumRadius", void 0);
+
+    _defineProperty(this, "centerToleranceSquared", _MathUtils.EPSILON1);
+
+    _defineProperty(this, "squaredXOverSquaredZ", void 0);
+
     assert(x >= 0.0);
     assert(y >= 0.0);
     assert(z >= 0.0);
@@ -28,7 +42,6 @@ export default class Ellipsoid {
     this.oneOverRadiiSquared = new Vector3(x === 0.0 ? 0.0 : 1.0 / (x * x), y === 0.0 ? 0.0 : 1.0 / (y * y), z === 0.0 ? 0.0 : 1.0 / (z * z));
     this.minimumRadius = Math.min(x, y, z);
     this.maximumRadius = Math.max(x, y, z);
-    this.centerToleranceSquared = _MathUtils.EPSILON1;
 
     if (this.radiiSquared.z !== 0) {
       this.squaredXOverSquaredZ = this.radiiSquared.x / this.radiiSquared.z;
@@ -122,7 +135,7 @@ export default class Ellipsoid {
     return scratchPosition.from(position).scale(this.radii).to(result);
   }
 
-  getSurfaceNormalIntersectionWithZAxis(position, buffer = 0.0, result = [0, 0, 0]) {
+  getSurfaceNormalIntersectionWithZAxis(position, buffer = 0, result = [0, 0, 0]) {
     assert(equals(this.radii.x, this.radii.y, _MathUtils.EPSILON15));
     assert(this.radii.z > 0);
     scratchPosition.from(position);
@@ -136,4 +149,6 @@ export default class Ellipsoid {
   }
 
 }
+
+_defineProperty(Ellipsoid, "WGS84", new Ellipsoid(WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z));
 //# sourceMappingURL=ellipsoid.js.map

package/dist/esm/ellipsoid/ellipsoid.js.map

@@ -1 +1 @@
-{"version":3,"sources":["../../../src/ellipsoid/ellipsoid.js"],"names":["Vector3","Matrix4","toRadians","toDegrees","assert","equals","_MathUtils","vec3","WGS84_RADIUS_X","WGS84_RADIUS_Y","WGS84_RADIUS_Z","fromCartographicToRadians","toCartographicFromRadians","scaleToGeodeticSurface","localFrameToFixedFrame","scratchVector","scratchNormal","scratchK","scratchPosition","scratchHeight","scratchCartesian","wgs84","Ellipsoid","WGS84","constructor","x","y","z","radii","radiiSquared","radiiToTheFourth","oneOverRadii","oneOverRadiiSquared","minimumRadius","Math","min","maximumRadius","max","centerToleranceSquared","EPSILON1","squaredXOverSquaredZ","Object","freeze","right","Boolean","toString","cartographicToCartesian","cartographic","result","normal","k","height","geodeticSurfaceNormalCartographic","copy","scale","gamma","sqrt","dot","add","to","cartesianToCartographic","cartesian","from","point","undefined","geodeticSurfaceNormal","h","subtract","longitude","atan2","latitude","asin","sign","length","eastNorthUpToFixedFrame","origin","firstAxis","secondAxis","thirdAxis","geocentricSurfaceNormal","normalize","cartographicVectorRadians","cosLatitude","cos","set","sin","scaleToGeocentricSurface","positionX","positionY","positionZ","beta","multiplyScalar","transformPositionToScaledSpace","position","transformPositionFromScaledSpace","getSurfaceNormalIntersectionWithZAxis","buffer","EPSILON15","abs"],"mappings":"AAIA,SAAQA,OAAR,EAAiBC,OAAjB,EAA0BC,SAA1B,EAAqCC,SAArC,EAAgDC,MAAhD,EAAwDC,MAAxD,EAAgEC,UAAhE,QAAiF,eAAjF;AACA,OAAO,KAAKC,IAAZ,MAAsB,gBAAtB;AAEA,SAAQC,cAAR,EAAwBC,cAAxB,EAAwCC,cAAxC,QAA6D,cAA7D;AACA,SAAQC,yBAAR,EAAmCC,yBAAnC,QAAmE,eAAnE;AAEA,OAAOC,sBAAP,MAAmC,qCAAnC;AACA,OAAOC,sBAAP,MAAmC,+BAAnC;AAEA,MAAMC,aAAa,GAAG,IAAIf,OAAJ,EAAtB;AACA,MAAMgB,aAAa,GAAG,IAAIhB,OAAJ,EAAtB;AACA,MAAMiB,QAAQ,GAAG,IAAIjB,OAAJ,EAAjB;AACA,MAAMkB,eAAe,GAAG,IAAIlB,OAAJ,EAAxB;AACA,MAAMmB,aAAa,GAAG,IAAInB,OAAJ,EAAtB;AACA,MAAMoB,gBAAgB,GAAG,IAAIpB,OAAJ,EAAzB;AAEA,IAAIqB,KAAJ;AAKA,eAAe,MAAMC,SAAN,CAAgB;AAEb,aAALC,KAAK,GAAG;AACjBF,IAAAA,KAAK,GAAGA,KAAK,IAAI,IAAIC,SAAJ,CAAcd,cAAd,EAA8BC,cAA9B,EAA8CC,cAA9C,CAAjB;AACA,WAAOW,KAAP;AACD;;AAGDG,EAAAA,WAAW,CAACC,CAAC,GAAG,GAAL,EAAUC,CAAC,GAAG,GAAd,EAAmBC,CAAC,GAAG,GAAvB,EAA4B;AACrCvB,IAAAA,MAAM,CAACqB,CAAC,IAAI,GAAN,CAAN;AACArB,IAAAA,MAAM,CAACsB,CAAC,IAAI,GAAN,CAAN;AACAtB,IAAAA,MAAM,CAACuB,CAAC,IAAI,GAAN,CAAN;AAEA,SAAKC,KAAL,GAAa,IAAI5B,OAAJ,CAAYyB,CAAZ,EAAeC,CAAf,EAAkBC,CAAlB,CAAb;AAEA,SAAKE,YAAL,GAAoB,IAAI7B,OAAJ,CAAYyB,CAAC,GAAGA,CAAhB,EAAmBC,CAAC,GAAGA,CAAvB,EAA0BC,CAAC,GAAGA,CAA9B,CAApB;AAEA,SAAKG,gBAAL,GAAwB,IAAI9B,OAAJ,CAAYyB,CAAC,GAAGA,CAAJ,GAAQA,CAAR,GAAYA,CAAxB,EAA2BC,CAAC,GAAGA,CAAJ,GAAQA,CAAR,GAAYA,CAAvC,EAA0CC,CAAC,GAAGA,CAAJ,GAAQA,CAAR,GAAYA,CAAtD,CAAxB;AAEA,SAAKI,YAAL,GAAoB,IAAI/B,OAAJ,CAClByB,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,MAAMA,CADN,EAElBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,MAAMA,CAFN,EAGlBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,MAAMA,CAHN,CAApB;AAMA,SAAKK,mBAAL,GAA2B,IAAIhC,OAAJ,CACzByB,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,OAAOA,CAAC,GAAGA,CAAX,CADO,EAEzBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,OAAOA,CAAC,GAAGA,CAAX,CAFO,EAGzBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,OAAOA,CAAC,GAAGA,CAAX,CAHO,CAA3B;AAMA,SAAKM,aAAL,GAAqBC,IAAI,CAACC,GAAL,CAASV,CAAT,EAAYC,CAAZ,EAAeC,CAAf,CAArB;AAEA,SAAKS,aAAL,GAAqBF,IAAI,CAACG,GAAL,CAASZ,CAAT,EAAYC,CAAZ,EAAeC,CAAf,CAArB;AAEA,SAAKW,sBAAL,GAA8BhC,UAAU,CAACiC,QAAzC;;AAEA,QAAI,KAAKV,YAAL,CAAkBF,CAAlB,KAAwB,CAA5B,EAA+B;AAC7B,WAAKa,oBAAL,GAA4B,KAAKX,YAAL,CAAkBJ,CAAlB,GAAsB,KAAKI,YAAL,CAAkBF,CAApE;AACD;;AAEDc,IAAAA,MAAM,CAACC,MAAP,CAAc,IAAd;AACD;;AAGDrC,EAAAA,MAAM,CAACsC,KAAD,EAAQ;AACZ,WAAO,SAASA,KAAT,IAAkBC,OAAO,CAACD,KAAK,IAAI,KAAKf,KAAL,CAAWvB,MAAX,CAAkBsC,KAAK,CAACf,KAAxB,CAAV,CAAhC;AACD;;AAGDiB,EAAAA,QAAQ,GAAG;AACT,WAAO,KAAKjB,KAAL,CAAWiB,QAAX,EAAP;AACD;;AAGDC,EAAAA,uBAAuB,CAACC,YAAD,EAAeC,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAxB,EAAmC;AACxD,UAAMC,MAAM,GAAGjC,aAAf;AACA,UAAMkC,CAAC,GAAGjC,QAAV;AAEA,UAAM,IAAKkC,MAAL,IAAeJ,YAArB;AACA,SAAKK,iCAAL,CAAuCL,YAAvC,EAAqDE,MAArD;AACAC,IAAAA,CAAC,CAACG,IAAF,CAAO,KAAKxB,YAAZ,EAA0ByB,KAA1B,CAAgCL,MAAhC;AAEA,UAAMM,KAAK,GAAGrB,IAAI,CAACsB,IAAL,CAAUP,MAAM,CAACQ,GAAP,CAAWP,CAAX,CAAV,CAAd;AACAA,IAAAA,CAAC,CAACI,KAAF,CAAQ,IAAIC,KAAZ;AAEAN,IAAAA,MAAM,CAACK,KAAP,CAAaH,MAAb;AAEAD,IAAAA,CAAC,CAACQ,GAAF,CAAMT,MAAN;AAEA,WAAOC,CAAC,CAACS,EAAF,CAAKX,MAAL,CAAP;AACD;;AAIDY,EAAAA,uBAAuB,CAACC,SAAD,EAAYb,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAArB,EAAgC;AACrD5B,IAAAA,gBAAgB,CAAC0C,IAAjB,CAAsBD,SAAtB;AACA,UAAME,KAAK,GAAG,KAAKlD,sBAAL,CAA4BO,gBAA5B,EAA8CF,eAA9C,CAAd;;AAEA,QAAI,CAAC6C,KAAL,EAAY;AACV,aAAOC,SAAP;AACD;;AAED,UAAMf,MAAM,GAAG,KAAKgB,qBAAL,CAA2BF,KAA3B,EAAkC/C,aAAlC,CAAf;AAEA,UAAMkD,CAAC,GAAG/C,aAAV;AACA+C,IAAAA,CAAC,CAACb,IAAF,CAAOjC,gBAAP,EAAyB+C,QAAzB,CAAkCJ,KAAlC;AAEA,UAAMK,SAAS,GAAGlC,IAAI,CAACmC,KAAL,CAAWpB,MAAM,CAACvB,CAAlB,EAAqBuB,MAAM,CAACxB,CAA5B,CAAlB;AACA,UAAM6C,QAAQ,GAAGpC,IAAI,CAACqC,IAAL,CAAUtB,MAAM,CAACtB,CAAjB,CAAjB;AACA,UAAMwB,MAAM,GAAGjB,IAAI,CAACsC,IAAL,CAAUjE,IAAI,CAACkD,GAAL,CAASS,CAAT,EAAY9C,gBAAZ,CAAV,IAA2Cb,IAAI,CAACkE,MAAL,CAAYP,CAAZ,CAA1D;AAEA,WAAOtD,yBAAyB,CAAC,CAACwD,SAAD,EAAYE,QAAZ,EAAsBnB,MAAtB,CAAD,EAAgCH,MAAhC,CAAhC;AACD;;AAID0B,EAAAA,uBAAuB,CAACC,MAAD,EAAS3B,MAAM,GAAG,IAAI/C,OAAJ,EAAlB,EAAiC;AACtD,WAAOa,sBAAsB,CAAC,IAAD,EAAO,MAAP,EAAe,OAAf,EAAwB,IAAxB,EAA8B6D,MAA9B,EAAsC3B,MAAtC,CAA7B;AACD;;AAIDlC,EAAAA,sBAAsB,CAAC8D,SAAD,EAAYC,UAAZ,EAAwBC,SAAxB,EAAmCH,MAAnC,EAA2C3B,MAAM,GAAG,IAAI/C,OAAJ,EAApD,EAAmE;AACvF,WAAOa,sBAAsB,CAAC,IAAD,EAAO8D,SAAP,EAAkBC,UAAlB,EAA8BC,SAA9B,EAAyCH,MAAzC,EAAiD3B,MAAjD,CAA7B;AACD;;AAID+B,EAAAA,uBAAuB,CAAClB,SAAD,EAAYb,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAArB,EAAgC;AACrD,WAAOjC,aAAa,CAAC+C,IAAd,CAAmBD,SAAnB,EAA8BmB,SAA9B,GAA0CrB,EAA1C,CAA6CX,MAA7C,CAAP;AACD;;AAGDI,EAAAA,iCAAiC,CAACL,YAAD,EAAeC,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAxB,EAAmC;AAClE,UAAMiC,yBAAyB,GAAGtE,yBAAyB,CAACoC,YAAD,CAA3D;AAEA,UAAMqB,SAAS,GAAGa,yBAAyB,CAAC,CAAD,CAA3C;AACA,UAAMX,QAAQ,GAAGW,yBAAyB,CAAC,CAAD,CAA1C;AAEA,UAAMC,WAAW,GAAGhD,IAAI,CAACiD,GAAL,CAASb,QAAT,CAApB;AAEAvD,IAAAA,aAAa,CACVqE,GADH,CACOF,WAAW,GAAGhD,IAAI,CAACiD,GAAL,CAASf,SAAT,CADrB,EAC0Cc,WAAW,GAAGhD,IAAI,CAACmD,GAAL,CAASjB,SAAT,CADxD,EAC6ElC,IAAI,CAACmD,GAAL,CAASf,QAAT,CAD7E,EAEGU,SAFH;AAIA,WAAOjE,aAAa,CAAC4C,EAAd,CAAiBX,MAAjB,CAAP;AACD;;AAGDiB,EAAAA,qBAAqB,CAACJ,SAAD,EAAYb,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAArB,EAAgC;AACnD,WAAOjC,aAAa,CAAC+C,IAAd,CAAmBD,SAAnB,EAA8BP,KAA9B,CAAoC,KAAKtB,mBAAzC,EAA8DgD,SAA9D,GAA0ErB,EAA1E,CAA6EX,MAA7E,CAAP;AACD;;AAKDnC,EAAAA,sBAAsB,CAACgD,SAAD,EAAYb,MAAZ,EAAoB;AACxC,WAAOnC,sBAAsB,CAACgD,SAAD,EAAY,IAAZ,EAAkBb,MAAlB,CAA7B;AACD;;AAIDsC,EAAAA,wBAAwB,CAACzB,SAAD,EAAYb,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAArB,EAAgC;AACtD9B,IAAAA,eAAe,CAAC4C,IAAhB,CAAqBD,SAArB;AAEA,UAAM0B,SAAS,GAAGrE,eAAe,CAACO,CAAlC;AACA,UAAM+D,SAAS,GAAGtE,eAAe,CAACQ,CAAlC;AACA,UAAM+D,SAAS,GAAGvE,eAAe,CAACS,CAAlC;AACA,UAAMK,mBAAmB,GAAG,KAAKA,mBAAjC;AAEA,UAAM0D,IAAI,GACR,MACAxD,IAAI,CAACsB,IAAL,CACE+B,SAAS,GAAGA,SAAZ,GAAwBvD,mBAAmB,CAACP,CAA5C,GACE+D,SAAS,GAAGA,SAAZ,GAAwBxD,mBAAmB,CAACN,CAD9C,GAEE+D,SAAS,GAAGA,SAAZ,GAAwBzD,mBAAmB,CAACL,CAHhD,CAFF;AAQA,WAAOT,eAAe,CAACyE,cAAhB,CAA+BD,IAA/B,EAAqC/B,EAArC,CAAwCX,MAAxC,CAAP;AACD;;AAID4C,EAAAA,8BAA8B,CAACC,QAAD,EAAW7C,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAApB,EAA+B;AAC3D,WAAO9B,eAAe,CAAC4C,IAAhB,CAAqB+B,QAArB,EAA+BvC,KAA/B,CAAqC,KAAKvB,YAA1C,EAAwD4B,EAAxD,CAA2DX,MAA3D,CAAP;AACD;;AAID8C,EAAAA,gCAAgC,CAACD,QAAD,EAAW7C,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAApB,EAA+B;AAC7D,WAAO9B,eAAe,CAAC4C,IAAhB,CAAqB+B,QAArB,EAA+BvC,KAA/B,CAAqC,KAAK1B,KAA1C,EAAiD+B,EAAjD,CAAoDX,MAApD,CAAP;AACD;;AAGD+C,EAAAA,qCAAqC,CAACF,QAAD,EAAWG,MAAM,GAAG,GAApB,EAAyBhD,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAlC,EAA6C;AAEhF5C,IAAAA,MAAM,CAACC,MAAM,CAAC,KAAKuB,KAAL,CAAWH,CAAZ,EAAe,KAAKG,KAAL,CAAWF,CAA1B,EAA6BpB,UAAU,CAAC2F,SAAxC,CAAP,CAAN;AACA7F,IAAAA,MAAM,CAAC,KAAKwB,KAAL,CAAWD,CAAX,GAAe,CAAhB,CAAN;AAEAT,IAAAA,eAAe,CAAC4C,IAAhB,CAAqB+B,QAArB;AACA,UAAMlE,CAAC,GAAGT,eAAe,CAACS,CAAhB,IAAqB,IAAI,KAAKa,oBAA9B,CAAV;;AAEA,QAAIN,IAAI,CAACgE,GAAL,CAASvE,CAAT,KAAe,KAAKC,KAAL,CAAWD,CAAX,GAAeqE,MAAlC,EAA0C;AACxC,aAAOhC,SAAP;AACD;;AAED,WAAO9C,eAAe,CAACkE,GAAhB,CAAoB,GAApB,EAAyB,GAAzB,EAA8BzD,CAA9B,EAAiCgC,EAAjC,CAAoCX,MAApC,CAAP;AACD;;AA5L4B","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\n/* eslint-disable */\nimport {Vector3, Matrix4, toRadians, toDegrees, assert, equals, _MathUtils} from '@math.gl/core';\nimport * as vec3 from 'gl-matrix/vec3';\n\nimport {WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z} from '../constants';\nimport {fromCartographicToRadians, toCartographicFromRadians} from '../type-utils';\n\nimport scaleToGeodeticSurface from './helpers/scale-to-geodetic-surface';\nimport localFrameToFixedFrame from './helpers/ellipsoid-transform';\n\nconst scratchVector = new Vector3();\nconst scratchNormal = new Vector3();\nconst scratchK = new Vector3();\nconst scratchPosition = new Vector3();\nconst scratchHeight = new Vector3();\nconst scratchCartesian = new Vector3();\n\nlet wgs84;\n\n// A quadratic surface defined in Cartesian coordinates by the equation\n// <code>(x / a)^2 + (y / b)^2 + (z / c)^2 = 1</code>.  Primarily used\n// to represent the shape of planetary bodies.\nexport default class Ellipsoid {\n  // An Ellipsoid instance initialized to the WGS84 standard.\n  static get WGS84() {\n    wgs84 = wgs84 || new Ellipsoid(WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z);\n    return wgs84;\n  }\n\n  // Creates an Ellipsoid from a Cartesian specifying the radii in x, y, and z directions.\n  constructor(x = 0.0, y = 0.0, z = 0.0) {\n    assert(x >= 0.0);\n    assert(y >= 0.0);\n    assert(z >= 0.0);\n\n    this.radii = new Vector3(x, y, z);\n\n    this.radiiSquared = new Vector3(x * x, y * y, z * z);\n\n    this.radiiToTheFourth = new Vector3(x * x * x * x, y * y * y * y, z * z * z * z);\n\n    this.oneOverRadii = new Vector3(\n      x === 0.0 ? 0.0 : 1.0 / x,\n      y === 0.0 ? 0.0 : 1.0 / y,\n      z === 0.0 ? 0.0 : 1.0 / z\n    );\n\n    this.oneOverRadiiSquared = new Vector3(\n      x === 0.0 ? 0.0 : 1.0 / (x * x),\n      y === 0.0 ? 0.0 : 1.0 / (y * y),\n      z === 0.0 ? 0.0 : 1.0 / (z * z)\n    );\n\n    this.minimumRadius = Math.min(x, y, z);\n\n    this.maximumRadius = Math.max(x, y, z);\n\n    this.centerToleranceSquared = _MathUtils.EPSILON1;\n\n    if (this.radiiSquared.z !== 0) {\n      this.squaredXOverSquaredZ = this.radiiSquared.x / this.radiiSquared.z;\n    }\n\n    Object.freeze(this);\n  }\n\n  // Compares this Ellipsoid against the provided Ellipsoid componentwise and returns\n  equals(right) {\n    return this === right || Boolean(right && this.radii.equals(right.radii));\n  }\n\n  // Creates a string representing this Ellipsoid in the format '(radii.x, radii.y, radii.z)'.\n  toString() {\n    return this.radii.toString();\n  }\n\n  // Converts the provided cartographic to Cartesian representation.\n  cartographicToCartesian(cartographic, result = [0, 0, 0]) {\n    const normal = scratchNormal;\n    const k = scratchK;\n\n    const [, , height] = cartographic;\n    this.geodeticSurfaceNormalCartographic(cartographic, normal);\n    k.copy(this.radiiSquared).scale(normal);\n\n    const gamma = Math.sqrt(normal.dot(k));\n    k.scale(1 / gamma);\n\n    normal.scale(height);\n\n    k.add(normal);\n\n    return k.to(result);\n  }\n\n  // Converts the provided cartesian to cartographic (lng/lat/z) representation.\n  // The cartesian is undefined at the center of the ellipsoid.\n  cartesianToCartographic(cartesian, result = [0, 0, 0]) {\n    scratchCartesian.from(cartesian);\n    const point = this.scaleToGeodeticSurface(scratchCartesian, scratchPosition);\n\n    if (!point) {\n      return undefined;\n    }\n\n    const normal = this.geodeticSurfaceNormal(point, scratchNormal);\n\n    const h = scratchHeight;\n    h.copy(scratchCartesian).subtract(point);\n\n    const longitude = Math.atan2(normal.y, normal.x);\n    const latitude = Math.asin(normal.z);\n    const height = Math.sign(vec3.dot(h, scratchCartesian)) * vec3.length(h);\n\n    return toCartographicFromRadians([longitude, latitude, height], result);\n  }\n\n  // Computes a 4x4 transformation matrix from a reference frame with an east-north-up axes\n  // centered at the provided origin to the provided ellipsoid's fixed reference frame.\n  eastNorthUpToFixedFrame(origin, result = new Matrix4()) {\n    return localFrameToFixedFrame(this, 'east', 'north', 'up', origin, result);\n  }\n\n  // Computes a 4x4 transformation matrix from a reference frame centered at\n  // the provided origin to the ellipsoid's fixed reference frame.\n  localFrameToFixedFrame(firstAxis, secondAxis, thirdAxis, origin, result = new Matrix4()) {\n    return localFrameToFixedFrame(this, firstAxis, secondAxis, thirdAxis, origin, result);\n  }\n\n  // Computes the unit vector directed from the center of this ellipsoid toward\n  // the provided Cartesian position.\n  geocentricSurfaceNormal(cartesian, result = [0, 0, 0]) {\n    return scratchVector.from(cartesian).normalize().to(result);\n  }\n\n  // Computes the normal of the plane tangent to the surface of the ellipsoid at provided position.\n  geodeticSurfaceNormalCartographic(cartographic, result = [0, 0, 0]) {\n    const cartographicVectorRadians = fromCartographicToRadians(cartographic);\n\n    const longitude = cartographicVectorRadians[0];\n    const latitude = cartographicVectorRadians[1];\n\n    const cosLatitude = Math.cos(latitude);\n\n    scratchVector\n      .set(cosLatitude * Math.cos(longitude), cosLatitude * Math.sin(longitude), Math.sin(latitude))\n      .normalize();\n\n    return scratchVector.to(result);\n  }\n\n  // Computes the normal of the plane tangent to the surface of the ellipsoid at the provided position.\n  geodeticSurfaceNormal(cartesian, result = [0, 0, 0]) {\n    return scratchVector.from(cartesian).scale(this.oneOverRadiiSquared).normalize().to(result);\n  }\n\n  // Scales the provided Cartesian position along the geodetic surface normal\n  // so that it is on the surface of this ellipsoid.  If the position is\n  // at the center of the ellipsoid, this function returns undefined.\n  scaleToGeodeticSurface(cartesian, result) {\n    return scaleToGeodeticSurface(cartesian, this, result);\n  }\n\n  // Scales the provided Cartesian position along the geocentric surface normal\n  // so that it is on the surface of this ellipsoid.\n  scaleToGeocentricSurface(cartesian, result = [0, 0, 0]) {\n    scratchPosition.from(cartesian);\n\n    const positionX = scratchPosition.x;\n    const positionY = scratchPosition.y;\n    const positionZ = scratchPosition.z;\n    const oneOverRadiiSquared = this.oneOverRadiiSquared;\n\n    const beta =\n      1.0 /\n      Math.sqrt(\n        positionX * positionX * oneOverRadiiSquared.x +\n          positionY * positionY * oneOverRadiiSquared.y +\n          positionZ * positionZ * oneOverRadiiSquared.z\n      );\n\n    return scratchPosition.multiplyScalar(beta).to(result);\n  }\n\n  // Transforms a Cartesian X, Y, Z position to the ellipsoid-scaled space by multiplying\n  // its components by the result of `Ellipsoid#oneOverRadii`\n  transformPositionToScaledSpace(position, result = [0, 0, 0]) {\n    return scratchPosition.from(position).scale(this.oneOverRadii).to(result);\n  }\n\n  // Transforms a Cartesian X, Y, Z position from the ellipsoid-scaled space by multiplying\n  // its components by the result of `Ellipsoid#radii`.\n  transformPositionFromScaledSpace(position, result = [0, 0, 0]) {\n    return scratchPosition.from(position).scale(this.radii).to(result);\n  }\n\n  // Computes a point which is the intersection of the surface normal with the z-axis.\n  getSurfaceNormalIntersectionWithZAxis(position, buffer = 0.0, result = [0, 0, 0]) {\n    // Ellipsoid must be an ellipsoid of revolution (radii.x == radii.y)\n    assert(equals(this.radii.x, this.radii.y, _MathUtils.EPSILON15));\n    assert(this.radii.z > 0);\n\n    scratchPosition.from(position);\n    const z = scratchPosition.z * (1 - this.squaredXOverSquaredZ);\n\n    if (Math.abs(z) >= this.radii.z - buffer) {\n      return undefined;\n    }\n\n    return scratchPosition.set(0.0, 0.0, z).to(result);\n  }\n}\n"],"file":"ellipsoid.js"}
+{"version":3,"sources":["../../../src/ellipsoid/ellipsoid.ts"],"names":["Vector3","Matrix4","assert","equals","_MathUtils","vec3","WGS84_RADIUS_X","WGS84_RADIUS_Y","WGS84_RADIUS_Z","fromCartographicToRadians","toCartographicFromRadians","scaleToGeodeticSurface","localFrameToFixedFrame","scratchVector","scratchNormal","scratchK","scratchPosition","scratchHeight","scratchCartesian","wgs84","Ellipsoid","constructor","x","y","z","EPSILON1","radii","radiiSquared","radiiToTheFourth","oneOverRadii","oneOverRadiiSquared","minimumRadius","Math","min","maximumRadius","max","squaredXOverSquaredZ","Object","freeze","right","Boolean","toString","cartographicToCartesian","cartographic","result","normal","k","height","geodeticSurfaceNormalCartographic","copy","scale","gamma","sqrt","dot","add","to","cartesianToCartographic","cartesian","from","point","undefined","geodeticSurfaceNormal","h","subtract","longitude","atan2","latitude","asin","sign","length","eastNorthUpToFixedFrame","origin","firstAxis","secondAxis","thirdAxis","geocentricSurfaceNormal","normalize","cartographicVectorRadians","cosLatitude","cos","set","sin","scaleToGeocentricSurface","positionX","positionY","positionZ","beta","multiplyScalar","transformPositionToScaledSpace","position","transformPositionFromScaledSpace","getSurfaceNormalIntersectionWithZAxis","buffer","EPSILON15","abs"],"mappings":";AAIA,SACEA,OADF,EAEEC,OAFF,EAKEC,MALF,EAMEC,MANF,EAOEC,UAPF,QASO,eATP;AAUA,OAAO,KAAKC,IAAZ,MAAsB,gBAAtB;AAEA,SAAQC,cAAR,EAAwBC,cAAxB,EAAwCC,cAAxC,QAA6D,cAA7D;AACA,SAAQC,yBAAR,EAAmCC,yBAAnC,QAAmE,eAAnE;AAEA,OAAOC,sBAAP,MAAmC,qCAAnC;AACA,OAAOC,sBAAP,MAAmC,+BAAnC;AAEA,MAAMC,aAAa,GAAG,IAAIb,OAAJ,EAAtB;AACA,MAAMc,aAAa,GAAG,IAAId,OAAJ,EAAtB;AACA,MAAMe,QAAQ,GAAG,IAAIf,OAAJ,EAAjB;AACA,MAAMgB,eAAe,GAAG,IAAIhB,OAAJ,EAAxB;AACA,MAAMiB,aAAa,GAAG,IAAIjB,OAAJ,EAAtB;AACA,MAAMkB,gBAAgB,GAAG,IAAIlB,OAAJ,EAAzB;AAEA,IAAImB,KAAJ;AAOA,eAAe,MAAMC,SAAN,CAAgB;AAkB7BC,EAAAA,WAAW,CAACC,CAAC,GAAG,GAAL,EAAUC,CAAC,GAAG,GAAd,EAAmBC,CAAC,GAAG,GAAvB,EAA4B;AAAA;;AAAA;;AAAA;;AAAA;;AAAA;;AAAA;;AAAA;;AAAA,oDAPGpB,UAAU,CAACqB,QAOd;;AAAA;;AACrCvB,IAAAA,MAAM,CAACoB,CAAC,IAAI,GAAN,CAAN;AACApB,IAAAA,MAAM,CAACqB,CAAC,IAAI,GAAN,CAAN;AACArB,IAAAA,MAAM,CAACsB,CAAC,IAAI,GAAN,CAAN;AAEA,SAAKE,KAAL,GAAa,IAAI1B,OAAJ,CAAYsB,CAAZ,EAAeC,CAAf,EAAkBC,CAAlB,CAAb;AAEA,SAAKG,YAAL,GAAoB,IAAI3B,OAAJ,CAAYsB,CAAC,GAAGA,CAAhB,EAAmBC,CAAC,GAAGA,CAAvB,EAA0BC,CAAC,GAAGA,CAA9B,CAApB;AAEA,SAAKI,gBAAL,GAAwB,IAAI5B,OAAJ,CAAYsB,CAAC,GAAGA,CAAJ,GAAQA,CAAR,GAAYA,CAAxB,EAA2BC,CAAC,GAAGA,CAAJ,GAAQA,CAAR,GAAYA,CAAvC,EAA0CC,CAAC,GAAGA,CAAJ,GAAQA,CAAR,GAAYA,CAAtD,CAAxB;AAEA,SAAKK,YAAL,GAAoB,IAAI7B,OAAJ,CAClBsB,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,MAAMA,CADN,EAElBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,MAAMA,CAFN,EAGlBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,MAAMA,CAHN,CAApB;AAMA,SAAKM,mBAAL,GAA2B,IAAI9B,OAAJ,CACzBsB,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,OAAOA,CAAC,GAAGA,CAAX,CADO,EAEzBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,OAAOA,CAAC,GAAGA,CAAX,CAFO,EAGzBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,OAAOA,CAAC,GAAGA,CAAX,CAHO,CAA3B;AAMA,SAAKO,aAAL,GAAqBC,IAAI,CAACC,GAAL,CAASX,CAAT,EAAYC,CAAZ,EAAeC,CAAf,CAArB;AAEA,SAAKU,aAAL,GAAqBF,IAAI,CAACG,GAAL,CAASb,CAAT,EAAYC,CAAZ,EAAeC,CAAf,CAArB;;AAEA,QAAI,KAAKG,YAAL,CAAkBH,CAAlB,KAAwB,CAA5B,EAA+B;AAC7B,WAAKY,oBAAL,GAA4B,KAAKT,YAAL,CAAkBL,CAAlB,GAAsB,KAAKK,YAAL,CAAkBH,CAApE;AACD;;AAEDa,IAAAA,MAAM,CAACC,MAAP,CAAc,IAAd;AACD;;AAGDnC,EAAAA,MAAM,CAACoC,KAAD,EAA4B;AAChC,WAAO,SAASA,KAAT,IAAkBC,OAAO,CAACD,KAAK,IAAI,KAAKb,KAAL,CAAWvB,MAAX,CAAkBoC,KAAK,CAACb,KAAxB,CAAV,CAAhC;AACD;;AAGDe,EAAAA,QAAQ,GAAW;AACjB,WAAO,KAAKf,KAAL,CAAWe,QAAX,EAAP;AACD;;AAMDC,EAAAA,uBAAuB,CAACC,YAAD,EAAeC,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAxB,EAAmC;AACxD,UAAMC,MAAM,GAAG/B,aAAf;AACA,UAAMgC,CAAC,GAAG/B,QAAV;AAEA,UAAM,IAAKgC,MAAL,IAAeJ,YAArB;AACA,SAAKK,iCAAL,CAAuCL,YAAvC,EAAqDE,MAArD;AACAC,IAAAA,CAAC,CAACG,IAAF,CAAO,KAAKtB,YAAZ,EAA0BuB,KAA1B,CAAgCL,MAAhC;AAEA,UAAMM,KAAK,GAAGnB,IAAI,CAACoB,IAAL,CAAUP,MAAM,CAACQ,GAAP,CAAWP,CAAX,CAAV,CAAd;AACAA,IAAAA,CAAC,CAACI,KAAF,CAAQ,IAAIC,KAAZ;AAEAN,IAAAA,MAAM,CAACK,KAAP,CAAaH,MAAb;AAEAD,IAAAA,CAAC,CAACQ,GAAF,CAAMT,MAAN;AAEA,WAAOC,CAAC,CAACS,EAAF,CAAKX,MAAL,CAAP;AACD;;AAODY,EAAAA,uBAAuB,CAACC,SAAD,EAAYb,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAArB,EAAgC;AACrD1B,IAAAA,gBAAgB,CAACwC,IAAjB,CAAsBD,SAAtB;AACA,UAAME,KAAK,GAAG,KAAKhD,sBAAL,CAA4BO,gBAA5B,EAA8CF,eAA9C,CAAd;;AAEA,QAAI,CAAC2C,KAAL,EAAY;AACV,aAAOC,SAAP;AACD;;AAED,UAAMf,MAAM,GAAG,KAAKgB,qBAAL,CAA2BF,KAA3B,EAAkC7C,aAAlC,CAAf;AAEA,UAAMgD,CAAC,GAAG7C,aAAV;AACA6C,IAAAA,CAAC,CAACb,IAAF,CAAO/B,gBAAP,EAAyB6C,QAAzB,CAAkCJ,KAAlC;AAEA,UAAMK,SAAS,GAAGhC,IAAI,CAACiC,KAAL,CAAWpB,MAAM,CAACtB,CAAlB,EAAqBsB,MAAM,CAACvB,CAA5B,CAAlB;AACA,UAAM4C,QAAQ,GAAGlC,IAAI,CAACmC,IAAL,CAAUtB,MAAM,CAACrB,CAAjB,CAAjB;AACA,UAAMuB,MAAM,GAAGf,IAAI,CAACoC,IAAL,CAAU/D,IAAI,CAACgD,GAAL,CAASS,CAAT,EAAY5C,gBAAZ,CAAV,IAA2Cb,IAAI,CAACgE,MAAL,CAAYP,CAAZ,CAA1D;AAEA,WAAOpD,yBAAyB,CAAC,CAACsD,SAAD,EAAYE,QAAZ,EAAsBnB,MAAtB,CAAD,EAAgCH,MAAhC,CAAhC;AACD;;AAOD0B,EAAAA,uBAAuB,CAACC,MAAD,EAAS3B,MAAM,GAAG,IAAI3C,OAAJ,EAAlB,EAAiC;AACtD,WAAOW,sBAAsB,CAAC,IAAD,EAAO,MAAP,EAAe,OAAf,EAAwB,IAAxB,EAA8B2D,MAA9B,EAAsC3B,MAAtC,CAA7B;AACD;;AAsBDhC,EAAAA,sBAAsB,CAAC4D,SAAD,EAAYC,UAAZ,EAAwBC,SAAxB,EAAmCH,MAAnC,EAA2C3B,MAAM,GAAG,IAAI3C,OAAJ,EAApD,EAAmE;AACvF,WAAOW,sBAAsB,CAAC,IAAD,EAAO4D,SAAP,EAAkBC,UAAlB,EAA8BC,SAA9B,EAAyCH,MAAzC,EAAiD3B,MAAjD,CAA7B;AACD;;AAMD+B,EAAAA,uBAAuB,CAAClB,SAAD,EAAYb,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAArB,EAAgC;AACrD,WAAO/B,aAAa,CAAC6C,IAAd,CAAmBD,SAAnB,EAA8BmB,SAA9B,GAA0CrB,EAA1C,CAA6CX,MAA7C,CAAP;AACD;;AAKDI,EAAAA,iCAAiC,CAACL,YAAD,EAAeC,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAxB,EAAmC;AAClE,UAAMiC,yBAAyB,GAAGpE,yBAAyB,CAACkC,YAAD,CAA3D;AAEA,UAAMqB,SAAS,GAAGa,yBAAyB,CAAC,CAAD,CAA3C;AACA,UAAMX,QAAQ,GAAGW,yBAAyB,CAAC,CAAD,CAA1C;AAEA,UAAMC,WAAW,GAAG9C,IAAI,CAAC+C,GAAL,CAASb,QAAT,CAApB;AAEArD,IAAAA,aAAa,CACVmE,GADH,CACOF,WAAW,GAAG9C,IAAI,CAAC+C,GAAL,CAASf,SAAT,CADrB,EAC0Cc,WAAW,GAAG9C,IAAI,CAACiD,GAAL,CAASjB,SAAT,CADxD,EAC6EhC,IAAI,CAACiD,GAAL,CAASf,QAAT,CAD7E,EAEGU,SAFH;AAIA,WAAO/D,aAAa,CAAC0C,EAAd,CAAiBX,MAAjB,CAAP;AACD;;AAKDiB,EAAAA,qBAAqB,CAACJ,SAAD,EAAYb,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAArB,EAAgC;AACnD,WAAO/B,aAAa,CAAC6C,IAAd,CAAmBD,SAAnB,EAA8BP,KAA9B,CAAoC,KAAKpB,mBAAzC,EAA8D8C,SAA9D,GAA0ErB,EAA1E,CAA6EX,MAA7E,CAAP;AACD;;AAKDjC,EAAAA,sBAAsB,CAAC8C,SAAD,EAAsBb,MAAtB,EAAmD;AACvE,WAAOjC,sBAAsB,CAAC8C,SAAD,EAAY,IAAZ,EAAkBb,MAAlB,CAA7B;AACD;;AAIDsC,EAAAA,wBAAwB,CAACzB,SAAD,EAAsBb,MAAgB,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAzC,EAA8D;AACpF5B,IAAAA,eAAe,CAAC0C,IAAhB,CAAqBD,SAArB;AAEA,UAAM0B,SAAS,GAAGnE,eAAe,CAACM,CAAlC;AACA,UAAM8D,SAAS,GAAGpE,eAAe,CAACO,CAAlC;AACA,UAAM8D,SAAS,GAAGrE,eAAe,CAACQ,CAAlC;AACA,UAAMM,mBAAmB,GAAG,KAAKA,mBAAjC;AAEA,UAAMwD,IAAI,GACR,MACAtD,IAAI,CAACoB,IAAL,CACE+B,SAAS,GAAGA,SAAZ,GAAwBrD,mBAAmB,CAACR,CAA5C,GACE8D,SAAS,GAAGA,SAAZ,GAAwBtD,mBAAmB,CAACP,CAD9C,GAEE8D,SAAS,GAAGA,SAAZ,GAAwBvD,mBAAmB,CAACN,CAHhD,CAFF;AAQA,WAAOR,eAAe,CAACuE,cAAhB,CAA+BD,IAA/B,EAAqC/B,EAArC,CAAwCX,MAAxC,CAAP;AACD;;AAID4C,EAAAA,8BAA8B,CAACC,QAAD,EAAqB7C,MAAgB,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAxC,EAA6D;AACzF,WAAO5B,eAAe,CAAC0C,IAAhB,CAAqB+B,QAArB,EAA+BvC,KAA/B,CAAqC,KAAKrB,YAA1C,EAAwD0B,EAAxD,CAA2DX,MAA3D,CAAP;AACD;;AAID8C,EAAAA,gCAAgC,CAACD,QAAD,EAAqB7C,MAAgB,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAxC,EAA6D;AAC3F,WAAO5B,eAAe,CAAC0C,IAAhB,CAAqB+B,QAArB,EAA+BvC,KAA/B,CAAqC,KAAKxB,KAA1C,EAAiD6B,EAAjD,CAAoDX,MAApD,CAAP;AACD;;AAGD+C,EAAAA,qCAAqC,CACnCF,QADmC,EAEnCG,MAAc,GAAG,CAFkB,EAGnChD,MAAgB,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAHgB,EAIzB;AAEV1C,IAAAA,MAAM,CAACC,MAAM,CAAC,KAAKuB,KAAL,CAAWJ,CAAZ,EAAe,KAAKI,KAAL,CAAWH,CAA1B,EAA6BnB,UAAU,CAACyF,SAAxC,CAAP,CAAN;AACA3F,IAAAA,MAAM,CAAC,KAAKwB,KAAL,CAAWF,CAAX,GAAe,CAAhB,CAAN;AAEAR,IAAAA,eAAe,CAAC0C,IAAhB,CAAqB+B,QAArB;AACA,UAAMjE,CAAC,GAAGR,eAAe,CAACQ,CAAhB,IAAqB,IAAI,KAAKY,oBAA9B,CAAV;;AAEA,QAAIJ,IAAI,CAAC8D,GAAL,CAAStE,CAAT,KAAe,KAAKE,KAAL,CAAWF,CAAX,GAAeoE,MAAlC,EAA0C;AACxC,aAAOhC,SAAP;AACD;;AAED,WAAO5C,eAAe,CAACgE,GAAhB,CAAoB,GAApB,EAAyB,GAAzB,EAA8BxD,CAA9B,EAAiC+B,EAAjC,CAAoCX,MAApC,CAAP;AACD;;AAzO4B;;gBAAVxB,S,WAEgB,IAAIA,SAAJ,CAAcd,cAAd,EAA8BC,cAA9B,EAA8CC,cAA9C,C","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\n/* eslint-disable */\nimport {\n  Vector3,\n  Matrix4,\n  toRadians,\n  toDegrees,\n  assert,\n  equals,\n  _MathUtils,\n  NumericArray\n} from '@math.gl/core';\nimport * as vec3 from 'gl-matrix/vec3';\n\nimport {WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z} from '../constants';\nimport {fromCartographicToRadians, toCartographicFromRadians} from '../type-utils';\n\nimport scaleToGeodeticSurface from './helpers/scale-to-geodetic-surface';\nimport localFrameToFixedFrame from './helpers/ellipsoid-transform';\n\nconst scratchVector = new Vector3();\nconst scratchNormal = new Vector3();\nconst scratchK = new Vector3();\nconst scratchPosition = new Vector3();\nconst scratchHeight = new Vector3();\nconst scratchCartesian = new Vector3();\n\nlet wgs84;\n\n/**\n * A quadratic surface defined in Cartesian coordinates by the equation\n * `(x / a)^2 + (y / b)^2 + (z / c)^2 = 1`.  Primarily used\n * to represent the shape of planetary bodies.\n */\nexport default class Ellipsoid {\n  /** An Ellipsoid instance initialized to the WGS84 standard. */\n  static readonly WGS84: Ellipsoid = new Ellipsoid(WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z);\n\n  readonly radii: Vector3;\n  readonly radiiSquared: Vector3;\n  readonly radiiToTheFourth: Vector3;\n  readonly oneOverRadii: Vector3;\n  readonly oneOverRadiiSquared: Vector3;\n  readonly minimumRadius: number;\n  readonly maximumRadius: number;\n  readonly centerToleranceSquared: number = _MathUtils.EPSILON1;\n  readonly squaredXOverSquaredZ: number;\n\n  /** Creates an Ellipsoid from a Cartesian specifying the radii in x, y, and z directions. */\n  constructor(x: number, y: number, z: number);\n  constructor();\n\n  constructor(x = 0.0, y = 0.0, z = 0.0) {\n    assert(x >= 0.0);\n    assert(y >= 0.0);\n    assert(z >= 0.0);\n\n    this.radii = new Vector3(x, y, z);\n\n    this.radiiSquared = new Vector3(x * x, y * y, z * z);\n\n    this.radiiToTheFourth = new Vector3(x * x * x * x, y * y * y * y, z * z * z * z);\n\n    this.oneOverRadii = new Vector3(\n      x === 0.0 ? 0.0 : 1.0 / x,\n      y === 0.0 ? 0.0 : 1.0 / y,\n      z === 0.0 ? 0.0 : 1.0 / z\n    );\n\n    this.oneOverRadiiSquared = new Vector3(\n      x === 0.0 ? 0.0 : 1.0 / (x * x),\n      y === 0.0 ? 0.0 : 1.0 / (y * y),\n      z === 0.0 ? 0.0 : 1.0 / (z * z)\n    );\n\n    this.minimumRadius = Math.min(x, y, z);\n\n    this.maximumRadius = Math.max(x, y, z);\n\n    if (this.radiiSquared.z !== 0) {\n      this.squaredXOverSquaredZ = this.radiiSquared.x / this.radiiSquared.z;\n    }\n\n    Object.freeze(this);\n  }\n\n  /** Compares this Ellipsoid against the provided Ellipsoid componentwise */\n  equals(right: Ellipsoid): boolean {\n    return this === right || Boolean(right && this.radii.equals(right.radii));\n  }\n\n  /** Creates a string representing this Ellipsoid in the format '(radii.x, radii.y, radii.z)'. */\n  toString(): string {\n    return this.radii.toString();\n  }\n\n  /** Converts the provided cartographic to Cartesian representation. */\n  cartographicToCartesian(cartographic: number[], result: Vector3): Vector3;\n  cartographicToCartesian(cartographic: number[], result?: number[]): number[];\n\n  cartographicToCartesian(cartographic, result = [0, 0, 0]) {\n    const normal = scratchNormal;\n    const k = scratchK;\n\n    const [, , height] = cartographic;\n    this.geodeticSurfaceNormalCartographic(cartographic, normal);\n    k.copy(this.radiiSquared).scale(normal);\n\n    const gamma = Math.sqrt(normal.dot(k));\n    k.scale(1 / gamma);\n\n    normal.scale(height);\n\n    k.add(normal);\n\n    return k.to(result);\n  }\n\n  /** Converts the provided cartesian to cartographic (lng/lat/z) representation.\n   * The cartesian is undefined at the center of the ellipsoid. */\n  cartesianToCartographic(cartesian: number[], result: Vector3): Vector3;\n  cartesianToCartographic(cartesian: number[], result?: number[]): number[];\n\n  cartesianToCartographic(cartesian, result = [0, 0, 0]) {\n    scratchCartesian.from(cartesian);\n    const point = this.scaleToGeodeticSurface(scratchCartesian, scratchPosition);\n\n    if (!point) {\n      return undefined;\n    }\n\n    const normal = this.geodeticSurfaceNormal(point, scratchNormal);\n\n    const h = scratchHeight;\n    h.copy(scratchCartesian).subtract(point);\n\n    const longitude = Math.atan2(normal.y, normal.x);\n    const latitude = Math.asin(normal.z);\n    const height = Math.sign(vec3.dot(h, scratchCartesian)) * vec3.length(h);\n\n    return toCartographicFromRadians([longitude, latitude, height], result);\n  }\n\n  /** Computes a 4x4 transformation matrix from a reference frame with an east-north-up axes\n   * centered at the provided origin to the provided ellipsoid's fixed reference frame. */\n  eastNorthUpToFixedFrame(origin: number[], result?: Matrix4): Matrix4;\n  eastNorthUpToFixedFrame(origin: number[], result: number[]): number[];\n\n  eastNorthUpToFixedFrame(origin, result = new Matrix4()) {\n    return localFrameToFixedFrame(this, 'east', 'north', 'up', origin, result);\n  }\n\n  /** Computes a 4x4 transformation matrix from a reference frame centered at\n   * the provided origin to the ellipsoid's fixed reference frame.\n   */\n  localFrameToFixedFrame(\n    firstAxis: string,\n    secondAxis: string,\n    thirdAxis: string,\n    origin: Readonly<NumericArray>,\n    result?: Matrix4\n  ): Matrix4;\n  localFrameToFixedFrame<Matrix4T>(\n    firstAxis: string,\n    secondAxis: string,\n    thirdAxis: string,\n    origin: Readonly<NumericArray>,\n    result: number[]\n  ): number[];\n\n  // Computes a 4x4 transformation matrix from a reference frame centered at\n  // the provided origin to the ellipsoid's fixed reference frame.\n  localFrameToFixedFrame(firstAxis, secondAxis, thirdAxis, origin, result = new Matrix4()) {\n    return localFrameToFixedFrame(this, firstAxis, secondAxis, thirdAxis, origin, result);\n  }\n\n  /** Computes the unit vector directed from the center of this ellipsoid toward\n   * the provided Cartesian position. */\n  geocentricSurfaceNormal(cartesian: number[], result?: number[]): number[];\n  geocentricSurfaceNormal<NumArray>(cartesian: number[], result: NumArray): NumArray;\n  geocentricSurfaceNormal(cartesian, result = [0, 0, 0]) {\n    return scratchVector.from(cartesian).normalize().to(result);\n  }\n\n  /** Computes the normal of the plane tangent to the surface of the ellipsoid at provided position. */\n  geodeticSurfaceNormalCartographic<NumArray>(cartographic: number[], result: NumArray): NumArray;\n  geodeticSurfaceNormalCartographic(cartographic: number[]): number[];\n  geodeticSurfaceNormalCartographic(cartographic, result = [0, 0, 0]) {\n    const cartographicVectorRadians = fromCartographicToRadians(cartographic);\n\n    const longitude = cartographicVectorRadians[0];\n    const latitude = cartographicVectorRadians[1];\n\n    const cosLatitude = Math.cos(latitude);\n\n    scratchVector\n      .set(cosLatitude * Math.cos(longitude), cosLatitude * Math.sin(longitude), Math.sin(latitude))\n      .normalize();\n\n    return scratchVector.to(result);\n  }\n\n  /** Computes the normal of the plane tangent to the surface of the ellipsoid at the provided position. */\n  geodeticSurfaceNormal<NumArrayT>(cartesian: number[], result: NumArrayT): NumArrayT;\n  geodeticSurfaceNormal(cartesian: number[]): number[];\n  geodeticSurfaceNormal(cartesian, result = [0, 0, 0]) {\n    return scratchVector.from(cartesian).scale(this.oneOverRadiiSquared).normalize().to(result);\n  }\n\n  /** Scales the provided Cartesian position along the geodetic surface normal\n   * so that it is on the surface of this ellipsoid.  If the position is\n   * at the center of the ellipsoid, this function returns undefined. */\n  scaleToGeodeticSurface(cartesian: number[], result?: number[]): number[] {\n    return scaleToGeodeticSurface(cartesian, this, result);\n  }\n\n  /** Scales the provided Cartesian position along the geocentric surface normal\n   * so that it is on the surface of this ellipsoid. */\n  scaleToGeocentricSurface(cartesian: number[], result: number[] = [0, 0, 0]): number[] {\n    scratchPosition.from(cartesian);\n\n    const positionX = scratchPosition.x;\n    const positionY = scratchPosition.y;\n    const positionZ = scratchPosition.z;\n    const oneOverRadiiSquared = this.oneOverRadiiSquared;\n\n    const beta =\n      1.0 /\n      Math.sqrt(\n        positionX * positionX * oneOverRadiiSquared.x +\n          positionY * positionY * oneOverRadiiSquared.y +\n          positionZ * positionZ * oneOverRadiiSquared.z\n      );\n\n    return scratchPosition.multiplyScalar(beta).to(result);\n  }\n\n  /** Transforms a Cartesian X, Y, Z position to the ellipsoid-scaled space by multiplying\n   * its components by the result of `Ellipsoid#oneOverRadii` */\n  transformPositionToScaledSpace(position: number[], result: number[] = [0, 0, 0]): number[] {\n    return scratchPosition.from(position).scale(this.oneOverRadii).to(result);\n  }\n\n  /** Transforms a Cartesian X, Y, Z position from the ellipsoid-scaled space by multiplying\n   * its components by the result of `Ellipsoid#radii`. */\n  transformPositionFromScaledSpace(position: number[], result: number[] = [0, 0, 0]): number[] {\n    return scratchPosition.from(position).scale(this.radii).to(result);\n  }\n\n  /** Computes a point which is the intersection of the surface normal with the z-axis. */\n  getSurfaceNormalIntersectionWithZAxis(\n    position: number[],\n    buffer: number = 0,\n    result: number[] = [0, 0, 0]\n  ): number[] {\n    // Ellipsoid must be an ellipsoid of revolution (radii.x == radii.y)\n    assert(equals(this.radii.x, this.radii.y, _MathUtils.EPSILON15));\n    assert(this.radii.z > 0);\n\n    scratchPosition.from(position);\n    const z = scratchPosition.z * (1 - this.squaredXOverSquaredZ);\n\n    if (Math.abs(z) >= this.radii.z - buffer) {\n      return undefined;\n    }\n\n    return scratchPosition.set(0.0, 0.0, z).to(result);\n  }\n}\n"],"file":"ellipsoid.js"}

package/dist/esm/ellipsoid/helpers/ellipsoid-transform.js.map

@@ -1 +1 @@
-{"version":3,"sources":["../../../../src/ellipsoid/helpers/ellipsoid-transform.js"],"names":["Vector3","assert","equals","equalsEpsilon","EPSILON14","scratchOrigin","VECTOR_PRODUCT_LOCAL_FRAME","up","south","north","west","east","down","degeneratePositionLocalFrame","scratchAxisVectors","scratchVector1","scratchVector2","scratchVector3","localFrameToFixedFrame","ellipsoid","firstAxis","secondAxis","thirdAxis","cartesianOrigin","result","thirdAxisInferred","firstAxisVector","secondAxisVector","thirdAxisVector","origin","copy","atPole","x","y","sign","Math","z","fromArray","scale","set","normalize","geodeticSurfaceNormal","cross"],"mappings":"AAAA,SAAQA,OAAR,EAAiBC,MAAjB,EAAyBC,MAAM,IAAIC,aAAnC,QAAuD,eAAvD;AAEA,MAAMC,SAAS,GAAG,KAAlB;AAEA,MAAMC,aAAa,GAAG,IAAIL,OAAJ,EAAtB;AAGA,MAAMM,0BAA0B,GAAG;AACjCC,EAAAA,EAAE,EAAE;AACFC,IAAAA,KAAK,EAAE,MADL;AAEFC,IAAAA,KAAK,EAAE,MAFL;AAGFC,IAAAA,IAAI,EAAE,OAHJ;AAIFC,IAAAA,IAAI,EAAE;AAJJ,GAD6B;AAOjCC,EAAAA,IAAI,EAAE;AACJJ,IAAAA,KAAK,EAAE,MADH;AAEJC,IAAAA,KAAK,EAAE,MAFH;AAGJC,IAAAA,IAAI,EAAE,OAHF;AAIJC,IAAAA,IAAI,EAAE;AAJF,GAP2B;AAajCH,EAAAA,KAAK,EAAE;AACLD,IAAAA,EAAE,EAAE,MADC;AAELK,IAAAA,IAAI,EAAE,MAFD;AAGLF,IAAAA,IAAI,EAAE,MAHD;AAILC,IAAAA,IAAI,EAAE;AAJD,GAb0B;AAmBjCF,EAAAA,KAAK,EAAE;AACLF,IAAAA,EAAE,EAAE,MADC;AAELK,IAAAA,IAAI,EAAE,MAFD;AAGLF,IAAAA,IAAI,EAAE,IAHD;AAILC,IAAAA,IAAI,EAAE;AAJD,GAnB0B;AAyBjCD,EAAAA,IAAI,EAAE;AACJH,IAAAA,EAAE,EAAE,OADA;AAEJK,IAAAA,IAAI,EAAE,OAFF;AAGJH,IAAAA,KAAK,EAAE,MAHH;AAIJD,IAAAA,KAAK,EAAE;AAJH,GAzB2B;AA+BjCG,EAAAA,IAAI,EAAE;AACJJ,IAAAA,EAAE,EAAE,OADA;AAEJK,IAAAA,IAAI,EAAE,OAFF;AAGJH,IAAAA,KAAK,EAAE,IAHH;AAIJD,IAAAA,KAAK,EAAE;AAJH;AA/B2B,CAAnC;AAuCA,MAAMK,4BAA4B,GAAG;AACnCJ,EAAAA,KAAK,EAAE,CAAC,CAAC,CAAF,EAAK,CAAL,EAAQ,CAAR,CAD4B;AAEnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAF6B;AAGnCJ,EAAAA,EAAE,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAH+B;AAInCC,EAAAA,KAAK,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAJ4B;AAKnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAC,CAAL,EAAQ,CAAR,CAL6B;AAMnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAC,CAAR;AAN6B,CAArC;AASA,MAAME,kBAAkB,GAAG;AACzBH,EAAAA,IAAI,EAAE,IAAIX,OAAJ,EADmB;AAEzBS,EAAAA,KAAK,EAAE,IAAIT,OAAJ,EAFkB;AAGzBO,EAAAA,EAAE,EAAE,IAAIP,OAAJ,EAHqB;AAIzBU,EAAAA,IAAI,EAAE,IAAIV,OAAJ,EAJmB;AAKzBQ,EAAAA,KAAK,EAAE,IAAIR,OAAJ,EALkB;AAMzBY,EAAAA,IAAI,EAAE,IAAIZ,OAAJ;AANmB,CAA3B;AASA,MAAMe,cAAc,GAAG,IAAIf,OAAJ,EAAvB;AACA,MAAMgB,cAAc,GAAG,IAAIhB,OAAJ,EAAvB;AACA,MAAMiB,cAAc,GAAG,IAAIjB,OAAJ,EAAvB;AAKA,eAAe,SAASkB,sBAAT,CACbC,SADa,EAEbC,SAFa,EAGbC,UAHa,EAIbC,SAJa,EAKbC,eALa,EAMbC,MANa,EAOb;AACA,QAAMC,iBAAiB,GACrBnB,0BAA0B,CAACc,SAAD,CAA1B,IAAyCd,0BAA0B,CAACc,SAAD,CAA1B,CAAsCC,UAAtC,CAD3C;AAGApB,EAAAA,MAAM,CAACwB,iBAAiB,KAAK,CAACH,SAAD,IAAcA,SAAS,KAAKG,iBAAjC,CAAlB,CAAN;AAEA,MAAIC,eAAJ;AACA,MAAIC,gBAAJ;AACA,MAAIC,eAAJ;AAEA,QAAMC,MAAM,GAAGxB,aAAa,CAACyB,IAAd,CAAmBP,eAAnB,CAAf;AAGA,QAAMQ,MAAM,GAAG5B,aAAa,CAAC0B,MAAM,CAACG,CAAR,EAAW,GAAX,EAAgB5B,SAAhB,CAAb,IAA2CD,aAAa,CAAC0B,MAAM,CAACI,CAAR,EAAW,GAAX,EAAgB7B,SAAhB,CAAvE;;AAEA,MAAI2B,MAAJ,EAAY;AAEV,UAAMG,IAAI,GAAGC,IAAI,CAACD,IAAL,CAAUL,MAAM,CAACO,CAAjB,CAAb;AAEAV,IAAAA,eAAe,GAAGX,cAAc,CAACsB,SAAf,CAAyBxB,4BAA4B,CAACO,SAAD,CAArD,CAAlB;;AACA,QAAIA,SAAS,KAAK,MAAd,IAAwBA,SAAS,KAAK,MAA1C,EAAkD;AAChDM,MAAAA,eAAe,CAACY,KAAhB,CAAsBJ,IAAtB;AACD;;AAEDP,IAAAA,gBAAgB,GAAGX,cAAc,CAACqB,SAAf,CAAyBxB,4BAA4B,CAACQ,UAAD,CAArD,CAAnB;;AACA,QAAIA,UAAU,KAAK,MAAf,IAAyBA,UAAU,KAAK,MAA5C,EAAoD;AAClDM,MAAAA,gBAAgB,CAACW,KAAjB,CAAuBJ,IAAvB;AACD;;AAEDN,IAAAA,eAAe,GAAGX,cAAc,CAACoB,SAAf,CAAyBxB,4BAA4B,CAACS,SAAD,CAArD,CAAlB;;AACA,QAAIA,SAAS,KAAK,MAAd,IAAwBA,SAAS,KAAK,MAA1C,EAAkD;AAChDM,MAAAA,eAAe,CAACU,KAAhB,CAAsBJ,IAAtB;AACD;AACF,GAlBD,MAkBO;AAEL,UAAM;AAAC3B,MAAAA,EAAD;AAAKI,MAAAA,IAAL;AAAWF,MAAAA;AAAX,QAAoBK,kBAA1B;AAEAH,IAAAA,IAAI,CAAC4B,GAAL,CAAS,CAACV,MAAM,CAACI,CAAjB,EAAoBJ,MAAM,CAACG,CAA3B,EAA8B,GAA9B,EAAmCQ,SAAnC;AACArB,IAAAA,SAAS,CAACsB,qBAAV,CAAgCZ,MAAhC,EAAwCtB,EAAxC;AACAE,IAAAA,KAAK,CAACqB,IAAN,CAAWvB,EAAX,EAAemC,KAAf,CAAqB/B,IAArB;AAEA,UAAM;AAACC,MAAAA,IAAD;AAAOF,MAAAA,IAAP;AAAaF,MAAAA;AAAb,QAAsBM,kBAA5B;AAEAF,IAAAA,IAAI,CAACkB,IAAL,CAAUvB,EAAV,EAAc+B,KAAd,CAAoB,CAAC,CAArB;AACA5B,IAAAA,IAAI,CAACoB,IAAL,CAAUnB,IAAV,EAAgB2B,KAAhB,CAAsB,CAAC,CAAvB;AACA9B,IAAAA,KAAK,CAACsB,IAAN,CAAWrB,KAAX,EAAkB6B,KAAlB,CAAwB,CAAC,CAAzB;AAGAZ,IAAAA,eAAe,GAAGZ,kBAAkB,CAACM,SAAD,CAApC;AACAO,IAAAA,gBAAgB,GAAGb,kBAAkB,CAACO,UAAD,CAArC;AACAO,IAAAA,eAAe,GAAGd,kBAAkB,CAACQ,SAAD,CAApC;AACD;;AAGDE,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACM,CAA5B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACO,CAA5B;AACAT,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACU,CAA5B;AACAZ,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAY,GAAZ;AACAA,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACK,CAA7B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACM,CAA7B;AACAT,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACS,CAA7B;AACAZ,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAY,GAAZ;AACAA,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYI,eAAe,CAACI,CAA5B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYI,eAAe,CAACK,CAA5B;AACAT,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaI,eAAe,CAACQ,CAA7B;AACAZ,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAa,GAAb;AACAA,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACG,CAApB;AACAR,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACI,CAApB;AACAT,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACO,CAApB;AACAZ,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAa,GAAb;AACA,SAAOA,MAAP;AACD","sourcesContent":["import {Vector3, assert, equals as equalsEpsilon} from '@math.gl/core';\n\nconst EPSILON14 = 1e-14;\n\nconst scratchOrigin = new Vector3();\n\n// Caclulate third axis from given two axii\nconst VECTOR_PRODUCT_LOCAL_FRAME = {\n  up: {\n    south: 'east',\n    north: 'west',\n    west: 'south',\n    east: 'north'\n  },\n  down: {\n    south: 'west',\n    north: 'east',\n    west: 'north',\n    east: 'south'\n  },\n  south: {\n    up: 'west',\n    down: 'east',\n    west: 'down',\n    east: 'up'\n  },\n  north: {\n    up: 'east',\n    down: 'west',\n    west: 'up',\n    east: 'down'\n  },\n  west: {\n    up: 'north',\n    down: 'south',\n    north: 'down',\n    south: 'up'\n  },\n  east: {\n    up: 'south',\n    down: 'north',\n    north: 'up',\n    south: 'down'\n  }\n};\n\nconst degeneratePositionLocalFrame = {\n  north: [-1, 0, 0],\n  east: [0, 1, 0],\n  up: [0, 0, 1],\n  south: [1, 0, 0],\n  west: [0, -1, 0],\n  down: [0, 0, -1]\n};\n\nconst scratchAxisVectors = {\n  east: new Vector3(),\n  north: new Vector3(),\n  up: new Vector3(),\n  west: new Vector3(),\n  south: new Vector3(),\n  down: new Vector3()\n};\n\nconst scratchVector1 = new Vector3();\nconst scratchVector2 = new Vector3();\nconst scratchVector3 = new Vector3();\n\n// Computes a 4x4 transformation matrix from a reference frame\n// centered at the provided origin to the provided ellipsoid's fixed reference frame.\n// eslint-disable-next-line max-statements, max-params, complexity\nexport default function localFrameToFixedFrame(\n  ellipsoid,\n  firstAxis,\n  secondAxis,\n  thirdAxis,\n  cartesianOrigin,\n  result\n) {\n  const thirdAxisInferred =\n    VECTOR_PRODUCT_LOCAL_FRAME[firstAxis] && VECTOR_PRODUCT_LOCAL_FRAME[firstAxis][secondAxis];\n  // firstAxis and secondAxis must be east, north, up, west, south or down.');\n  assert(thirdAxisInferred && (!thirdAxis || thirdAxis === thirdAxisInferred));\n\n  let firstAxisVector;\n  let secondAxisVector;\n  let thirdAxisVector;\n\n  const origin = scratchOrigin.copy(cartesianOrigin);\n\n  // If x and y are zero, assume origin is at a pole, which is a special case.\n  const atPole = equalsEpsilon(origin.x, 0.0, EPSILON14) && equalsEpsilon(origin.y, 0.0, EPSILON14);\n\n  if (atPole) {\n    // Look up axis value and adjust\n    const sign = Math.sign(origin.z);\n\n    firstAxisVector = scratchVector1.fromArray(degeneratePositionLocalFrame[firstAxis]);\n    if (firstAxis !== 'east' && firstAxis !== 'west') {\n      firstAxisVector.scale(sign);\n    }\n\n    secondAxisVector = scratchVector2.fromArray(degeneratePositionLocalFrame[secondAxis]);\n    if (secondAxis !== 'east' && secondAxis !== 'west') {\n      secondAxisVector.scale(sign);\n    }\n\n    thirdAxisVector = scratchVector3.fromArray(degeneratePositionLocalFrame[thirdAxis]);\n    if (thirdAxis !== 'east' && thirdAxis !== 'west') {\n      thirdAxisVector.scale(sign);\n    }\n  } else {\n    // Calculate all axis\n    const {up, east, north} = scratchAxisVectors;\n\n    east.set(-origin.y, origin.x, 0.0).normalize();\n    ellipsoid.geodeticSurfaceNormal(origin, up);\n    north.copy(up).cross(east);\n\n    const {down, west, south} = scratchAxisVectors;\n\n    down.copy(up).scale(-1);\n    west.copy(east).scale(-1);\n    south.copy(north).scale(-1);\n\n    // Pick three axis based on desired orientation\n    firstAxisVector = scratchAxisVectors[firstAxis];\n    secondAxisVector = scratchAxisVectors[secondAxis];\n    thirdAxisVector = scratchAxisVectors[thirdAxis];\n  }\n\n  // TODO - assuming the result is column-major\n  result[0] = firstAxisVector.x;\n  result[1] = firstAxisVector.y;\n  result[2] = firstAxisVector.z;\n  result[3] = 0.0;\n  result[4] = secondAxisVector.x;\n  result[5] = secondAxisVector.y;\n  result[6] = secondAxisVector.z;\n  result[7] = 0.0;\n  result[8] = thirdAxisVector.x;\n  result[9] = thirdAxisVector.y;\n  result[10] = thirdAxisVector.z;\n  result[11] = 0.0;\n  result[12] = origin.x;\n  result[13] = origin.y;\n  result[14] = origin.z;\n  result[15] = 1.0;\n  return result;\n}\n"],"file":"ellipsoid-transform.js"}
+{"version":3,"sources":["../../../../src/ellipsoid/helpers/ellipsoid-transform.ts"],"names":["Vector3","assert","equals","equalsEpsilon","EPSILON14","scratchOrigin","VECTOR_PRODUCT_LOCAL_FRAME","up","south","north","west","east","down","degeneratePositionLocalFrame","scratchAxisVectors","scratchVector1","scratchVector2","scratchVector3","localFrameToFixedFrame","ellipsoid","firstAxis","secondAxis","thirdAxis","cartesianOrigin","result","thirdAxisInferred","firstAxisVector","secondAxisVector","thirdAxisVector","origin","copy","atPole","x","y","sign","Math","z","fromArray","scale","set","normalize","geodeticSurfaceNormal","cross"],"mappings":"AAAA,SAAQA,OAAR,EAAiBC,MAAjB,EAAyBC,MAAM,IAAIC,aAAnC,QAAuD,eAAvD;AAIA,MAAMC,SAAS,GAAG,KAAlB;AAEA,MAAMC,aAAa,GAAG,IAAIL,OAAJ,EAAtB;AAGA,MAAMM,0BAA0B,GAAG;AACjCC,EAAAA,EAAE,EAAE;AACFC,IAAAA,KAAK,EAAE,MADL;AAEFC,IAAAA,KAAK,EAAE,MAFL;AAGFC,IAAAA,IAAI,EAAE,OAHJ;AAIFC,IAAAA,IAAI,EAAE;AAJJ,GAD6B;AAOjCC,EAAAA,IAAI,EAAE;AACJJ,IAAAA,KAAK,EAAE,MADH;AAEJC,IAAAA,KAAK,EAAE,MAFH;AAGJC,IAAAA,IAAI,EAAE,OAHF;AAIJC,IAAAA,IAAI,EAAE;AAJF,GAP2B;AAajCH,EAAAA,KAAK,EAAE;AACLD,IAAAA,EAAE,EAAE,MADC;AAELK,IAAAA,IAAI,EAAE,MAFD;AAGLF,IAAAA,IAAI,EAAE,MAHD;AAILC,IAAAA,IAAI,EAAE;AAJD,GAb0B;AAmBjCF,EAAAA,KAAK,EAAE;AACLF,IAAAA,EAAE,EAAE,MADC;AAELK,IAAAA,IAAI,EAAE,MAFD;AAGLF,IAAAA,IAAI,EAAE,IAHD;AAILC,IAAAA,IAAI,EAAE;AAJD,GAnB0B;AAyBjCD,EAAAA,IAAI,EAAE;AACJH,IAAAA,EAAE,EAAE,OADA;AAEJK,IAAAA,IAAI,EAAE,OAFF;AAGJH,IAAAA,KAAK,EAAE,MAHH;AAIJD,IAAAA,KAAK,EAAE;AAJH,GAzB2B;AA+BjCG,EAAAA,IAAI,EAAE;AACJJ,IAAAA,EAAE,EAAE,OADA;AAEJK,IAAAA,IAAI,EAAE,OAFF;AAGJH,IAAAA,KAAK,EAAE,IAHH;AAIJD,IAAAA,KAAK,EAAE;AAJH;AA/B2B,CAAnC;AAuCA,MAAMK,4BAA4B,GAAG;AACnCJ,EAAAA,KAAK,EAAE,CAAC,CAAC,CAAF,EAAK,CAAL,EAAQ,CAAR,CAD4B;AAEnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAF6B;AAGnCJ,EAAAA,EAAE,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAH+B;AAInCC,EAAAA,KAAK,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAJ4B;AAKnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAC,CAAL,EAAQ,CAAR,CAL6B;AAMnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAC,CAAR;AAN6B,CAArC;AASA,MAAME,kBAAkB,GAAG;AACzBH,EAAAA,IAAI,EAAE,IAAIX,OAAJ,EADmB;AAEzBS,EAAAA,KAAK,EAAE,IAAIT,OAAJ,EAFkB;AAGzBO,EAAAA,EAAE,EAAE,IAAIP,OAAJ,EAHqB;AAIzBU,EAAAA,IAAI,EAAE,IAAIV,OAAJ,EAJmB;AAKzBQ,EAAAA,KAAK,EAAE,IAAIR,OAAJ,EALkB;AAMzBY,EAAAA,IAAI,EAAE,IAAIZ,OAAJ;AANmB,CAA3B;AASA,MAAMe,cAAc,GAAG,IAAIf,OAAJ,EAAvB;AACA,MAAMgB,cAAc,GAAG,IAAIhB,OAAJ,EAAvB;AACA,MAAMiB,cAAc,GAAG,IAAIjB,OAAJ,EAAvB;AAOA,eAAe,SAASkB,sBAAT,CACbC,SADa,EAEbC,SAFa,EAGbC,UAHa,EAIbC,SAJa,EAKbC,eALa,EAMbC,MANa,EAOH;AACV,QAAMC,iBAAiB,GACrBnB,0BAA0B,CAACc,SAAD,CAA1B,IACCd,0BAA0B,CAACc,SAAD,CAA1B,CAAsCC,UAAtC,CAFH;AAIApB,EAAAA,MAAM,CAACwB,iBAAiB,KAAK,CAACH,SAAD,IAAcA,SAAS,KAAKG,iBAAjC,CAAlB,CAAN;AAEA,MAAIC,eAAJ;AACA,MAAIC,gBAAJ;AACA,MAAIC,eAAJ;AAEA,QAAMC,MAAM,GAAGxB,aAAa,CAACyB,IAAd,CAAmBP,eAAnB,CAAf;AAGA,QAAMQ,MAAM,GAAG5B,aAAa,CAAC0B,MAAM,CAACG,CAAR,EAAW,GAAX,EAAgB5B,SAAhB,CAAb,IAA2CD,aAAa,CAAC0B,MAAM,CAACI,CAAR,EAAW,GAAX,EAAgB7B,SAAhB,CAAvE;;AAEA,MAAI2B,MAAJ,EAAY;AAEV,UAAMG,IAAI,GAAGC,IAAI,CAACD,IAAL,CAAUL,MAAM,CAACO,CAAjB,CAAb;AAEAV,IAAAA,eAAe,GAAGX,cAAc,CAACsB,SAAf,CAAyBxB,4BAA4B,CAACO,SAAD,CAArD,CAAlB;;AACA,QAAIA,SAAS,KAAK,MAAd,IAAwBA,SAAS,KAAK,MAA1C,EAAkD;AAChDM,MAAAA,eAAe,CAACY,KAAhB,CAAsBJ,IAAtB;AACD;;AAEDP,IAAAA,gBAAgB,GAAGX,cAAc,CAACqB,SAAf,CAAyBxB,4BAA4B,CAACQ,UAAD,CAArD,CAAnB;;AACA,QAAIA,UAAU,KAAK,MAAf,IAAyBA,UAAU,KAAK,MAA5C,EAAoD;AAClDM,MAAAA,gBAAgB,CAACW,KAAjB,CAAuBJ,IAAvB;AACD;;AAEDN,IAAAA,eAAe,GAAGX,cAAc,CAACoB,SAAf,CAAyBxB,4BAA4B,CAACS,SAAD,CAArD,CAAlB;;AACA,QAAIA,SAAS,KAAK,MAAd,IAAwBA,SAAS,KAAK,MAA1C,EAAkD;AAChDM,MAAAA,eAAe,CAACU,KAAhB,CAAsBJ,IAAtB;AACD;AACF,GAlBD,MAkBO;AAEL,UAAM;AAAC3B,MAAAA,EAAD;AAAKI,MAAAA,IAAL;AAAWF,MAAAA;AAAX,QAAoBK,kBAA1B;AAEAH,IAAAA,IAAI,CAAC4B,GAAL,CAAS,CAACV,MAAM,CAACI,CAAjB,EAAoBJ,MAAM,CAACG,CAA3B,EAA8B,GAA9B,EAAmCQ,SAAnC;AACArB,IAAAA,SAAS,CAACsB,qBAAV,CAAgCZ,MAAhC,EAAwCtB,EAAxC;AACAE,IAAAA,KAAK,CAACqB,IAAN,CAAWvB,EAAX,EAAemC,KAAf,CAAqB/B,IAArB;AAEA,UAAM;AAACC,MAAAA,IAAD;AAAOF,MAAAA,IAAP;AAAaF,MAAAA;AAAb,QAAsBM,kBAA5B;AAEAF,IAAAA,IAAI,CAACkB,IAAL,CAAUvB,EAAV,EAAc+B,KAAd,CAAoB,CAAC,CAArB;AACA5B,IAAAA,IAAI,CAACoB,IAAL,CAAUnB,IAAV,EAAgB2B,KAAhB,CAAsB,CAAC,CAAvB;AACA9B,IAAAA,KAAK,CAACsB,IAAN,CAAWrB,KAAX,EAAkB6B,KAAlB,CAAwB,CAAC,CAAzB;AAGAZ,IAAAA,eAAe,GAAGZ,kBAAkB,CAACM,SAAD,CAApC;AACAO,IAAAA,gBAAgB,GAAGb,kBAAkB,CAACO,UAAD,CAArC;AACAO,IAAAA,eAAe,GAAGd,kBAAkB,CAACQ,SAAD,CAApC;AACD;;AAGDE,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACM,CAA5B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACO,CAA5B;AACAT,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACU,CAA5B;AACAZ,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAY,GAAZ;AACAA,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACK,CAA7B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACM,CAA7B;AACAT,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACS,CAA7B;AACAZ,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAY,GAAZ;AACAA,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYI,eAAe,CAACI,CAA5B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYI,eAAe,CAACK,CAA5B;AACAT,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaI,eAAe,CAACQ,CAA7B;AACAZ,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAa,GAAb;AACAA,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACG,CAApB;AACAR,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACI,CAApB;AACAT,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACO,CAApB;AACAZ,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAa,GAAb;AACA,SAAOA,MAAP;AACD","sourcesContent":["import {Vector3, assert, equals as equalsEpsilon} from '@math.gl/core';\n\nimport type Ellipsoid from '../ellipsoid';\n\nconst EPSILON14 = 1e-14;\n\nconst scratchOrigin = new Vector3();\n\n// Caclulate third axis from given two axii\nconst VECTOR_PRODUCT_LOCAL_FRAME = {\n  up: {\n    south: 'east',\n    north: 'west',\n    west: 'south',\n    east: 'north'\n  },\n  down: {\n    south: 'west',\n    north: 'east',\n    west: 'north',\n    east: 'south'\n  },\n  south: {\n    up: 'west',\n    down: 'east',\n    west: 'down',\n    east: 'up'\n  },\n  north: {\n    up: 'east',\n    down: 'west',\n    west: 'up',\n    east: 'down'\n  },\n  west: {\n    up: 'north',\n    down: 'south',\n    north: 'down',\n    south: 'up'\n  },\n  east: {\n    up: 'south',\n    down: 'north',\n    north: 'up',\n    south: 'down'\n  }\n} as const;\n\nconst degeneratePositionLocalFrame = {\n  north: [-1, 0, 0],\n  east: [0, 1, 0],\n  up: [0, 0, 1],\n  south: [1, 0, 0],\n  west: [0, -1, 0],\n  down: [0, 0, -1]\n} as const;\n\nconst scratchAxisVectors = {\n  east: new Vector3(),\n  north: new Vector3(),\n  up: new Vector3(),\n  west: new Vector3(),\n  south: new Vector3(),\n  down: new Vector3()\n};\n\nconst scratchVector1 = new Vector3();\nconst scratchVector2 = new Vector3();\nconst scratchVector3 = new Vector3();\n\ntype Axis = 'up' | 'down' | 'north' | 'east' | 'south' | 'west';\n\n// Computes a 4x4 transformation matrix from a reference frame\n// centered at the provided origin to the provided ellipsoid's fixed reference frame.\n// eslint-disable-next-line max-statements, max-params, complexity\nexport default function localFrameToFixedFrame(\n  ellipsoid: Ellipsoid,\n  firstAxis: Axis,\n  secondAxis: Axis,\n  thirdAxis: Axis,\n  cartesianOrigin: number[],\n  result: number[]\n): number[] {\n  const thirdAxisInferred =\n    VECTOR_PRODUCT_LOCAL_FRAME[firstAxis] &&\n    (VECTOR_PRODUCT_LOCAL_FRAME[firstAxis][secondAxis] as Axis);\n  // firstAxis and secondAxis must be east, north, up, west, south or down.');\n  assert(thirdAxisInferred && (!thirdAxis || thirdAxis === thirdAxisInferred));\n\n  let firstAxisVector: Vector3;\n  let secondAxisVector: Vector3;\n  let thirdAxisVector: Vector3;\n\n  const origin = scratchOrigin.copy(cartesianOrigin);\n\n  // If x and y are zero, assume origin is at a pole, which is a special case.\n  const atPole = equalsEpsilon(origin.x, 0.0, EPSILON14) && equalsEpsilon(origin.y, 0.0, EPSILON14);\n\n  if (atPole) {\n    // Look up axis value and adjust\n    const sign = Math.sign(origin.z);\n\n    firstAxisVector = scratchVector1.fromArray(degeneratePositionLocalFrame[firstAxis]);\n    if (firstAxis !== 'east' && firstAxis !== 'west') {\n      firstAxisVector.scale(sign);\n    }\n\n    secondAxisVector = scratchVector2.fromArray(degeneratePositionLocalFrame[secondAxis]);\n    if (secondAxis !== 'east' && secondAxis !== 'west') {\n      secondAxisVector.scale(sign);\n    }\n\n    thirdAxisVector = scratchVector3.fromArray(degeneratePositionLocalFrame[thirdAxis]);\n    if (thirdAxis !== 'east' && thirdAxis !== 'west') {\n      thirdAxisVector.scale(sign);\n    }\n  } else {\n    // Calculate all axis\n    const {up, east, north} = scratchAxisVectors;\n\n    east.set(-origin.y, origin.x, 0.0).normalize();\n    ellipsoid.geodeticSurfaceNormal(origin, up);\n    north.copy(up).cross(east);\n\n    const {down, west, south} = scratchAxisVectors;\n\n    down.copy(up).scale(-1);\n    west.copy(east).scale(-1);\n    south.copy(north).scale(-1);\n\n    // Pick three axis based on desired orientation\n    firstAxisVector = scratchAxisVectors[firstAxis];\n    secondAxisVector = scratchAxisVectors[secondAxis];\n    thirdAxisVector = scratchAxisVectors[thirdAxis];\n  }\n\n  // TODO - assuming the result is column-major\n  result[0] = firstAxisVector.x;\n  result[1] = firstAxisVector.y;\n  result[2] = firstAxisVector.z;\n  result[3] = 0.0;\n  result[4] = secondAxisVector.x;\n  result[5] = secondAxisVector.y;\n  result[6] = secondAxisVector.z;\n  result[7] = 0.0;\n  result[8] = thirdAxisVector.x;\n  result[9] = thirdAxisVector.y;\n  result[10] = thirdAxisVector.z;\n  result[11] = 0.0;\n  result[12] = origin.x;\n  result[13] = origin.y;\n  result[14] = origin.z;\n  result[15] = 1.0;\n  return result;\n}\n"],"file":"ellipsoid-transform.js"}

package/dist/esm/ellipsoid/helpers/scale-to-geodetic-surface.js

@@ -1,18 +1,17 @@
-import { Vector3, assert, _MathUtils } from '@math.gl/core';
-import * as vec3 from 'gl-matrix/vec3';
+import { Vector3, _MathUtils } from '@math.gl/core';
 const scratchVector = new Vector3();
 const scaleToGeodeticSurfaceIntersection = new Vector3();
 const scaleToGeodeticSurfaceGradient = new Vector3();
-export default function scaleToGeodeticSurface(cartesian, ellipsoid, result = new Vector3()) {
+export default function scaleToGeodeticSurface(cartesian, ellipsoid, result = []) {
   const {
     oneOverRadii,
     oneOverRadiiSquared,
     centerToleranceSquared
   } = ellipsoid;
   scratchVector.from(cartesian);
-  const positionX = cartesian.x;
-  const positionY = cartesian.y;
-  const positionZ = cartesian.z;
+  const positionX = scratchVector.x;
+  const positionY = scratchVector.y;
+  const positionZ = scratchVector.z;
   const oneOverRadiiX = oneOverRadii.x;
   const oneOverRadiiY = oneOverRadii.y;
   const oneOverRadiiZ = oneOverRadii.z;
@@ -38,7 +37,7 @@ export default function scaleToGeodeticSurface(cartesian, ellipsoid, result = ne
   const oneOverRadiiSquaredZ = oneOverRadiiSquared.z;
   const gradient = scaleToGeodeticSurfaceGradient;
   gradient.set(intersection.x * oneOverRadiiSquaredX * 2.0, intersection.y * oneOverRadiiSquaredY * 2.0, intersection.z * oneOverRadiiSquaredZ * 2.0);
-  let lambda = (1.0 - ratio) * cartesian.len() / (0.5 * gradient.len());
+  let lambda = (1.0 - ratio) * scratchVector.len() / (0.5 * gradient.len());
   let correction = 0.0;
   let xMultiplier;
   let yMultiplier;

package/dist/esm/ellipsoid/helpers/scale-to-geodetic-surface.js.map

@@ -1 +1 @@
-{"version":3,"sources":["../../../../src/ellipsoid/helpers/scale-to-geodetic-surface.js"],"names":["Vector3","assert","_MathUtils","vec3","scratchVector","scaleToGeodeticSurfaceIntersection","scaleToGeodeticSurfaceGradient","scaleToGeodeticSurface","cartesian","ellipsoid","result","oneOverRadii","oneOverRadiiSquared","centerToleranceSquared","from","positionX","x","positionY","y","positionZ","z","oneOverRadiiX","oneOverRadiiY","oneOverRadiiZ","x2","y2","z2","squaredNorm","ratio","Math","sqrt","Number","isFinite","undefined","intersection","copy","scale","to","oneOverRadiiSquaredX","oneOverRadiiSquaredY","oneOverRadiiSquaredZ","gradient","set","lambda","len","correction","xMultiplier","yMultiplier","zMultiplier","func","xMultiplier2","yMultiplier2","zMultiplier2","xMultiplier3","yMultiplier3","zMultiplier3","denominator","derivative","abs","EPSILON12"],"mappings":"AACA,SAAQA,OAAR,EAAiBC,MAAjB,EAAyBC,UAAzB,QAA0C,eAA1C;AACA,OAAO,KAAKC,IAAZ,MAAsB,gBAAtB;AAEA,MAAMC,aAAa,GAAG,IAAIJ,OAAJ,EAAtB;AACA,MAAMK,kCAAkC,GAAG,IAAIL,OAAJ,EAA3C;AACA,MAAMM,8BAA8B,GAAG,IAAIN,OAAJ,EAAvC;AAKA,eAAe,SAASO,sBAAT,CAAgCC,SAAhC,EAA2CC,SAA3C,EAAsDC,MAAM,GAAG,IAAIV,OAAJ,EAA/D,EAA8E;AAC3F,QAAM;AAACW,IAAAA,YAAD;AAAeC,IAAAA,mBAAf;AAAoCC,IAAAA;AAApC,MAA8DJ,SAApE;AAEAL,EAAAA,aAAa,CAACU,IAAd,CAAmBN,SAAnB;AAEA,QAAMO,SAAS,GAAGP,SAAS,CAACQ,CAA5B;AACA,QAAMC,SAAS,GAAGT,SAAS,CAACU,CAA5B;AACA,QAAMC,SAAS,GAAGX,SAAS,CAACY,CAA5B;AAEA,QAAMC,aAAa,GAAGV,YAAY,CAACK,CAAnC;AACA,QAAMM,aAAa,GAAGX,YAAY,CAACO,CAAnC;AACA,QAAMK,aAAa,GAAGZ,YAAY,CAACS,CAAnC;AAEA,QAAMI,EAAE,GAAGT,SAAS,GAAGA,SAAZ,GAAwBM,aAAxB,GAAwCA,aAAnD;AACA,QAAMI,EAAE,GAAGR,SAAS,GAAGA,SAAZ,GAAwBK,aAAxB,GAAwCA,aAAnD;AACA,QAAMI,EAAE,GAAGP,SAAS,GAAGA,SAAZ,GAAwBI,aAAxB,GAAwCA,aAAnD;AAGA,QAAMI,WAAW,GAAGH,EAAE,GAAGC,EAAL,GAAUC,EAA9B;AACA,QAAME,KAAK,GAAGC,IAAI,CAACC,IAAL,CAAU,MAAMH,WAAhB,CAAd;;AAGA,MAAI,CAACI,MAAM,CAACC,QAAP,CAAgBJ,KAAhB,CAAL,EAA6B;AAC3B,WAAOK,SAAP;AACD;;AAGD,QAAMC,YAAY,GAAG7B,kCAArB;AACA6B,EAAAA,YAAY,CAACC,IAAb,CAAkB3B,SAAlB,EAA6B4B,KAA7B,CAAmCR,KAAnC;;AAGA,MAAID,WAAW,GAAGd,sBAAlB,EAA0C;AACxC,WAAOqB,YAAY,CAACG,EAAb,CAAgB3B,MAAhB,CAAP;AACD;;AAED,QAAM4B,oBAAoB,GAAG1B,mBAAmB,CAACI,CAAjD;AACA,QAAMuB,oBAAoB,GAAG3B,mBAAmB,CAACM,CAAjD;AACA,QAAMsB,oBAAoB,GAAG5B,mBAAmB,CAACQ,CAAjD;AAIA,QAAMqB,QAAQ,GAAGnC,8BAAjB;AACAmC,EAAAA,QAAQ,CAACC,GAAT,CACER,YAAY,CAAClB,CAAb,GAAiBsB,oBAAjB,GAAwC,GAD1C,EAEEJ,YAAY,CAAChB,CAAb,GAAiBqB,oBAAjB,GAAwC,GAF1C,EAGEL,YAAY,CAACd,CAAb,GAAiBoB,oBAAjB,GAAwC,GAH1C;AAOA,MAAIG,MAAM,GAAI,CAAC,MAAMf,KAAP,IAAgBpB,SAAS,CAACoC,GAAV,EAAjB,IAAqC,MAAMH,QAAQ,CAACG,GAAT,EAA3C,CAAb;AACA,MAAIC,UAAU,GAAG,GAAjB;AAEA,MAAIC,WAAJ;AACA,MAAIC,WAAJ;AACA,MAAIC,WAAJ;AACA,MAAIC,IAAJ;;AAEA,KAAG;AACDN,IAAAA,MAAM,IAAIE,UAAV;AAEAC,IAAAA,WAAW,GAAG,OAAO,MAAMH,MAAM,GAAGL,oBAAtB,CAAd;AACAS,IAAAA,WAAW,GAAG,OAAO,MAAMJ,MAAM,GAAGJ,oBAAtB,CAAd;AACAS,IAAAA,WAAW,GAAG,OAAO,MAAML,MAAM,GAAGH,oBAAtB,CAAd;AAEA,UAAMU,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AACA,UAAMK,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AACA,UAAMK,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AAEA,UAAMK,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AACA,UAAMQ,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AACA,UAAMQ,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AAEAC,IAAAA,IAAI,GAAGzB,EAAE,GAAG0B,YAAL,GAAoBzB,EAAE,GAAG0B,YAAzB,GAAwCzB,EAAE,GAAG0B,YAA7C,GAA4D,GAAnE;AAIA,UAAMI,WAAW,GACfhC,EAAE,GAAG6B,YAAL,GAAoBf,oBAApB,GACAb,EAAE,GAAG6B,YAAL,GAAoBf,oBADpB,GAEAb,EAAE,GAAG6B,YAAL,GAAoBf,oBAHtB;AAKA,UAAMiB,UAAU,GAAG,CAAC,GAAD,GAAOD,WAA1B;AAEAX,IAAAA,UAAU,GAAGI,IAAI,GAAGQ,UAApB;AACD,GA3BD,QA2BS5B,IAAI,CAAC6B,GAAL,CAAST,IAAT,IAAiB/C,UAAU,CAACyD,SA3BrC;;AA6BA,SAAOvD,aAAa,CAACgC,KAAd,CAAoB,CAACU,WAAD,EAAcC,WAAd,EAA2BC,WAA3B,CAApB,EAA6DX,EAA7D,CAAgE3B,MAAhE,CAAP;AACD","sourcesContent":["/* eslint-disable */\nimport {Vector3, assert, _MathUtils} from '@math.gl/core';\nimport * as vec3 from 'gl-matrix/vec3';\n\nconst scratchVector = new Vector3();\nconst scaleToGeodeticSurfaceIntersection = new Vector3();\nconst scaleToGeodeticSurfaceGradient = new Vector3();\n\n// Scales the provided Cartesian position along the geodetic surface normal\n// so that it is on the surface of this ellipsoid.  If the position is\n// at the center of the ellipsoid, this function returns undefined.\nexport default function scaleToGeodeticSurface(cartesian, ellipsoid, result = new Vector3()) {\n  const {oneOverRadii, oneOverRadiiSquared, centerToleranceSquared} = ellipsoid;\n\n  scratchVector.from(cartesian);\n\n  const positionX = cartesian.x;\n  const positionY = cartesian.y;\n  const positionZ = cartesian.z;\n\n  const oneOverRadiiX = oneOverRadii.x;\n  const oneOverRadiiY = oneOverRadii.y;\n  const oneOverRadiiZ = oneOverRadii.z;\n\n  const x2 = positionX * positionX * oneOverRadiiX * oneOverRadiiX;\n  const y2 = positionY * positionY * oneOverRadiiY * oneOverRadiiY;\n  const z2 = positionZ * positionZ * oneOverRadiiZ * oneOverRadiiZ;\n\n  // Compute the squared ellipsoid norm.\n  const squaredNorm = x2 + y2 + z2;\n  const ratio = Math.sqrt(1.0 / squaredNorm);\n\n  // When very close to center or at center\n  if (!Number.isFinite(ratio)) {\n    return undefined;\n  }\n\n  // As an initial approximation, assume that the radial intersection is the projection point.\n  const intersection = scaleToGeodeticSurfaceIntersection;\n  intersection.copy(cartesian).scale(ratio);\n\n  // If the position is near the center, the iteration will not converge.\n  if (squaredNorm < centerToleranceSquared) {\n    return intersection.to(result);\n  }\n\n  const oneOverRadiiSquaredX = oneOverRadiiSquared.x;\n  const oneOverRadiiSquaredY = oneOverRadiiSquared.y;\n  const oneOverRadiiSquaredZ = oneOverRadiiSquared.z;\n\n  // Use the gradient at the intersection point in place of the true unit normal.\n  // The difference in magnitude will be absorbed in the multiplier.\n  const gradient = scaleToGeodeticSurfaceGradient;\n  gradient.set(\n    intersection.x * oneOverRadiiSquaredX * 2.0,\n    intersection.y * oneOverRadiiSquaredY * 2.0,\n    intersection.z * oneOverRadiiSquaredZ * 2.0\n  );\n\n  // Compute the initial guess at the normal vector multiplier, lambda.\n  let lambda = ((1.0 - ratio) * cartesian.len()) / (0.5 * gradient.len());\n  let correction = 0.0;\n\n  let xMultiplier;\n  let yMultiplier;\n  let zMultiplier;\n  let func;\n\n  do {\n    lambda -= correction;\n\n    xMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredX);\n    yMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredY);\n    zMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredZ);\n\n    const xMultiplier2 = xMultiplier * xMultiplier;\n    const yMultiplier2 = yMultiplier * yMultiplier;\n    const zMultiplier2 = zMultiplier * zMultiplier;\n\n    const xMultiplier3 = xMultiplier2 * xMultiplier;\n    const yMultiplier3 = yMultiplier2 * yMultiplier;\n    const zMultiplier3 = zMultiplier2 * zMultiplier;\n\n    func = x2 * xMultiplier2 + y2 * yMultiplier2 + z2 * zMultiplier2 - 1.0;\n\n    // \"denominator\" here refers to the use of this expression in the velocity and acceleration\n    // computations in the sections to follow.\n    const denominator =\n      x2 * xMultiplier3 * oneOverRadiiSquaredX +\n      y2 * yMultiplier3 * oneOverRadiiSquaredY +\n      z2 * zMultiplier3 * oneOverRadiiSquaredZ;\n\n    const derivative = -2.0 * denominator;\n\n    correction = func / derivative;\n  } while (Math.abs(func) > _MathUtils.EPSILON12);\n\n  return scratchVector.scale([xMultiplier, yMultiplier, zMultiplier]).to(result);\n}\n"],"file":"scale-to-geodetic-surface.js"}
+{"version":3,"sources":["../../../../src/ellipsoid/helpers/scale-to-geodetic-surface.ts"],"names":["Vector3","_MathUtils","scratchVector","scaleToGeodeticSurfaceIntersection","scaleToGeodeticSurfaceGradient","scaleToGeodeticSurface","cartesian","ellipsoid","result","oneOverRadii","oneOverRadiiSquared","centerToleranceSquared","from","positionX","x","positionY","y","positionZ","z","oneOverRadiiX","oneOverRadiiY","oneOverRadiiZ","x2","y2","z2","squaredNorm","ratio","Math","sqrt","Number","isFinite","undefined","intersection","copy","scale","to","oneOverRadiiSquaredX","oneOverRadiiSquaredY","oneOverRadiiSquaredZ","gradient","set","lambda","len","correction","xMultiplier","yMultiplier","zMultiplier","func","xMultiplier2","yMultiplier2","zMultiplier2","xMultiplier3","yMultiplier3","zMultiplier3","denominator","derivative","abs","EPSILON12"],"mappings":"AACA,SAAQA,OAAR,EAAiBC,UAAjB,QAAkC,eAAlC;AAGA,MAAMC,aAAa,GAAG,IAAIF,OAAJ,EAAtB;AACA,MAAMG,kCAAkC,GAAG,IAAIH,OAAJ,EAA3C;AACA,MAAMI,8BAA8B,GAAG,IAAIJ,OAAJ,EAAvC;AAKA,eAAe,SAASK,sBAAT,CACbC,SADa,EAEbC,SAFa,EAGbC,MAAgB,GAAG,EAHN,EAIH;AACV,QAAM;AAACC,IAAAA,YAAD;AAAeC,IAAAA,mBAAf;AAAoCC,IAAAA;AAApC,MAA8DJ,SAApE;AAEAL,EAAAA,aAAa,CAACU,IAAd,CAAmBN,SAAnB;AAEA,QAAMO,SAAS,GAAGX,aAAa,CAACY,CAAhC;AACA,QAAMC,SAAS,GAAGb,aAAa,CAACc,CAAhC;AACA,QAAMC,SAAS,GAAGf,aAAa,CAACgB,CAAhC;AAEA,QAAMC,aAAa,GAAGV,YAAY,CAACK,CAAnC;AACA,QAAMM,aAAa,GAAGX,YAAY,CAACO,CAAnC;AACA,QAAMK,aAAa,GAAGZ,YAAY,CAACS,CAAnC;AAEA,QAAMI,EAAE,GAAGT,SAAS,GAAGA,SAAZ,GAAwBM,aAAxB,GAAwCA,aAAnD;AACA,QAAMI,EAAE,GAAGR,SAAS,GAAGA,SAAZ,GAAwBK,aAAxB,GAAwCA,aAAnD;AACA,QAAMI,EAAE,GAAGP,SAAS,GAAGA,SAAZ,GAAwBI,aAAxB,GAAwCA,aAAnD;AAGA,QAAMI,WAAW,GAAGH,EAAE,GAAGC,EAAL,GAAUC,EAA9B;AACA,QAAME,KAAK,GAAGC,IAAI,CAACC,IAAL,CAAU,MAAMH,WAAhB,CAAd;;AAGA,MAAI,CAACI,MAAM,CAACC,QAAP,CAAgBJ,KAAhB,CAAL,EAA6B;AAC3B,WAAOK,SAAP;AACD;;AAGD,QAAMC,YAAY,GAAG7B,kCAArB;AACA6B,EAAAA,YAAY,CAACC,IAAb,CAAkB3B,SAAlB,EAA6B4B,KAA7B,CAAmCR,KAAnC;;AAGA,MAAID,WAAW,GAAGd,sBAAlB,EAA0C;AACxC,WAAOqB,YAAY,CAACG,EAAb,CAAgB3B,MAAhB,CAAP;AACD;;AAED,QAAM4B,oBAAoB,GAAG1B,mBAAmB,CAACI,CAAjD;AACA,QAAMuB,oBAAoB,GAAG3B,mBAAmB,CAACM,CAAjD;AACA,QAAMsB,oBAAoB,GAAG5B,mBAAmB,CAACQ,CAAjD;AAIA,QAAMqB,QAAQ,GAAGnC,8BAAjB;AACAmC,EAAAA,QAAQ,CAACC,GAAT,CACER,YAAY,CAAClB,CAAb,GAAiBsB,oBAAjB,GAAwC,GAD1C,EAEEJ,YAAY,CAAChB,CAAb,GAAiBqB,oBAAjB,GAAwC,GAF1C,EAGEL,YAAY,CAACd,CAAb,GAAiBoB,oBAAjB,GAAwC,GAH1C;AAOA,MAAIG,MAAM,GAAI,CAAC,MAAMf,KAAP,IAAgBxB,aAAa,CAACwC,GAAd,EAAjB,IAAyC,MAAMH,QAAQ,CAACG,GAAT,EAA/C,CAAb;AACA,MAAIC,UAAU,GAAG,GAAjB;AAEA,MAAIC,WAAJ;AACA,MAAIC,WAAJ;AACA,MAAIC,WAAJ;AACA,MAAIC,IAAJ;;AAEA,KAAG;AACDN,IAAAA,MAAM,IAAIE,UAAV;AAEAC,IAAAA,WAAW,GAAG,OAAO,MAAMH,MAAM,GAAGL,oBAAtB,CAAd;AACAS,IAAAA,WAAW,GAAG,OAAO,MAAMJ,MAAM,GAAGJ,oBAAtB,CAAd;AACAS,IAAAA,WAAW,GAAG,OAAO,MAAML,MAAM,GAAGH,oBAAtB,CAAd;AAEA,UAAMU,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AACA,UAAMK,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AACA,UAAMK,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AAEA,UAAMK,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AACA,UAAMQ,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AACA,UAAMQ,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AAEAC,IAAAA,IAAI,GAAGzB,EAAE,GAAG0B,YAAL,GAAoBzB,EAAE,GAAG0B,YAAzB,GAAwCzB,EAAE,GAAG0B,YAA7C,GAA4D,GAAnE;AAIA,UAAMI,WAAW,GACfhC,EAAE,GAAG6B,YAAL,GAAoBf,oBAApB,GACAb,EAAE,GAAG6B,YAAL,GAAoBf,oBADpB,GAEAb,EAAE,GAAG6B,YAAL,GAAoBf,oBAHtB;AAKA,UAAMiB,UAAU,GAAG,CAAC,GAAD,GAAOD,WAA1B;AAEAX,IAAAA,UAAU,GAAGI,IAAI,GAAGQ,UAApB;AACD,GA3BD,QA2BS5B,IAAI,CAAC6B,GAAL,CAAST,IAAT,IAAiB9C,UAAU,CAACwD,SA3BrC;;AA6BA,SAAOvD,aAAa,CAACgC,KAAd,CAAoB,CAACU,WAAD,EAAcC,WAAd,EAA2BC,WAA3B,CAApB,EAA6DX,EAA7D,CAAgE3B,MAAhE,CAAP;AACD","sourcesContent":["/* eslint-disable */\nimport {Vector3, _MathUtils} from '@math.gl/core';\nimport type Ellipsoid from '../ellipsoid';\n\nconst scratchVector = new Vector3();\nconst scaleToGeodeticSurfaceIntersection = new Vector3();\nconst scaleToGeodeticSurfaceGradient = new Vector3();\n\n// Scales the provided Cartesian position along the geodetic surface normal\n// so that it is on the surface of this ellipsoid.  If the position is\n// at the center of the ellipsoid, this function returns undefined.\nexport default function scaleToGeodeticSurface(\n  cartesian: number[],\n  ellipsoid: Ellipsoid,\n  result: number[] = []\n): number[] {\n  const {oneOverRadii, oneOverRadiiSquared, centerToleranceSquared} = ellipsoid;\n\n  scratchVector.from(cartesian);\n\n  const positionX = scratchVector.x;\n  const positionY = scratchVector.y;\n  const positionZ = scratchVector.z;\n\n  const oneOverRadiiX = oneOverRadii.x;\n  const oneOverRadiiY = oneOverRadii.y;\n  const oneOverRadiiZ = oneOverRadii.z;\n\n  const x2 = positionX * positionX * oneOverRadiiX * oneOverRadiiX;\n  const y2 = positionY * positionY * oneOverRadiiY * oneOverRadiiY;\n  const z2 = positionZ * positionZ * oneOverRadiiZ * oneOverRadiiZ;\n\n  // Compute the squared ellipsoid norm.\n  const squaredNorm = x2 + y2 + z2;\n  const ratio = Math.sqrt(1.0 / squaredNorm);\n\n  // When very close to center or at center\n  if (!Number.isFinite(ratio)) {\n    return undefined;\n  }\n\n  // As an initial approximation, assume that the radial intersection is the projection point.\n  const intersection = scaleToGeodeticSurfaceIntersection;\n  intersection.copy(cartesian).scale(ratio);\n\n  // If the position is near the center, the iteration will not converge.\n  if (squaredNorm < centerToleranceSquared) {\n    return intersection.to(result);\n  }\n\n  const oneOverRadiiSquaredX = oneOverRadiiSquared.x;\n  const oneOverRadiiSquaredY = oneOverRadiiSquared.y;\n  const oneOverRadiiSquaredZ = oneOverRadiiSquared.z;\n\n  // Use the gradient at the intersection point in place of the true unit normal.\n  // The difference in magnitude will be absorbed in the multiplier.\n  const gradient = scaleToGeodeticSurfaceGradient;\n  gradient.set(\n    intersection.x * oneOverRadiiSquaredX * 2.0,\n    intersection.y * oneOverRadiiSquaredY * 2.0,\n    intersection.z * oneOverRadiiSquaredZ * 2.0\n  );\n\n  // Compute the initial guess at the normal vector multiplier, lambda.\n  let lambda = ((1.0 - ratio) * scratchVector.len()) / (0.5 * gradient.len());\n  let correction = 0.0;\n\n  let xMultiplier;\n  let yMultiplier;\n  let zMultiplier;\n  let func;\n\n  do {\n    lambda -= correction;\n\n    xMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredX);\n    yMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredY);\n    zMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredZ);\n\n    const xMultiplier2 = xMultiplier * xMultiplier;\n    const yMultiplier2 = yMultiplier * yMultiplier;\n    const zMultiplier2 = zMultiplier * zMultiplier;\n\n    const xMultiplier3 = xMultiplier2 * xMultiplier;\n    const yMultiplier3 = yMultiplier2 * yMultiplier;\n    const zMultiplier3 = zMultiplier2 * zMultiplier;\n\n    func = x2 * xMultiplier2 + y2 * yMultiplier2 + z2 * zMultiplier2 - 1.0;\n\n    // \"denominator\" here refers to the use of this expression in the velocity and acceleration\n    // computations in the sections to follow.\n    const denominator =\n      x2 * xMultiplier3 * oneOverRadiiSquaredX +\n      y2 * yMultiplier3 * oneOverRadiiSquaredY +\n      z2 * zMultiplier3 * oneOverRadiiSquaredZ;\n\n    const derivative = -2.0 * denominator;\n\n    correction = func / derivative;\n  } while (Math.abs(func) > _MathUtils.EPSILON12);\n\n  return scratchVector.scale([xMultiplier, yMultiplier, zMultiplier]).to(result);\n}\n"],"file":"scale-to-geodetic-surface.js"}

package/dist/esm/index.js.map

@@ -1 +1 @@
-{"version":3,"sources":["../../src/index.js"],"names":["default","Ellipsoid","isWGS84"],"mappings":"AAAA,SAAQA,OAAO,IAAIC,SAAnB,QAAmC,uBAAnC;AACA,SAAQC,OAAR,QAAsB,cAAtB","sourcesContent":["export {default as Ellipsoid} from './ellipsoid/ellipsoid';\nexport {isWGS84} from './type-utils';\n"],"file":"index.js"}
+{"version":3,"sources":["../../src/index.ts"],"names":["default","Ellipsoid","isWGS84"],"mappings":"AAAA,SAAQA,OAAO,IAAIC,SAAnB,QAAmC,uBAAnC;AACA,SAAQC,OAAR,QAAsB,cAAtB","sourcesContent":["export {default as Ellipsoid} from './ellipsoid/ellipsoid';\nexport {isWGS84} from './type-utils';\n"],"file":"index.js"}