@itwin/core-geometry 4.0.0-dev.70 → 4.0.0-dev.73

package/lib/cjs/geometry3d/Ray3d.js.map

@@ -1 +1 @@
-{"version":3,"file":"Ray3d.js","sourceRoot":"","sources":["../../../src/geometry3d/Ray3d.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AACH,sEAAoH;AACpH,0CAAmE;AACnE,yDAAsD;AACtD,yCAAsC;AAEtC,uDAA6C;AAC7C,uDAAsD;AACtD,mCAA2C;AAC3C,2CAAwC;AAGxC;;;;;;GAMG;AACH,MAAa,KAAK;IAOhB,oCAAoC;IACpC,YAAoB,MAAe,EAAE,SAAmB;QACtD,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;QACrB,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC;QAC3B,IAAI,CAAC,CAAC,GAAG,SAAS,CAAC;IACrB,CAAC;IACO,MAAM,CAAC,OAAO,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS;QACrF,OAAO,IAAI,KAAK,CAAC,yBAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,0BAAQ,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACtE,CAAC;IACD,kCAAkC;IAC3B,MAAM,CAAC,WAAW;QACvB,OAAO,KAAK,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACzC,CAAC;IACD,kCAAkC;IAC3B,MAAM,CAAC,WAAW;QACvB,OAAO,KAAK,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACzC,CAAC;IACD,kCAAkC;IAC3B,MAAM,CAAC,WAAW;QACvB,OAAO,KAAK,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACzC,CAAC;IACD,mCAAmC;IAC5B,MAAM,CAAC,UAAU,CAAC,MAAc;QACrC,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,MAAM,CAAC,OAAO,EAAE,CAAC;YACxB,MAAM,CAAC,SAAS,CAAC,OAAO,EAAE,CAAC;YAC3B,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,KAAK,CAAC,yBAAO,CAAC,UAAU,EAAE,EAAE,0BAAQ,CAAC,UAAU,EAAE,CAAC,CAAC;IAChE,CAAC;IACD;;;;;OAKG;IACI,aAAa,CAAC,KAAY;QAC/B,OAAO,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,MAAM,CAAC,IAAI,IAAI,CAAC,SAAS,CAAC,aAAa,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC;IAClG,CAAC;IACD,kHAAkH;IAC3G,iBAAiB,CAAC,UAAmB;QAC1C,OAAO,IAAI,CAAC,SAAS,CAAC,kBAAkB,CAAC,IAAI,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;IACpE,CAAC;IACD,8GAA8G;IACvG,eAAe,CAAC,UAAmB;QACxC,OAAO,mBAAQ,CAAC,kBAAkB,CAChC,IAAI,CAAC,iBAAiB,CAAC,UAAU,CAAC,EAClC,IAAI,CAAC,SAAS,CAAC,gBAAgB,EAAE,EACjC,CAAC,CACF,CAAC;IACJ,CAAC;IACD,sDAAsD;IAC/C,iBAAiB,CAAC,UAAmB;QAC1C;;;;;;WAMG;QACH,OAAO,IAAI,CAAC,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,SAAS,EAAE,IAAI,CAAC,eAAe,CAAC,UAAU,CAAC,CAAC,CAAC;IAClF,CAAC;IACD;;;OAGG;IACI,qBAAqB,CAAC,KAAY;QACvC;;;;WAIG;QACH,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,YAAY,CAAC,KAAK,CAAC,SAAS,EAAE,IAAI,CAAC;YACrD,OAAO,KAAK,CAAC;QACf;;;;WAIG;QACH,IAAI,SAAS,GAAG,IAAI,CAAC,iBAAiB,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QACrD,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,mBAAmB,CAAC,SAAS,CAAC;YAC9C,OAAO,KAAK,CAAC;QACf,SAAS,GAAG,KAAK,CAAC,iBAAiB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACjD,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,mBAAmB,CAAC,SAAS,CAAC;YAC7C,OAAO,KAAK,CAAC;QACf,OAAO,IAAI,CAAC;IACd,CAAC;IACD,8CAA8C;IACvC,MAAM,CAAC,MAAM,CAAC,MAAe,EAAE,SAAmB,EAAE,MAAc;QACvE,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,GAAG,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;YAC9B,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,KAAK,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC;IACtD,CAAC;IACD;;;;;;OAMG;IACI,MAAM,CAAC,wBAAwB,CACpC,aAA2B,EAAE,kBAAgC,EAAE,MAAc;QAE7E,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,EAAE,GAAG,kBAAkB,CAAC,CAAC,CAAC,CAAC;QACjC,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,EAAE,GAAG,kBAAkB,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAC7D,MAAM,EAAE,GAAG,kBAAkB,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAC7D,MAAM,EAAE,GAAG,kBAAkB,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAC7D,IAAI,mBAAQ,CAAC,qBAAqB,CAAC,CAAC,CAAC;YACnC,OAAO,SAAS,CAAC;QACnB,MAAM,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC;QACrB,MAAM,KAAK,GAAG,IAAI,GAAG,IAAI,CAAC;QAC1B,OAAO,KAAK,CAAC,YAAY,CACvB,CAAC,GAAG,IAAI,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,GAAG,IAAI,EAC5B,EAAE,GAAG,KAAK,EAAE,EAAE,GAAG,KAAK,EAAE,EAAE,GAAG,KAAK,EAClC,MAAM,CACP,CAAC;IACJ,CAAC;IACD,2DAA2D;IACpD,MAAM,CAAC,YAAY,CACxB,OAAe,EAAE,OAAe,EAAE,OAAe,EACjD,UAAkB,EAAE,UAAkB,EAAE,UAAkB,EAC1D,MAAc;QAEd,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,YAAY,EAAE,CAAC,GAAG,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;YACrD,MAAM,CAAC,eAAe,EAAE,CAAC,GAAG,CAAC,UAAU,EAAE,UAAU,EAAE,UAAU,CAAC,CAAC;YACjE,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,KAAK,CAAC,yBAAO,CAAC,MAAM,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,EAAE,0BAAQ,CAAC,MAAM,CAAC,UAAU,EAAE,UAAU,EAAE,UAAU,CAAC,CAAC,CAAC;IACnH,CAAC;IACD,mDAAmD;IAC5C,MAAM,CAAC,aAAa,CAAC,MAAe,EAAE,SAAmB;QAC9D,OAAO,IAAI,KAAK,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACtC,CAAC;IACD,qEAAqE;IAC9D,MAAM,CAAC,uBAAuB,CAAC,MAAe,EAAE,SAAmB,EAAE,CAAS,EAAE,MAAc;QACnG,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;YAC9B,MAAM,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;YACpC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;YACb,OAAO,MAAM,CAAC;SACf;QACD,MAAM,GAAG,IAAI,KAAK,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC;QACtD,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;QACb,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,sHAAsH;IAC/G,MAAM,CAAC,cAAc,CAAC,MAAe,EAAE,MAAe,EAAE,MAAc;QAC3E,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;YAC9B,MAAM,CAAC,SAAS,CAAC,WAAW,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;YAC7C,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,KAAK,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,0BAAQ,CAAC,cAAc,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,CAAC;IAC5E,CAAC;IACD,8CAA8C;IACvC,YAAY;QACjB,OAAO,IAAI,CAAC,MAAM,CAAC;IACrB,CAAC;IACD,wDAAwD;IACjD,eAAe;QACpB,OAAO,IAAI,CAAC,SAAS,CAAC;IACxB,CAAC;IACD,kDAAkD;IAC3C,GAAG,CAAC,MAAe,EAAE,SAAmB;QAC7C,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;QAC5B,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;IACpC,CAAC;IACD,qBAAqB;IACd,KAAK,CAAC,MAAc;QACzB,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC;YACxD,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC;IAChE,CAAC;IACD,4DAA4D;IACrD,gBAAgB,CAAC,SAAoB;QAC1C,OAAO,IAAI,KAAK,CAAC,SAAS,CAAC,eAAe,CAAC,IAAI,CAAC,MAAM,CAAC,EAAE,SAAS,CAAC,cAAc,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC;IACrG,CAAC;IACD,oEAAoE;IAC7D,uBAAuB,CAAC,SAAoB;QACjD,MAAM,MAAM,GAAG,SAAS,CAAC,sBAAsB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QAC7D,MAAM,SAAS,GAAG,SAAS,CAAC,MAAM,CAAC,4BAA4B,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QACtH,IAAI,SAAS,KAAK,MAAM,IAAI,SAAS,KAAK,SAAS;YACjD,OAAO,IAAI,KAAK,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;QACtC,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,kCAAkC;IAC3B,gBAAgB,CAAC,SAAoB;QAC1C,SAAS,CAAC,eAAe,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QACpD,SAAS,CAAC,cAAc,CAAC,IAAI,CAAC,SAAS,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;IAC3D,CAAC;IACD,kCAAkC;IAC3B,OAAO,CAAC,MAAa;QAC1B,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,SAAS,CAAC,CAAC;IAC5C,CAAC;IACD;;;;OAIG;IACI,eAAe,CAAC,QAAgB,EAAE,MAAgB;QACvD,OAAO,IAAI,CAAC,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,SAAS,EAAE,QAAQ,EAAE,MAAM,CAAC,CAAC;IAClE,CAAC;IACD;;;OAGG;IACI,aAAa;QAClB,MAAM,IAAI,GAAG,mBAAQ,CAAC,kBAAkB,CAAC,IAAI,CAAC,SAAS,EAAE,oBAAS,CAAC,GAAG,CAAC,CAAC;QACxE,OAAO,qBAAS,CAAC,qBAAqB,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;IAC5D,CAAC;IACD,8DAA8D;IACvD,WAAW,CAAC,IAAU;QAC3B,IAAI,CAAC,IAAI,EAAE;YACT,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACzB,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAC5B,OAAO;SACR;QACD,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACrC,IAAI,CAAC,SAAS,CAAC,WAAW,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;IAC7C,CAAC;IACD;;;OAGG;IACI,+BAA+B,CAAC,YAAoB,GAAG;QAC5D,IAAI,IAAI,CAAC,SAAS,CAAC,mBAAmB,EAAE,EAAE;YACxC,IAAI,CAAC,SAAS,CAAC,YAAY,CAAC,SAAS,CAAC,CAAC;YACvC,OAAO,IAAI,CAAC;SACb;QACD,IAAI,CAAC,SAAS,CAAC,OAAO,EAAE,CAAC;QACzB,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC;QACb,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;;;;;;OAQG;IACH,6CAA6C;IAC7C,0GAA0G;IAC1G,+EAA+E;IACxE,iCAAiC,CAAC,CAAS;QAChD,MAAM,SAAS,GAAG,mBAAQ,CAAC,0BAA0B,CAAC;QACtD,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;QACX,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,SAAS,IAAI,IAAI,CAAC,SAAS,CAAC,mBAAmB,CAAC,SAAS,CAAC;YAC1E,OAAO,IAAI,CAAC;QACd,IAAI,CAAC,SAAS,CAAC,OAAO,EAAE,CAAC;QACzB,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC;QACb,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;OAGG;IACI,MAAM;QACX,OAAO,EAAE,MAAM,EAAE,IAAI,CAAC,MAAM,CAAC,MAAM,EAAE,EAAE,SAAS,EAAE,IAAI,CAAC,SAAS,CAAC,MAAM,EAAE,EAAE,CAAC;IAC9E,CAAC;IACD,+EAA+E;IACxE,MAAM,CAAC,QAAQ,CAAC,IAAU;QAC/B,MAAM,MAAM,GAAG,KAAK,CAAC,WAAW,EAAE,CAAC;QACnC,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,CAAC;QACzB,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,qDAAqD;IAC9C,QAAQ,CAAC,UAAmB;QACjC,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,CAAC,gBAAgB,EAAE,CAAC;QAC7C,MAAM,EAAE,GAAG,IAAI,CAAC,iBAAiB,CAAC,UAAU,CAAC,CAAC;QAC9C,MAAM,EAAE,GAAG,mBAAQ,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;QACrD,IAAI,EAAE;YACJ,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,eAAe,CAAC,UAAU,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;;YAEzE,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,eAAe,CAAC,UAAU,CAAC,CAAC,CAAC;IAC9D,CAAC;IACD;;;;;OAKG;IACI,qBAAqB,CAAC,KAAmC,EAAE,MAAgB;QAChF,MAAM,OAAO,GAAG,0BAAQ,CAAC,cAAc,CAAC,KAAK,CAAC,YAAY,EAAE,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QAC3E,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,CAAC,UAAU,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC;QAC9D,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,CAAC,gBAAgB,EAAE,CAAC;QAChD,MAAM,KAAK,GAAG,OAAO,CAAC,UAAU,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC;QACvD,MAAM,QAAQ,GAAG,mBAAQ,CAAC,yBAAyB,CAAC,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;QACnE,IAAI,SAAS,KAAK,QAAQ;YACxB,OAAO,SAAS,CAAC;QACnB,MAAM,SAAS,GAAG,mBAAQ,CAAC,yBAAyB,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;QACnE,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,SAAS,CAAC;QACnB,IAAI,MAAM,EAAE;YACV,IAAI,CAAC,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,SAAS,EAAE,QAAQ,EAAE,MAAM,CAAC,CAAC;SAC1D;QACD,OAAO,QAAQ,CAAC;IAClB,CAAC;IACD;;;;OAIG;IACI,uBAAuB,CAAC,KAAc,EAAE,MAAgB;QAC7D,IAAI,KAAK,CAAC,MAAM;YACd,OAAO,eAAO,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC;QACpC,MAAM,QAAQ,GAAG,eAAO,CAAC,QAAQ,CAAC,CAAC,mBAAQ,CAAC,qBAAqB,EAAE,mBAAQ,CAAC,qBAAqB,EAAE,MAAM,CAAC,CAAC;QAC3G,IAAI,QAAQ,CAAC,uBAAuB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC;eAC3F,QAAQ,CAAC,uBAAuB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC;eAC5F,QAAQ,CAAC,uBAAuB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC;YAC/F,OAAO,QAAQ,CAAC;QAClB,OAAO,QAAQ,CAAC;IAClB,CAAC;IAED;;;;OAIG;IACI,iCAAiC,CAAC,WAAmB,EAAE,MAAiB;QAC7E,MAAM,OAAO,GAAG,0BAAQ,CAAC,cAAc,CAAC,IAAI,CAAC,MAAM,EAAE,WAAW,CAAC,CAAC;QAClE,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,CAAC,gBAAgB,EAAE,CAAC;QAC7C,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,CAAC,kBAAkB,CAAC,IAAI,CAAC,MAAM,EAAE,WAAW,CAAC,CAAC;QACvE,MAAM,QAAQ,GAAG,mBAAQ,CAAC,kBAAkB,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;QAC1D,OAAO,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,SAAS,EAAE,CAAC,QAAQ,EAAE,MAAM,CAAC,CAAC;IAC/D,CAAC;IACD;;;;;;;;;;;;;;OAcG;IACI,MAAM,CAAC,yBAAyB,CAAC,IAAW,EAAE,IAAW;QAC9D,MAAM,qBAAqB,GAAG,0BAAQ,CAAC,MAAM,EAAE,CAAC;QAChD,IAAI,SAAS,EAAE,SAAS,CAAC;QACzB,IAAI,MAAM,EAAE,MAAM,CAAC;QACnB,IAAI,QAAQ,CAAC;QACb,IAAI,yBAAW,CAAC,mCAAmC,CACjD,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EACjG,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,qBAAqB,CAAC,EAAE;YAC3H,SAAS,GAAG,qBAAqB,CAAC,CAAC,CAAC;YACpC,SAAS,GAAG,qBAAqB,CAAC,CAAC,CAAC;YACpC,MAAM,GAAG,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;YACzC,MAAM,GAAG,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;YACzC,QAAQ,GAAG,MAAM,CAAC,mBAAmB,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,4CAAsB,CAAC,YAAY,CAAC,CAAC,CAAC,4CAAsB,CAAC,kBAAkB,CAAC;SACjI;aAAM;YACL,SAAS,GAAG,GAAG,CAAC;YAChB,SAAS,GAAG,IAAI,CAAC,eAAe,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;YAC9C,MAAM,GAAG,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;YACzC,MAAM,GAAG,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;YACzC,QAAQ,GAAG,MAAM,CAAC,mBAAmB,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,4CAAsB,CAAC,kBAAkB,CAAC,CAAC,CAAC,4CAAsB,CAAC,gBAAgB,CAAC;SACrI;QACD,MAAM,IAAI,GAAG,6CAAuB,CAAC,aAAa,CAChD,yCAAmB,CAAC,sBAAsB,CAAC,IAAI,EAAE,SAAS,EAAE,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC,EAC5F,yCAAmB,CAAC,sBAAsB,CAAC,IAAI,EAAE,SAAS,EAAE,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QAChG,IAAI,CAAC,YAAY,GAAG,QAAQ,CAAC;QAC7B,OAAO,IAAI,CAAC;IACd,CAAC;IACD;;;;;;;;OAQG;IACI,MAAM,CAAC,0BAA0B,CAAC,GAAW,EAAE,QAAgB,EAAE,GAAW,EAAE,YAAoB,EAAE,MAAc;QACvH,MAAM,GAAG,MAAM,IAAI,KAAK,CAAC,UAAU,EAAE,CAAC;QACtC,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QACzB,MAAM,CAAC,SAAS,CAAC,GAAG,CAAC,YAAY,GAAG,EAAE,EAAE,YAAY,GAAG,EAAE,EAAE,YAAY,GAAG,EAAE,CAAC,CAAC;QAC9E,IAAI,QAAQ,IAAI,GAAG;YACjB,MAAM,CAAC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,GAAG,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,QAAQ,GAAG,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,QAAQ,GAAG,EAAE,CAAC,CAAC;aACpF;YACH,MAAM,CAAC,GAAW,QAAQ,GAAG,GAAG,CAAC;YACjC,MAAM,CAAC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;SACnE;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;CACF;AAvZD,sBAuZC","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n\r\n/** @packageDocumentation\r\n * @module CartesianGeometry\r\n */\r\nimport { CurveCurveApproachType, CurveLocationDetail, CurveLocationDetailPair } from \"../curve/CurveLocationDetail\";\r\nimport { AxisOrder, BeJSONFunctions, Geometry } from \"../Geometry\";\r\nimport { SmallSystem } from \"../numerics/Polynomials\";\r\nimport { Matrix3d } from \"./Matrix3d\";\r\nimport { Plane3dByOriginAndUnitNormal } from \"./Plane3dByOriginAndUnitNormal\";\r\nimport { Vector2d } from \"./Point2dVector2d\";\r\nimport { Point3d, Vector3d } from \"./Point3dVector3d\";\r\nimport { Range1d, Range3d } from \"./Range\";\r\nimport { Transform } from \"./Transform\";\r\nimport { XYAndZ } from \"./XYZProps\";\r\n\r\n/**\r\n * A Ray3d contains\r\n * * an `origin` point.\r\n * * a `direction` vector (The vector is not required to be normalized).\r\n * * an optional weight (number).\r\n * @public\r\n */\r\nexport class Ray3d implements BeJSONFunctions {\r\n  /** The ray origin */\r\n  public origin: Point3d;\r\n  /** The ray direction. This is commonly (but not always) a unit vector. */\r\n  public direction: Vector3d;\r\n  /** Numeric annotation. */\r\n  public a?: number; // optional (e.g. weight)\r\n  // constructor (captures references)\r\n  private constructor(origin: Point3d, direction: Vector3d) {\r\n    this.origin = origin;\r\n    this.direction = direction;\r\n    this.a = undefined;\r\n  }\r\n  private static _create(x: number, y: number, z: number, u: number, v: number, w: number) {\r\n    return new Ray3d(Point3d.create(x, y, z), Vector3d.create(u, v, w));\r\n  }\r\n  /** Create a ray on the x axis. */\r\n  public static createXAxis(): Ray3d {\r\n    return Ray3d._create(0, 0, 0, 1, 0, 0);\r\n  }\r\n  /** Create a ray on the y axis. */\r\n  public static createYAxis(): Ray3d {\r\n    return Ray3d._create(0, 0, 0, 0, 1, 0);\r\n  }\r\n  /** Create a ray on the z axis. */\r\n  public static createZAxis(): Ray3d {\r\n    return Ray3d._create(0, 0, 0, 0, 0, 1);\r\n  }\r\n  /** Create a ray with all zeros. */\r\n  public static createZero(result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.origin.setZero();\r\n      result.direction.setZero();\r\n      return result;\r\n    }\r\n    return new Ray3d(Point3d.createZero(), Vector3d.createZero());\r\n  }\r\n  /**\r\n   * Test for nearly equal Ray3d objects.\r\n   * * This tests for near equality of origin and direction -- i.e. member-by-member comparison.\r\n   * * Use [[isAlmostEqualPointSet]] to allow origins to be anywhere along the common ray and to have to allow the\r\n   * directions to be scaled or opposing.\r\n   */\r\n  public isAlmostEqual(other: Ray3d): boolean {\r\n    return this.origin.isAlmostEqual(other.origin) && this.direction.isAlmostEqual(other.direction);\r\n  }\r\n  /** Return the dot product of the ray's direction vector with a vector from the ray origin to the `spacePoint`. */\r\n  public dotProductToPoint(spacePoint: Point3d): number {\r\n    return this.direction.dotProductStartEnd(this.origin, spacePoint);\r\n  }\r\n  /** Return the fractional coordinate (along the direction vector) of the `spacePoint` projected to the ray. */\r\n  public pointToFraction(spacePoint: Point3d): number {\r\n    return Geometry.safeDivideFraction(\r\n      this.dotProductToPoint(spacePoint),\r\n      this.direction.magnitudeSquared(),\r\n      0\r\n    );\r\n  }\r\n  /** Return the `spacePoint` projected onto the ray. */\r\n  public projectPointToRay(spacePoint: Point3d): Point3d {\r\n    /**\r\n     * To project a point to the ray, we can create a vector called \"v\" from ray origin to the spacePoint and project\r\n     * that vector to the ray direction vector \"r\". The projection is \"((v.r)/||r||^2) r\" where \"v.r\" is the dot\r\n     * product. Note that pointToFraction returns \"(v.r)/||r||^2\".\r\n     * Note: If r is the normal of a plane, then projection length \"(v.r)/||r||\" is the signed altitude of the\r\n     * spacePoint with respect to the plane.\r\n     */\r\n    return this.origin.plusScaled(this.direction, this.pointToFraction(spacePoint));\r\n  }\r\n  /**\r\n   * Test for nearly equal rays, allowing origin float and direction scaling.\r\n   * * Use [[isAlmostEqual]] to require member-by-member comparison.\r\n   */\r\n  public isAlmostEqualPointSet(other: Ray3d): boolean {\r\n    /**\r\n     * This function tests two rays to determine if they define the same infinite lines.\r\n     * So the origins can be different as long as they are on the infinite line (they can\r\n     * \"float\") but the directions must be parallel or antiparallel.\r\n     */\r\n    if (!this.direction.isParallelTo(other.direction, true))\r\n      return false;\r\n    /**\r\n     * In exact math, we consider a ray to have an infinite line as direction (not a finite vector).\r\n     * Therefore, in exact math it is not possible for one origin to be on the other ray but not vice\r\n     * versa. However, we test both ways because first check may pass due to round-off errors.\r\n     */\r\n    let workPoint = this.projectPointToRay(other.origin);\r\n    if (!other.origin.isAlmostEqualMetric(workPoint))\r\n      return false;\r\n    workPoint = other.projectPointToRay(this.origin);\r\n    if (!this.origin.isAlmostEqualMetric(workPoint))\r\n      return false;\r\n    return true;\r\n  }\r\n  /** Create a ray from origin and direction. */\r\n  public static create(origin: Point3d, direction: Vector3d, result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.set(origin, direction);\r\n      return result;\r\n    }\r\n    return new Ray3d(origin.clone(), direction.clone());\r\n  }\r\n  /**\r\n   * Given a homogeneous point and its derivative components, construct a Ray3d with cartesian\r\n   * coordinates and derivatives.\r\n   * @param weightedPoint `[x,y,z,w]` parts of weighted point.\r\n   * @param weightedDerivative `[x,y,z,w]` derivatives\r\n   * @param result\r\n   */\r\n  public static createWeightedDerivative(\r\n    weightedPoint: Float64Array, weightedDerivative: Float64Array, result?: Ray3d\r\n  ): Ray3d | undefined {\r\n    const w = weightedPoint[3];\r\n    const dw = weightedDerivative[3];\r\n    const x = weightedPoint[0];\r\n    const y = weightedPoint[1];\r\n    const z = weightedPoint[2];\r\n    const dx = weightedDerivative[0] * w - weightedPoint[0] * dw;\r\n    const dy = weightedDerivative[1] * w - weightedPoint[1] * dw;\r\n    const dz = weightedDerivative[2] * w - weightedPoint[2] * dw;\r\n    if (Geometry.isSmallMetricDistance(w))\r\n      return undefined;\r\n    const divW = 1.0 / w;\r\n    const divWW = divW * divW;\r\n    return Ray3d.createXYZUVW(\r\n      x * divW, y * divW, z * divW,\r\n      dx * divWW, dy * divWW, dz * divWW,\r\n      result\r\n    );\r\n  }\r\n  /** Create from coordinates of the origin and direction. */\r\n  public static createXYZUVW(\r\n    originX: number, originY: number, originZ: number,\r\n    directionX: number, directionY: number, directionZ: number,\r\n    result?: Ray3d\r\n  ): Ray3d {\r\n    if (result) {\r\n      result.getOriginRef().set(originX, originY, originZ);\r\n      result.getDirectionRef().set(directionX, directionY, directionZ);\r\n      return result;\r\n    }\r\n    return new Ray3d(Point3d.create(originX, originY, originZ), Vector3d.create(directionX, directionY, directionZ));\r\n  }\r\n  /** Capture origin and direction in a new Ray3d. */\r\n  public static createCapture(origin: Point3d, direction: Vector3d): Ray3d {\r\n    return new Ray3d(origin, direction);\r\n  }\r\n  /** Create from (clones of) origin, direction, and numeric weight. */\r\n  public static createPointVectorNumber(origin: Point3d, direction: Vector3d, a: number, result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.origin.setFrom(origin);\r\n      result.direction.setFrom(direction);\r\n      result.a = a;\r\n      return result;\r\n    }\r\n    result = new Ray3d(origin.clone(), direction.clone());\r\n    result.a = a;\r\n    return result;\r\n  }\r\n  /** Create from origin and target. The direction vector is the full length (non-unit) vector from origin to target. */\r\n  public static createStartEnd(origin: Point3d, target: Point3d, result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.origin.setFrom(origin);\r\n      result.direction.setStartEnd(origin, target);\r\n      return result;\r\n    }\r\n    return new Ray3d(origin.clone(), Vector3d.createStartEnd(origin, target));\r\n  }\r\n  /** Return a reference to the ray's origin. */\r\n  public getOriginRef(): Point3d {\r\n    return this.origin;\r\n  }\r\n  /** Return a reference to the ray's direction vector. */\r\n  public getDirectionRef(): Vector3d {\r\n    return this.direction;\r\n  }\r\n  /** Copy coordinates from origin and direction. */\r\n  public set(origin: Point3d, direction: Vector3d): void {\r\n    this.origin.setFrom(origin);\r\n    this.direction.setFrom(direction);\r\n  }\r\n  /** Clone the ray. */\r\n  public clone(result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.set(this.origin.clone(), this.direction.clone());\r\n      return result;\r\n    }\r\n    return new Ray3d(this.origin.clone(), this.direction.clone());\r\n  }\r\n  /** Create a clone and return the transform of the clone. */\r\n  public cloneTransformed(transform: Transform): Ray3d {\r\n    return new Ray3d(transform.multiplyPoint3d(this.origin), transform.multiplyVector(this.direction));\r\n  }\r\n  /** Create a clone and return the inverse transform of the clone. */\r\n  public cloneInverseTransformed(transform: Transform): Ray3d | undefined {\r\n    const origin = transform.multiplyInversePoint3d(this.origin);\r\n    const direction = transform.matrix.multiplyInverseXYZAsVector3d(this.direction.x, this.direction.y, this.direction.z);\r\n    if (undefined !== origin && undefined !== direction)\r\n      return new Ray3d(origin, direction);\r\n    return undefined;\r\n  }\r\n  /** Apply a transform in place. */\r\n  public transformInPlace(transform: Transform) {\r\n    transform.multiplyPoint3d(this.origin, this.origin);\r\n    transform.multiplyVector(this.direction, this.direction);\r\n  }\r\n  /** Copy data from another ray. */\r\n  public setFrom(source: Ray3d): void {\r\n    this.set(source.origin, source.direction);\r\n  }\r\n  /**\r\n   * Return a point at fractional position along the ray.\r\n   * * fraction 0 is the ray origin.\r\n   * * fraction 1 is at the end of the direction vector when placed at the origin.\r\n   */\r\n  public fractionToPoint(fraction: number, result?: Point3d): Point3d {\r\n    return this.origin.plusScaled(this.direction, fraction, result);\r\n  }\r\n  /**\r\n   * Return a transform for rigid axes at ray origin with z in ray direction.\r\n   * * If the direction vector is zero, axes default to identity (from [[Matrix3d.createRigidHeadsUp]])\r\n   */\r\n  public toRigidZFrame(): Transform | undefined {\r\n    const axes = Matrix3d.createRigidHeadsUp(this.direction, AxisOrder.ZXY);\r\n    return Transform.createOriginAndMatrix(this.origin, axes);\r\n  }\r\n  /** Convert {origin:[x,y,z], direction:[u,v,w]} to a Ray3d. */\r\n  public setFromJSON(json?: any) {\r\n    if (!json) {\r\n      this.origin.set(0, 0, 0);\r\n      this.direction.set(0, 0, 1);\r\n      return;\r\n    }\r\n    this.origin.setFromJSON(json.origin);\r\n    this.direction.setFromJSON(json.direction);\r\n  }\r\n  /**\r\n   * Try to scale the direction vector to a given magnitude.\r\n   * * Returns false if ray direction is a zero vector.\r\n   */\r\n  public trySetDirectionMagnitudeInPlace(magnitude: number = 1.0): boolean {\r\n    if (this.direction.tryNormalizeInPlace()) {\r\n      this.direction.scaleInPlace(magnitude);\r\n      return true;\r\n    }\r\n    this.direction.setZero();\r\n    this.a = 0.0;\r\n    return false;\r\n  }\r\n  /**\r\n   * * If parameter `a` is clearly nonzero and the direction vector can be normalized,\r\n   *    * save the parameter `a` as the optional `a` member of the ray.\r\n   *    * normalize the ray's direction vector\r\n   * * If parameter `a` is nearly zero,\r\n   *   * Set the `a` member to zero\r\n   *   * Set the ray's direction vector to zero.\r\n   * @param a area to be saved.\r\n   */\r\n  // input a ray and \"a\" understood as an area.\r\n  // if a is clearly nonzero metric squared and the vector can be normalized, install those and return true.\r\n  // otherwise set ray.z to zero and zero the vector of the ray and return false.\r\n  public tryNormalizeInPlaceWithAreaWeight(a: number): boolean {\r\n    const tolerance = Geometry.smallMetricDistanceSquared;\r\n    this.a = a;\r\n    if (Math.abs(a) > tolerance && this.direction.tryNormalizeInPlace(tolerance))\r\n      return true;\r\n    this.direction.setZero();\r\n    this.a = 0.0;\r\n    return false;\r\n  }\r\n  /**\r\n   * Construct a JSON object from this Ray3d.\r\n   * @return {*} [origin,normal]\r\n   */\r\n  public toJSON(): any {\r\n    return { origin: this.origin.toJSON(), direction: this.direction.toJSON() };\r\n  }\r\n  /** Create a new ray from json object. See `setFromJSON` for json structure; */\r\n  public static fromJSON(json?: any) {\r\n    const result = Ray3d.createXAxis();\r\n    result.setFromJSON(json);\r\n    return result;\r\n  }\r\n  /** return distance from the ray to point in space */\r\n  public distance(spacePoint: Point3d): number {\r\n    const uu = this.direction.magnitudeSquared();\r\n    const uv = this.dotProductToPoint(spacePoint);\r\n    const aa = Geometry.inverseMetricDistanceSquared(uu);\r\n    if (aa)\r\n      return Math.sqrt(this.origin.distanceSquared(spacePoint) - uv * uv * aa);\r\n    else\r\n      return Math.sqrt(this.origin.distanceSquared(spacePoint));\r\n  }\r\n  /**\r\n   * Return the intersection of the unbounded ray with a plane.\r\n   * Stores the point of intersection in the result point given as a parameter,\r\n   * and returns the parameter along the ray where the intersection occurs.\r\n   * Returns undefined if the ray and plane are parallel or coplanar.\r\n   */\r\n  public intersectionWithPlane(plane: Plane3dByOriginAndUnitNormal, result?: Point3d): number | undefined {\r\n    const vectorA = Vector3d.createStartEnd(plane.getOriginRef(), this.origin);\r\n    const uDotN = this.direction.dotProduct(plane.getNormalRef());\r\n    const nDotN = this.direction.magnitudeSquared();\r\n    const aDotN = vectorA.dotProduct(plane.getNormalRef());\r\n    const division = Geometry.conditionalDivideFraction(-aDotN, uDotN);\r\n    if (undefined === division)\r\n      return undefined;\r\n    const division1 = Geometry.conditionalDivideFraction(nDotN, uDotN);\r\n    if (undefined === division1)\r\n      return undefined;\r\n    if (result) {\r\n      this.origin.plusScaled(this.direction, division, result);\r\n    }\r\n    return division;\r\n  }\r\n  /**\r\n   * Find intersection of the ray with a Range3d.\r\n   * * return the range of fractions (on the ray) which are \"inside\" the range.\r\n   * * Note that a range is always returned;  if there is no intersection it is indicated by the test `result.sNull`\r\n   */\r\n  public intersectionWithRange3d(range: Range3d, result?: Range1d): Range1d {\r\n    if (range.isNull)\r\n      return Range1d.createNull(result);\r\n    const interval = Range1d.createXX(-Geometry.largeCoordinateResult, Geometry.largeCoordinateResult, result);\r\n    if (interval.clipLinearMapToInterval(this.origin.x, this.direction.x, range.low.x, range.high.x)\r\n      && interval.clipLinearMapToInterval(this.origin.y, this.direction.y, range.low.y, range.high.y)\r\n      && interval.clipLinearMapToInterval(this.origin.z, this.direction.z, range.low.z, range.high.z))\r\n      return interval;\r\n    return interval;\r\n  }\r\n\r\n  /**\r\n   * Construct a vector from `ray.origin` to target point.\r\n   * * return the part of the vector that is perpendicular to `ray.direction`.\r\n   * * i.e., return the shortest vector from the ray to the point.\r\n   */\r\n  public perpendicularPartOfVectorToTarget(targetPoint: XYAndZ, result?: Vector3d): Vector3d {\r\n    const vectorV = Vector3d.createStartEnd(this.origin, targetPoint);\r\n    const uu = this.direction.magnitudeSquared();\r\n    const uv = this.direction.dotProductStartEnd(this.origin, targetPoint);\r\n    const fraction = Geometry.safeDivideFraction(uv, uu, 0.0);\r\n    return vectorV.plusScaled(this.direction, -fraction, result);\r\n  }\r\n  /**\r\n   * Determine if two rays intersect, are fully overlapped, parallel but no coincident, or skew\r\n   * * Return a CurveLocationDetailPair which\r\n   * * contains fraction and point on each ray.\r\n   * * has (in the CurveLocationDetailPair structure, as member approachType) annotation indicating one of\r\n   * these relationships\r\n   *   * CurveCurveApproachType.Intersection -- the rays have a simple intersection, at fractions indicated\r\n   * in detailA and detailB\r\n   *   * CurveCurveApproachType.PerpendicularChord -- there is pair of where the rays have closest approach.\r\n   * The rays are skew in space.\r\n   *   * CurveCurveApproachType.CoincidentGeometry -- the rays are the same unbounded line in space. The\r\n   * fractions and points are a representative single common point.\r\n   *   * CurveCurveApproachType.Parallel -- the rays are parallel (and not coincident). The two points are\r\n   * at the minimum distance\r\n   */\r\n  public static closestApproachRay3dRay3d(rayA: Ray3d, rayB: Ray3d): CurveLocationDetailPair {\r\n    const intersectionFractions = Vector2d.create();\r\n    let fractionA, fractionB;\r\n    let pointA, pointB;\r\n    let pairType;\r\n    if (SmallSystem.ray3dXYZUVWClosestApproachUnbounded(\r\n      rayA.origin.x, rayA.origin.y, rayA.origin.z, rayA.direction.x, rayA.direction.y, rayA.direction.z,\r\n      rayB.origin.x, rayB.origin.y, rayB.origin.z, rayB.direction.x, rayB.direction.y, rayB.direction.z, intersectionFractions)) {\r\n      fractionA = intersectionFractions.x;\r\n      fractionB = intersectionFractions.y;\r\n      pointA = rayA.fractionToPoint(fractionA);\r\n      pointB = rayB.fractionToPoint(fractionB);\r\n      pairType = pointA.isAlmostEqualMetric(pointB) ? CurveCurveApproachType.Intersection : CurveCurveApproachType.PerpendicularChord;\r\n    } else {\r\n      fractionB = 0.0;\r\n      fractionA = rayA.pointToFraction(rayB.origin);\r\n      pointA = rayA.fractionToPoint(fractionA);\r\n      pointB = rayB.fractionToPoint(fractionB);\r\n      pairType = pointA.isAlmostEqualMetric(pointB) ? CurveCurveApproachType.CoincidentGeometry : CurveCurveApproachType.ParallelGeometry;\r\n    }\r\n    const pair = CurveLocationDetailPair.createCapture(\r\n      CurveLocationDetail.createRayFractionPoint(rayA, fractionA, rayA.fractionToPoint(fractionA)),\r\n      CurveLocationDetail.createRayFractionPoint(rayB, fractionB, rayB.fractionToPoint(fractionB)));\r\n    pair.approachType = pairType;\r\n    return pair;\r\n  }\r\n  /**\r\n   * Return a ray with `ray.origin` interpolated between pt1 and pt2 at the given fraction\r\n   * and `ray.direction` set to the vector from pt1 to pt2 multiplied by the given scale factor.\r\n   * @param pt1 start point of interpolation.\r\n   * @param fraction fractional position between points.\r\n   * @param pt2 endpoint of interpolation.\r\n   * @param tangentScale scale factor to apply to the startToEnd vector.\r\n   * @param result optional receiver.\r\n   */\r\n  public static interpolatePointAndTangent(pt1: XYAndZ, fraction: number, pt2: XYAndZ, tangentScale: number, result?: Ray3d): Ray3d {\r\n    result = result ?? Ray3d.createZero();\r\n    const dx = pt2.x - pt1.x;\r\n    const dy = pt2.y - pt1.y;\r\n    const dz = pt2.z - pt1.z;\r\n    result.direction.set(tangentScale * dx, tangentScale * dy, tangentScale * dz);\r\n    if (fraction <= 0.5)\r\n      result.origin.set(pt1.x + fraction * dx, pt1.y + fraction * dy, pt1.z + fraction * dz);\r\n    else {\r\n      const t: number = fraction - 1.0;\r\n      result.origin.set(pt2.x + t * dx, pt2.y + t * dy, pt2.z + t * dz);\r\n    }\r\n    return result;\r\n  }\r\n}\r\n"]}
+{"version":3,"file":"Ray3d.js","sourceRoot":"","sources":["../../../src/geometry3d/Ray3d.ts"],"names":[],"mappings":";AAAA;;;+FAG+F;;;AAE/F;;GAEG;AACH,sEAAoH;AACpH,0CAAmE;AACnE,yDAAsD;AACtD,yCAAsC;AAEtC,uDAA6C;AAC7C,uDAAsD;AACtD,mCAA2C;AAC3C,2CAAwC;AAGxC;;;;;;GAMG;AACH,MAAa,KAAK;IAOhB,oCAAoC;IACpC,YAAoB,MAAe,EAAE,SAAmB;QACtD,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;QACrB,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC;QAC3B,IAAI,CAAC,CAAC,GAAG,SAAS,CAAC;IACrB,CAAC;IACO,MAAM,CAAC,OAAO,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS;QACrF,OAAO,IAAI,KAAK,CAAC,yBAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,0BAAQ,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACtE,CAAC;IACD,kCAAkC;IAC3B,MAAM,CAAC,WAAW;QACvB,OAAO,KAAK,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACzC,CAAC;IACD,kCAAkC;IAC3B,MAAM,CAAC,WAAW;QACvB,OAAO,KAAK,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACzC,CAAC;IACD,kCAAkC;IAC3B,MAAM,CAAC,WAAW;QACvB,OAAO,KAAK,CAAC,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACzC,CAAC;IACD,mCAAmC;IAC5B,MAAM,CAAC,UAAU,CAAC,MAAc;QACrC,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,MAAM,CAAC,OAAO,EAAE,CAAC;YACxB,MAAM,CAAC,SAAS,CAAC,OAAO,EAAE,CAAC;YAC3B,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,KAAK,CAAC,yBAAO,CAAC,UAAU,EAAE,EAAE,0BAAQ,CAAC,UAAU,EAAE,CAAC,CAAC;IAChE,CAAC;IACD;;;;;OAKG;IACI,aAAa,CAAC,KAAY;QAC/B,OAAO,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,KAAK,CAAC,MAAM,CAAC,IAAI,IAAI,CAAC,SAAS,CAAC,aAAa,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC;IAClG,CAAC;IACD,kHAAkH;IAC3G,iBAAiB,CAAC,UAAmB;QAC1C,OAAO,IAAI,CAAC,SAAS,CAAC,kBAAkB,CAAC,IAAI,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;IACpE,CAAC;IACD,8GAA8G;IACvG,eAAe,CAAC,UAAmB;QACxC,OAAO,mBAAQ,CAAC,kBAAkB,CAChC,IAAI,CAAC,iBAAiB,CAAC,UAAU,CAAC,EAClC,IAAI,CAAC,SAAS,CAAC,gBAAgB,EAAE,EACjC,CAAC,CACF,CAAC;IACJ,CAAC;IACD,sDAAsD;IAC/C,iBAAiB,CAAC,UAAmB;QAC1C;;;;;;WAMG;QACH,OAAO,IAAI,CAAC,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,SAAS,EAAE,IAAI,CAAC,eAAe,CAAC,UAAU,CAAC,CAAC,CAAC;IAClF,CAAC;IACD;;;OAGG;IACI,qBAAqB,CAAC,KAAY;QACvC;;;;WAIG;QACH,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,YAAY,CAAC,KAAK,CAAC,SAAS,EAAE,IAAI,CAAC;YACrD,OAAO,KAAK,CAAC;QACf;;;;WAIG;QACH,IAAI,SAAS,GAAG,IAAI,CAAC,iBAAiB,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QACrD,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,mBAAmB,CAAC,SAAS,CAAC;YAC9C,OAAO,KAAK,CAAC;QACf,SAAS,GAAG,KAAK,CAAC,iBAAiB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACjD,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,mBAAmB,CAAC,SAAS,CAAC;YAC7C,OAAO,KAAK,CAAC;QACf,OAAO,IAAI,CAAC;IACd,CAAC;IACD,8CAA8C;IACvC,MAAM,CAAC,MAAM,CAAC,MAAe,EAAE,SAAmB,EAAE,MAAc;QACvE,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,GAAG,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;YAC9B,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,KAAK,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC;IACtD,CAAC;IACD;;;;;;OAMG;IACI,MAAM,CAAC,wBAAwB,CACpC,aAA2B,EAAE,kBAAgC,EAAE,MAAc;QAE7E,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,EAAE,GAAG,kBAAkB,CAAC,CAAC,CAAC,CAAC;QACjC,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAAC;QAC3B,MAAM,EAAE,GAAG,kBAAkB,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAC7D,MAAM,EAAE,GAAG,kBAAkB,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAC7D,MAAM,EAAE,GAAG,kBAAkB,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QAC7D,IAAI,mBAAQ,CAAC,qBAAqB,CAAC,CAAC,CAAC;YACnC,OAAO,SAAS,CAAC;QACnB,MAAM,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC;QACrB,MAAM,KAAK,GAAG,IAAI,GAAG,IAAI,CAAC;QAC1B,OAAO,KAAK,CAAC,YAAY,CACvB,CAAC,GAAG,IAAI,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,GAAG,IAAI,EAC5B,EAAE,GAAG,KAAK,EAAE,EAAE,GAAG,KAAK,EAAE,EAAE,GAAG,KAAK,EAClC,MAAM,CACP,CAAC;IACJ,CAAC;IACD,2DAA2D;IACpD,MAAM,CAAC,YAAY,CACxB,OAAe,EAAE,OAAe,EAAE,OAAe,EACjD,UAAkB,EAAE,UAAkB,EAAE,UAAkB,EAC1D,MAAc;QAEd,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,YAAY,EAAE,CAAC,GAAG,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;YACrD,MAAM,CAAC,eAAe,EAAE,CAAC,GAAG,CAAC,UAAU,EAAE,UAAU,EAAE,UAAU,CAAC,CAAC;YACjE,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,KAAK,CAAC,yBAAO,CAAC,MAAM,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,EAAE,0BAAQ,CAAC,MAAM,CAAC,UAAU,EAAE,UAAU,EAAE,UAAU,CAAC,CAAC,CAAC;IACnH,CAAC;IACD,mDAAmD;IAC5C,MAAM,CAAC,aAAa,CAAC,MAAe,EAAE,SAAmB;QAC9D,OAAO,IAAI,KAAK,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;IACtC,CAAC;IACD,qEAAqE;IAC9D,MAAM,CAAC,uBAAuB,CAAC,MAAe,EAAE,SAAmB,EAAE,CAAS,EAAE,MAAc;QACnG,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;YAC9B,MAAM,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;YACpC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;YACb,OAAO,MAAM,CAAC;SACf;QACD,MAAM,GAAG,IAAI,KAAK,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC;QACtD,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC;QACb,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,sHAAsH;IAC/G,MAAM,CAAC,cAAc,CAAC,MAAe,EAAE,MAAe,EAAE,MAAc;QAC3E,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;YAC9B,MAAM,CAAC,SAAS,CAAC,WAAW,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;YAC7C,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,KAAK,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,0BAAQ,CAAC,cAAc,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,CAAC;IAC5E,CAAC;IACD,8CAA8C;IACvC,YAAY;QACjB,OAAO,IAAI,CAAC,MAAM,CAAC;IACrB,CAAC;IACD,wDAAwD;IACjD,eAAe;QACpB,OAAO,IAAI,CAAC,SAAS,CAAC;IACxB,CAAC;IACD,kDAAkD;IAC3C,GAAG,CAAC,MAAe,EAAE,SAAmB;QAC7C,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;QAC5B,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;IACpC,CAAC;IACD,qBAAqB;IACd,KAAK,CAAC,MAAc;QACzB,IAAI,MAAM,EAAE;YACV,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC;YACxD,OAAO,MAAM,CAAC;SACf;QACD,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC;IAChE,CAAC;IACD,4DAA4D;IACrD,gBAAgB,CAAC,SAAoB;QAC1C,OAAO,IAAI,KAAK,CAAC,SAAS,CAAC,eAAe,CAAC,IAAI,CAAC,MAAM,CAAC,EAAE,SAAS,CAAC,cAAc,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC;IACrG,CAAC;IACD,oEAAoE;IAC7D,uBAAuB,CAAC,SAAoB;QACjD,MAAM,MAAM,GAAG,SAAS,CAAC,sBAAsB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QAC7D,MAAM,SAAS,GAAG,SAAS,CAAC,MAAM,CAAC,4BAA4B,CAC7D,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,CACrD,CAAC;QACF,IAAI,SAAS,KAAK,MAAM,IAAI,SAAS,KAAK,SAAS;YACjD,OAAO,IAAI,KAAK,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;QACtC,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,kCAAkC;IAC3B,gBAAgB,CAAC,SAAoB;QAC1C,SAAS,CAAC,eAAe,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QACpD,SAAS,CAAC,cAAc,CAAC,IAAI,CAAC,SAAS,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC;IAC3D,CAAC;IACD,kCAAkC;IAC3B,OAAO,CAAC,MAAa;QAC1B,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,SAAS,CAAC,CAAC;IAC5C,CAAC;IACD;;;;OAIG;IACI,eAAe,CAAC,QAAgB,EAAE,MAAgB;QACvD,OAAO,IAAI,CAAC,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,SAAS,EAAE,QAAQ,EAAE,MAAM,CAAC,CAAC;IAClE,CAAC;IACD;;;OAGG;IACI,aAAa;QAClB,MAAM,IAAI,GAAG,mBAAQ,CAAC,kBAAkB,CAAC,IAAI,CAAC,SAAS,EAAE,oBAAS,CAAC,GAAG,CAAC,CAAC;QACxE,OAAO,qBAAS,CAAC,qBAAqB,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;IAC5D,CAAC;IACD,8DAA8D;IACvD,WAAW,CAAC,IAAU;QAC3B,IAAI,CAAC,IAAI,EAAE;YACT,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACzB,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAC5B,OAAO;SACR;QACD,IAAI,CAAC,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QACrC,IAAI,CAAC,SAAS,CAAC,WAAW,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;IAC7C,CAAC;IACD;;;OAGG;IACI,MAAM;QACX,OAAO,EAAE,MAAM,EAAE,IAAI,CAAC,MAAM,CAAC,MAAM,EAAE,EAAE,SAAS,EAAE,IAAI,CAAC,SAAS,CAAC,MAAM,EAAE,EAAE,CAAC;IAC9E,CAAC;IACD,+EAA+E;IACxE,MAAM,CAAC,QAAQ,CAAC,IAAU;QAC/B,MAAM,MAAM,GAAG,KAAK,CAAC,WAAW,EAAE,CAAC;QACnC,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,CAAC;QACzB,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;OAGG;IACI,+BAA+B,CAAC,YAAoB,GAAG;QAC5D,IAAI,IAAI,CAAC,SAAS,CAAC,mBAAmB,EAAE,EAAE;YACxC,IAAI,CAAC,SAAS,CAAC,YAAY,CAAC,SAAS,CAAC,CAAC;YACvC,OAAO,IAAI,CAAC;SACb;QACD,IAAI,CAAC,SAAS,CAAC,OAAO,EAAE,CAAC;QACzB,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC;QACb,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;;;;;;;;OAUG;IACI,iCAAiC,CAAC,CAAS;QAChD,MAAM,SAAS,GAAG,mBAAQ,CAAC,0BAA0B,CAAC;QACtD,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;QACX,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,SAAS,IAAI,IAAI,CAAC,SAAS,CAAC,mBAAmB,CAAC,SAAS,CAAC;YAC1E,OAAO,IAAI,CAAC;QACd,IAAI,CAAC,SAAS,CAAC,OAAO,EAAE,CAAC;QACzB,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC;QACb,OAAO,KAAK,CAAC;IACf,CAAC;IACD,sDAAsD;IAC/C,QAAQ,CAAC,UAAmB;QACjC,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,CAAC,gBAAgB,EAAE,CAAC;QAC7C,MAAM,EAAE,GAAG,IAAI,CAAC,iBAAiB,CAAC,UAAU,CAAC,CAAC;QAC9C,MAAM,EAAE,GAAG,mBAAQ,CAAC,4BAA4B,CAAC,EAAE,CAAC,CAAC;QACrD,IAAI,EAAE;YACJ,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,eAAe,CAAC,UAAU,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;;YAEzE,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,eAAe,CAAC,UAAU,CAAC,CAAC,CAAC;IAC9D,CAAC;IACD;;;;;;;;;;OAUG;IACI,qBAAqB,CAAC,KAAmC,EAAE,MAAgB;QAChF,MAAM,OAAO,GAAG,0BAAQ,CAAC,cAAc,CAAC,KAAK,CAAC,YAAY,EAAE,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QAC3E,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,CAAC,UAAU,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC;QAC9D,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,CAAC,gBAAgB,EAAE,CAAC;QAChD,MAAM,KAAK,GAAG,OAAO,CAAC,UAAU,CAAC,KAAK,CAAC,YAAY,EAAE,CAAC,CAAC;QACvD,MAAM,QAAQ,GAAG,mBAAQ,CAAC,yBAAyB,CAAC,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;QACnE,IAAI,SAAS,KAAK,QAAQ;YACxB,OAAO,SAAS,CAAC;QACnB,MAAM,SAAS,GAAG,mBAAQ,CAAC,yBAAyB,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;QACnE,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,SAAS,CAAC;QACnB,IAAI,MAAM,EAAE;YACV,IAAI,CAAC,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,SAAS,EAAE,QAAQ,EAAE,MAAM,CAAC,CAAC;SAC1D;QACD,OAAO,QAAQ,CAAC;IAClB,CAAC;IACD;;;;OAIG;IACI,uBAAuB,CAAC,KAAc,EAAE,MAAgB;QAC7D,IAAI,KAAK,CAAC,MAAM;YACd,OAAO,eAAO,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC;QACpC,MAAM,QAAQ,GAAG,eAAO,CAAC,QAAQ,CAAC,CAAC,mBAAQ,CAAC,qBAAqB,EAAE,mBAAQ,CAAC,qBAAqB,EAAE,MAAM,CAAC,CAAC;QAC3G,IAAI,QAAQ,CAAC,uBAAuB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC;eAC3F,QAAQ,CAAC,uBAAuB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC;eAC5F,QAAQ,CAAC,uBAAuB,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,KAAK,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC;YAE/F,OAAO,QAAQ,CAAC;QAClB,OAAO,QAAQ,CAAC;IAClB,CAAC;IACD;;;OAGG;IACI,iCAAiC,CAAC,WAAmB,EAAE,MAAiB;QAC7E,MAAM,OAAO,GAAG,0BAAQ,CAAC,cAAc,CAAC,IAAI,CAAC,MAAM,EAAE,WAAW,CAAC,CAAC;QAClE,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,CAAC,gBAAgB,EAAE,CAAC;QAC7C,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,CAAC,kBAAkB,CAAC,IAAI,CAAC,MAAM,EAAE,WAAW,CAAC,CAAC;QACvE,MAAM,QAAQ,GAAG,mBAAQ,CAAC,kBAAkB,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,CAAC,CAAC;QAC1D,OAAO,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,SAAS,EAAE,CAAC,QAAQ,EAAE,MAAM,CAAC,CAAC;IAC/D,CAAC;IACD;;;;;;;;;;;;OAYG;IACI,MAAM,CAAC,yBAAyB,CAAC,IAAW,EAAE,IAAW;QAC9D,MAAM,qBAAqB,GAAG,0BAAQ,CAAC,MAAM,EAAE,CAAC;QAChD,IAAI,SAAS,EAAE,SAAS,CAAC;QACzB,IAAI,MAAM,EAAE,MAAM,CAAC;QACnB,IAAI,QAAQ,CAAC;QACb,IAAI,yBAAW,CAAC,mCAAmC,CACjD,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EACjG,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,EACjG,qBAAqB,CAAC,EACtB;YACA,SAAS,GAAG,qBAAqB,CAAC,CAAC,CAAC;YACpC,SAAS,GAAG,qBAAqB,CAAC,CAAC,CAAC;YACpC,MAAM,GAAG,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;YACzC,MAAM,GAAG,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;YACzC,QAAQ,GAAG,MAAM,CAAC,mBAAmB,CAAC,MAAM,CAAC,CAAC,CAAC;gBAC7C,4CAAsB,CAAC,YAAY,CAAC,CAAC,CAAC,4CAAsB,CAAC,kBAAkB,CAAC;SACnF;aAAM;YACL,SAAS,GAAG,GAAG,CAAC;YAChB,SAAS,GAAG,IAAI,CAAC,eAAe,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;YAC9C,MAAM,GAAG,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;YACzC,MAAM,GAAG,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC;YACzC,QAAQ,GAAG,MAAM,CAAC,mBAAmB,CAAC,MAAM,CAAC,CAAC,CAAC;gBAC7C,4CAAsB,CAAC,kBAAkB,CAAC,CAAC,CAAC,4CAAsB,CAAC,gBAAgB,CAAC;SACvF;QACD,MAAM,IAAI,GAAG,6CAAuB,CAAC,aAAa,CAChD,yCAAmB,CAAC,sBAAsB,CAAC,IAAI,EAAE,SAAS,EAAE,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC,EAC5F,yCAAmB,CAAC,sBAAsB,CAAC,IAAI,EAAE,SAAS,EAAE,IAAI,CAAC,eAAe,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;QAChG,IAAI,CAAC,YAAY,GAAG,QAAQ,CAAC;QAC7B,OAAO,IAAI,CAAC;IACd,CAAC;IACD;;;;;;;;OAQG;IACI,MAAM,CAAC,0BAA0B,CACtC,GAAW,EAAE,QAAgB,EAAE,GAAW,EAAE,YAAoB,EAAE,MAAc;QAEhF,MAAM,GAAG,MAAM,IAAI,KAAK,CAAC,UAAU,EAAE,CAAC;QACtC,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QACzB,MAAM,CAAC,SAAS,CAAC,GAAG,CAAC,YAAY,GAAG,EAAE,EAAE,YAAY,GAAG,EAAE,EAAE,YAAY,GAAG,EAAE,CAAC,CAAC;QAC9E,IAAI,QAAQ,IAAI,GAAG;YACjB,MAAM,CAAC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,GAAG,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,QAAQ,GAAG,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,QAAQ,GAAG,EAAE,CAAC,CAAC;aACpF;YACH,MAAM,CAAC,GAAW,QAAQ,GAAG,GAAG,CAAC;YACjC,MAAM,CAAC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;SACnE;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;CACF;AAhaD,sBAgaC","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n\r\n/** @packageDocumentation\r\n * @module CartesianGeometry\r\n */\r\nimport { CurveCurveApproachType, CurveLocationDetail, CurveLocationDetailPair } from \"../curve/CurveLocationDetail\";\r\nimport { AxisOrder, BeJSONFunctions, Geometry } from \"../Geometry\";\r\nimport { SmallSystem } from \"../numerics/Polynomials\";\r\nimport { Matrix3d } from \"./Matrix3d\";\r\nimport { Plane3dByOriginAndUnitNormal } from \"./Plane3dByOriginAndUnitNormal\";\r\nimport { Vector2d } from \"./Point2dVector2d\";\r\nimport { Point3d, Vector3d } from \"./Point3dVector3d\";\r\nimport { Range1d, Range3d } from \"./Range\";\r\nimport { Transform } from \"./Transform\";\r\nimport { XYAndZ } from \"./XYZProps\";\r\n\r\n/**\r\n * A Ray3d contains\r\n * * an `origin` point.\r\n * * a `direction` vector (The vector is not required to be normalized).\r\n * * an optional weight (number).\r\n * @public\r\n */\r\nexport class Ray3d implements BeJSONFunctions {\r\n  /** The ray origin */\r\n  public origin: Point3d;\r\n  /** The ray direction. This is commonly (but not always) a unit vector. */\r\n  public direction: Vector3d;\r\n  /** Numeric annotation. */\r\n  public a?: number; // optional (e.g. weight)\r\n  // constructor (captures references)\r\n  private constructor(origin: Point3d, direction: Vector3d) {\r\n    this.origin = origin;\r\n    this.direction = direction;\r\n    this.a = undefined;\r\n  }\r\n  private static _create(x: number, y: number, z: number, u: number, v: number, w: number) {\r\n    return new Ray3d(Point3d.create(x, y, z), Vector3d.create(u, v, w));\r\n  }\r\n  /** Create a ray on the x axis. */\r\n  public static createXAxis(): Ray3d {\r\n    return Ray3d._create(0, 0, 0, 1, 0, 0);\r\n  }\r\n  /** Create a ray on the y axis. */\r\n  public static createYAxis(): Ray3d {\r\n    return Ray3d._create(0, 0, 0, 0, 1, 0);\r\n  }\r\n  /** Create a ray on the z axis. */\r\n  public static createZAxis(): Ray3d {\r\n    return Ray3d._create(0, 0, 0, 0, 0, 1);\r\n  }\r\n  /** Create a ray with all zeros. */\r\n  public static createZero(result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.origin.setZero();\r\n      result.direction.setZero();\r\n      return result;\r\n    }\r\n    return new Ray3d(Point3d.createZero(), Vector3d.createZero());\r\n  }\r\n  /**\r\n   * Test for nearly equal Ray3d objects.\r\n   * * This tests for near equality of origin and direction -- i.e. member-by-member comparison.\r\n   * * Use [[isAlmostEqualPointSet]] to allow origins to be anywhere along the common ray and to have to allow the\r\n   * directions to be scaled or opposing.\r\n   */\r\n  public isAlmostEqual(other: Ray3d): boolean {\r\n    return this.origin.isAlmostEqual(other.origin) && this.direction.isAlmostEqual(other.direction);\r\n  }\r\n  /** Return the dot product of the ray's direction vector with a vector from the ray origin to the `spacePoint`. */\r\n  public dotProductToPoint(spacePoint: Point3d): number {\r\n    return this.direction.dotProductStartEnd(this.origin, spacePoint);\r\n  }\r\n  /** Return the fractional coordinate (along the direction vector) of the `spacePoint` projected to the ray. */\r\n  public pointToFraction(spacePoint: Point3d): number {\r\n    return Geometry.safeDivideFraction(\r\n      this.dotProductToPoint(spacePoint),\r\n      this.direction.magnitudeSquared(),\r\n      0\r\n    );\r\n  }\r\n  /** Return the `spacePoint` projected onto the ray. */\r\n  public projectPointToRay(spacePoint: Point3d): Point3d {\r\n    /**\r\n     * To project a point to the ray, we can create a vector called \"v\" from ray origin to the spacePoint and project\r\n     * that vector to the ray direction vector \"r\". The projection is \"((v.r)/||r||^2) r\" where \"v.r\" is the dot\r\n     * product. Note that pointToFraction returns \"(v.r)/||r||^2\".\r\n     * Note: If r is the normal of a plane, then projection length \"(v.r)/||r||\" is the signed altitude of the\r\n     * spacePoint with respect to the plane.\r\n     */\r\n    return this.origin.plusScaled(this.direction, this.pointToFraction(spacePoint));\r\n  }\r\n  /**\r\n   * Test for nearly equal rays, allowing origin float and direction scaling.\r\n   * * Use [[isAlmostEqual]] to require member-by-member comparison.\r\n   */\r\n  public isAlmostEqualPointSet(other: Ray3d): boolean {\r\n    /**\r\n     * This function tests two rays to determine if they define the same infinite lines.\r\n     * So the origins can be different as long as they are on the infinite line (they can\r\n     * \"float\") but the directions must be parallel or antiparallel.\r\n     */\r\n    if (!this.direction.isParallelTo(other.direction, true))\r\n      return false;\r\n    /**\r\n     * In exact math, we consider a ray to have an infinite line as direction (not a finite vector).\r\n     * Therefore, in exact math it is not possible for one origin to be on the other ray but not vice\r\n     * versa. However, we test both ways because first check may pass due to round-off errors.\r\n     */\r\n    let workPoint = this.projectPointToRay(other.origin);\r\n    if (!other.origin.isAlmostEqualMetric(workPoint))\r\n      return false;\r\n    workPoint = other.projectPointToRay(this.origin);\r\n    if (!this.origin.isAlmostEqualMetric(workPoint))\r\n      return false;\r\n    return true;\r\n  }\r\n  /** Create a ray from origin and direction. */\r\n  public static create(origin: Point3d, direction: Vector3d, result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.set(origin, direction);\r\n      return result;\r\n    }\r\n    return new Ray3d(origin.clone(), direction.clone());\r\n  }\r\n  /**\r\n   * Given a homogeneous point and its derivative components, construct a Ray3d with cartesian\r\n   * coordinates and derivatives.\r\n   * @param weightedPoint `[x,y,z,w]` parts of weighted point.\r\n   * @param weightedDerivative `[x,y,z,w]` derivatives\r\n   * @param result\r\n   */\r\n  public static createWeightedDerivative(\r\n    weightedPoint: Float64Array, weightedDerivative: Float64Array, result?: Ray3d\r\n  ): Ray3d | undefined {\r\n    const w = weightedPoint[3];\r\n    const dw = weightedDerivative[3];\r\n    const x = weightedPoint[0];\r\n    const y = weightedPoint[1];\r\n    const z = weightedPoint[2];\r\n    const dx = weightedDerivative[0] * w - weightedPoint[0] * dw;\r\n    const dy = weightedDerivative[1] * w - weightedPoint[1] * dw;\r\n    const dz = weightedDerivative[2] * w - weightedPoint[2] * dw;\r\n    if (Geometry.isSmallMetricDistance(w))\r\n      return undefined;\r\n    const divW = 1.0 / w;\r\n    const divWW = divW * divW;\r\n    return Ray3d.createXYZUVW(\r\n      x * divW, y * divW, z * divW,\r\n      dx * divWW, dy * divWW, dz * divWW,\r\n      result\r\n    );\r\n  }\r\n  /** Create from coordinates of the origin and direction. */\r\n  public static createXYZUVW(\r\n    originX: number, originY: number, originZ: number,\r\n    directionX: number, directionY: number, directionZ: number,\r\n    result?: Ray3d\r\n  ): Ray3d {\r\n    if (result) {\r\n      result.getOriginRef().set(originX, originY, originZ);\r\n      result.getDirectionRef().set(directionX, directionY, directionZ);\r\n      return result;\r\n    }\r\n    return new Ray3d(Point3d.create(originX, originY, originZ), Vector3d.create(directionX, directionY, directionZ));\r\n  }\r\n  /** Capture origin and direction in a new Ray3d. */\r\n  public static createCapture(origin: Point3d, direction: Vector3d): Ray3d {\r\n    return new Ray3d(origin, direction);\r\n  }\r\n  /** Create from (clones of) origin, direction, and numeric weight. */\r\n  public static createPointVectorNumber(origin: Point3d, direction: Vector3d, a: number, result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.origin.setFrom(origin);\r\n      result.direction.setFrom(direction);\r\n      result.a = a;\r\n      return result;\r\n    }\r\n    result = new Ray3d(origin.clone(), direction.clone());\r\n    result.a = a;\r\n    return result;\r\n  }\r\n  /** Create from origin and target. The direction vector is the full length (non-unit) vector from origin to target. */\r\n  public static createStartEnd(origin: Point3d, target: Point3d, result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.origin.setFrom(origin);\r\n      result.direction.setStartEnd(origin, target);\r\n      return result;\r\n    }\r\n    return new Ray3d(origin.clone(), Vector3d.createStartEnd(origin, target));\r\n  }\r\n  /** Return a reference to the ray's origin. */\r\n  public getOriginRef(): Point3d {\r\n    return this.origin;\r\n  }\r\n  /** Return a reference to the ray's direction vector. */\r\n  public getDirectionRef(): Vector3d {\r\n    return this.direction;\r\n  }\r\n  /** Copy coordinates from origin and direction. */\r\n  public set(origin: Point3d, direction: Vector3d): void {\r\n    this.origin.setFrom(origin);\r\n    this.direction.setFrom(direction);\r\n  }\r\n  /** Clone the ray. */\r\n  public clone(result?: Ray3d): Ray3d {\r\n    if (result) {\r\n      result.set(this.origin.clone(), this.direction.clone());\r\n      return result;\r\n    }\r\n    return new Ray3d(this.origin.clone(), this.direction.clone());\r\n  }\r\n  /** Create a clone and return the transform of the clone. */\r\n  public cloneTransformed(transform: Transform): Ray3d {\r\n    return new Ray3d(transform.multiplyPoint3d(this.origin), transform.multiplyVector(this.direction));\r\n  }\r\n  /** Create a clone and return the inverse transform of the clone. */\r\n  public cloneInverseTransformed(transform: Transform): Ray3d | undefined {\r\n    const origin = transform.multiplyInversePoint3d(this.origin);\r\n    const direction = transform.matrix.multiplyInverseXYZAsVector3d(\r\n      this.direction.x, this.direction.y, this.direction.z\r\n    );\r\n    if (undefined !== origin && undefined !== direction)\r\n      return new Ray3d(origin, direction);\r\n    return undefined;\r\n  }\r\n  /** Apply a transform in place. */\r\n  public transformInPlace(transform: Transform) {\r\n    transform.multiplyPoint3d(this.origin, this.origin);\r\n    transform.multiplyVector(this.direction, this.direction);\r\n  }\r\n  /** Copy data from another ray. */\r\n  public setFrom(source: Ray3d): void {\r\n    this.set(source.origin, source.direction);\r\n  }\r\n  /**\r\n   * Return a point at fractional position along the ray.\r\n   * * fraction 0 is the ray origin.\r\n   * * fraction 1 is at the end of the direction vector when placed at the origin.\r\n   */\r\n  public fractionToPoint(fraction: number, result?: Point3d): Point3d {\r\n    return this.origin.plusScaled(this.direction, fraction, result);\r\n  }\r\n  /**\r\n   * Return a transform for rigid axes at ray origin with z in ray direction.\r\n   * * If the direction vector is zero, axes default to identity (from [[Matrix3d.createRigidHeadsUp]])\r\n   */\r\n  public toRigidZFrame(): Transform | undefined {\r\n    const axes = Matrix3d.createRigidHeadsUp(this.direction, AxisOrder.ZXY);\r\n    return Transform.createOriginAndMatrix(this.origin, axes);\r\n  }\r\n  /** Convert {origin:[x,y,z], direction:[u,v,w]} to a Ray3d. */\r\n  public setFromJSON(json?: any) {\r\n    if (!json) {\r\n      this.origin.set(0, 0, 0);\r\n      this.direction.set(0, 0, 1);\r\n      return;\r\n    }\r\n    this.origin.setFromJSON(json.origin);\r\n    this.direction.setFromJSON(json.direction);\r\n  }\r\n  /**\r\n   * Construct a JSON object from this Ray3d.\r\n   * @return {*} [origin,normal]\r\n   */\r\n  public toJSON(): any {\r\n    return { origin: this.origin.toJSON(), direction: this.direction.toJSON() };\r\n  }\r\n  /** Create a new ray from json object. See `setFromJSON` for json structure; */\r\n  public static fromJSON(json?: any) {\r\n    const result = Ray3d.createXAxis();\r\n    result.setFromJSON(json);\r\n    return result;\r\n  }\r\n  /**\r\n   * Try to scale the direction vector to a given `magnitude`.\r\n   * * Returns `false` if the ray direction is a zero vector.\r\n   */\r\n  public trySetDirectionMagnitudeInPlace(magnitude: number = 1.0): boolean {\r\n    if (this.direction.tryNormalizeInPlace()) {\r\n      this.direction.scaleInPlace(magnitude);\r\n      return true;\r\n    }\r\n    this.direction.setZero();\r\n    this.a = 0.0;\r\n    return false;\r\n  }\r\n  /**\r\n   * Normalize the ray direction in place.\r\n   * * If parameter `a` is clearly nonzero and the direction vector can be normalized,\r\n   *    * Save the parameter `a` as the optional `a` member of the ray.\r\n   *    * Normalize the ray's direction vector.\r\n   * * If parameter `a` is nearly zero,\r\n   *    * Set the `a` member to zero.\r\n   *    * Set the ray's direction vector to zero.\r\n   * @param a value to be saved (e.g,. area).\r\n   * @returns `true` if `a` is nonzero and normalization was successful. Otherwise, return `false`.\r\n   */\r\n  public tryNormalizeInPlaceWithAreaWeight(a: number): boolean {\r\n    const tolerance = Geometry.smallMetricDistanceSquared;\r\n    this.a = a;\r\n    if (Math.abs(a) > tolerance && this.direction.tryNormalizeInPlace(tolerance))\r\n      return true;\r\n    this.direction.setZero();\r\n    this.a = 0.0;\r\n    return false;\r\n  }\r\n  /** Return distance from the ray to point in space. */\r\n  public distance(spacePoint: Point3d): number {\r\n    const uu = this.direction.magnitudeSquared();\r\n    const uv = this.dotProductToPoint(spacePoint);\r\n    const aa = Geometry.inverseMetricDistanceSquared(uu);\r\n    if (aa)\r\n      return Math.sqrt(this.origin.distanceSquared(spacePoint) - uv * uv * aa);\r\n    else\r\n      return Math.sqrt(this.origin.distanceSquared(spacePoint));\r\n  }\r\n  /**\r\n   * Return the intersection parameter of the line defined by the ray with a `plane`.\r\n   * * Stores the point of intersection in the `result` point (if passed as a parameter) and returns the parameter\r\n   * along the ray where the intersection occurs. If we call the parameter 'f' then the point of intersection would\r\n   * be `ray.origin + f * ray.direction`. Therefore:\r\n   *    * if ray intersects the plane at its origin, the function returns f = 0.\r\n   *    * if intersects at `ray.origin + ray.direction`, the function returns f = 1.\r\n   *    * if intersects behind the ray origin, the function returns f < 0.\r\n   *    * if intersects after `ray.origin + ray.direction`, the function returns f > 1.\r\n   * * Returns `undefined` if the ray and plane are parallel or coplanar.\r\n   */\r\n  public intersectionWithPlane(plane: Plane3dByOriginAndUnitNormal, result?: Point3d): number | undefined {\r\n    const vectorA = Vector3d.createStartEnd(plane.getOriginRef(), this.origin);\r\n    const uDotN = this.direction.dotProduct(plane.getNormalRef());\r\n    const nDotN = this.direction.magnitudeSquared();\r\n    const aDotN = vectorA.dotProduct(plane.getNormalRef());\r\n    const division = Geometry.conditionalDivideFraction(-aDotN, uDotN);\r\n    if (undefined === division)\r\n      return undefined;\r\n    const division1 = Geometry.conditionalDivideFraction(nDotN, uDotN);\r\n    if (undefined === division1)\r\n      return undefined;\r\n    if (result) {\r\n      this.origin.plusScaled(this.direction, division, result);\r\n    }\r\n    return division;\r\n  }\r\n  /**\r\n   * Find the intersection of the line defined by the ray with a Range3d.\r\n   * * Return the range of parameters (on the ray) which are \"inside\" the range.\r\n   * * Note that a range is always returned; if there is no intersection it is indicated by the test `result.isNull`.\r\n   */\r\n  public intersectionWithRange3d(range: Range3d, result?: Range1d): Range1d {\r\n    if (range.isNull)\r\n      return Range1d.createNull(result);\r\n    const interval = Range1d.createXX(-Geometry.largeCoordinateResult, Geometry.largeCoordinateResult, result);\r\n    if (interval.clipLinearMapToInterval(this.origin.x, this.direction.x, range.low.x, range.high.x)\r\n      && interval.clipLinearMapToInterval(this.origin.y, this.direction.y, range.low.y, range.high.y)\r\n      && interval.clipLinearMapToInterval(this.origin.z, this.direction.z, range.low.z, range.high.z)\r\n    )\r\n      return interval;\r\n    return interval;\r\n  }\r\n  /**\r\n   * Return the shortest vector `v` to `targetPoint` from the line defined by this ray.\r\n   * * If the projection of `targetPoint` onto the line defined by this ray is q, then `v  = targetPoint - q`.\r\n   */\r\n  public perpendicularPartOfVectorToTarget(targetPoint: XYAndZ, result?: Vector3d): Vector3d {\r\n    const vectorV = Vector3d.createStartEnd(this.origin, targetPoint);\r\n    const uu = this.direction.magnitudeSquared();\r\n    const uv = this.direction.dotProductStartEnd(this.origin, targetPoint);\r\n    const fraction = Geometry.safeDivideFraction(uv, uu, 0.0);\r\n    return vectorV.plusScaled(this.direction, -fraction, result);\r\n  }\r\n  /**\r\n   * Determine if two rays intersect, or are fully overlapped, or parallel but not coincident, or skew.\r\n   * * Return a CurveLocationDetailPair which contains fraction and point on each ray and has\r\n   * annotation (in member `approachType`) indicating one of these relationships:\r\n   *   * CurveCurveApproachType.Intersection -- the rays have a simple intersection, at fractions indicated\r\n   * in detailA and detailB\r\n   *   * CurveCurveApproachType.PerpendicularChord -- there is pair of where the rays have closest approach.\r\n   * The rays are skew in space.\r\n   *   * CurveCurveApproachType.CoincidentGeometry -- the rays are the same unbounded line in space. The\r\n   * fractions and points are a representative single common point.\r\n   *   * CurveCurveApproachType.Parallel -- the rays are parallel (and not coincident). The two points are\r\n   * at the minimum distance\r\n   */\r\n  public static closestApproachRay3dRay3d(rayA: Ray3d, rayB: Ray3d): CurveLocationDetailPair {\r\n    const intersectionFractions = Vector2d.create();\r\n    let fractionA, fractionB;\r\n    let pointA, pointB;\r\n    let pairType;\r\n    if (SmallSystem.ray3dXYZUVWClosestApproachUnbounded(\r\n      rayA.origin.x, rayA.origin.y, rayA.origin.z, rayA.direction.x, rayA.direction.y, rayA.direction.z,\r\n      rayB.origin.x, rayB.origin.y, rayB.origin.z, rayB.direction.x, rayB.direction.y, rayB.direction.z,\r\n      intersectionFractions)\r\n    ) {\r\n      fractionA = intersectionFractions.x;\r\n      fractionB = intersectionFractions.y;\r\n      pointA = rayA.fractionToPoint(fractionA);\r\n      pointB = rayB.fractionToPoint(fractionB);\r\n      pairType = pointA.isAlmostEqualMetric(pointB) ?\r\n        CurveCurveApproachType.Intersection : CurveCurveApproachType.PerpendicularChord;\r\n    } else {\r\n      fractionB = 0.0;\r\n      fractionA = rayA.pointToFraction(rayB.origin);\r\n      pointA = rayA.fractionToPoint(fractionA);\r\n      pointB = rayB.fractionToPoint(fractionB);\r\n      pairType = pointA.isAlmostEqualMetric(pointB) ?\r\n        CurveCurveApproachType.CoincidentGeometry : CurveCurveApproachType.ParallelGeometry;\r\n    }\r\n    const pair = CurveLocationDetailPair.createCapture(\r\n      CurveLocationDetail.createRayFractionPoint(rayA, fractionA, rayA.fractionToPoint(fractionA)),\r\n      CurveLocationDetail.createRayFractionPoint(rayB, fractionB, rayB.fractionToPoint(fractionB)));\r\n    pair.approachType = pairType;\r\n    return pair;\r\n  }\r\n  /**\r\n   * Return a ray with `ray.origin` interpolated between `pt1` and `pt2` at the given `fraction`\r\n   * and `ray.direction` set to the vector from `pt1` to `pt2` multiplied by the given `tangentScale`.\r\n   * @param pt1 start point of the interpolation.\r\n   * @param fraction fractional position between points.\r\n   * @param pt2 end point of the interpolation.\r\n   * @param tangentScale scale factor to apply to the startToEnd vector.\r\n   * @param result optional receiver.\r\n   */\r\n  public static interpolatePointAndTangent(\r\n    pt1: XYAndZ, fraction: number, pt2: XYAndZ, tangentScale: number, result?: Ray3d\r\n  ): Ray3d {\r\n    result = result ?? Ray3d.createZero();\r\n    const dx = pt2.x - pt1.x;\r\n    const dy = pt2.y - pt1.y;\r\n    const dz = pt2.z - pt1.z;\r\n    result.direction.set(tangentScale * dx, tangentScale * dy, tangentScale * dz);\r\n    if (fraction <= 0.5)\r\n      result.origin.set(pt1.x + fraction * dx, pt1.y + fraction * dy, pt1.z + fraction * dz);\r\n    else {\r\n      const t: number = fraction - 1.0;\r\n      result.origin.set(pt2.x + t * dx, pt2.y + t * dy, pt2.z + t * dz);\r\n    }\r\n    return result;\r\n  }\r\n}\r\n"]}

package/lib/esm/curve/Query/PlanarSubdivision.d.ts

@@ -19,7 +19,13 @@ export declare class PlanarSubdivision {
      * @param point1 end point
      */
     private static addHalfEdge;
-    static collectSignedLoop(loop: Loop, signedAreas: SignedLoops, zeroAreaTolerance?: number): number;
+    /**
+     * Based on computed (and toleranced) area, push the loop (pointer) onto the appropriate array of positive, negative, or sliver loops.
+     * @param zeroAreaTolerance absolute area tolerance for sliver face detection
+     * @param isSliverFace whether the loop is known a priori (e.g., via topology) to have zero area
+     * @returns the area (forced to zero if within tolerance)
+     */
+    static collectSignedLoop(loop: Loop, outLoops: SignedLoops, zeroAreaTolerance?: number, isSliverFace?: boolean): number;
     static createLoopInFace(faceSeed: HalfEdge, announce?: (he: HalfEdge, curve: CurvePrimitive, loop: Loop) => void): Loop;
     private static isNullFace;
     private static nonNullEdgeMate;

package/lib/esm/curve/Query/PlanarSubdivision.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"PlanarSubdivision.d.ts","sourceRoot":"","sources":["../../../../src/curve/Query/PlanarSubdivision.ts"],"names":[],"mappings":"AAMA,OAAO,EAAE,QAAQ,EAAE,aAAa,EAAE,MAAM,sBAAsB,CAAC;AAG/D,OAAO,EAAuB,uBAAuB,EAAE,MAAM,wBAAwB,CAAC;AACtF,OAAO,EAAE,cAAc,EAAE,MAAM,mBAAmB,CAAC;AACnD,OAAO,EAAE,IAAI,EAAsB,WAAW,EAAE,MAAM,SAAS,CAAC;AAqDhE;;GAEG;AACH,qBAAa,iBAAiB;IAC5B,qGAAqG;WACvF,qBAAqB,CAAC,UAAU,EAAE,cAAc,EAAE,EAAE,QAAQ,EAAE,uBAAuB,EAAE,GAAG,aAAa;IA+BvH;;;;;;;;;OASG;IACD,OAAO,CAAC,MAAM,CAAC,WAAW;WAeZ,iBAAiB,CAAC,IAAI,EAAE,IAAI,EAAE,WAAW,EAAE,WAAW,EAAE,iBAAiB,GAAE,MAAgB,GAAG,MAAM;WAepG,gBAAgB,CAAC,QAAQ,EAAE,QAAQ,EAC/C,QAAQ,CAAC,EAAE,CAAC,EAAE,EAAE,QAAQ,EAAE,KAAK,EAAE,cAAc,EAAE,IAAI,EAAE,IAAI,KAAK,IAAI,GAAG,IAAI;IAsB7E,OAAO,CAAC,MAAM,CAAC,UAAU;IAYzB,OAAO,CAAC,MAAM,CAAC,eAAe;WAWhB,oCAAoC,CAAC,KAAK,EAAE,aAAa,EAAE,iBAAiB,GAAE,MAAgB,GAAG,WAAW,EAAE;CAmC7H"}
+{"version":3,"file":"PlanarSubdivision.d.ts","sourceRoot":"","sources":["../../../../src/curve/Query/PlanarSubdivision.ts"],"names":[],"mappings":"AAMA,OAAO,EAAE,QAAQ,EAAE,aAAa,EAAE,MAAM,sBAAsB,CAAC;AAG/D,OAAO,EAAuB,uBAAuB,EAAE,MAAM,wBAAwB,CAAC;AACtF,OAAO,EAAE,cAAc,EAAE,MAAM,mBAAmB,CAAC;AACnD,OAAO,EAAE,IAAI,EAAsB,WAAW,EAAE,MAAM,SAAS,CAAC;AAqDhE;;GAEG;AACH,qBAAa,iBAAiB;IAC5B,qGAAqG;WACvF,qBAAqB,CAAC,UAAU,EAAE,cAAc,EAAE,EAAE,QAAQ,EAAE,uBAAuB,EAAE,GAAG,aAAa;IA+BrH;;;;;;;;;OASG;IACH,OAAO,CAAC,MAAM,CAAC,WAAW;IAa1B;;;;;OAKG;WACW,iBAAiB,CAAC,IAAI,EAAE,IAAI,EAAE,QAAQ,EAAE,WAAW,EAAE,iBAAiB,GAAE,MAAgB,EAAE,YAAY,CAAC,EAAE,OAAO,GAAG,MAAM;WAezH,gBAAgB,CAAC,QAAQ,EAAE,QAAQ,EAC/C,QAAQ,CAAC,EAAE,CAAC,EAAE,EAAE,QAAQ,EAAE,KAAK,EAAE,cAAc,EAAE,IAAI,EAAE,IAAI,KAAK,IAAI,GAAG,IAAI;IAsB7E,OAAO,CAAC,MAAM,CAAC,UAAU;IAYzB,OAAO,CAAC,MAAM,CAAC,eAAe;WAWhB,oCAAoC,CAAC,KAAK,EAAE,aAAa,EAAE,iBAAiB,GAAE,MAAgB,GAAG,WAAW,EAAE;CAkC7H"}

package/lib/esm/curve/Query/PlanarSubdivision.js

@@ -120,21 +120,25 @@ export class PlanarSubdivision {
             mate.sortAngle = sortAngle(detail01.curve, detail01.fraction1, true);
         }
     }
-    // based on computed (and toleranced) area, push the loop (pointer) onto the appropriate array of positive, negative, or sliver loops.
-    // return the area (forced to zero if within tolerance)
-    static collectSignedLoop(loop, signedAreas, zeroAreaTolerance = 1.0e-10) {
-        let area = RegionOps.computeXYArea(loop);
+    /**
+     * Based on computed (and toleranced) area, push the loop (pointer) onto the appropriate array of positive, negative, or sliver loops.
+     * @param zeroAreaTolerance absolute area tolerance for sliver face detection
+     * @param isSliverFace whether the loop is known a priori (e.g., via topology) to have zero area
+     * @returns the area (forced to zero if within tolerance)
+     */
+    static collectSignedLoop(loop, outLoops, zeroAreaTolerance = 1.0e-10, isSliverFace) {
+        let area = isSliverFace ? 0.0 : RegionOps.computeXYArea(loop);
         if (area === undefined)
             area = 0;
         if (Math.abs(area) < zeroAreaTolerance)
             area = 0.0;
         loop.computedAreaInPlanarSubdivision = area;
         if (area > 0)
-            signedAreas.positiveAreaLoops.push(loop);
+            outLoops.positiveAreaLoops.push(loop);
         else if (area < 0)
-            signedAreas.negativeAreaLoops.push(loop);
+            outLoops.negativeAreaLoops.push(loop);
         else
-            signedAreas.slivers.push(loop);
+            outLoops.slivers.push(loop);
         return area;
     }
     static createLoopInFace(faceSeed, announce) {
@@ -190,11 +194,9 @@ export class PlanarSubdivision {
             const componentAreas = { positiveAreaLoops: [], negativeAreaLoops: [], slivers: [] };
             const edges = [];
             for (const faceSeed of faceSeeds) {
+                const isNullFace = this.isNullFace(faceSeed);
                 const loop = this.createLoopInFace(faceSeed, (he, curveC, loopC) => {
-                    if (this.isNullFace(he)) {
-                        // Ignore all edges of null faces.
-                    }
-                    else {
+                    if (!isNullFace) {
                         const mate = this.nonNullEdgeMate(graph, he);
                         if (mate !== undefined) {
                             const e = edgeMap.get(mate);
@@ -211,7 +213,7 @@ export class PlanarSubdivision {
                         }
                     }
                 });
-                this.collectSignedLoop(loop, componentAreas, zeroAreaTolerance);
+                this.collectSignedLoop(loop, componentAreas, zeroAreaTolerance, isNullFace);
             }
             componentAreas.edges = edges;
             result.push(componentAreas);

package/lib/esm/curve/Query/PlanarSubdivision.js.map

@@ -1 +1 @@
-{"version":3,"file":"PlanarSubdivision.js","sourceRoot":"","sources":["../../../../src/curve/Query/PlanarSubdivision.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F,OAAO,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AAE1C,OAAO,EAAY,aAAa,EAAE,MAAM,sBAAsB,CAAC;AAC/D,OAAO,EAAE,mBAAmB,EAAE,MAAM,oCAAoC,CAAC;AACzE,OAAO,EAAE,kBAAkB,EAAE,MAAM,wBAAwB,CAAC;AAC5D,OAAO,EAAE,mBAAmB,EAA2B,MAAM,wBAAwB,CAAC;AAEtF,OAAO,EAAE,IAAI,EAAE,kBAAkB,EAAe,MAAM,SAAS,CAAC;AAChE,OAAO,EAAE,SAAS,EAAE,MAAM,cAAc,CAAC;AAEzC;;GAEG;AAEH,MAAM,+CAA+C;IAArD;QACS,oBAAe,GAAG,IAAI,GAAG,EAA6C,CAAC;QAC9E,yDAAyD;QAClD,qBAAgB,GAAG,IAAI,GAAG,EAA0B,CAAC;QACpD,0BAAqB,GAAW,CAAC,CAAC;IA4B5C,CAAC;IA3BQ,oBAAoB,CAAC,SAAqC;QAC/D,IAAI,SAAS,KAAK,SAAS,EAAE;YAC3B,MAAM,KAAK,GAAG,IAAI,CAAC,gBAAgB,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;YACnD,IAAI,KAAK,KAAK,SAAS,EAAE;gBACvB,IAAI,CAAC,gBAAgB,CAAC,GAAG,CAAC,SAAS,EAAE,IAAI,CAAC,qBAAqB,EAAE,CAAC,CAAC;aACpE;SACF;IACH,CAAC;IACM,qBAAqB,CAAC,SAAqC,EAAE,IAA6B;QAC/F,IAAI,SAAS,EAAE;YACb,MAAM,UAAU,GAAG,IAAI,CAAC,eAAe,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;YACvD,IAAI,CAAC,oBAAoB,CAAC,SAAS,CAAC,CAAC;YACrC,IAAI,CAAC,UAAU,EAAE;gBACf,IAAI,CAAC,eAAe,CAAC,GAAG,CAAC,SAAS,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC;aAC7C;iBAAM;gBACL,UAAU,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;aACvB;SACF;IACH,CAAC;IACM,UAAU,CAAC,IAA6B;QAC7C,MAAM,UAAU,GAAG,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC;QACtC,IAAI,UAAU;YACZ,IAAI,CAAC,qBAAqB,CAAC,UAAU,EAAE,IAAI,CAAC,CAAC;QAC/C,MAAM,UAAU,GAAG,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC;QACtC,IAAI,UAAU;YACZ,IAAI,CAAC,qBAAqB,CAAC,UAAU,EAAE,IAAI,CAAC,CAAC;IACjD,CAAC;CACF;AACD;;;;;;;;;;;;EAYE;AACF;;GAEG;AACH,MAAM,OAAO,iBAAiB;IAC5B,qGAAqG;IAC9F,MAAM,CAAC,qBAAqB,CAAC,UAA4B,EAAE,QAAmC;QACnG,MAAM,iBAAiB,GAAG,IAAI,+CAA+C,EAAE,CAAC,CAAG,0DAA0D;QAC7I,KAAK,MAAM,CAAC,IAAI,UAAU;YACxB,iBAAiB,CAAC,oBAAoB,CAAC,CAAC,CAAC,CAAC;QAC5C,KAAK,MAAM,IAAI,IAAI,QAAQ,EAAE;YAC3B,iBAAiB,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC;SACpC;QACD,MAAM,KAAK,GAAG,IAAI,aAAa,EAAE,CAAC;QAClC,KAAK,MAAM,KAAK,IAAI,iBAAiB,CAAC,eAAe,CAAC,OAAO,EAAE,EAAE;YAC/D,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;YACnB,MAAM,OAAO,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;YACzB,OAAO,CAAC,IAAI,CAAC,CAAC,KAA8B,EAAE,KAA8B,EAAE,EAAE;gBAC9E,MAAM,SAAS,GAAG,kBAAkB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;gBAC/C,MAAM,SAAS,GAAG,kBAAkB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;gBAC/C,IAAI,SAAS,KAAK,SAAS,IAAI,SAAS,KAAK,SAAS;oBACpD,OAAO,CAAC,MAAM,CAAC;gBACjB,OAAO,SAAS,GAAG,SAAS,CAAC;YAC/B,CAAC,CAAC,CAAC;YACH,IAAI,OAAO,GAAG,gBAAgB,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAE,CAAC;YAC/C,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,CAAC,EAAE,CAAC,CAAC,UAAU,EAAG,EAAE,GAAG,EAAE,OAAO,CAAC,KAAK,EAAE,OAAO,CAAC,QAAQ,CAAC,CAAC;YAClF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBACvC,4CAA4C;gBAC5C,MAAM,OAAO,GAAG,gBAAgB,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAE,CAAC;gBACjD,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,CAAC,EAAE,OAAO,CAAC,KAAK,EAAE,OAAO,CAAC,QAAQ,EAAE,OAAO,CAAC,KAAK,EAAE,OAAO,CAAC,QAAQ,CAAC,CAAC;gBAC7F,OAAO,GAAG,OAAO,CAAC;aACnB;YACD,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,CAAC,EAAE,OAAO,CAAC,KAAK,EAAE,OAAO,CAAC,QAAQ,EAAE,CAAC,CAAC,QAAQ,EAAE,EAAE,GAAG,CAAC,CAAC;SAChF;QACD,kBAAkB,CAAC,sBAAsB,CAAC,KAAK,EAAE,CAAC,EAAY,EAAE,EAAE,CAAC,EAAE,CAAC,SAAU,CAAC,CAAC;QAClF,OAAO,KAAK,CAAC;IACf,CAAC;IACH;;;;;;;;;OASG;IACO,MAAM,CAAC,WAAW,CAAC,KAAoB,EAAE,CAAiB,EAAE,MAAe,EAAE,SAAiB,EAAE,MAAe,EAAE,SAAiB;QACxI,IAAI,CAAC,MAAM,CAAC,aAAa,CAAE,MAAM,CAAC,EAAC;YACjC,MAAM,QAAQ,GAAG,KAAK,CAAC,gBAAgB,CAAC,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC;YAC9D,MAAM,QAAQ,GAAG,mBAAmB,CAAC,oCAAoC,CAAC,CAAC,EAAE,SAAS,EAAE,SAAS,CAAC,CAAC;YACnG,MAAM,IAAI,GAAG,QAAQ,CAAC,QAAQ,CAAC;YAC/B,QAAQ,CAAC,OAAO,GAAG,QAAQ,CAAC;YAC5B,QAAQ,CAAC,QAAQ,GAAG,GAAG,CAAC;YACxB,IAAI,CAAC,OAAO,GAAG,QAAQ,CAAC;YACxB,IAAI,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC;YACrB,QAAQ,CAAC,SAAS,GAAG,SAAS,CAAC,QAAQ,CAAC,KAAM,EAAE,QAAQ,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;YAC1E,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC,QAAQ,CAAC,KAAM,EAAE,QAAQ,CAAC,SAAU,EAAE,IAAI,CAAC,CAAC;SACtE;IACH,CAAC;IACL,sIAAsI;IACpI,uDAAuD;IAChD,MAAM,CAAC,iBAAiB,CAAC,IAAU,EAAE,WAAwB,EAAE,oBAA4B,OAAO;QACvG,IAAI,IAAI,GAAG,SAAS,CAAC,aAAa,CAAC,IAAI,CAAC,CAAC;QACzC,IAAI,IAAI,KAAK,SAAS;YACpB,IAAI,GAAG,CAAC,CAAC;QACX,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,iBAAiB;YACpC,IAAI,GAAG,GAAG,CAAC;QACZ,IAAY,CAAC,+BAA+B,GAAG,IAAI,CAAC;QACrD,IAAI,IAAI,GAAG,CAAC;YACV,WAAW,CAAC,iBAAiB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;aACtC,IAAI,IAAI,GAAG,CAAC;YACf,WAAW,CAAC,iBAAiB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;;YAEzC,WAAW,CAAC,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QACjC,OAAO,IAAI,CAAC;IACd,CAAC;IACM,MAAM,CAAC,gBAAgB,CAAC,QAAkB,EAC/C,QAAoE;QACpE,IAAI,EAAE,GAAG,QAAQ,CAAC;QAClB,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC;QAC3B,GAAG;YACD,MAAM,MAAM,GAAG,EAAE,CAAC,OAA8B,CAAC;YACjD,IAAI,MAAM,EAAE;gBACV,IAAI,KAAK,CAAC;gBACV,IAAI,EAAE,CAAC,QAAS,GAAG,CAAC;oBAClB,KAAK,GAAG,MAAM,CAAC,KAAM,CAAC,iBAAiB,CAAC,MAAM,CAAC,QAAQ,EAAE,MAAM,CAAC,SAAU,CAAC,CAAC;;oBAE5E,KAAK,GAAG,MAAM,CAAC,KAAM,CAAC,iBAAiB,CAAC,MAAM,CAAC,SAAU,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC;gBAC9E,IAAI,KAAK,EAAE;oBACT,IAAI,QAAQ,KAAK,SAAS;wBACxB,QAAQ,CAAC,EAAE,EAAE,KAAK,EAAE,IAAI,CAAC,CAAC;oBAC5B,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC;iBACzB;aACF;YACD,EAAE,GAAG,EAAE,CAAC,aAAa,CAAC;SACvB,QAAQ,EAAE,KAAK,QAAQ,EAAE;QAC1B,OAAO,IAAI,CAAC;IACd,CAAC;IACD,qGAAqG;IAC7F,MAAM,CAAC,UAAU,CAAC,EAAY;QACpC,MAAM,eAAe,GAAG,CAAC,EAAE,CAAC,aAAa,CAAC,aAAa,KAAK,EAAE,CAAC,CAAC;QAChE,IAAI,YAAY,GAAG,KAAK,CAAC;QACzB,IAAI,eAAe,EAAE;YACnB,MAAM,EAAE,GAAG,kBAAkB,CAAC,gBAAgB,CAAC,EAAE,CAAC,CAAC;YACnD,MAAM,EAAE,GAAG,kBAAkB,CAAC,gBAAgB,CAAC,EAAE,CAAC,aAAa,CAAC,QAAQ,CAAC,CAAC;YAC1E,IAAI,CAAC,QAAQ,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,eAAe;gBACrD,YAAY,GAAG,IAAI,CAAC,CAAE,8EAA8E;SACvG;QACD,OAAO,eAAe,IAAI,CAAC,YAAY,CAAC;IAC1C,CAAC;IACD,qEAAqE;IAC7D,MAAM,CAAC,eAAe,CAAC,MAAqB,EAAE,CAAW;QAC/D,IAAI,IAAI,CAAC,UAAU,CAAE,CAAC,CAAC;YACrB,OAAO,SAAS,CAAC;QACnB,IAAI,EAAE,GAAG,CAAC,CAAC,QAAQ,CAAC;QACpB,OAAO,IAAI,CAAC,UAAU,CAAC,EAAE,CAAC,EAAC;YACzB,EAAE,GAAG,EAAE,CAAC,aAAa,CAAC,QAAQ,CAAC;YAC/B,IAAI,EAAE,KAAK,CAAC;gBACV,OAAO,SAAS,CAAC;SACpB;QACD,OAAO,EAAE,CAAC;IACZ,CAAC;IACM,MAAM,CAAC,oCAAoC,CAAC,KAAoB,EAAE,oBAA4B,OAAO;QAC1G,MAAM,CAAC,GAAG,mBAAmB,CAAC,iDAAiD,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;QAClG,MAAM,MAAM,GAAkB,EAAE,CAAC;QACjC,MAAM,OAAO,GAAG,IAAI,GAAG,EAAgC,CAAC;QACxD,KAAK,MAAM,SAAS,IAAI,CAAC,EAAE;YACzB,MAAM,cAAc,GAAgB,EAAE,iBAAiB,EAAE,EAAE,EAAE,iBAAiB,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,CAAC;YAClG,MAAM,KAAK,GAAyB,EAAE,CAAC;YACvC,KAAK,MAAM,QAAQ,IAAI,SAAS,EAAE;gBAChC,MAAM,IAAI,GAAG,IAAI,CAAC,gBAAgB,CAAC,QAAQ,EAAE,CAAC,EAAY,EAAE,MAAsB,EAAE,KAAW,EAAE,EAAE;oBACjG,IAAI,IAAI,CAAC,UAAU,CAAC,EAAE,CAAC,EAAE;wBACvB,kCAAkC;qBACnC;yBAAM;wBACL,MAAM,IAAI,GAAG,IAAI,CAAC,eAAe,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;wBAC3C,IAAI,IAAI,KAAK,SAAS,EAAC;4BACrB,MAAM,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;4BAC5B,IAAI,CAAC,KAAK,SAAS,EAAE;gCACnB,oGAAoG;gCACpG,MAAM,EAAE,GAAG,IAAI,kBAAkB,CAAC,KAAK,EAAE,MAAM,EAAE,SAAS,EAAE,SAAS,CAAC,CAAC;gCACvE,OAAO,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;6BACrB;iCAAM,IAAI,CAAC,YAAY,kBAAkB,EAAE;gCAC1C,CAAC,CAAC,IAAI,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC;gCACtB,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;gCACd,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;6BACtB;yBACJ;qBACF;gBACH,CAAC,CAAC,CAAC;gBACH,IAAI,CAAC,iBAAiB,CAAC,IAAI,EAAE,cAAc,EAAE,iBAAiB,CAAC,CAAC;aACjE;YACD,cAAc,CAAC,KAAK,GAAG,KAAK,CAAC;YAC7B,MAAM,CAAC,IAAI,CAAC,cAAc,CAAC,CAAC;YAC5B,OAAO,CAAC,KAAK,EAAE,CAAC;SACjB;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;CACF;AAED,SAAS,SAAS,CAAC,KAAqB,EAAE,QAAgB,EAAE,OAAgB;IAC1E,MAAM,GAAG,GAAG,KAAK,CAAC,4BAA4B,CAAC,QAAQ,CAAC,CAAC;IACzD,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;IAC/B,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,CAAC,SAAS,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;AAC9D,CAAC;AAED,SAAS,kBAAkB,CAAC,IAA6B,EAAE,KAAqB;IAC9E,IAAI,IAAI,CAAC,OAAO,CAAC,KAAK,KAAK,KAAK;QAC9B,OAAO,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC;IAC/B,IAAI,IAAI,CAAC,OAAO,CAAC,KAAK,KAAK,KAAK;QAC9B,OAAO,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC;IAC/B,OAAO,SAAS,CAAC;AACnB,CAAC;AACD,SAAS,gBAAgB,CAAC,IAA6B,EAAE,KAAqB;IAC5E,IAAI,IAAI,CAAC,OAAO,CAAC,KAAK,KAAK,KAAK;QAC9B,OAAO,IAAI,CAAC,OAAO,CAAC;IACtB,IAAI,IAAI,CAAC,OAAO,CAAC,KAAK,KAAK,KAAK;QAC9B,OAAO,IAAI,CAAC,OAAO,CAAC;IACtB,OAAO,SAAS,CAAC;AACnB,CAAC","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\nimport { Geometry } from \"../../Geometry\";\r\nimport { Point3d } from \"../../geometry3d/Point3dVector3d\";\r\nimport { HalfEdge, HalfEdgeGraph } from \"../../topology/Graph\";\r\nimport { HalfEdgeGraphSearch } from \"../../topology/HalfEdgeGraphSearch\";\r\nimport { HalfEdgeGraphMerge } from \"../../topology/Merging\";\r\nimport { CurveLocationDetail, CurveLocationDetailPair } from \"../CurveLocationDetail\";\r\nimport { CurvePrimitive } from \"../CurvePrimitive\";\r\nimport { Loop, LoopCurveLoopCurve, SignedLoops } from \"../Loop\";\r\nimport { RegionOps } from \"../RegionOps\";\r\n\r\n/** @packageDocumentation\r\n * @module Curve\r\n */\r\n\r\nclass MapCurvePrimitiveToCurveLocationDetailPairArray {\r\n  public primitiveToPair = new Map<CurvePrimitive, CurveLocationDetailPair[]>();\r\n  // index assigned to this primitive for this calculation.\r\n  public primitiveToIndex = new Map<CurvePrimitive, number>();\r\n  private _numIndexedPrimitives: number = 0;\r\n  public assignPrimitiveIndex(primitive: CurvePrimitive | undefined) {\r\n    if (primitive !== undefined) {\r\n      const index = this.primitiveToIndex.get(primitive);\r\n      if (index === undefined) {\r\n        this.primitiveToIndex.set(primitive, this._numIndexedPrimitives++);\r\n      }\r\n    }\r\n  }\r\n  public insertPrimitiveToPair(primitive: CurvePrimitive | undefined, pair: CurveLocationDetailPair) {\r\n    if (primitive) {\r\n      const priorPairs = this.primitiveToPair.get(primitive);\r\n      this.assignPrimitiveIndex(primitive);\r\n      if (!priorPairs) {\r\n        this.primitiveToPair.set(primitive, [pair]);\r\n      } else {\r\n        priorPairs.push(pair);\r\n      }\r\n    }\r\n  }\r\n  public insertPair(pair: CurveLocationDetailPair) {\r\n    const primitiveA = pair.detailA.curve;\r\n    if (primitiveA)\r\n      this.insertPrimitiveToPair(primitiveA, pair);\r\n    const primitiveB = pair.detailB.curve;\r\n    if (primitiveB)\r\n      this.insertPrimitiveToPair(primitiveB, pair);\r\n  }\r\n}\r\n/*\r\n  function getDetailString(detail: CurveLocationDetail | undefined): string {\r\n    if (!detail)\r\n      return \"{}\";\r\n    else return tagString(\"primitive\", this.primitiveToIndex.get(detail.curve!)) + tagString(\"f0\", detail.fraction) + tagString(\"f1\", detail.fraction1);\r\n  }\r\n}\r\nfunction tagString(name: string, value: number | undefined): string {\r\n  if (value !== undefined)\r\n    return \"(\" + name + \" \" + value + \")\";\r\n  return \"\";\r\n}\r\n*/\r\n/**\r\n * @internal\r\n */\r\nexport class PlanarSubdivision {\r\n  /** Create a graph from an array of curves, and an array of the curves' precomputed intersections. */\r\n  public static assembleHalfEdgeGraph(primitives: CurvePrimitive[], allPairs: CurveLocationDetailPair[]): HalfEdgeGraph {\r\n    const detailByPrimitive = new MapCurvePrimitiveToCurveLocationDetailPairArray();   // map from key CurvePrimitive to CurveLocationDetailPair.\r\n    for (const p of primitives)\r\n      detailByPrimitive.assignPrimitiveIndex(p);\r\n    for (const pair of allPairs) {\r\n      detailByPrimitive.insertPair(pair);\r\n    }\r\n    const graph = new HalfEdgeGraph();\r\n    for (const entry of detailByPrimitive.primitiveToPair.entries()) {\r\n      const p = entry[0];\r\n      const details = entry[1];\r\n      details.sort((pairA: CurveLocationDetailPair, pairB: CurveLocationDetailPair) => {\r\n        const fractionA = getFractionOnCurve(pairA, p);\r\n        const fractionB = getFractionOnCurve(pairB, p);\r\n        if (fractionA === undefined || fractionB === undefined)\r\n          return -1000.0;\r\n        return fractionA - fractionB;\r\n      });\r\n      let detail0 = getDetailOnCurve(details[0], p)!;\r\n      this.addHalfEdge(graph, p, p.startPoint (), 0.0, detail0.point, detail0.fraction);\r\n      for (let i = 1; i < details.length; i++) {\r\n        // create (both sides of) a graph edge . . .\r\n        const detail1 = getDetailOnCurve(details[i], p)!;\r\n        this.addHalfEdge(graph, p, detail0.point, detail0.fraction, detail1.point, detail1.fraction);\r\n        detail0 = detail1;\r\n      }\r\n      this.addHalfEdge(graph, p, detail0.point, detail0.fraction, p.endPoint(), 1.0);\r\n    }\r\n    HalfEdgeGraphMerge.clusterAndMergeXYTheta(graph, (he: HalfEdge) => he.sortAngle!);\r\n    return graph;\r\n  }\r\n/**\r\n * Create a pair of mated half edges referencing an interval of a primitive\r\n *   * no action if start and end points are identical.\r\n * @param graph containing graph.\r\n * @param p the curve\r\n * @param fraction0 starting fraction\r\n * @param point0 start point\r\n * @param fraction1 end fraction\r\n * @param point1 end point\r\n */\r\n  private static addHalfEdge(graph: HalfEdgeGraph, p: CurvePrimitive, point0: Point3d, fraction0: number, point1: Point3d, fraction1: number) {\r\n    if (!point0.isAlmostEqual (point1)){\r\n      const halfEdge = graph.createEdgeXYAndZ(point0, 0, point1, 0);\r\n      const detail01 = CurveLocationDetail.createCurveEvaluatedFractionFraction(p, fraction0, fraction1);\r\n      const mate = halfEdge.edgeMate;\r\n      halfEdge.edgeTag = detail01;\r\n      halfEdge.sortData = 1.0;\r\n      mate.edgeTag = detail01;\r\n      mate.sortData = -1.0;\r\n      halfEdge.sortAngle = sortAngle(detail01.curve!, detail01.fraction, false);\r\n      mate.sortAngle = sortAngle(detail01.curve!, detail01.fraction1!, true);\r\n      }\r\n    }\r\n// based on computed (and toleranced) area, push the loop (pointer) onto the appropriate array of positive, negative, or sliver loops.\r\n  // return the area (forced to zero if within tolerance)\r\n  public static collectSignedLoop(loop: Loop, signedAreas: SignedLoops, zeroAreaTolerance: number = 1.0e-10): number{\r\n    let area = RegionOps.computeXYArea(loop);\r\n    if (area === undefined)\r\n      area = 0;\r\n    if (Math.abs(area) < zeroAreaTolerance)\r\n      area = 0.0;\r\n    (loop as any).computedAreaInPlanarSubdivision = area;\r\n    if (area > 0)\r\n      signedAreas.positiveAreaLoops.push(loop);\r\n    else if (area < 0)\r\n      signedAreas.negativeAreaLoops.push(loop);\r\n    else\r\n      signedAreas.slivers.push(loop);\r\n    return area;\r\n  }\r\n  public static createLoopInFace(faceSeed: HalfEdge,\r\n    announce?: (he: HalfEdge, curve: CurvePrimitive, loop: Loop) => void): Loop {\r\n    let he = faceSeed;\r\n    const loop = Loop.create();\r\n    do {\r\n      const detail = he.edgeTag as CurveLocationDetail;\r\n      if (detail) {\r\n        let curve;\r\n        if (he.sortData! > 0)\r\n          curve = detail.curve!.clonePartialCurve(detail.fraction, detail.fraction1!);\r\n        else\r\n          curve = detail.curve!.clonePartialCurve(detail.fraction1!, detail.fraction);\r\n        if (curve) {\r\n          if (announce !== undefined)\r\n            announce(he, curve, loop);\r\n          loop.tryAddChild(curve);\r\n        }\r\n      }\r\n      he = he.faceSuccessor;\r\n    } while (he !== faceSeed);\r\n    return loop;\r\n  }\r\n  // Return true if there are only two edges in the face loop, and their start curvatures are the same.\r\n  private static isNullFace(he: HalfEdge): boolean {\r\n    const faceHasTwoEdges = (he.faceSuccessor.faceSuccessor === he);\r\n    let faceIsBanana = false;\r\n    if (faceHasTwoEdges) {\r\n      const c0 = HalfEdgeGraphMerge.curvatureSortKey(he);\r\n      const c1 = HalfEdgeGraphMerge.curvatureSortKey(he.faceSuccessor.edgeMate);\r\n      if (!Geometry.isSameCoordinate(c0, c1)) // default tol!\r\n        faceIsBanana = true;  // heuristic: we could also check end curvatures, and/or higher derivatives...\r\n    }\r\n    return faceHasTwoEdges && !faceIsBanana;\r\n  }\r\n  // Look across edge mates (possibly several) for a nonnull mate face.\r\n  private static nonNullEdgeMate(_graph: HalfEdgeGraph, e: HalfEdge): HalfEdge | undefined {\r\n    if (this.isNullFace (e))\r\n      return undefined;\r\n    let e1 = e.edgeMate;\r\n    while (this.isNullFace(e1)){\r\n      e1 = e1.faceSuccessor.edgeMate;\r\n      if (e1 === e)\r\n        return undefined;\r\n    }\r\n    return e1;\r\n  }\r\n  public static collectSignedLoopSetsInHalfEdgeGraph(graph: HalfEdgeGraph, zeroAreaTolerance: number = 1.0e-10): SignedLoops[] {\r\n    const q = HalfEdgeGraphSearch.collectConnectedComponentsWithExteriorParityMasks(graph, undefined);\r\n    const result: SignedLoops[] = [];\r\n    const edgeMap = new Map<HalfEdge, LoopCurveLoopCurve>();\r\n    for (const faceSeeds of q) {\r\n      const componentAreas: SignedLoops = { positiveAreaLoops: [], negativeAreaLoops: [], slivers: [] };\r\n      const edges: LoopCurveLoopCurve[] = [];\r\n      for (const faceSeed of faceSeeds) {\r\n        const loop = this.createLoopInFace(faceSeed, (he: HalfEdge, curveC: CurvePrimitive, loopC: Loop) => {\r\n          if (this.isNullFace(he)) {\r\n            // Ignore all edges of null faces.\r\n          } else {\r\n            const mate = this.nonNullEdgeMate(graph, he);\r\n              if (mate !== undefined){\r\n                const e = edgeMap.get(mate);\r\n                if (e === undefined) {\r\n                  // Record this as loopA,edgeA of a shared edge to be completed later from the other side of the edge\r\n                  const e1 = new LoopCurveLoopCurve(loopC, curveC, undefined, undefined);\r\n                  edgeMap.set(he, e1);\r\n                } else if (e instanceof LoopCurveLoopCurve) {\r\n                  e.setB(loopC, curveC);\r\n                  edges.push(e);\r\n                  edgeMap.delete(mate);\r\n                }\r\n            }\r\n          }\r\n        });\r\n        this.collectSignedLoop(loop, componentAreas, zeroAreaTolerance);\r\n      }\r\n      componentAreas.edges = edges;\r\n      result.push(componentAreas);\r\n      edgeMap.clear();\r\n    }\r\n    return result;\r\n  }\r\n}\r\n\r\nfunction sortAngle(curve: CurvePrimitive, fraction: number, reverse: boolean): number {\r\n  const ray = curve.fractionToPointAndDerivative(fraction);\r\n  const s = reverse ? -1.0 : 1.0;\r\n  return Math.atan2(s * ray.direction.y, s * ray.direction.x);\r\n}\r\n\r\nfunction getFractionOnCurve(pair: CurveLocationDetailPair, curve: CurvePrimitive): number | undefined {\r\n  if (pair.detailA.curve === curve)\r\n    return pair.detailA.fraction;\r\n  if (pair.detailB.curve === curve)\r\n    return pair.detailB.fraction;\r\n  return undefined;\r\n}\r\nfunction getDetailOnCurve(pair: CurveLocationDetailPair, curve: CurvePrimitive): CurveLocationDetail | undefined {\r\n  if (pair.detailA.curve === curve)\r\n    return pair.detailA;\r\n  if (pair.detailB.curve === curve)\r\n    return pair.detailB;\r\n  return undefined;\r\n}\r\n"]}
+{"version":3,"file":"PlanarSubdivision.js","sourceRoot":"","sources":["../../../../src/curve/Query/PlanarSubdivision.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F,OAAO,EAAE,QAAQ,EAAE,MAAM,gBAAgB,CAAC;AAE1C,OAAO,EAAY,aAAa,EAAE,MAAM,sBAAsB,CAAC;AAC/D,OAAO,EAAE,mBAAmB,EAAE,MAAM,oCAAoC,CAAC;AACzE,OAAO,EAAE,kBAAkB,EAAE,MAAM,wBAAwB,CAAC;AAC5D,OAAO,EAAE,mBAAmB,EAA2B,MAAM,wBAAwB,CAAC;AAEtF,OAAO,EAAE,IAAI,EAAE,kBAAkB,EAAe,MAAM,SAAS,CAAC;AAChE,OAAO,EAAE,SAAS,EAAE,MAAM,cAAc,CAAC;AAEzC;;GAEG;AAEH,MAAM,+CAA+C;IAArD;QACS,oBAAe,GAAG,IAAI,GAAG,EAA6C,CAAC;QAC9E,yDAAyD;QAClD,qBAAgB,GAAG,IAAI,GAAG,EAA0B,CAAC;QACpD,0BAAqB,GAAW,CAAC,CAAC;IA4B5C,CAAC;IA3BQ,oBAAoB,CAAC,SAAqC;QAC/D,IAAI,SAAS,KAAK,SAAS,EAAE;YAC3B,MAAM,KAAK,GAAG,IAAI,CAAC,gBAAgB,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;YACnD,IAAI,KAAK,KAAK,SAAS,EAAE;gBACvB,IAAI,CAAC,gBAAgB,CAAC,GAAG,CAAC,SAAS,EAAE,IAAI,CAAC,qBAAqB,EAAE,CAAC,CAAC;aACpE;SACF;IACH,CAAC;IACM,qBAAqB,CAAC,SAAqC,EAAE,IAA6B;QAC/F,IAAI,SAAS,EAAE;YACb,MAAM,UAAU,GAAG,IAAI,CAAC,eAAe,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;YACvD,IAAI,CAAC,oBAAoB,CAAC,SAAS,CAAC,CAAC;YACrC,IAAI,CAAC,UAAU,EAAE;gBACf,IAAI,CAAC,eAAe,CAAC,GAAG,CAAC,SAAS,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC;aAC7C;iBAAM;gBACL,UAAU,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;aACvB;SACF;IACH,CAAC;IACM,UAAU,CAAC,IAA6B;QAC7C,MAAM,UAAU,GAAG,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC;QACtC,IAAI,UAAU;YACZ,IAAI,CAAC,qBAAqB,CAAC,UAAU,EAAE,IAAI,CAAC,CAAC;QAC/C,MAAM,UAAU,GAAG,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC;QACtC,IAAI,UAAU;YACZ,IAAI,CAAC,qBAAqB,CAAC,UAAU,EAAE,IAAI,CAAC,CAAC;IACjD,CAAC;CACF;AACD;;;;;;;;;;;;EAYE;AACF;;GAEG;AACH,MAAM,OAAO,iBAAiB;IAC5B,qGAAqG;IAC9F,MAAM,CAAC,qBAAqB,CAAC,UAA4B,EAAE,QAAmC;QACnG,MAAM,iBAAiB,GAAG,IAAI,+CAA+C,EAAE,CAAC,CAAG,0DAA0D;QAC7I,KAAK,MAAM,CAAC,IAAI,UAAU;YACxB,iBAAiB,CAAC,oBAAoB,CAAC,CAAC,CAAC,CAAC;QAC5C,KAAK,MAAM,IAAI,IAAI,QAAQ,EAAE;YAC3B,iBAAiB,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC;SACpC;QACD,MAAM,KAAK,GAAG,IAAI,aAAa,EAAE,CAAC;QAClC,KAAK,MAAM,KAAK,IAAI,iBAAiB,CAAC,eAAe,CAAC,OAAO,EAAE,EAAE;YAC/D,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;YACnB,MAAM,OAAO,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;YACzB,OAAO,CAAC,IAAI,CAAC,CAAC,KAA8B,EAAE,KAA8B,EAAE,EAAE;gBAC9E,MAAM,SAAS,GAAG,kBAAkB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;gBAC/C,MAAM,SAAS,GAAG,kBAAkB,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;gBAC/C,IAAI,SAAS,KAAK,SAAS,IAAI,SAAS,KAAK,SAAS;oBACpD,OAAO,CAAC,MAAM,CAAC;gBACjB,OAAO,SAAS,GAAG,SAAS,CAAC;YAC/B,CAAC,CAAC,CAAC;YACH,IAAI,OAAO,GAAG,gBAAgB,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAE,CAAC;YAC/C,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,CAAC,EAAE,CAAC,CAAC,UAAU,EAAG,EAAE,GAAG,EAAE,OAAO,CAAC,KAAK,EAAE,OAAO,CAAC,QAAQ,CAAC,CAAC;YAClF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBACvC,4CAA4C;gBAC5C,MAAM,OAAO,GAAG,gBAAgB,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAE,CAAC;gBACjD,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,CAAC,EAAE,OAAO,CAAC,KAAK,EAAE,OAAO,CAAC,QAAQ,EAAE,OAAO,CAAC,KAAK,EAAE,OAAO,CAAC,QAAQ,CAAC,CAAC;gBAC7F,OAAO,GAAG,OAAO,CAAC;aACnB;YACD,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,CAAC,EAAE,OAAO,CAAC,KAAK,EAAE,OAAO,CAAC,QAAQ,EAAE,CAAC,CAAC,QAAQ,EAAE,EAAE,GAAG,CAAC,CAAC;SAChF;QACD,kBAAkB,CAAC,sBAAsB,CAAC,KAAK,EAAE,CAAC,EAAY,EAAE,EAAE,CAAC,EAAE,CAAC,SAAU,CAAC,CAAC;QAClF,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;;;;;;;OASG;IACK,MAAM,CAAC,WAAW,CAAC,KAAoB,EAAE,CAAiB,EAAE,MAAe,EAAE,SAAiB,EAAE,MAAe,EAAE,SAAiB;QACxI,IAAI,CAAC,MAAM,CAAC,aAAa,CAAE,MAAM,CAAC,EAAC;YACjC,MAAM,QAAQ,GAAG,KAAK,CAAC,gBAAgB,CAAC,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC;YAC9D,MAAM,QAAQ,GAAG,mBAAmB,CAAC,oCAAoC,CAAC,CAAC,EAAE,SAAS,EAAE,SAAS,CAAC,CAAC;YACnG,MAAM,IAAI,GAAG,QAAQ,CAAC,QAAQ,CAAC;YAC/B,QAAQ,CAAC,OAAO,GAAG,QAAQ,CAAC;YAC5B,QAAQ,CAAC,QAAQ,GAAG,GAAG,CAAC;YACxB,IAAI,CAAC,OAAO,GAAG,QAAQ,CAAC;YACxB,IAAI,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC;YACrB,QAAQ,CAAC,SAAS,GAAG,SAAS,CAAC,QAAQ,CAAC,KAAM,EAAE,QAAQ,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;YAC1E,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC,QAAQ,CAAC,KAAM,EAAE,QAAQ,CAAC,SAAU,EAAE,IAAI,CAAC,CAAC;SACtE;IACH,CAAC;IACH;;;;;OAKG;IACI,MAAM,CAAC,iBAAiB,CAAC,IAAU,EAAE,QAAqB,EAAE,oBAA4B,OAAO,EAAE,YAAsB;QAC5H,IAAI,IAAI,GAAG,YAAY,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,SAAS,CAAC,aAAa,CAAC,IAAI,CAAC,CAAC;QAC9D,IAAI,IAAI,KAAK,SAAS;YACpB,IAAI,GAAG,CAAC,CAAC;QACX,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,iBAAiB;YACpC,IAAI,GAAG,GAAG,CAAC;QACZ,IAAY,CAAC,+BAA+B,GAAG,IAAI,CAAC;QACrD,IAAI,IAAI,GAAG,CAAC;YACV,QAAQ,CAAC,iBAAiB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;aACnC,IAAI,IAAI,GAAG,CAAC;YACf,QAAQ,CAAC,iBAAiB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;;YAEtC,QAAQ,CAAC,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAC9B,OAAO,IAAI,CAAC;IACd,CAAC;IACM,MAAM,CAAC,gBAAgB,CAAC,QAAkB,EAC/C,QAAoE;QACpE,IAAI,EAAE,GAAG,QAAQ,CAAC;QAClB,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC;QAC3B,GAAG;YACD,MAAM,MAAM,GAAG,EAAE,CAAC,OAA8B,CAAC;YACjD,IAAI,MAAM,EAAE;gBACV,IAAI,KAAK,CAAC;gBACV,IAAI,EAAE,CAAC,QAAS,GAAG,CAAC;oBAClB,KAAK,GAAG,MAAM,CAAC,KAAM,CAAC,iBAAiB,CAAC,MAAM,CAAC,QAAQ,EAAE,MAAM,CAAC,SAAU,CAAC,CAAC;;oBAE5E,KAAK,GAAG,MAAM,CAAC,KAAM,CAAC,iBAAiB,CAAC,MAAM,CAAC,SAAU,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC;gBAC9E,IAAI,KAAK,EAAE;oBACT,IAAI,QAAQ,KAAK,SAAS;wBACxB,QAAQ,CAAC,EAAE,EAAE,KAAK,EAAE,IAAI,CAAC,CAAC;oBAC5B,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC;iBACzB;aACF;YACD,EAAE,GAAG,EAAE,CAAC,aAAa,CAAC;SACvB,QAAQ,EAAE,KAAK,QAAQ,EAAE;QAC1B,OAAO,IAAI,CAAC;IACd,CAAC;IACD,qGAAqG;IAC7F,MAAM,CAAC,UAAU,CAAC,EAAY;QACpC,MAAM,eAAe,GAAG,CAAC,EAAE,CAAC,aAAa,CAAC,aAAa,KAAK,EAAE,CAAC,CAAC;QAChE,IAAI,YAAY,GAAG,KAAK,CAAC;QACzB,IAAI,eAAe,EAAE;YACnB,MAAM,EAAE,GAAG,kBAAkB,CAAC,gBAAgB,CAAC,EAAE,CAAC,CAAC;YACnD,MAAM,EAAE,GAAG,kBAAkB,CAAC,gBAAgB,CAAC,EAAE,CAAC,aAAa,CAAC,QAAQ,CAAC,CAAC;YAC1E,IAAI,CAAC,QAAQ,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,eAAe;gBACrD,YAAY,GAAG,IAAI,CAAC,CAAE,8EAA8E;SACvG;QACD,OAAO,eAAe,IAAI,CAAC,YAAY,CAAC;IAC1C,CAAC;IACD,qEAAqE;IAC7D,MAAM,CAAC,eAAe,CAAC,MAAqB,EAAE,CAAW;QAC/D,IAAI,IAAI,CAAC,UAAU,CAAE,CAAC,CAAC;YACrB,OAAO,SAAS,CAAC;QACnB,IAAI,EAAE,GAAG,CAAC,CAAC,QAAQ,CAAC;QACpB,OAAO,IAAI,CAAC,UAAU,CAAC,EAAE,CAAC,EAAC;YACzB,EAAE,GAAG,EAAE,CAAC,aAAa,CAAC,QAAQ,CAAC;YAC/B,IAAI,EAAE,KAAK,CAAC;gBACV,OAAO,SAAS,CAAC;SACpB;QACD,OAAO,EAAE,CAAC;IACZ,CAAC;IACM,MAAM,CAAC,oCAAoC,CAAC,KAAoB,EAAE,oBAA4B,OAAO;QAC1G,MAAM,CAAC,GAAG,mBAAmB,CAAC,iDAAiD,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;QAClG,MAAM,MAAM,GAAkB,EAAE,CAAC;QACjC,MAAM,OAAO,GAAG,IAAI,GAAG,EAAgC,CAAC;QACxD,KAAK,MAAM,SAAS,IAAI,CAAC,EAAE;YACzB,MAAM,cAAc,GAAgB,EAAE,iBAAiB,EAAE,EAAE,EAAE,iBAAiB,EAAE,EAAE,EAAE,OAAO,EAAE,EAAE,EAAE,CAAC;YAClG,MAAM,KAAK,GAAyB,EAAE,CAAC;YACvC,KAAK,MAAM,QAAQ,IAAI,SAAS,EAAE;gBAChC,MAAM,UAAU,GAAG,IAAI,CAAC,UAAU,CAAC,QAAQ,CAAC,CAAC;gBAC7C,MAAM,IAAI,GAAG,IAAI,CAAC,gBAAgB,CAAC,QAAQ,EAAE,CAAC,EAAY,EAAE,MAAsB,EAAE,KAAW,EAAE,EAAE;oBACjG,IAAI,CAAC,UAAU,EAAE;wBACf,MAAM,IAAI,GAAG,IAAI,CAAC,eAAe,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;wBAC7C,IAAI,IAAI,KAAK,SAAS,EAAE;4BACtB,MAAM,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;4BAC5B,IAAI,CAAC,KAAK,SAAS,EAAE;gCACnB,oGAAoG;gCACpG,MAAM,EAAE,GAAG,IAAI,kBAAkB,CAAC,KAAK,EAAE,MAAM,EAAE,SAAS,EAAE,SAAS,CAAC,CAAC;gCACvE,OAAO,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;6BACrB;iCAAM,IAAI,CAAC,YAAY,kBAAkB,EAAE;gCAC1C,CAAC,CAAC,IAAI,CAAC,KAAK,EAAE,MAAM,CAAC,CAAC;gCACtB,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;gCACd,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;6BACtB;yBACF;qBACF;gBACH,CAAC,CAAC,CAAC;gBACH,IAAI,CAAC,iBAAiB,CAAC,IAAI,EAAE,cAAc,EAAE,iBAAiB,EAAE,UAAU,CAAC,CAAC;aAC7E;YACD,cAAc,CAAC,KAAK,GAAG,KAAK,CAAC;YAC7B,MAAM,CAAC,IAAI,CAAC,cAAc,CAAC,CAAC;YAC5B,OAAO,CAAC,KAAK,EAAE,CAAC;SACjB;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;CACF;AAED,SAAS,SAAS,CAAC,KAAqB,EAAE,QAAgB,EAAE,OAAgB;IAC1E,MAAM,GAAG,GAAG,KAAK,CAAC,4BAA4B,CAAC,QAAQ,CAAC,CAAC;IACzD,MAAM,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;IAC/B,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,CAAC,SAAS,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;AAC9D,CAAC;AAED,SAAS,kBAAkB,CAAC,IAA6B,EAAE,KAAqB;IAC9E,IAAI,IAAI,CAAC,OAAO,CAAC,KAAK,KAAK,KAAK;QAC9B,OAAO,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC;IAC/B,IAAI,IAAI,CAAC,OAAO,CAAC,KAAK,KAAK,KAAK;QAC9B,OAAO,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC;IAC/B,OAAO,SAAS,CAAC;AACnB,CAAC;AACD,SAAS,gBAAgB,CAAC,IAA6B,EAAE,KAAqB;IAC5E,IAAI,IAAI,CAAC,OAAO,CAAC,KAAK,KAAK,KAAK;QAC9B,OAAO,IAAI,CAAC,OAAO,CAAC;IACtB,IAAI,IAAI,CAAC,OAAO,CAAC,KAAK,KAAK,KAAK;QAC9B,OAAO,IAAI,CAAC,OAAO,CAAC;IACtB,OAAO,SAAS,CAAC;AACnB,CAAC","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\nimport { Geometry } from \"../../Geometry\";\r\nimport { Point3d } from \"../../geometry3d/Point3dVector3d\";\r\nimport { HalfEdge, HalfEdgeGraph } from \"../../topology/Graph\";\r\nimport { HalfEdgeGraphSearch } from \"../../topology/HalfEdgeGraphSearch\";\r\nimport { HalfEdgeGraphMerge } from \"../../topology/Merging\";\r\nimport { CurveLocationDetail, CurveLocationDetailPair } from \"../CurveLocationDetail\";\r\nimport { CurvePrimitive } from \"../CurvePrimitive\";\r\nimport { Loop, LoopCurveLoopCurve, SignedLoops } from \"../Loop\";\r\nimport { RegionOps } from \"../RegionOps\";\r\n\r\n/** @packageDocumentation\r\n * @module Curve\r\n */\r\n\r\nclass MapCurvePrimitiveToCurveLocationDetailPairArray {\r\n  public primitiveToPair = new Map<CurvePrimitive, CurveLocationDetailPair[]>();\r\n  // index assigned to this primitive for this calculation.\r\n  public primitiveToIndex = new Map<CurvePrimitive, number>();\r\n  private _numIndexedPrimitives: number = 0;\r\n  public assignPrimitiveIndex(primitive: CurvePrimitive | undefined) {\r\n    if (primitive !== undefined) {\r\n      const index = this.primitiveToIndex.get(primitive);\r\n      if (index === undefined) {\r\n        this.primitiveToIndex.set(primitive, this._numIndexedPrimitives++);\r\n      }\r\n    }\r\n  }\r\n  public insertPrimitiveToPair(primitive: CurvePrimitive | undefined, pair: CurveLocationDetailPair) {\r\n    if (primitive) {\r\n      const priorPairs = this.primitiveToPair.get(primitive);\r\n      this.assignPrimitiveIndex(primitive);\r\n      if (!priorPairs) {\r\n        this.primitiveToPair.set(primitive, [pair]);\r\n      } else {\r\n        priorPairs.push(pair);\r\n      }\r\n    }\r\n  }\r\n  public insertPair(pair: CurveLocationDetailPair) {\r\n    const primitiveA = pair.detailA.curve;\r\n    if (primitiveA)\r\n      this.insertPrimitiveToPair(primitiveA, pair);\r\n    const primitiveB = pair.detailB.curve;\r\n    if (primitiveB)\r\n      this.insertPrimitiveToPair(primitiveB, pair);\r\n  }\r\n}\r\n/*\r\n  function getDetailString(detail: CurveLocationDetail | undefined): string {\r\n    if (!detail)\r\n      return \"{}\";\r\n    else return tagString(\"primitive\", this.primitiveToIndex.get(detail.curve!)) + tagString(\"f0\", detail.fraction) + tagString(\"f1\", detail.fraction1);\r\n  }\r\n}\r\nfunction tagString(name: string, value: number | undefined): string {\r\n  if (value !== undefined)\r\n    return \"(\" + name + \" \" + value + \")\";\r\n  return \"\";\r\n}\r\n*/\r\n/**\r\n * @internal\r\n */\r\nexport class PlanarSubdivision {\r\n  /** Create a graph from an array of curves, and an array of the curves' precomputed intersections. */\r\n  public static assembleHalfEdgeGraph(primitives: CurvePrimitive[], allPairs: CurveLocationDetailPair[]): HalfEdgeGraph {\r\n    const detailByPrimitive = new MapCurvePrimitiveToCurveLocationDetailPairArray();   // map from key CurvePrimitive to CurveLocationDetailPair.\r\n    for (const p of primitives)\r\n      detailByPrimitive.assignPrimitiveIndex(p);\r\n    for (const pair of allPairs) {\r\n      detailByPrimitive.insertPair(pair);\r\n    }\r\n    const graph = new HalfEdgeGraph();\r\n    for (const entry of detailByPrimitive.primitiveToPair.entries()) {\r\n      const p = entry[0];\r\n      const details = entry[1];\r\n      details.sort((pairA: CurveLocationDetailPair, pairB: CurveLocationDetailPair) => {\r\n        const fractionA = getFractionOnCurve(pairA, p);\r\n        const fractionB = getFractionOnCurve(pairB, p);\r\n        if (fractionA === undefined || fractionB === undefined)\r\n          return -1000.0;\r\n        return fractionA - fractionB;\r\n      });\r\n      let detail0 = getDetailOnCurve(details[0], p)!;\r\n      this.addHalfEdge(graph, p, p.startPoint (), 0.0, detail0.point, detail0.fraction);\r\n      for (let i = 1; i < details.length; i++) {\r\n        // create (both sides of) a graph edge . . .\r\n        const detail1 = getDetailOnCurve(details[i], p)!;\r\n        this.addHalfEdge(graph, p, detail0.point, detail0.fraction, detail1.point, detail1.fraction);\r\n        detail0 = detail1;\r\n      }\r\n      this.addHalfEdge(graph, p, detail0.point, detail0.fraction, p.endPoint(), 1.0);\r\n    }\r\n    HalfEdgeGraphMerge.clusterAndMergeXYTheta(graph, (he: HalfEdge) => he.sortAngle!);\r\n    return graph;\r\n  }\r\n  /**\r\n   * Create a pair of mated half edges referencing an interval of a primitive\r\n   *   * no action if start and end points are identical.\r\n   * @param graph containing graph.\r\n   * @param p the curve\r\n   * @param fraction0 starting fraction\r\n   * @param point0 start point\r\n   * @param fraction1 end fraction\r\n   * @param point1 end point\r\n   */\r\n  private static addHalfEdge(graph: HalfEdgeGraph, p: CurvePrimitive, point0: Point3d, fraction0: number, point1: Point3d, fraction1: number) {\r\n    if (!point0.isAlmostEqual (point1)){\r\n      const halfEdge = graph.createEdgeXYAndZ(point0, 0, point1, 0);\r\n      const detail01 = CurveLocationDetail.createCurveEvaluatedFractionFraction(p, fraction0, fraction1);\r\n      const mate = halfEdge.edgeMate;\r\n      halfEdge.edgeTag = detail01;\r\n      halfEdge.sortData = 1.0;\r\n      mate.edgeTag = detail01;\r\n      mate.sortData = -1.0;\r\n      halfEdge.sortAngle = sortAngle(detail01.curve!, detail01.fraction, false);\r\n      mate.sortAngle = sortAngle(detail01.curve!, detail01.fraction1!, true);\r\n      }\r\n    }\r\n  /**\r\n   * Based on computed (and toleranced) area, push the loop (pointer) onto the appropriate array of positive, negative, or sliver loops.\r\n   * @param zeroAreaTolerance absolute area tolerance for sliver face detection\r\n   * @param isSliverFace whether the loop is known a priori (e.g., via topology) to have zero area\r\n   * @returns the area (forced to zero if within tolerance)\r\n   */\r\n  public static collectSignedLoop(loop: Loop, outLoops: SignedLoops, zeroAreaTolerance: number = 1.0e-10, isSliverFace?: boolean): number{\r\n    let area = isSliverFace ? 0.0 : RegionOps.computeXYArea(loop);\r\n    if (area === undefined)\r\n      area = 0;\r\n    if (Math.abs(area) < zeroAreaTolerance)\r\n      area = 0.0;\r\n    (loop as any).computedAreaInPlanarSubdivision = area;\r\n    if (area > 0)\r\n      outLoops.positiveAreaLoops.push(loop);\r\n    else if (area < 0)\r\n      outLoops.negativeAreaLoops.push(loop);\r\n    else\r\n      outLoops.slivers.push(loop);\r\n    return area;\r\n  }\r\n  public static createLoopInFace(faceSeed: HalfEdge,\r\n    announce?: (he: HalfEdge, curve: CurvePrimitive, loop: Loop) => void): Loop {\r\n    let he = faceSeed;\r\n    const loop = Loop.create();\r\n    do {\r\n      const detail = he.edgeTag as CurveLocationDetail;\r\n      if (detail) {\r\n        let curve;\r\n        if (he.sortData! > 0)\r\n          curve = detail.curve!.clonePartialCurve(detail.fraction, detail.fraction1!);\r\n        else\r\n          curve = detail.curve!.clonePartialCurve(detail.fraction1!, detail.fraction);\r\n        if (curve) {\r\n          if (announce !== undefined)\r\n            announce(he, curve, loop);\r\n          loop.tryAddChild(curve);\r\n        }\r\n      }\r\n      he = he.faceSuccessor;\r\n    } while (he !== faceSeed);\r\n    return loop;\r\n  }\r\n  // Return true if there are only two edges in the face loop, and their start curvatures are the same.\r\n  private static isNullFace(he: HalfEdge): boolean {\r\n    const faceHasTwoEdges = (he.faceSuccessor.faceSuccessor === he);\r\n    let faceIsBanana = false;\r\n    if (faceHasTwoEdges) {\r\n      const c0 = HalfEdgeGraphMerge.curvatureSortKey(he);\r\n      const c1 = HalfEdgeGraphMerge.curvatureSortKey(he.faceSuccessor.edgeMate);\r\n      if (!Geometry.isSameCoordinate(c0, c1)) // default tol!\r\n        faceIsBanana = true;  // heuristic: we could also check end curvatures, and/or higher derivatives...\r\n    }\r\n    return faceHasTwoEdges && !faceIsBanana;\r\n  }\r\n  // Look across edge mates (possibly several) for a nonnull mate face.\r\n  private static nonNullEdgeMate(_graph: HalfEdgeGraph, e: HalfEdge): HalfEdge | undefined {\r\n    if (this.isNullFace (e))\r\n      return undefined;\r\n    let e1 = e.edgeMate;\r\n    while (this.isNullFace(e1)){\r\n      e1 = e1.faceSuccessor.edgeMate;\r\n      if (e1 === e)\r\n        return undefined;\r\n    }\r\n    return e1;\r\n  }\r\n  public static collectSignedLoopSetsInHalfEdgeGraph(graph: HalfEdgeGraph, zeroAreaTolerance: number = 1.0e-10): SignedLoops[] {\r\n    const q = HalfEdgeGraphSearch.collectConnectedComponentsWithExteriorParityMasks(graph, undefined);\r\n    const result: SignedLoops[] = [];\r\n    const edgeMap = new Map<HalfEdge, LoopCurveLoopCurve>();\r\n    for (const faceSeeds of q) {\r\n      const componentAreas: SignedLoops = { positiveAreaLoops: [], negativeAreaLoops: [], slivers: [] };\r\n      const edges: LoopCurveLoopCurve[] = [];\r\n      for (const faceSeed of faceSeeds) {\r\n        const isNullFace = this.isNullFace(faceSeed);\r\n        const loop = this.createLoopInFace(faceSeed, (he: HalfEdge, curveC: CurvePrimitive, loopC: Loop) => {\r\n          if (!isNullFace) {\r\n            const mate = this.nonNullEdgeMate(graph, he);\r\n            if (mate !== undefined) {\r\n              const e = edgeMap.get(mate);\r\n              if (e === undefined) {\r\n                // Record this as loopA,edgeA of a shared edge to be completed later from the other side of the edge\r\n                const e1 = new LoopCurveLoopCurve(loopC, curveC, undefined, undefined);\r\n                edgeMap.set(he, e1);\r\n              } else if (e instanceof LoopCurveLoopCurve) {\r\n                e.setB(loopC, curveC);\r\n                edges.push(e);\r\n                edgeMap.delete(mate);\r\n              }\r\n            }\r\n          }\r\n        });\r\n        this.collectSignedLoop(loop, componentAreas, zeroAreaTolerance, isNullFace);\r\n      }\r\n      componentAreas.edges = edges;\r\n      result.push(componentAreas);\r\n      edgeMap.clear();\r\n    }\r\n    return result;\r\n  }\r\n}\r\n\r\nfunction sortAngle(curve: CurvePrimitive, fraction: number, reverse: boolean): number {\r\n  const ray = curve.fractionToPointAndDerivative(fraction);\r\n  const s = reverse ? -1.0 : 1.0;\r\n  return Math.atan2(s * ray.direction.y, s * ray.direction.x);\r\n}\r\n\r\nfunction getFractionOnCurve(pair: CurveLocationDetailPair, curve: CurvePrimitive): number | undefined {\r\n  if (pair.detailA.curve === curve)\r\n    return pair.detailA.fraction;\r\n  if (pair.detailB.curve === curve)\r\n    return pair.detailB.fraction;\r\n  return undefined;\r\n}\r\nfunction getDetailOnCurve(pair: CurveLocationDetailPair, curve: CurvePrimitive): CurveLocationDetail | undefined {\r\n  if (pair.detailA.curve === curve)\r\n    return pair.detailA;\r\n  if (pair.detailB.curve === curve)\r\n    return pair.detailB;\r\n  return undefined;\r\n}\r\n"]}

package/lib/esm/curve/RegionOps.d.ts

@@ -291,8 +291,8 @@ export declare class RegionOps {
     /**
      * Find all areas bounded by the unstructured, possibly intersecting curves.
      * * This method performs no merging of nearly coincident edges and vertices, which can lead to unexpected results
-     * given sufficiently imprecise input. Input geometry consisting of regions can be merged for better results: call
-     * [[regionBooleanXY]](regions, undefined, RegionBinaryOpType.Union) and pass the merged region into [[constructAllXYRegionLoops]].
+     * given sufficiently imprecise input. Input geometry consisting of regions can be merged for better results by pre-processing with
+     * [[regionBooleanXY]].
      * @param curvesAndRegions Any collection of curves. Each Loop/ParityRegion/UnionRegion contributes its curve primitives.
      * @returns array of [[SignedLoops]], each entry of which describes the faces in a single connected component:
      *    * `positiveAreaLoops` contains "interior" loops, _including holes in ParityRegion input_. These loops have positive area and counterclockwise orientation.

package/lib/esm/curve/RegionOps.js

@@ -563,8 +563,8 @@ export class RegionOps {
     /**
      * Find all areas bounded by the unstructured, possibly intersecting curves.
      * * This method performs no merging of nearly coincident edges and vertices, which can lead to unexpected results
-     * given sufficiently imprecise input. Input geometry consisting of regions can be merged for better results: call
-     * [[regionBooleanXY]](regions, undefined, RegionBinaryOpType.Union) and pass the merged region into [[constructAllXYRegionLoops]].
+     * given sufficiently imprecise input. Input geometry consisting of regions can be merged for better results by pre-processing with
+     * [[regionBooleanXY]].
      * @param curvesAndRegions Any collection of curves. Each Loop/ParityRegion/UnionRegion contributes its curve primitives.
      * @returns array of [[SignedLoops]], each entry of which describes the faces in a single connected component:
      *    * `positiveAreaLoops` contains "interior" loops, _including holes in ParityRegion input_. These loops have positive area and counterclockwise orientation.