@itwin/core-geometry 4.0.0-dev.46 → 4.0.0-dev.50

package/lib/esm/curve/CurveCurveIntersectXYZ.js.map

@@ -1 +1 @@
-{"version":3,"file":"CurveCurveIntersectXYZ.js","sourceRoot":"","sources":["../../../src/curve/CurveCurveIntersectXYZ.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AAEH,OAAO,EAAE,cAAc,EAAE,kBAAkB,EAAE,MAAM,yBAAyB,CAAC;AAK7E,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AACvC,OAAO,EAAE,mBAAmB,EAAE,MAAM,+BAA+B,CAAC;AACpE,OAAO,EAAE,4BAA4B,EAAE,MAAM,4CAA4C,CAAC;AAC1F,mCAAmC;AACnC,OAAO,EAAE,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AACzD,mDAAmD;AACnD,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AAClE,iDAAiD;AACjD,OAAO,EAAE,WAAW,EAAE,cAAc,EAAE,MAAM,yBAAyB,CAAC;AACtE,mDAAmD;AACnD,uDAAuD;AACvD,qDAAqD;AACrD,OAAO,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC;AAChC,iDAAiD;AACjD,OAAO,EAAE,4BAA4B,EAAE,MAAM,yBAAyB,CAAC;AACvE,OAAO,EAAE,iBAAiB,EAAE,mBAAmB,EAAE,MAAM,uBAAuB,CAAC;AAG/E,OAAO,EAAE,aAAa,EAAE,MAAM,iBAAiB,CAAC;AAChD,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAE9C,mBAAmB;AACnB;;;;;;;;GAQG;AACH,MAAM,OAAO,sBAAuB,SAAQ,mBAAmB;IAS7D;;;;;;OAMG;IACH,YAAmB,UAAyB,EAAE,OAAgB,EAAE,SAAwB,EAAE,OAAgB;QACxG,KAAK,EAAE,CAAC;QACR,+BAA+B;QAC/B,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;QACxB,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;QAC5B,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;QACxB,IAAI,CAAC,YAAY,EAAE,CAAC;IACtB,CAAC;IAjBO,YAAY;QAClB,IAAI,CAAC,QAAQ,GAAG,IAAI,4BAA4B,EAAE,CAAC;IACrD,CAAC;IAgBD;;;;OAIG;IACI,WAAW,CAAC,eAAwB,KAAK;QAC9C,MAAM,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC;QAC7B,IAAI,YAAY;YACd,IAAI,CAAC,YAAY,EAAE,CAAC;QACtB,OAAO,MAAM,CAAC;IAChB,CAAC;IAIO,cAAc,CAAC,OAAgB,EAAE,QAAgB,EAAE,OAAgB;QACzE,IAAI,CAAC,OAAO,IAAI,QAAQ,GAAG,GAAG;YAC5B,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,OAAO,IAAI,QAAQ,GAAG,GAAG;YAC5B,OAAO,KAAK,CAAC;QACf,OAAO,IAAI,CAAC;IACd,CAAC;IAED;;;;OAIG;IACK,6BAA6B,CACnC,cAAsB,EACtB,GAAmB,EACnB,UAAkB,EAClB,UAAkB,EAClB,cAAsB,EAAI,iCAAiC;IAC3D,GAAmB,EACnB,UAAkB,EAClB,UAAkB,EAClB,QAAiB;QAEjB,MAAM,eAAe,GAAG,QAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,cAAc,EAAE,UAAU,CAAC,CAAC;QACrF,MAAM,eAAe,GAAG,QAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,cAAc,EAAE,UAAU,CAAC,CAAC;QACrF,8CAA8C;QAC9C,MAAM,WAAW,GAAG,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,MAAM,CAAC;QAC/C,IAAI,WAAW,GAAG,CAAC,EAAE;YACnB,MAAM,YAAY,GAAG,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC;YACnE,MAAM,YAAY,GAAG,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC;YACnE,IAAI,QAAQ,EAAE;gBACZ,IAAI,QAAQ,CAAC,mBAAmB,CAAC,YAAY,EAAE,eAAe,CAAC,IAAI,QAAQ,CAAC,mBAAmB,CAAC,YAAY,EAAE,eAAe,CAAC;oBAC5H,OAAO;aACV;iBAAM;gBACL,IAAI,QAAQ,CAAC,mBAAmB,CAAC,YAAY,EAAE,eAAe,CAAC,IAAI,QAAQ,CAAC,mBAAmB,CAAC,YAAY,EAAE,eAAe,CAAC;oBAC5H,OAAO;aACV;SACF;QACD,MAAM,MAAM,GAAG,GAAG,CAAC,eAAe,CAAC,eAAe,CAAC,CAAC;QACpD,MAAM,MAAM,GAAG,GAAG,CAAC,eAAe,CAAC,eAAe,CAAC,CAAC;QACpD,IAAI,CAAC,MAAM,CAAC,mBAAmB,CAAC,MAAM,CAAC;YACrC,OAAO;QAET,MAAM,OAAO,GAAG,mBAAmB,CAAC,wBAAwB,CAAC,GAAG,EAC9D,eAAe,EAAE,MAAM,CAAC,CAAC;QAC3B,OAAO,CAAC,eAAe,CAAC,iBAAiB,CAAC,QAAQ,CAAC,CAAC;QACpD,MAAM,OAAO,GAAG,mBAAmB,CAAC,wBAAwB,CAAC,GAAG,EAC9D,eAAe,EAAE,MAAM,CAAC,CAAC;QAC3B,OAAO,CAAC,eAAe,CAAC,iBAAiB,CAAC,QAAQ,CAAC,CAAC;QACpD,IAAI,QAAQ,EAAE;YACZ,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;YAClC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;SACnC;aAAM;YACL,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;YAClC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;SACnC;IACH,CAAC;IACD;;;OAGG;IACK,uBAAuB,CAC7B,GAAmB,EACnB,QAAiB,EACjB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,QAAiB,EACjB,GAAmB,EACnB,QAAiB,EACjB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,QAAiB,EACjB,QAAiB;QAEjB,MAAM,EAAE,GAAG,sBAAsB,CAAC,cAAc,CAAC;QACjD,IAAI,WAAW,CAAC,qCAAqC,CACnD,OAAO,EAAE,OAAO,EAChB,OAAO,EAAE,OAAO,EAAE,EAAE,CAAC;eAClB,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,EAAE,CAAC,CAAC,EAAE,QAAQ,CAAC;eAC7C,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,EAAE,CAAC,CAAC,EAAE,QAAQ,CAAC,EAChD;YACA,IAAI,CAAC,6BAA6B,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,UAAU,EAAE,UAAU,EAAE,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,UAAU,EAAE,UAAU,EAAE,QAAQ,CAAC,CAAC;SACpH;IACH,CAAC;IACD,+DAA+D;IAC/D,4IAA4I;IAC5I,0CAA0C;IAC1C,iCAAiC;IACzB,sBAAsB,CAC5B,GAAmB,EACnB,QAAiB,EACjB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,QAAiB,EACjB,GAAmB,EACnB,QAAiB,EACjB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,QAAiB,EACjB,QAAiB;QAEjB,IAAI,CAAC,uBAAuB,CAC1B,GAAG,EAAE,QAAQ,EAAE,OAAO,EAAE,UAAU,EAAE,OAAO,EAAE,UAAU,EAAE,QAAQ,EACjE,GAAG,EAAE,QAAQ,EAAE,OAAO,EAAE,UAAU,EAAE,OAAO,EAAE,UAAU,EAAE,QAAQ,EACjE,QAAQ,CAAC,CAAC;IACd,CAAC;IACD;;;;;;;;;;;OAWG;IACI,qCAAqC,CAAC,MAAe,EAAE,OAAiB,EAAE,WAAmB,EAAE,OAAiB,EAAE,OAAiB;QACxI,WAAW,GAAG,QAAQ,CAAC,kBAAkB,CAAC,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,EAAE,GAAG,EAAE,GAAG,GAAG,QAAQ,CAAC,aAAa,CAAC,CAAC;QACpG,MAAM,KAAK,GAAG,OAAO,CAAC,gBAAgB,EAAE,CAAC;QACzC,MAAM,KAAK,GAAG,OAAO,CAAC,gBAAgB,EAAE,CAAC;QACzC,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,UAAU,CAAC,OAAO,CAAC,CAAC,CAAC;QACpD,MAAM,KAAK,GAAG,OAAO,CAAC,gBAAgB,CAAC,KAAK,GAAG,KAAK,IAAI,WAAW,GAAG,WAAW,GAAG,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC;QACvH,IAAI,KAAK;YACP,OAAO,4BAA4B,CAAC,MAAM,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;QAC5D,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,wEAAwE;IACxE,4IAA4I;IAC5I,0CAA0C;IAClC,kBAAkB,CACxB,GAAmB,EACnB,QAAiB,EACjB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,QAAiB,EACjB,GAAU,EACV,QAAiB,EACjB,QAAiB,EACjB,QAAiB;QAEjB,MAAM,UAAU,GAAG,QAAQ,CAAC,cAAc,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAC7D,MAAM,KAAK,GAAG,IAAI,CAAC,qCAAqC,CAAC,OAAO,EAAE,UAAU,EAAE,IAAI,EAAE,GAAG,CAAC,mBAAmB,EAAE,GAAG,CAAC,OAAO,CAAC,CAAC;QAC1H,IAAI,KAAK,KAAK,SAAS,EAAE;YACvB,MAAM,UAAU,GAA0B,EAAE,CAAC;YAC7C,GAAG,CAAC,6BAA6B,CAAC,KAAK,EAAE,UAAU,CAAC,CAAC;YACrD,IAAI,YAAY,CAAC;YACjB,IAAI,SAA8B,CAAC;YACnC,KAAK,MAAM,CAAC,IAAI,UAAU,EAAE;gBAC1B,IAAI,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,CAAC,CAAC,QAAQ,EAAE,QAAQ,CAAC,EAAE;oBACvD,YAAY,GAAG,WAAW,CAAC,kCAAkC,CAAC,OAAO,EAAE,OAAO,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC;oBACzF,IAAI,YAAY,KAAK,SAAS,EAAE;wBAC9B,SAAS,GAAG,OAAO,CAAC,WAAW,CAAC,YAAY,EAAE,OAAO,EAAE,SAAS,CAAC,CAAC;wBAClE,IAAI,SAAS,CAAC,mBAAmB,CAAC,CAAC,CAAC,KAAK,CAAC;+BACrC,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,YAAY,EAAE,QAAQ,CAAC,EAAE;4BAC1D,IAAI,CAAC,6BAA6B,CAAC,YAAY,EAAE,GAAG,EAAE,UAAU,EAAE,UAAU,EAC1E,CAAC,CAAC,QAAQ,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAAC,CAAC;yBAEpC;qBACF;iBACF;aACF;SACF;IACH,CAAC;IAED,qCAAqC;IACrC,yEAAyE;IACzE,wEAAwE;IACxE,oDAAoD;IAC5C,qBAAqB,CAC3B,GAAU,EACV,OAAgB,EAChB,GAAU,EACV,OAAgB,EAChB,QAAiB;QAEjB,MAAM,YAAY,GAAG,GAAG,CAAC,sBAAsB,CAAC,GAAG,CAAC,CAAC;QACrD,IAAI,YAAY,KAAK,SAAS,EAAE;YAC9B,MAAM,cAAc,GAAa,EAAE,CAAC;YACpC,MAAM,aAAa,GAAa,EAAE,CAAC;YACnC,cAAc,CAAC,6CAA6C,CAC1D,YAAY,CAAC,MAAM,CAAC,CAAC,EAAE,YAAY,CAAC,MAAM,CAAC,CAAC,EAAE,GAAG,EACjD,YAAY,CAAC,OAAO,CAAC,CAAC,EAAE,YAAY,CAAC,OAAO,CAAC,CAAC,EAAE,GAAG,EACnD,YAAY,CAAC,QAAQ,CAAC,CAAC,EAAE,YAAY,CAAC,QAAQ,CAAC,CAAC,EAAE,GAAG,EACrD,cAAc,EAAE,aAAa,CAAC,CAAC;YACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,cAAc,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBAC9C,MAAM,SAAS,GAAG,GAAG,CAAC,KAAK,CAAC,+BAA+B,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC9E,MAAM,SAAS,GAAG,GAAG,CAAC,KAAK,CAAC,+BAA+B,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC/E,4FAA4F;gBAC5F,IAAI,IAAI,CAAC,cAAc,CAAC,OAAO,EAAE,SAAS,EAAE,OAAO,CAAC,IAAI,IAAI,CAAC,cAAc,CAAC,OAAO,EAAE,SAAS,EAAE,OAAO,CAAC,EAAE;oBACxG,IAAI,CAAC,6BAA6B,CAAC,SAAS,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EACrD,SAAS,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAAC,CAAC;iBACnC;aACF;SACF;IACH,CAAC;IACD,sCAAsC;IACtC,6EAA6E;IAC7E,+BAA+B;IACvB,cAAc,CACpB,GAAU,EACV,OAAgB,EAChB,GAAU,EACV,OAAgB,EAChB,QAAiB;QAEjB,mCAAmC;QACnC,yDAAyD;QACzD,0DAA0D;QAC1D,qDAAqD;QACrD,kEAAkE;QAClE,MAAM,MAAM,GAAG,4BAA4B,CAAC,MAAM,CAAC,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,mBAAmB,CAAC,CAAC;QACxF,MAAM,MAAM,GAAG,4BAA4B,CAAC,MAAM,CAAC,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,mBAAmB,CAAC,CAAC;QACxF,IAAI,MAAM,KAAK,SAAS,IAAI,MAAM,KAAK,SAAS;YAC9C,OAAO;QAET,IAAI,MAAM,CAAC,YAAY,EAAE,CAAC,YAAY,CAAC,MAAM,CAAC,YAAY,EAAE,CAAC,EAAE;YAC7D,IAAI,MAAM,CAAC,cAAc,CAAC,MAAM,CAAC,YAAY,EAAE,CAAC,IAAI,MAAM,CAAC,cAAc,CAAC,MAAM,CAAC,YAAY,EAAE,CAAC,EAAE;gBAChG,eAAe;gBACf,IAAI,CAAC,qBAAqB,CAAC,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC;aAClE;SACF;aAAM;YACL,MAAM,UAAU,GAA0B,EAAE,CAAC;YAC7C,GAAG,CAAC,6BAA6B,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;YACtD,MAAM,UAAU,GAA0B,EAAE,CAAC;YAC7C,GAAG,CAAC,6BAA6B,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;YACtD,KAAK,MAAM,OAAO,IAAI,UAAU,EAAE;gBAChC,KAAK,MAAM,OAAO,IAAI,UAAU,EAAE;oBAChC,IAAI,OAAO,CAAC,KAAK,CAAC,aAAa,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE;wBAC9C,IAAI,IAAI,CAAC,cAAc,CAAC,OAAO,EAAE,OAAO,CAAC,QAAQ,EAAE,OAAO,CAAC;+BACtD,IAAI,CAAC,cAAc,CAAC,OAAO,EAAE,OAAO,CAAC,QAAQ,EAAE,OAAO,CAAC,EAAE;4BAC5D,IAAI,CAAC,6BAA6B,CAAC,OAAO,CAAC,QAAQ,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAC5D,OAAO,CAAC,QAAQ,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAAC,CAAC;yBAC1C;qBACF;iBACF;aACF;SACF;IACH,CAAC;IACD,sCAAsC;IACtC,6EAA6E;IAC7E,+BAA+B;IACvB,yBAAyB,CAC/B,CAAQ,EACR,QAAiB,EACjB,IAAoB,EACpB,QAAiB,EACjB,SAAkB;QAElB;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;WA2DG;IACL,CAAC;IACD;;;;;;;;;OASG;IACH;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;iBAiGa;IAEb,sCAAsC;IACtC,6EAA6E;IAC7E,+BAA+B;IACvB,oCAAoC,CAC1C,QAA4B,EAC5B,QAA4B,EAC5B,SAAkB;QAClB;;;;;;;;;;;;;;;;;;;;;;;;;qBAyBa;IACf,CAAC;IAED;;;;OAIG;IACH;;;;;;;;;;;;;;;;;;;;;;OAsBG;IACH,qDAAqD;IACrD,6EAA6E;IAC7E,+BAA+B;IACvB,2BAA2B,CACjC,IAAoB,EACpB,SAAkB,EAClB,QAAiB,EACjB,WAAmB,EACnB,QAAiB,EACjB,WAAmB,EACnB,SAAkB,EAClB,OAAuB,EACvB,QAAiB,EACjB,SAAkB;QAElB;;;;;;;;;;;;;;;;;;;;aAoBK;IACP,CAAC;IAMD,0CAA0C;IACnC,8BAA8B,CAAC,IAAkB,EAAE,QAAiB,EAAE,OAAuB,EAAE,QAAiB,EAAE,SAAkB;QACzI;;;;;;;;;;;;;;;;;;;WAmBG;IACL,CAAC;IAED,gDAAgD;IACzC,wBAAwB,CAAC,GAAkB,EAAE,OAAgB,EAAE,GAAiB,EAAE,OAAgB,EAAE,QAAiB;QAC1H,MAAM,OAAO,GAAG,GAAG,CAAC,SAAS,CAAC;QAC9B,MAAM,OAAO,GAAG,GAAG,CAAC,SAAS,CAAC;QAC9B,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;QACrD,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;QACrD,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAC7B,IAAI,IAAI,GAAG,CAAC,EAAE;YACZ,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YAC7B,IAAI,GAAG,CAAC;YACR,IAAI,GAAG,CAAC;YACR,GAAG,GAAG,GAAG,CAAC;YACV,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;YACxB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE;gBACrE,GAAG,CAAC,OAAO,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;gBACzB,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;gBACf,IAAI,CAAC,sBAAsB,CACzB,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EACjD,GAAG,EAAE,EAAE,KAAK,CAAC,IAAI,OAAO,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,KAAK,IAAI,IAAI,OAAO,EAClF,QAAQ,CAAC,CAAC;aACb;SACF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IAED,4CAA4C;IACrC,oBAAoB,CAAC,IAAW,EAAE,OAAgB,EAAE,GAAiB,EAAE,OAAgB,EAAE,QAAiB;QAC/G,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;QACrD,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;QACrD,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAC7B,IAAI,IAAI,GAAG,CAAC,EAAE;YACZ,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YAC7B,IAAI,GAAG,CAAC;YACR,IAAI,GAAG,CAAC;YACR,GAAG,GAAG,GAAG,CAAC;YACV,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;YACxB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE;gBACrE,GAAG,CAAC,OAAO,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;gBACzB,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;gBACf,IAAI,CAAC,kBAAkB,CACrB,GAAG,EAAE,EAAE,KAAK,CAAC,IAAI,OAAO,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,KAAK,IAAI,IAAI,OAAO,EAClF,IAAI,EAAE,OAAO,EAAE,OAAO,EACtB,CAAC,QAAQ,CAAC,CAAC;aACd;SACF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IAED,2DAA2D;IAC3C,mBAAmB,CAAC,QAAuB;QACzD,IAAI,IAAI,CAAC,UAAU,YAAY,aAAa,EAAE;YAC5C,MAAM,QAAQ,GAAG,IAAI,CAAC,UAAU,CAAC;YACjC,IAAI,CAAC,sBAAsB,CACzB,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EACxF,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EACxF,KAAK,CAAC,CAAC;SACV;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,YAAY,EAAE;YAClD,IAAI,CAAC,wBAAwB,CAAC,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;SAC/F;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,KAAK,EAAE;YAC3C,IAAI,CAAC,kBAAkB,CACrB,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EACxF,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;SACzD;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,cAAc,EAAE;YACpD,IAAI,CAAC,2BAA2B,CAC9B,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EACxF,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;SAC1C;IACH,CAAC;IACD,+DAA+D;IAC/C,kBAAkB,CAAC,GAAiB;QAClD,IAAI,IAAI,CAAC,UAAU,YAAY,YAAY,EAAE;YAC3C,MAAM,GAAG,GAAG,IAAI,CAAC,UAAU,CAAC;YAC5B,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;YACrD,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;YACrD,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;YACrD,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;YACrD,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;YAC7B,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;YAC7B,IAAI,IAAI,GAAG,CAAC,IAAI,IAAI,GAAG,CAAC,EAAE;gBACxB,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;gBACxB,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;gBAC7B,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;gBAC7B,IAAI,GAAG,GAAG,GAAG,CAAC;gBACd,IAAI,GAAG,CAAC;gBACR,IAAI,GAAG,CAAC;gBACR,IAAI,GAAG,CAAC;gBACR,MAAM,OAAO,GAAG,IAAI,CAAC,QAAQ,CAAC;gBAC9B,MAAM,OAAO,GAAG,IAAI,CAAC,QAAQ,CAAC;gBAC9B,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;gBACxB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE;oBACrE,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;oBACf,GAAG,GAAG,GAAG,CAAC;oBACV,GAAG,CAAC,OAAO,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;oBACzB,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;oBACxB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE;wBACrE,GAAG,CAAC,OAAO,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;wBACzB,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;wBACf,IAAI,CAAC,sBAAsB,CACzB,GAAG,EAAE,EAAE,KAAK,CAAC,IAAI,OAAO,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,KAAK,IAAI,IAAI,OAAO,EAClF,GAAG,EAAE,EAAE,KAAK,CAAC,IAAI,OAAO,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,KAAK,IAAI,IAAI,OAAO,EAClF,KAAK,CAAC,CAAC;qBACV;iBACF;aACF;SACF;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,aAAa,EAAE;YACnD,IAAI,CAAC,wBAAwB,CAAC,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC;SACzF;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,KAAK,EAAE;YAC3C,IAAI,CAAC,oBAAoB,CAAC,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC;SACrF;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,cAAc,EAAE;YACpD,IAAI,CAAC,8BAA8B,CAAC,GAAG,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;SAChG;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,wDAAwD;IACxC,WAAW,CAAC,IAAW;QACrC,IAAI,IAAI,CAAC,UAAU,YAAY,aAAa,EAAE;YAC5C,IAAI,CAAC,kBAAkB,CACrB,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,UAAU,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EAC7G,IAAI,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC;SAC7C;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,YAAY,EAAE;YAClD,IAAI,CAAC,oBAAoB,CAAC,IAAI,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;SACvF;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,KAAK,EAAE;YAC3C,IAAI,CAAC,cAAc,CAAC,IAAI,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;SACjF;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,cAAc,EAAE;YACpD,IAAI,CAAC,yBAAyB,CAAC,IAAI,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;SAC5F;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,iEAAiE;IACjD,oBAAoB,CAAC,KAAqB;QACxD,IAAI,IAAI,CAAC,UAAU,YAAY,aAAa,EAAE;YAC5C,IAAI,CAAC,2BAA2B,CAC9B,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,UAAU,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EAC7G,KAAK,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC;SAC/B;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,YAAY,EAAE;YAClD,IAAI,CAAC,8BAA8B,CAAC,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC;SACjG;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,KAAK,EAAE;YAC3C,IAAI,CAAC,yBAAyB,CAAC,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC;SAC5F;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,kBAAkB,EAAE;YACxD,IAAI,CAAC,oCAAoC,CAAC,KAAK,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SAC1E;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,4EAA4E;IAC5D,qBAAqB,CAAC,MAAuB;QAC3D;;;;;;;;;;WAUG;QACH,OAAO,SAAS,CAAC;IACnB,CAAC;;AA/sBc,qCAAc,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC;AAshBnC,oCAAa,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;AACjC,oCAAa,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;AACjC,oCAAa,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;AACjC,oCAAa,GAAG,OAAO,CAAC,MAAM,EAAE,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\n\r\n/** @packageDocumentation\r\n * @module Curve\r\n */\r\n\r\nimport { BSplineCurve3d, BSplineCurve3dBase } from \"../bspline/BSplineCurve\";\r\n// import { BezierCurveBase } from \"../bspline/BezierCurveBase\";\r\n// import { BezierCurve3dH } from \"../bspline/BezierCurve3dH\";\r\n// import { UnivariateBezier } from \"../numerics/BezierPolynomials\";\r\nimport { BSplineCurve3dH } from \"../bspline/BSplineCurve3dH\";\r\nimport { Geometry } from \"../Geometry\";\r\nimport { NullGeometryHandler } from \"../geometry3d/GeometryHandler\";\r\nimport { Plane3dByOriginAndUnitNormal } from \"../geometry3d/Plane3dByOriginAndUnitNormal\";\r\n// import { Arc3d } from \"./Arc3d\";\r\nimport { Vector2d } from \"../geometry3d/Point2dVector2d\";\r\n// import { XYAndZ } from \"../geometry3d/XYZProps\";\r\nimport { Point3d, Vector3d } from \"../geometry3d/Point3dVector3d\";\r\n// import { LineString3d } from \"./LineString3d\";\r\nimport { SmallSystem, TrigPolynomial } from \"../numerics/Polynomials\";\r\n// import { Point4d } from \"../geometry4d/Point4d\";\r\n// import { Transform } from \"../geometry3d/Transform\";\r\n// import { Matrix3d } from \"../geometry3d/Matrix3d\";\r\nimport { Arc3d } from \"./Arc3d\";\r\n// import { Range3d } from \"../geometry3d/Range\";\r\nimport { CurveLocationDetailArrayPair } from \"./CurveCurveIntersectXY\";\r\nimport { CurveIntervalRole, CurveLocationDetail } from \"./CurveLocationDetail\";\r\nimport { CurvePrimitive } from \"./CurvePrimitive\";\r\nimport { GeometryQuery } from \"./GeometryQuery\";\r\nimport { LineSegment3d } from \"./LineSegment3d\";\r\nimport { LineString3d } from \"./LineString3d\";\r\n\r\n// cspell:word XYRR\r\n/**\r\n * * Handler class for XYZ intersections.\r\n * * Instances are initialized and called from CurveCurve.\r\n * * Constructor is told two geometry items A and B\r\n *   * geometryB is saved for later reference\r\n *   * type-specific handler methods will \"see\" geometry A repeatedly.\r\n *   * Hence geometryA is NOT saved by the constructor.\r\n * @internal\r\n */\r\nexport class CurveCurveIntersectXYZ extends NullGeometryHandler {\r\n  // private geometryA: GeometryQuery;  // nb never used -- passed through handlers.\r\n  private _extendA: boolean;\r\n  private _geometryB: GeometryQuery;\r\n  private _extendB: boolean;\r\n  private _results!: CurveLocationDetailArrayPair;\r\n  private reinitialize() {\r\n    this._results = new CurveLocationDetailArrayPair();\r\n  }\r\n  /**\r\n   *\r\n   * @param _geometryA first curve for intersection.  This is NOT saved.\r\n   * @param extendA flag to enable using extension of geometryA.\r\n   * @param geometryB second curve for intersection.  Saved for reference by specific handler methods.\r\n   * @param extendB flag for extension of geometryB.\r\n   */\r\n  public constructor(_geometryA: GeometryQuery, extendA: boolean, geometryB: GeometryQuery, extendB: boolean) {\r\n    super();\r\n    // this.geometryA = _geometryA;\r\n    this._extendA = extendA;\r\n    this._geometryB = geometryB;\r\n    this._extendB = extendB;\r\n    this.reinitialize();\r\n  }\r\n  /**\r\n   * * Return the results structure for the intersection calculation.\r\n   * @param reinitialize if true, a new results structure is created for use by later calls.\r\n   *\r\n   */\r\n  public grabResults(reinitialize: boolean = false): CurveLocationDetailArrayPair {\r\n    const result = this._results;\r\n    if (reinitialize)\r\n      this.reinitialize();\r\n    return result;\r\n  }\r\n\r\n  private static _workVector2dA = Vector2d.create();\r\n\r\n  private acceptFraction(extend0: boolean, fraction: number, extend1: boolean) {\r\n    if (!extend0 && fraction < 0.0)\r\n      return false;\r\n    if (!extend1 && fraction > 1.0)\r\n      return false;\r\n    return true;\r\n  }\r\n\r\n  /** compute intersection of two line segments.\r\n   * filter by extension rules.\r\n   * reject if evaluated points do not match coordinates (e.g. close approach point)\r\n   * record with fraction mapping.\r\n   */\r\n  private recordPointWithLocalFractions(\r\n    localFractionA: number,\r\n    cpA: CurvePrimitive,\r\n    fractionA0: number,\r\n    fractionA1: number,\r\n    localFractionB: number,   // Computed intersection fraction\r\n    cpB: CurvePrimitive,\r\n    fractionB0: number,\r\n    fractionB1: number,\r\n    reversed: boolean,\r\n  ) {\r\n    const globalFractionA = Geometry.interpolate(fractionA0, localFractionA, fractionA1);\r\n    const globalFractionB = Geometry.interpolate(fractionB0, localFractionB, fractionB1);\r\n    // ignore duplicate of most recent point .  ..\r\n    const numPrevious = this._results.dataA.length;\r\n    if (numPrevious > 0) {\r\n      const topFractionA = this._results.dataA[numPrevious - 1].fraction;\r\n      const topFractionB = this._results.dataB[numPrevious - 1].fraction;\r\n      if (reversed) {\r\n        if (Geometry.isAlmostEqualNumber(topFractionA, globalFractionB) && Geometry.isAlmostEqualNumber(topFractionB, globalFractionA))\r\n          return;\r\n      } else {\r\n        if (Geometry.isAlmostEqualNumber(topFractionA, globalFractionA) && Geometry.isAlmostEqualNumber(topFractionB, globalFractionB))\r\n          return;\r\n      }\r\n    }\r\n    const pointA = cpA.fractionToPoint(globalFractionA);\r\n    const pointB = cpB.fractionToPoint(globalFractionB);\r\n    if (!pointA.isAlmostEqualMetric(pointB))\r\n      return;\r\n\r\n    const detailA = CurveLocationDetail.createCurveFractionPoint(cpA,\r\n      globalFractionA, pointA);\r\n    detailA.setIntervalRole(CurveIntervalRole.isolated);\r\n    const detailB = CurveLocationDetail.createCurveFractionPoint(cpB,\r\n      globalFractionB, pointB);\r\n    detailB.setIntervalRole(CurveIntervalRole.isolated);\r\n    if (reversed) {\r\n      this._results.dataA.push(detailB);\r\n      this._results.dataB.push(detailA);\r\n    } else {\r\n      this._results.dataA.push(detailA);\r\n      this._results.dataB.push(detailB);\r\n    }\r\n  }\r\n  /** compute intersection of two line segments.\r\n   * filter by extension rules.\r\n   * record with fraction mapping.\r\n   */\r\n  private computeSegmentSegment3D(\r\n    cpA: CurvePrimitive,\r\n    extendA0: boolean,\r\n    pointA0: Point3d,\r\n    fractionA0: number,\r\n    pointA1: Point3d,\r\n    fractionA1: number,\r\n    extendA1: boolean,\r\n    cpB: CurvePrimitive,\r\n    extendB0: boolean,\r\n    pointB0: Point3d,\r\n    fractionB0: number,\r\n    pointB1: Point3d,\r\n    fractionB1: number,\r\n    extendB1: boolean,\r\n    reversed: boolean,\r\n  ) {\r\n    const uv = CurveCurveIntersectXYZ._workVector2dA;\r\n    if (SmallSystem.lineSegment3dClosestApproachUnbounded(\r\n      pointA0, pointA1,\r\n      pointB0, pointB1, uv)\r\n      && this.acceptFraction(extendA0, uv.x, extendA1)\r\n      && this.acceptFraction(extendB0, uv.y, extendB1)\r\n    ) {\r\n      this.recordPointWithLocalFractions(uv.x, cpA, fractionA0, fractionA1, uv.y, cpB, fractionB0, fractionB1, reversed);\r\n    }\r\n  }\r\n  // Caller accesses data from a line segment and passes to here.\r\n  // (The line segment in question might be (a) a full line segment or (b) a fragment within a linestring.  The fraction and extend parameters\r\n  // allow all combinations to be passed in)\r\n  // This method applies transform.\r\n  private dispatchSegmentSegment(\r\n    cpA: CurvePrimitive,\r\n    extendA0: boolean,\r\n    pointA0: Point3d,\r\n    fractionA0: number,\r\n    pointA1: Point3d,\r\n    fractionA1: number,\r\n    extendA1: boolean,\r\n    cpB: CurvePrimitive,\r\n    extendB0: boolean,\r\n    pointB0: Point3d,\r\n    fractionB0: number,\r\n    pointB1: Point3d,\r\n    fractionB1: number,\r\n    extendB1: boolean,\r\n    reversed: boolean,\r\n  ) {\r\n    this.computeSegmentSegment3D(\r\n      cpA, extendA0, pointA0, fractionA0, pointA1, fractionA1, extendA1,\r\n      cpB, extendB0, pointB0, fractionB0, pointB1, fractionB1, extendB1,\r\n      reversed);\r\n  }\r\n  /**\r\n   * Create a plane whose normal is the \"better\" cross product: `vectorA.crossProduct(vectorB)` or `vectorA.crossProduct(vectorC)`\r\n   * * The heuristic for \"better\" is:\r\n   *   * first choice is cross product with `vectorB`, if `vectorA` and `vectorB` are sufficiently far from parallel (or anti-parallel).\r\n   *   * otherwise use vectorC\r\n   * @param origin plane origin\r\n   * @param vectorA vector which must be in the plane.\r\n   * @param cosineValue largest cosine of the angle theta between vectorA and vectorB to prefer their cross product,\r\n   *                    e.g. passing 0.94 ~ cos(20deg) will switch to using vectorC in the cross product if theta < ~20deg or theta > ~160deg.\r\n   * @param vectorB first candidate for additional in-plane vector\r\n   * @param vectorC second candidate for additional in-plane vector\r\n   */\r\n  public createPlaneWithPreferredPerpendicular(origin: Point3d, vectorA: Vector3d, cosineValue: number, vectorB: Vector3d, vectorC: Vector3d): Plane3dByOriginAndUnitNormal | undefined {\r\n    cosineValue = Geometry.restrictToInterval(Math.abs(cosineValue), 0.0, 1.0 - Geometry.smallFraction);\r\n    const dotAA = vectorA.magnitudeSquared();\r\n    const dotBB = vectorB.magnitudeSquared();\r\n    const dotAB = Math.abs(vectorA.dotProduct(vectorB));\r\n    const cross = vectorA.unitCrossProduct(dotAB * dotAB <= cosineValue * cosineValue * dotAA * dotBB ? vectorB : vectorC);\r\n    if (cross)\r\n      return Plane3dByOriginAndUnitNormal.create(origin, cross);\r\n    return undefined;\r\n  }\r\n  // Caller accesses data from a linestring or segment and passes it here.\r\n  // (The line segment in question might be (a) a full line segment or (b) a fragment within a linestring.  The fraction and extend parameters\r\n  // allow all combinations to be passed in)\r\n  private dispatchSegmentArc(\r\n    cpA: CurvePrimitive,\r\n    extendA0: boolean,\r\n    pointA0: Point3d,\r\n    fractionA0: number,\r\n    pointA1: Point3d,\r\n    fractionA1: number,\r\n    extendA1: boolean,\r\n    arc: Arc3d,\r\n    extendB0: boolean,\r\n    extendB1: boolean,\r\n    reversed: boolean,\r\n  ) {\r\n    const lineVector = Vector3d.createStartEnd(pointA0, pointA1);\r\n    const plane = this.createPlaneWithPreferredPerpendicular(pointA0, lineVector, 0.94, arc.perpendicularVector, arc.vector0);\r\n    if (plane !== undefined) {\r\n      const candidates: CurveLocationDetail[] = [];\r\n      arc.appendPlaneIntersectionPoints(plane, candidates);\r\n      let lineFraction;\r\n      let linePoint: Point3d | undefined;\r\n      for (const c of candidates) {\r\n        if (this.acceptFraction(extendB0, c.fraction, extendB1)) {\r\n          lineFraction = SmallSystem.lineSegment3dClosestPointUnbounded(pointA0, pointA1, c.point);\r\n          if (lineFraction !== undefined) {\r\n            linePoint = pointA0.interpolate(lineFraction, pointA1, linePoint);\r\n            if (linePoint.isAlmostEqualMetric(c.point)\r\n              && this.acceptFraction(extendA0, lineFraction, extendA1)) {\r\n              this.recordPointWithLocalFractions(lineFraction, cpA, fractionA0, fractionA1,\r\n                c.fraction, arc, 0, 1, reversed);\r\n\r\n            }\r\n          }\r\n        }\r\n      }\r\n    }\r\n  }\r\n\r\n  // Caller promises arcs are coplanar.\r\n  // Passes \"other\" as {center, vector0, vector90} in local xy space of cpA\r\n  // Solves the arc-arc equations for that local ellipse with unit circle.\r\n  // Solution fractions map directly to original arcs.\r\n  private dispatchArcArcInPlane(\r\n    cpA: Arc3d,\r\n    extendA: boolean,\r\n    cpB: Arc3d,\r\n    extendB: boolean,\r\n    reversed: boolean,\r\n  ) {\r\n    const otherVectors = cpA.otherArcAsLocalVectors(cpB);\r\n    if (otherVectors !== undefined) {\r\n      const ellipseRadians: number[] = [];\r\n      const circleRadians: number[] = [];\r\n      TrigPolynomial.solveUnitCircleHomogeneousEllipseIntersection(\r\n        otherVectors.center.x, otherVectors.center.y, 1.0,\r\n        otherVectors.vector0.x, otherVectors.vector0.y, 0.0,\r\n        otherVectors.vector90.x, otherVectors.vector90.y, 0.0,\r\n        ellipseRadians, circleRadians);\r\n      for (let i = 0; i < ellipseRadians.length; i++) {\r\n        const fractionA = cpA.sweep.radiansToSignedPeriodicFraction(circleRadians[i]);\r\n        const fractionB = cpA.sweep.radiansToSignedPeriodicFraction(ellipseRadians[i]);\r\n        // hm .. do we really need to check the fractions?  We know they are internal to the beziers\r\n        if (this.acceptFraction(extendA, fractionA, extendA) && this.acceptFraction(extendB, fractionB, extendB)) {\r\n          this.recordPointWithLocalFractions(fractionA, cpA, 0, 1,\r\n            fractionB, cpB, 0, 1, reversed);\r\n        }\r\n      }\r\n    }\r\n  }\r\n  // Caller accesses data from two arcs.\r\n  // Selects the best conditioned arc (in xy parts) as \"circle after inversion\"\r\n  // Solves the arc-arc equations\r\n  private dispatchArcArc(\r\n    cpA: Arc3d,\r\n    extendA: boolean,\r\n    cpB: Arc3d,\r\n    extendB: boolean,\r\n    reversed: boolean,\r\n  ) {\r\n    // If arcs are in different planes:\r\n    // 1) Intersect each plane with the other arc (quadratic)\r\n    // 2) accept points that appear in both intersection sets.\r\n    // If arcs are in parallel planes -- no intersections\r\n    // If arcs are in the same plane -- xy intersection in that plane.\r\n    const planeA = Plane3dByOriginAndUnitNormal.create(cpA.center, cpA.perpendicularVector);\r\n    const planeB = Plane3dByOriginAndUnitNormal.create(cpB.center, cpB.perpendicularVector);\r\n    if (planeA === undefined || planeB === undefined)\r\n      return;\r\n\r\n    if (planeA.getNormalRef().isParallelTo(planeB.getNormalRef())) {\r\n      if (planeA.isPointInPlane(planeB.getOriginRef()) && planeB.isPointInPlane(planeA.getOriginRef())) {\r\n        // coplanar !!!\r\n        this.dispatchArcArcInPlane(cpA, extendA, cpB, extendB, reversed);\r\n      }\r\n    } else {\r\n      const arcBPoints: CurveLocationDetail[] = [];\r\n      cpB.appendPlaneIntersectionPoints(planeA, arcBPoints);\r\n      const arcAPoints: CurveLocationDetail[] = [];\r\n      cpA.appendPlaneIntersectionPoints(planeB, arcAPoints);\r\n      for (const detailB of arcBPoints) {\r\n        for (const detailA of arcAPoints) {\r\n          if (detailA.point.isAlmostEqual(detailB.point)) {\r\n            if (this.acceptFraction(extendA, detailA.fraction, extendA)\r\n              && this.acceptFraction(extendB, detailB.fraction, extendB)) {\r\n              this.recordPointWithLocalFractions(detailA.fraction, cpA, 0, 1,\r\n                detailB.fraction, cpB, 0, 1, reversed);\r\n            }\r\n          }\r\n        }\r\n      }\r\n    }\r\n  }\r\n  // Caller accesses data from two arcs.\r\n  // Selects the best conditioned arc (in xy parts) as \"circle after inversion\"\r\n  // Solves the arc-arc equations\r\n  private dispatchArcBsplineCurve3d(\r\n    _: Arc3d,\r\n    _extendA: boolean,\r\n    _cpB: BSplineCurve3d,\r\n    _extendB: boolean,\r\n    _reversed: boolean,\r\n  ) {\r\n    /*\r\n    // Arc: X = C + cU + sV\r\n    // implicitize the arc as viewed.  This \"3d\" matrix is homogeneous \"XYW\" not \"xyz\"\r\n    let matrixA: Matrix3d;\r\n    if (this._worldToLocalPerspective) {\r\n      const dataA = cpA.toTransformedPoint4d(this._worldToLocalPerspective);\r\n      matrixA = Matrix3d.createColumnsXYW(dataA.vector0, dataA.vector0.w, dataA.vector90, dataA.vector90.w, dataA.center, dataA.center.w);\r\n    } else {\r\n      const dataA = cpA.toTransformedVectors(this._worldToLocalAffine);\r\n      matrixA = Matrix3d.createColumnsXYW(dataA.vector0, 0, dataA.vector90, 0, dataA.center, 1);\r\n    }\r\n    // The worldToLocal has moved the arc vectors into screen space.\r\n    // matrixA captures the xyw parts (ignoring z)\r\n    // for any point in world space,\r\n    // THIS CODE ONLY WORKS FOR\r\n    const matrixAInverse = matrixA.inverse();\r\n    if (matrixAInverse) {\r\n      const orderF = cpB.order; // order of the beziers for simple coordinates\r\n      const orderG = 2 * orderF - 1;  // order of the (single) bezier for squared coordinates.\r\n      const coffF = new Float64Array(orderF);\r\n      const univariateBezierG = new UnivariateBezier(orderG);\r\n      const axx = matrixAInverse.at(0, 0); const axy = matrixAInverse.at(0, 1); const axz = 0.0; const axw = matrixAInverse.at(0, 2);\r\n      const ayx = matrixAInverse.at(1, 0); const ayy = matrixAInverse.at(1, 1); const ayz = 0.0; const ayw = matrixAInverse.at(1, 2);\r\n      const awx = matrixAInverse.at(2, 0); const awy = matrixAInverse.at(2, 1); const awz = 0.0; const aww = matrixAInverse.at(2, 2);\r\n\r\n      if (matrixAInverse) {\r\n        let bezier: BezierCurve3dH | undefined;\r\n        for (let spanIndex = 0; ; spanIndex++) {\r\n          bezier = cpB.getSaturatedBezierSpan3dH(spanIndex, bezier);\r\n          if (!bezier) break;\r\n          if (this._worldToLocalPerspective)\r\n            bezier.tryMultiplyMatrix4dInPlace(this._worldToLocalPerspective);\r\n          else if (this._worldToLocalAffine)\r\n            bezier.tryTransformInPlace(this._worldToLocalAffine);\r\n          univariateBezierG.zero();\r\n          bezier.poleProductsXYZW(coffF, axx, axy, axz, axw);\r\n          univariateBezierG.addSquaredSquaredBezier(coffF, 1.0);\r\n          bezier.poleProductsXYZW(coffF, ayx, ayy, ayz, ayw);\r\n          univariateBezierG.addSquaredSquaredBezier(coffF, 1.0);\r\n          bezier.poleProductsXYZW(coffF, awx, awy, awz, aww);\r\n          univariateBezierG.addSquaredSquaredBezier(coffF, -1.0);\r\n          const roots = univariateBezierG.roots(0.0, true);\r\n          if (roots) {\r\n            for (const root of roots) {\r\n              const fractionB = bezier.fractionToParentFraction(root);\r\n              // The univariate bezier (which has been transformed by the view transform) evaluates into xyw space\r\n              const bcurvePoint4d = bezier.fractionToPoint4d(root);\r\n              const c = bcurvePoint4d.dotProductXYZW(axx, axy, axz, axw);\r\n              const s = bcurvePoint4d.dotProductXYZW(ayx, ayy, ayz, ayw);\r\n              const arcFraction = cpA.sweep.radiansToSignedPeriodicFraction(Math.atan2(s, c));\r\n              if (this.acceptFraction(extendA, arcFraction, extendA) && this.acceptFraction(extendB, fractionB, extendB)) {\r\n                this.recordPointWithLocalFractions(arcFraction, cpA, 0, 1,\r\n                  fractionB, cpB, 0, 1, reversed);\r\n              }\r\n            }\r\n          }\r\n        }\r\n      }\r\n    }\r\n     */\r\n  }\r\n  /*\r\n  // apply the transformation to bezier curves. optionally construct ranges.\r\n  private transformBeziers(beziers: BezierCurve3dH[]) {\r\n    if (this._worldToLocalAffine) {\r\n      for (const bezier of beziers) bezier.tryTransformInPlace(this._worldToLocalAffine);\r\n    } else if (this._worldToLocalPerspective) {\r\n      for (const bezier of beziers) bezier.tryMultiplyMatrix4dInPlace(this._worldToLocalPerspective);\r\n    }\r\n  }\r\n   */\r\n  /*\r\n  private getRanges(beziers: BezierCurveBase[]): Range3d[] {\r\n    const ranges: Range3d[] = [];\r\n    ranges.length = 0;\r\n    for (const b of beziers) {\r\n      ranges.push(b.range());\r\n    }\r\n    return ranges;\r\n  }\r\n  private _xyzwA0?: Point4d;\r\n  private _xyzwA1?: Point4d;\r\n  private _xyzwPlane?: Point4d;\r\n  private _xyzwB?: Point4d;\r\n\r\n  private dispatchBezierBezierStrokeFirst(\r\n    bezierA: BezierCurve3dH,\r\n    bcurveA: BSplineCurve3dBase,\r\n    strokeCountA: number,\r\n    bezierB: BezierCurve3dH,\r\n    bcurveB: BSplineCurve3dBase,\r\n    _strokeCountB: number,\r\n    univariateBezierB: UnivariateBezier,  // caller-allocated for univariate coefficients.\r\n    reversed: boolean) {\r\n    if (!this._xyzwA0) this._xyzwA0 = Point4d.create();\r\n    if (!this._xyzwA1) this._xyzwA1 = Point4d.create();\r\n    if (!this._xyzwPlane) this._xyzwPlane = Point4d.create();\r\n    if (!this._xyzwB) this._xyzwB = Point4d.create();\r\n    /-*\r\n\r\n              const roots = univariateBezierG.roots(0.0, true);\r\n              if (roots) {\r\n                for (const root of roots) {\r\n                  const fractionB = bezier.fractionToParentFraction(root);\r\n                  // The univariate bezier (which has been transformed by the view transform) evaluates into xyw space\r\n                  const bcurvePoint4d = bezier.fractionToPoint4d(root);\r\n                  const c = bcurvePoint4d.dotProductXYZW(axx, axy, axz, axw);\r\n                  const s = bcurvePoint4d.dotProductXYZW(ayx, ayy, ayz, ayw);\r\n                  const arcFraction = cpA.sweep.radiansToSignedPeriodicFraction(Math.atan2(s, c));\r\n                  if (this.acceptFraction(extendA, arcFraction, extendA) && this.acceptFraction(extendB, fractionB, extendB)) {\r\n                    this.recordPointWithLocalFractions(arcFraction, cpA, 0, 1,\r\n                      fractionB, cpB, 0, 1, reversed);\r\n                  }\r\n                }\r\n   *-/\r\n    bezierA.fractionToPoint4d(0.0, this._xyzwA0);\r\n    let f0 = 0.0;\r\n    let f1 = 1.0;\r\n    const intervalTolerance = 1.0e-5;\r\n    const df = 1.0 / strokeCountA;\r\n    for (let i = 1; i <= strokeCountA; i++ , f0 = f1, this._xyzwA0.setFrom(this._xyzwA1)) {\r\n      f1 = i * df;\r\n      bezierA.fractionToPoint4d(f1, this._xyzwA1);\r\n      Point4d.createPlanePointPointZ(this._xyzwA0, this._xyzwA1, this._xyzwPlane);\r\n      bezierB.poleProductsXYZW(univariateBezierB.coffs, this._xyzwPlane.x, this._xyzwPlane.y, this._xyzwPlane.z, this._xyzwPlane.w);\r\n      let errors = 0;\r\n      const roots = univariateBezierB.roots(0.0, true);\r\n      if (roots)\r\n        for (const r of roots) {\r\n          const bezierBFraction = r;\r\n          bezierB.fractionToPoint4d(bezierBFraction, this._xyzwB);\r\n          const segmentAFraction = SmallSystem.lineSegment3dHXYClosestPointUnbounded(this._xyzwA0, this._xyzwA1, this._xyzwB);\r\n          if (segmentAFraction && Geometry.isIn01WithTolerance(segmentAFraction, intervalTolerance)) {\r\n            const bezierAFraction = Geometry.interpolate(f0, segmentAFraction, f1);\r\n            /*- TODO implement newton search\r\n            const xyMatchingFunction = new BezierBezierIntersectionXYRRToRRD(bezierA, bezierB);\r\n            const newtonSearcher = new Newton2dUnboundedWithDerivative(xyMatchingFunction);\r\n            newtonSearcher.setUV(bezierAFraction, bezierBFraction);\r\n            if (newtonSearcher.runIterations()) {\r\n              bezierAFraction = newtonSearcher.getU();\r\n              bezierBFraction = newtonSearcher.getV();\r\n            }\r\n             *-/\r\n              // We have a near intersection at fractions on the two beziers !!!\r\n              // Iterate on the curves for a true intersection ....\r\n              // NEEDS WORK -- just accept . . .\r\n            const bcurveAFraction = bezierA.fractionToParentFraction(bezierAFraction);\r\n            const bcurveBFraction = bezierB.fractionToParentFraction(bezierBFraction);\r\n            const xyzA0 = bezierA.fractionToPoint(bezierAFraction);\r\n            const xyzA1 = bcurveA.fractionToPoint(bcurveAFraction);\r\n            const xyzB0 = bezierB.fractionToPoint(bezierBFraction);\r\n            const xyzB1 = bcurveB.fractionToPoint(bcurveBFraction);\r\n            if (!xyzA0.isAlmostEqualXY(xyzA1))\r\n              errors++;\r\n            if (!xyzB0.isAlmostEqualXY(xyzB1))\r\n              errors++;\r\n            if (errors > 0 && !xyzA0.isAlmostEqual(xyzB0))\r\n              errors++;\r\n            if (errors > 0 && !xyzA1.isAlmostEqual(xyzB1))\r\n              errors++;\r\n            if (this.acceptFraction(false, bcurveAFraction, false) && this.acceptFraction(false, bcurveBFraction, false)) {\r\n              this.recordPointWithLocalFractions(bcurveAFraction, bcurveA, 0, 1,\r\n                bcurveBFraction, bcurveB, 0, 1, reversed);\r\n            }\r\n          }\r\n        }\r\n    }\r\n  }\r\n             */\r\n\r\n  // Caller accesses data from two arcs.\r\n  // Selects the best conditioned arc (in xy parts) as \"circle after inversion\"\r\n  // Solves the arc-arc equations\r\n  private dispatchBSplineCurve3dBSplineCurve3d(\r\n    _bcurveA: BSplineCurve3dBase,\r\n    _bcurveB: BSplineCurve3dBase,\r\n    _reversed: boolean) {\r\n    /*\r\n  const bezierSpanA = bcurveA.collectBezierSpans(true) as BezierCurve3dH[];\r\n  const bezierSpanB = bcurveB.collectBezierSpans(true) as BezierCurve3dH[];\r\n  const numA = bezierSpanA.length;\r\n  const numB = bezierSpanB.length;\r\n  this.transformBeziers(bezierSpanA);\r\n  this.transformBeziers(bezierSpanB);\r\n  const rangeA = this.getRanges(bezierSpanA);\r\n  const rangeB = this.getRanges(bezierSpanB);\r\n  const orderA = bcurveA.order;\r\n  const orderB = bcurveB.order;\r\n  const univariateCoffsA = new UnivariateBezier(orderA);\r\n  const univariateCoffsB = new UnivariateBezier(orderB);\r\n  for (let a = 0; a < numA; a++) {\r\n    for (let b = 0; b < numB; b++) {\r\n      if (rangeA[a].intersectsRangeXY(rangeB[b])) {\r\n        const strokeCountA = bezierSpanA[a].computeStrokeCountForOptions();\r\n        const strokeCountB = bezierSpanB[b].computeStrokeCountForOptions();\r\n        if (strokeCountA < strokeCountB)\r\n          this.dispatchBezierBezierStrokeFirst(bezierSpanA[a], bcurveA, strokeCountA, bezierSpanB[b], bcurveB, strokeCountB, univariateCoffsB, !_reversed);\r\n        else\r\n          this.dispatchBezierBezierStrokeFirst(bezierSpanB[b], bcurveB, strokeCountB, bezierSpanA[a], bcurveA, strokeCountA, univariateCoffsA, _reversed);\r\n      }\r\n    }\r\n  }\r\n               */\r\n  }\r\n\r\n  /**\r\n   * Apply the projection transform (if any) to (xyz, w)\r\n   * @param xyz xyz parts of input point.\r\n   * @param w   weight to use for homogeneous effects\r\n   */\r\n  /*\r\n  private projectPoint(xyz: XYAndZ, w: number = 1.0): Point4d {\r\n    if (this._worldToLocalPerspective)\r\n      return this._worldToLocalPerspective.multiplyPoint3d(xyz, w);\r\n    if (this._worldToLocalAffine)\r\n      return this._worldToLocalAffine.multiplyXYZW(xyz.x, xyz.y, xyz.z, w);\r\n    return Point4d.createFromPointAndWeight(xyz, w);\r\n  }\r\n  private mapNPCPlaneToWorld(npcPlane: Point4d, worldPlane: Point4d) {\r\n    // for NPC pointY, Y^ * H = 0 is \"on\" plane H.  (Hat is transpose)\r\n    // NPC Y is A*X for our transform A and worldPointX.\r\n    // hence (A X)^ * H = 0\r\n    // hence X^ * A^ * H = 0\r\n    // hence K = A^ * H\r\n    if (this._worldToLocalAffine) {\r\n      this._worldToLocalAffine.multiplyTransposeXYZW(npcPlane.x, npcPlane.y, npcPlane.z, npcPlane.w, worldPlane);\r\n    } else if (this._worldToLocalPerspective) {\r\n      this._worldToLocalPerspective.multiplyTransposePoint4d(npcPlane, worldPlane);\r\n    } else {\r\n      npcPlane.clone(worldPlane);\r\n    }\r\n  }\r\n   */\r\n  // Caller accesses data from segment and bsplineCurve\r\n  // Selects the best conditioned arc (in xy parts) as \"circle after inversion\"\r\n  // Solves the arc-arc equations\r\n  private dispatchSegmentBsplineCurve(\r\n    _cpA: CurvePrimitive,\r\n    _extendA0: boolean,\r\n    _pointA0: Point3d,\r\n    _fractionA0: number,\r\n    _pointA1: Point3d,\r\n    _fractionA1: number,\r\n    _extendA1: boolean,\r\n    _bcurve: BSplineCurve3d,\r\n    _extendB: boolean,\r\n    _reversed: boolean,\r\n  ) {\r\n    /*\r\n    const pointA0H = this.projectPoint(pointA0);\r\n    const pointA1H = this.projectPoint(pointA1);\r\n    const planeCoffs = Point4d.createPlanePointPointZ(pointA0H, pointA1H);\r\n    this.mapNPCPlaneToWorld(planeCoffs, planeCoffs);\r\n      // NOW .. we have a plane in world space.  Intersect it with the bspline:\r\n    const intersections: CurveLocationDetail[] = [];\r\n    bcurve.appendPlaneIntersectionPoints(planeCoffs, intersections);\r\n      // intersections has WORLD points with bspline fractions.   (The bspline fractions are all good 0..1 fractions within the spline.)\r\n      // accept those that are within the segment range.\r\n    for (const detail of intersections) {\r\n      const fractionB = detail.fraction;\r\n      const curvePoint = detail.point;\r\n      const curvePointH = this.projectPoint(curvePoint);\r\n      const lineFraction = SmallSystem.lineSegment3dHXYClosestPointUnbounded(pointA0H, pointA1H, curvePointH);\r\n      if (lineFraction !== undefined && this.acceptFraction(extendA0, lineFraction, extendA1) && this.acceptFraction(extendB, fractionB, extendB)) {\r\n        this.recordPointWithLocalFractions(lineFraction, cpA, fractionA0, fractionA1,\r\n          fractionB, bcurve, 0, 1, reversed);\r\n      }\r\n    }\r\n       */\r\n  }\r\n\r\n  private static _workPointAA0 = Point3d.create();\r\n  private static _workPointAA1 = Point3d.create();\r\n  private static _workPointBB0 = Point3d.create();\r\n  private static _workPointBB1 = Point3d.create();\r\n  /** low lever bspline curve -- STUB  .. */\r\n  public dispatchLineStringBSplineCurve(_lsA: LineString3d, _extendA: boolean, _curveB: BSplineCurve3d, _extendB: boolean, _reversed: boolean): any {\r\n    /*\r\n    const numA = lsA.numPoints();\r\n    if (numA > 1) {\r\n      const dfA = 1.0 / (numA - 1);\r\n      let fA0;\r\n      let fA1;\r\n      fA0 = 0.0;\r\n      const pointA0 = CurveCurveIntersectXYZ._workPointA0;\r\n      const pointA1 = CurveCurveIntersectXYZ._workPointA1;\r\n      lsA.pointAt(0, pointA0);\r\n      for (let iA = 1; iA < numA; iA++ , pointA0.setFrom(pointA1), fA0 = fA1) {\r\n        lsA.pointAt(iA, pointA1);\r\n        fA1 = iA * dfA;\r\n        this.dispatchSegmentBsplineCurve(\r\n          lsA, iA === 1 && extendA, pointA0, fA0, pointA1, fA1, (iA + 1) === numA && extendA,\r\n          curveB, extendB, reversed);\r\n      }\r\n    }\r\n    return undefined;\r\n     */\r\n  }\r\n\r\n  /** low lever segment intersect linestring .. */\r\n  public computeSegmentLineString(lsA: LineSegment3d, extendA: boolean, lsB: LineString3d, extendB: boolean, reversed: boolean): any {\r\n    const pointA0 = lsA.point0Ref;\r\n    const pointA1 = lsA.point1Ref;\r\n    const pointB0 = CurveCurveIntersectXYZ._workPointBB0;\r\n    const pointB1 = CurveCurveIntersectXYZ._workPointBB1;\r\n    const numB = lsB.numPoints();\r\n    if (numB > 1) {\r\n      const dfB = 1.0 / (numB - 1);\r\n      let fB0;\r\n      let fB1;\r\n      fB0 = 0.0;\r\n      lsB.pointAt(0, pointB0);\r\n      for (let ib = 1; ib < numB; ib++, pointB0.setFrom(pointB1), fB0 = fB1) {\r\n        lsB.pointAt(ib, pointB1);\r\n        fB1 = ib * dfB;\r\n        this.dispatchSegmentSegment(\r\n          lsA, extendA, pointA0, 0.0, pointA1, 1.0, extendA,\r\n          lsB, ib === 1 && extendB, pointB0, fB0, pointB1, fB1, (ib + 1) === numB && extendB,\r\n          reversed);\r\n      }\r\n    }\r\n    return undefined;\r\n  }\r\n\r\n  /** low lever arc intersect linestring .. */\r\n  public computeArcLineString(arcA: Arc3d, extendA: boolean, lsB: LineString3d, extendB: boolean, reversed: boolean): any {\r\n    const pointB0 = CurveCurveIntersectXYZ._workPointBB0;\r\n    const pointB1 = CurveCurveIntersectXYZ._workPointBB1;\r\n    const numB = lsB.numPoints();\r\n    if (numB > 1) {\r\n      const dfB = 1.0 / (numB - 1);\r\n      let fB0;\r\n      let fB1;\r\n      fB0 = 0.0;\r\n      lsB.pointAt(0, pointB0);\r\n      for (let ib = 1; ib < numB; ib++, pointB0.setFrom(pointB1), fB0 = fB1) {\r\n        lsB.pointAt(ib, pointB1);\r\n        fB1 = ib * dfB;\r\n        this.dispatchSegmentArc(\r\n          lsB, ib === 1 && extendB, pointB0, fB0, pointB1, fB1, (ib + 1) === numB && extendB,\r\n          arcA, extendA, extendA,\r\n          !reversed);\r\n      }\r\n    }\r\n    return undefined;\r\n  }\r\n\r\n  /** double dispatch handler for strongly typed segment.. */\r\n  public override handleLineSegment3d(segmentA: LineSegment3d): any {\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      const segmentB = this._geometryB;\r\n      this.dispatchSegmentSegment(\r\n        segmentA, this._extendA, segmentA.point0Ref, 0.0, segmentA.point1Ref, 1.0, this._extendA,\r\n        segmentB, this._extendB, segmentB.point0Ref, 0.0, segmentB.point1Ref, 1.0, this._extendB,\r\n        false);\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.computeSegmentLineString(segmentA, this._extendA, this._geometryB, this._extendB, false);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchSegmentArc(\r\n        segmentA, this._extendA, segmentA.point0Ref, 0.0, segmentA.point1Ref, 1.0, this._extendA,\r\n        this._geometryB, this._extendB, this._extendB, false);\r\n    } else if (this._geometryB instanceof BSplineCurve3d) {\r\n      this.dispatchSegmentBsplineCurve(\r\n        segmentA, this._extendA, segmentA.point0Ref, 0.0, segmentA.point1Ref, 1.0, this._extendA,\r\n        this._geometryB, this._extendB, false);\r\n    }\r\n  }\r\n  /** double dispatch handler for strongly typed linestring .. */\r\n  public override handleLineString3d(lsA: LineString3d): any {\r\n    if (this._geometryB instanceof LineString3d) {\r\n      const lsB = this._geometryB;\r\n      const pointA0 = CurveCurveIntersectXYZ._workPointAA0;\r\n      const pointA1 = CurveCurveIntersectXYZ._workPointAA1;\r\n      const pointB0 = CurveCurveIntersectXYZ._workPointBB0;\r\n      const pointB1 = CurveCurveIntersectXYZ._workPointBB1;\r\n      const numA = lsA.numPoints();\r\n      const numB = lsB.numPoints();\r\n      if (numA > 1 && numB > 1) {\r\n        lsA.pointAt(0, pointA0);\r\n        const dfA = 1.0 / (numA - 1);\r\n        const dfB = 1.0 / (numB - 1);\r\n        let fA0 = 0.0;\r\n        let fB0;\r\n        let fA1;\r\n        let fB1;\r\n        const extendA = this._extendA;\r\n        const extendB = this._extendB;\r\n        lsA.pointAt(0, pointA0);\r\n        for (let ia = 1; ia < numA; ia++, pointA0.setFrom(pointA1), fA0 = fA1) {\r\n          fA1 = ia * dfA;\r\n          fB0 = 0.0;\r\n          lsA.pointAt(ia, pointA1);\r\n          lsB.pointAt(0, pointB0);\r\n          for (let ib = 1; ib < numB; ib++, pointB0.setFrom(pointB1), fB0 = fB1) {\r\n            lsB.pointAt(ib, pointB1);\r\n            fB1 = ib * dfB;\r\n            this.dispatchSegmentSegment(\r\n              lsA, ia === 1 && extendA, pointA0, fA0, pointA1, fA1, (ia + 1) === numA && extendA,\r\n              lsB, ib === 1 && extendB, pointB0, fB0, pointB1, fB1, (ib + 1) === numB && extendB,\r\n              false);\r\n          }\r\n        }\r\n      }\r\n    } else if (this._geometryB instanceof LineSegment3d) {\r\n      this.computeSegmentLineString(this._geometryB, this._extendB, lsA, this._extendA, true);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.computeArcLineString(this._geometryB, this._extendB, lsA, this._extendA, true);\r\n    } else if (this._geometryB instanceof BSplineCurve3d) {\r\n      this.dispatchLineStringBSplineCurve(lsA, this._extendA, this._geometryB, this._extendB, false);\r\n    }\r\n    return undefined;\r\n  }\r\n  /** double dispatch handler for strongly typed arc .. */\r\n  public override handleArc3d(arc0: Arc3d): any {\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      this.dispatchSegmentArc(\r\n        this._geometryB, this._extendB, this._geometryB.point0Ref, 0.0, this._geometryB.point1Ref, 1.0, this._extendB,\r\n        arc0, this._extendA, this._extendA, true);\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.computeArcLineString(arc0, this._extendA, this._geometryB, this._extendB, false);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchArcArc(arc0, this._extendA, this._geometryB, this._extendB, false);\r\n    } else if (this._geometryB instanceof BSplineCurve3d) {\r\n      this.dispatchArcBsplineCurve3d(arc0, this._extendA, this._geometryB, this._extendB, false);\r\n    }\r\n    return undefined;\r\n  }\r\n  /** double dispatch handler for strongly typed bspline curve.. */\r\n  public override handleBSplineCurve3d(curve: BSplineCurve3d): any {\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      this.dispatchSegmentBsplineCurve(\r\n        this._geometryB, this._extendB, this._geometryB.point0Ref, 0.0, this._geometryB.point1Ref, 1.0, this._extendB,\r\n        curve, this._extendA, true);\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.dispatchLineStringBSplineCurve(this._geometryB, this._extendB, curve, this._extendA, true);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchArcBsplineCurve3d(this._geometryB, this._extendB, curve, this._extendA, true);\r\n    } else if (this._geometryB instanceof BSplineCurve3dBase) {\r\n      this.dispatchBSplineCurve3dBSplineCurve3d(curve, this._geometryB, false);\r\n    }\r\n    return undefined;\r\n  }\r\n  /** double dispatch handler for strongly typed homogeneous bspline curve. */\r\n  public override handleBSplineCurve3dH(_curve: BSplineCurve3dH): any {\r\n    /* NEEDS WORK -- make \"dispatch\" methods tolerant of both 3d and 3dH ...\"easy\" if both present BezierCurve3dH span loaders\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      this.dispatchSegmentBsplineCurve(\r\n        this._geometryB, this._extendB, this._geometryB.point0Ref, 0.0, this._geometryB.point1Ref, 1.0, this._extendB,\r\n        curve, this._extendA, true);\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.dispatchLineStringBSplineCurve(this._geometryB, this._extendB, curve, this._extendA, true);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchArcBsplineCurve3d(this._geometryB, this._extendB, curve, this._extendA, true);\r\n    }\r\n     */\r\n    return undefined;\r\n  }\r\n}\r\n"]}
+{"version":3,"file":"CurveCurveIntersectXYZ.js","sourceRoot":"","sources":["../../../src/curve/CurveCurveIntersectXYZ.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AAEH,OAAO,EAAE,cAAc,EAAE,kBAAkB,EAAE,MAAM,yBAAyB,CAAC;AAK7E,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AACvC,OAAO,EAAE,mBAAmB,EAAE,MAAM,+BAA+B,CAAC;AACpE,OAAO,EAAE,4BAA4B,EAAE,MAAM,4CAA4C,CAAC;AAC1F,mCAAmC;AACnC,OAAO,EAAE,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AACzD,mDAAmD;AACnD,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AAClE,iDAAiD;AACjD,OAAO,EAAE,WAAW,EAAE,cAAc,EAAE,MAAM,yBAAyB,CAAC;AACtE,mDAAmD;AACnD,uDAAuD;AACvD,qDAAqD;AACrD,OAAO,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC;AAChC,iDAAiD;AACjD,OAAO,EAAE,4BAA4B,EAAE,MAAM,yBAAyB,CAAC;AACvE,OAAO,EAAE,iBAAiB,EAAE,mBAAmB,EAAE,MAAM,uBAAuB,CAAC;AAG/E,OAAO,EAAE,aAAa,EAAE,MAAM,iBAAiB,CAAC;AAChD,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAE9C,mBAAmB;AACnB;;;;;;;;GAQG;AACH,MAAM,OAAO,sBAAuB,SAAQ,mBAAmB;IAS7D;;;;;;OAMG;IACH,YAAmB,UAAyB,EAAE,OAAgB,EAAE,SAAwB,EAAE,OAAgB;QACxG,KAAK,EAAE,CAAC;QACR,+BAA+B;QAC/B,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;QACxB,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;QAC5B,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;QACxB,IAAI,CAAC,YAAY,EAAE,CAAC;IACtB,CAAC;IAjBO,YAAY;QAClB,IAAI,CAAC,QAAQ,GAAG,IAAI,4BAA4B,EAAE,CAAC;IACrD,CAAC;IAgBD;;;;OAIG;IACI,WAAW,CAAC,eAAwB,KAAK;QAC9C,MAAM,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC;QAC7B,IAAI,YAAY;YACd,IAAI,CAAC,YAAY,EAAE,CAAC;QACtB,OAAO,MAAM,CAAC;IAChB,CAAC;IAIO,cAAc,CAAC,OAAgB,EAAE,QAAgB,EAAE,OAAgB;QACzE,IAAI,CAAC,OAAO,IAAI,QAAQ,GAAG,GAAG;YAC5B,OAAO,KAAK,CAAC;QACf,IAAI,CAAC,OAAO,IAAI,QAAQ,GAAG,GAAG;YAC5B,OAAO,KAAK,CAAC;QACf,OAAO,IAAI,CAAC;IACd,CAAC;IAED;;;;OAIG;IACK,6BAA6B,CACnC,cAAsB,EACtB,GAAmB,EACnB,UAAkB,EAClB,UAAkB,EAClB,cAAsB,EAAI,iCAAiC;IAC3D,GAAmB,EACnB,UAAkB,EAClB,UAAkB,EAClB,QAAiB;QAEjB,MAAM,eAAe,GAAG,QAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,cAAc,EAAE,UAAU,CAAC,CAAC;QACrF,MAAM,eAAe,GAAG,QAAQ,CAAC,WAAW,CAAC,UAAU,EAAE,cAAc,EAAE,UAAU,CAAC,CAAC;QACrF,8CAA8C;QAC9C,MAAM,WAAW,GAAG,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,MAAM,CAAC;QAC/C,IAAI,WAAW,GAAG,CAAC,EAAE;YACnB,MAAM,YAAY,GAAG,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC;YACnE,MAAM,YAAY,GAAG,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC;YACnE,IAAI,QAAQ,EAAE;gBACZ,IAAI,QAAQ,CAAC,mBAAmB,CAAC,YAAY,EAAE,eAAe,CAAC,IAAI,QAAQ,CAAC,mBAAmB,CAAC,YAAY,EAAE,eAAe,CAAC;oBAC5H,OAAO;aACV;iBAAM;gBACL,IAAI,QAAQ,CAAC,mBAAmB,CAAC,YAAY,EAAE,eAAe,CAAC,IAAI,QAAQ,CAAC,mBAAmB,CAAC,YAAY,EAAE,eAAe,CAAC;oBAC5H,OAAO;aACV;SACF;QACD,MAAM,MAAM,GAAG,GAAG,CAAC,eAAe,CAAC,eAAe,CAAC,CAAC;QACpD,MAAM,MAAM,GAAG,GAAG,CAAC,eAAe,CAAC,eAAe,CAAC,CAAC;QACpD,IAAI,CAAC,MAAM,CAAC,mBAAmB,CAAC,MAAM,CAAC;YACrC,OAAO;QAET,MAAM,OAAO,GAAG,mBAAmB,CAAC,wBAAwB,CAAC,GAAG,EAC9D,eAAe,EAAE,MAAM,CAAC,CAAC;QAC3B,OAAO,CAAC,eAAe,CAAC,iBAAiB,CAAC,QAAQ,CAAC,CAAC;QACpD,MAAM,OAAO,GAAG,mBAAmB,CAAC,wBAAwB,CAAC,GAAG,EAC9D,eAAe,EAAE,MAAM,CAAC,CAAC;QAC3B,OAAO,CAAC,eAAe,CAAC,iBAAiB,CAAC,QAAQ,CAAC,CAAC;QACpD,IAAI,QAAQ,EAAE;YACZ,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;YAClC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;SACnC;aAAM;YACL,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;YAClC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;SACnC;IACH,CAAC;IACD;;;OAGG;IACK,uBAAuB,CAC7B,GAAmB,EACnB,QAAiB,EACjB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,QAAiB,EACjB,GAAmB,EACnB,QAAiB,EACjB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,QAAiB,EACjB,QAAiB;QAEjB,MAAM,EAAE,GAAG,sBAAsB,CAAC,cAAc,CAAC;QACjD,IAAI,WAAW,CAAC,qCAAqC,CACnD,OAAO,EAAE,OAAO,EAChB,OAAO,EAAE,OAAO,EAAE,EAAE,CAAC;eAClB,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,EAAE,CAAC,CAAC,EAAE,QAAQ,CAAC;eAC7C,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,EAAE,CAAC,CAAC,EAAE,QAAQ,CAAC,EAChD;YACA,IAAI,CAAC,6BAA6B,CAAC,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,UAAU,EAAE,UAAU,EAAE,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,UAAU,EAAE,UAAU,EAAE,QAAQ,CAAC,CAAC;SACpH;IACH,CAAC;IACD,+DAA+D;IAC/D,4IAA4I;IAC5I,0CAA0C;IAC1C,iCAAiC;IACzB,sBAAsB,CAC5B,GAAmB,EACnB,QAAiB,EACjB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,QAAiB,EACjB,GAAmB,EACnB,QAAiB,EACjB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,QAAiB,EACjB,QAAiB;QAEjB,IAAI,CAAC,uBAAuB,CAC1B,GAAG,EAAE,QAAQ,EAAE,OAAO,EAAE,UAAU,EAAE,OAAO,EAAE,UAAU,EAAE,QAAQ,EACjE,GAAG,EAAE,QAAQ,EAAE,OAAO,EAAE,UAAU,EAAE,OAAO,EAAE,UAAU,EAAE,QAAQ,EACjE,QAAQ,CAAC,CAAC;IACd,CAAC;IACD;;;;;;;;;;;OAWG;IACI,qCAAqC,CAAC,MAAe,EAAE,OAAiB,EAAE,WAAmB,EAAE,OAAiB,EAAE,OAAiB;QACxI,WAAW,GAAG,QAAQ,CAAC,kBAAkB,CAAC,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,EAAE,GAAG,EAAE,GAAG,GAAG,QAAQ,CAAC,aAAa,CAAC,CAAC;QACpG,MAAM,KAAK,GAAG,OAAO,CAAC,gBAAgB,EAAE,CAAC;QACzC,MAAM,KAAK,GAAG,OAAO,CAAC,gBAAgB,EAAE,CAAC;QACzC,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,UAAU,CAAC,OAAO,CAAC,CAAC,CAAC;QACpD,MAAM,KAAK,GAAG,OAAO,CAAC,gBAAgB,CAAC,KAAK,GAAG,KAAK,IAAI,WAAW,GAAG,WAAW,GAAG,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC;QACvH,IAAI,KAAK;YACP,OAAO,4BAA4B,CAAC,MAAM,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;QAC5D,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,wEAAwE;IACxE,4IAA4I;IAC5I,0CAA0C;IAClC,kBAAkB,CACxB,GAAmB,EACnB,QAAiB,EACjB,OAAgB,EAChB,UAAkB,EAClB,OAAgB,EAChB,UAAkB,EAClB,QAAiB,EACjB,GAAU,EACV,QAAiB,EACjB,QAAiB,EACjB,QAAiB;QAEjB,MAAM,UAAU,GAAG,QAAQ,CAAC,cAAc,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;QAC7D,MAAM,KAAK,GAAG,IAAI,CAAC,qCAAqC,CAAC,OAAO,EAAE,UAAU,EAAE,IAAI,EAAE,GAAG,CAAC,mBAAmB,EAAE,GAAG,CAAC,OAAO,CAAC,CAAC;QAC1H,IAAI,KAAK,KAAK,SAAS,EAAE;YACvB,MAAM,UAAU,GAA0B,EAAE,CAAC;YAC7C,GAAG,CAAC,6BAA6B,CAAC,KAAK,EAAE,UAAU,CAAC,CAAC;YACrD,IAAI,YAAY,CAAC;YACjB,IAAI,SAA8B,CAAC;YACnC,KAAK,MAAM,CAAC,IAAI,UAAU,EAAE;gBAC1B,IAAI,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,CAAC,CAAC,QAAQ,EAAE,QAAQ,CAAC,EAAE;oBACvD,YAAY,GAAG,WAAW,CAAC,kCAAkC,CAAC,OAAO,EAAE,OAAO,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC;oBACzF,IAAI,YAAY,KAAK,SAAS,EAAE;wBAC9B,SAAS,GAAG,OAAO,CAAC,WAAW,CAAC,YAAY,EAAE,OAAO,EAAE,SAAS,CAAC,CAAC;wBAClE,IAAI,SAAS,CAAC,mBAAmB,CAAC,CAAC,CAAC,KAAK,CAAC;+BACrC,IAAI,CAAC,cAAc,CAAC,QAAQ,EAAE,YAAY,EAAE,QAAQ,CAAC,EAAE;4BAC1D,IAAI,CAAC,6BAA6B,CAAC,YAAY,EAAE,GAAG,EAAE,UAAU,EAAE,UAAU,EAC1E,CAAC,CAAC,QAAQ,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAAC,CAAC;yBAEpC;qBACF;iBACF;aACF;SACF;IACH,CAAC;IAED,qCAAqC;IACrC,yEAAyE;IACzE,wEAAwE;IACxE,oDAAoD;IAC5C,qBAAqB,CAC3B,GAAU,EACV,OAAgB,EAChB,GAAU,EACV,OAAgB,EAChB,QAAiB;QAEjB,MAAM,YAAY,GAAG,GAAG,CAAC,sBAAsB,CAAC,GAAG,CAAC,CAAC;QACrD,IAAI,YAAY,KAAK,SAAS,EAAE;YAC9B,MAAM,cAAc,GAAa,EAAE,CAAC;YACpC,MAAM,aAAa,GAAa,EAAE,CAAC;YACnC,cAAc,CAAC,6CAA6C,CAC1D,YAAY,CAAC,MAAM,CAAC,CAAC,EAAE,YAAY,CAAC,MAAM,CAAC,CAAC,EAAE,GAAG,EACjD,YAAY,CAAC,OAAO,CAAC,CAAC,EAAE,YAAY,CAAC,OAAO,CAAC,CAAC,EAAE,GAAG,EACnD,YAAY,CAAC,QAAQ,CAAC,CAAC,EAAE,YAAY,CAAC,QAAQ,CAAC,CAAC,EAAE,GAAG,EACrD,cAAc,EAAE,aAAa,CAAC,CAAC;YACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,cAAc,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBAC9C,MAAM,SAAS,GAAG,GAAG,CAAC,KAAK,CAAC,+BAA+B,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC9E,MAAM,SAAS,GAAG,GAAG,CAAC,KAAK,CAAC,+BAA+B,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC/E,4FAA4F;gBAC5F,IAAI,IAAI,CAAC,cAAc,CAAC,OAAO,EAAE,SAAS,EAAE,OAAO,CAAC,IAAI,IAAI,CAAC,cAAc,CAAC,OAAO,EAAE,SAAS,EAAE,OAAO,CAAC,EAAE;oBACxG,IAAI,CAAC,6BAA6B,CAAC,SAAS,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EACrD,SAAS,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAAC,CAAC;iBACnC;aACF;SACF;IACH,CAAC;IACD,sCAAsC;IACtC,6EAA6E;IAC7E,+BAA+B;IACvB,cAAc,CACpB,GAAU,EACV,OAAgB,EAChB,GAAU,EACV,OAAgB,EAChB,QAAiB;QAEjB,mCAAmC;QACnC,yDAAyD;QACzD,0DAA0D;QAC1D,qDAAqD;QACrD,kEAAkE;QAClE,MAAM,MAAM,GAAG,4BAA4B,CAAC,MAAM,CAAC,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,mBAAmB,CAAC,CAAC;QACxF,MAAM,MAAM,GAAG,4BAA4B,CAAC,MAAM,CAAC,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,mBAAmB,CAAC,CAAC;QACxF,IAAI,MAAM,KAAK,SAAS,IAAI,MAAM,KAAK,SAAS;YAC9C,OAAO;QAET,IAAI,MAAM,CAAC,YAAY,EAAE,CAAC,YAAY,CAAC,MAAM,CAAC,YAAY,EAAE,CAAC,EAAE;YAC7D,IAAI,MAAM,CAAC,cAAc,CAAC,MAAM,CAAC,YAAY,EAAE,CAAC,IAAI,MAAM,CAAC,cAAc,CAAC,MAAM,CAAC,YAAY,EAAE,CAAC,EAAE;gBAChG,eAAe;gBACf,IAAI,CAAC,qBAAqB,CAAC,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,QAAQ,CAAC,CAAC;aAClE;SACF;aAAM;YACL,MAAM,UAAU,GAA0B,EAAE,CAAC;YAC7C,GAAG,CAAC,6BAA6B,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;YACtD,MAAM,UAAU,GAA0B,EAAE,CAAC;YAC7C,GAAG,CAAC,6BAA6B,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;YACtD,KAAK,MAAM,OAAO,IAAI,UAAU,EAAE;gBAChC,KAAK,MAAM,OAAO,IAAI,UAAU,EAAE;oBAChC,IAAI,OAAO,CAAC,KAAK,CAAC,aAAa,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE;wBAC9C,IAAI,IAAI,CAAC,cAAc,CAAC,OAAO,EAAE,OAAO,CAAC,QAAQ,EAAE,OAAO,CAAC;+BACtD,IAAI,CAAC,cAAc,CAAC,OAAO,EAAE,OAAO,CAAC,QAAQ,EAAE,OAAO,CAAC,EAAE;4BAC5D,IAAI,CAAC,6BAA6B,CAAC,OAAO,CAAC,QAAQ,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAC5D,OAAO,CAAC,QAAQ,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,QAAQ,CAAC,CAAC;yBAC1C;qBACF;iBACF;aACF;SACF;IACH,CAAC;IACD,sCAAsC;IACtC,6EAA6E;IAC7E,+BAA+B;IACvB,yBAAyB,CAC/B,CAAQ,EACR,QAAiB,EACjB,IAAoB,EACpB,QAAiB,EACjB,SAAkB;QAElB;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;WA2DG;IACL,CAAC;IACD;;;;;;;;;OASG;IACH;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;iBAiGa;IAEb,sCAAsC;IACtC,6EAA6E;IAC7E,+BAA+B;IACvB,oCAAoC,CAC1C,QAA4B,EAC5B,QAA4B,EAC5B,SAAkB;QAClB;;;;;;;;;;;;;;;;;;;;;;;;;qBAyBa;IACf,CAAC;IAED;;;;OAIG;IACH;;;;;;;;;;;;;;;;;;;;;;OAsBG;IACH,qDAAqD;IACrD,6EAA6E;IAC7E,+BAA+B;IACvB,2BAA2B,CACjC,IAAoB,EACpB,SAAkB,EAClB,QAAiB,EACjB,WAAmB,EACnB,QAAiB,EACjB,WAAmB,EACnB,SAAkB,EAClB,OAAuB,EACvB,QAAiB,EACjB,SAAkB;QAElB;;;;;;;;;;;;;;;;;;;;aAoBK;IACP,CAAC;IAMD,0CAA0C;IACnC,8BAA8B,CAAC,IAAkB,EAAE,QAAiB,EAAE,OAAuB,EAAE,QAAiB,EAAE,SAAkB;QACzI;;;;;;;;;;;;;;;;;;;WAmBG;IACL,CAAC;IAED,gDAAgD;IACzC,wBAAwB,CAAC,GAAkB,EAAE,OAAgB,EAAE,GAAiB,EAAE,OAAgB,EAAE,QAAiB;QAC1H,MAAM,OAAO,GAAG,GAAG,CAAC,SAAS,CAAC;QAC9B,MAAM,OAAO,GAAG,GAAG,CAAC,SAAS,CAAC;QAC9B,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;QACrD,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;QACrD,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAC7B,IAAI,IAAI,GAAG,CAAC,EAAE;YACZ,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YAC7B,IAAI,GAAG,CAAC;YACR,IAAI,GAAG,CAAC;YACR,GAAG,GAAG,GAAG,CAAC;YACV,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;YACxB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE;gBACrE,GAAG,CAAC,OAAO,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;gBACzB,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;gBACf,IAAI,CAAC,sBAAsB,CACzB,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EACjD,GAAG,EAAE,EAAE,KAAK,CAAC,IAAI,OAAO,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,KAAK,IAAI,IAAI,OAAO,EAClF,QAAQ,CAAC,CAAC;aACb;SACF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IAED,4CAA4C;IACrC,oBAAoB,CAAC,IAAW,EAAE,OAAgB,EAAE,GAAiB,EAAE,OAAgB,EAAE,QAAiB;QAC/G,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;QACrD,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;QACrD,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAC7B,IAAI,IAAI,GAAG,CAAC,EAAE;YACZ,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YAC7B,IAAI,GAAG,CAAC;YACR,IAAI,GAAG,CAAC;YACR,GAAG,GAAG,GAAG,CAAC;YACV,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;YACxB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE;gBACrE,GAAG,CAAC,OAAO,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;gBACzB,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;gBACf,IAAI,CAAC,kBAAkB,CACrB,GAAG,EAAE,EAAE,KAAK,CAAC,IAAI,OAAO,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,KAAK,IAAI,IAAI,OAAO,EAClF,IAAI,EAAE,OAAO,EAAE,OAAO,EACtB,CAAC,QAAQ,CAAC,CAAC;aACd;SACF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IAED,2DAA2D;IAC3C,mBAAmB,CAAC,QAAuB;QACzD,IAAI,IAAI,CAAC,UAAU,YAAY,aAAa,EAAE;YAC5C,MAAM,QAAQ,GAAG,IAAI,CAAC,UAAU,CAAC;YACjC,IAAI,CAAC,sBAAsB,CACzB,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EACxF,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EACxF,KAAK,CAAC,CAAC;SACV;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,YAAY,EAAE;YAClD,IAAI,CAAC,wBAAwB,CAAC,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;SAC/F;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,KAAK,EAAE;YAC3C,IAAI,CAAC,kBAAkB,CACrB,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EACxF,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;SACzD;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,cAAc,EAAE;YACpD,IAAI,CAAC,2BAA2B,CAC9B,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,QAAQ,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EACxF,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;SAC1C;IACH,CAAC;IACD,+DAA+D;IAC/C,kBAAkB,CAAC,GAAiB;QAClD,IAAI,IAAI,CAAC,UAAU,YAAY,YAAY,EAAE;YAC3C,MAAM,GAAG,GAAG,IAAI,CAAC,UAAU,CAAC;YAC5B,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;YACrD,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;YACrD,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;YACrD,MAAM,OAAO,GAAG,sBAAsB,CAAC,aAAa,CAAC;YACrD,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;YAC7B,MAAM,IAAI,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;YAC7B,IAAI,IAAI,GAAG,CAAC,IAAI,IAAI,GAAG,CAAC,EAAE;gBACxB,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;gBACxB,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;gBAC7B,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;gBAC7B,IAAI,GAAG,GAAG,GAAG,CAAC;gBACd,IAAI,GAAG,CAAC;gBACR,IAAI,GAAG,CAAC;gBACR,IAAI,GAAG,CAAC;gBACR,MAAM,OAAO,GAAG,IAAI,CAAC,QAAQ,CAAC;gBAC9B,MAAM,OAAO,GAAG,IAAI,CAAC,QAAQ,CAAC;gBAC9B,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;gBACxB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE;oBACrE,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;oBACf,GAAG,GAAG,GAAG,CAAC;oBACV,GAAG,CAAC,OAAO,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;oBACzB,GAAG,CAAC,OAAO,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;oBACxB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,IAAI,EAAE,EAAE,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE;wBACrE,GAAG,CAAC,OAAO,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;wBACzB,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;wBACf,IAAI,CAAC,sBAAsB,CACzB,GAAG,EAAE,EAAE,KAAK,CAAC,IAAI,OAAO,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,KAAK,IAAI,IAAI,OAAO,EAClF,GAAG,EAAE,EAAE,KAAK,CAAC,IAAI,OAAO,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,KAAK,IAAI,IAAI,OAAO,EAClF,KAAK,CAAC,CAAC;qBACV;iBACF;aACF;SACF;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,aAAa,EAAE;YACnD,IAAI,CAAC,wBAAwB,CAAC,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC;SACzF;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,KAAK,EAAE;YAC3C,IAAI,CAAC,oBAAoB,CAAC,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC;SACrF;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,cAAc,EAAE;YACpD,IAAI,CAAC,8BAA8B,CAAC,GAAG,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;SAChG;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,wDAAwD;IACxC,WAAW,CAAC,IAAW;QACrC,IAAI,IAAI,CAAC,UAAU,YAAY,aAAa,EAAE;YAC5C,IAAI,CAAC,kBAAkB,CACrB,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,UAAU,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EAC7G,IAAI,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC;SAC7C;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,YAAY,EAAE;YAClD,IAAI,CAAC,oBAAoB,CAAC,IAAI,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;SACvF;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,KAAK,EAAE;YAC3C,IAAI,CAAC,cAAc,CAAC,IAAI,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;SACjF;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,cAAc,EAAE;YACpD,IAAI,CAAC,yBAAyB,CAAC,IAAI,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;SAC5F;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,iEAAiE;IACjD,oBAAoB,CAAC,KAAqB;QACxD,IAAI,IAAI,CAAC,UAAU,YAAY,aAAa,EAAE;YAC5C,IAAI,CAAC,2BAA2B,CAC9B,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,UAAU,CAAC,SAAS,EAAE,GAAG,EAAE,IAAI,CAAC,QAAQ,EAC7G,KAAK,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC;SAC/B;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,YAAY,EAAE;YAClD,IAAI,CAAC,8BAA8B,CAAC,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC;SACjG;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,KAAK,EAAE;YAC3C,IAAI,CAAC,yBAAyB,CAAC,IAAI,CAAC,UAAU,EAAE,IAAI,CAAC,QAAQ,EAAE,KAAK,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC;SAC5F;aAAM,IAAI,IAAI,CAAC,UAAU,YAAY,kBAAkB,EAAE;YACxD,IAAI,CAAC,oCAAoC,CAAC,KAAK,EAAE,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC;SAC1E;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,4EAA4E;IAC5D,qBAAqB,CAAC,MAAuB;QAC3D;;;;;;;;;;WAUG;QACH,OAAO,SAAS,CAAC;IACnB,CAAC;;AA/sBc,qCAAc,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC;AAshBnC,oCAAa,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;AACjC,oCAAa,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;AACjC,oCAAa,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;AACjC,oCAAa,GAAG,OAAO,CAAC,MAAM,EAAE,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\n\r\n/** @packageDocumentation\r\n * @module Curve\r\n */\r\n\r\nimport { BSplineCurve3d, BSplineCurve3dBase } from \"../bspline/BSplineCurve\";\r\n// import { BezierCurveBase } from \"../bspline/BezierCurveBase\";\r\n// import { BezierCurve3dH } from \"../bspline/BezierCurve3dH\";\r\n// import { UnivariateBezier } from \"../numerics/BezierPolynomials\";\r\nimport { BSplineCurve3dH } from \"../bspline/BSplineCurve3dH\";\r\nimport { Geometry } from \"../Geometry\";\r\nimport { NullGeometryHandler } from \"../geometry3d/GeometryHandler\";\r\nimport { Plane3dByOriginAndUnitNormal } from \"../geometry3d/Plane3dByOriginAndUnitNormal\";\r\n// import { Arc3d } from \"./Arc3d\";\r\nimport { Vector2d } from \"../geometry3d/Point2dVector2d\";\r\n// import { XYAndZ } from \"../geometry3d/XYZProps\";\r\nimport { Point3d, Vector3d } from \"../geometry3d/Point3dVector3d\";\r\n// import { LineString3d } from \"./LineString3d\";\r\nimport { SmallSystem, TrigPolynomial } from \"../numerics/Polynomials\";\r\n// import { Point4d } from \"../geometry4d/Point4d\";\r\n// import { Transform } from \"../geometry3d/Transform\";\r\n// import { Matrix3d } from \"../geometry3d/Matrix3d\";\r\nimport { Arc3d } from \"./Arc3d\";\r\n// import { Range3d } from \"../geometry3d/Range\";\r\nimport { CurveLocationDetailArrayPair } from \"./CurveCurveIntersectXY\";\r\nimport { CurveIntervalRole, CurveLocationDetail } from \"./CurveLocationDetail\";\r\nimport { CurvePrimitive } from \"./CurvePrimitive\";\r\nimport { GeometryQuery } from \"./GeometryQuery\";\r\nimport { LineSegment3d } from \"./LineSegment3d\";\r\nimport { LineString3d } from \"./LineString3d\";\r\n\r\n// cspell:word XYRR\r\n/**\r\n * * Handler class for XYZ intersections.\r\n * * Instances are initialized and called from CurveCurve.\r\n * * Constructor is told two geometry items A and B\r\n *   * geometryB is saved for later reference\r\n *   * type-specific handler methods will \"see\" geometry A repeatedly.\r\n *   * Hence geometryA is NOT saved by the constructor.\r\n * @internal\r\n */\r\nexport class CurveCurveIntersectXYZ extends NullGeometryHandler {\r\n  // private geometryA: GeometryQuery;  // nb never used -- passed through handlers.\r\n  private _extendA: boolean;\r\n  private _geometryB: GeometryQuery;\r\n  private _extendB: boolean;\r\n  private _results!: CurveLocationDetailArrayPair;\r\n  private reinitialize() {\r\n    this._results = new CurveLocationDetailArrayPair();\r\n  }\r\n  /**\r\n   *\r\n   * @param _geometryA first curve for intersection.  This is NOT saved.\r\n   * @param extendA flag to enable using extension of geometryA.\r\n   * @param geometryB second curve for intersection.  Saved for reference by specific handler methods.\r\n   * @param extendB flag for extension of geometryB.\r\n   */\r\n  public constructor(_geometryA: GeometryQuery, extendA: boolean, geometryB: GeometryQuery, extendB: boolean) {\r\n    super();\r\n    // this.geometryA = _geometryA;\r\n    this._extendA = extendA;\r\n    this._geometryB = geometryB;\r\n    this._extendB = extendB;\r\n    this.reinitialize();\r\n  }\r\n  /**\r\n   * * Return the results structure for the intersection calculation.\r\n   * @param reinitialize if true, a new results structure is created for use by later calls.\r\n   *\r\n   */\r\n  public grabResults(reinitialize: boolean = false): CurveLocationDetailArrayPair {\r\n    const result = this._results;\r\n    if (reinitialize)\r\n      this.reinitialize();\r\n    return result;\r\n  }\r\n\r\n  private static _workVector2dA = Vector2d.create();\r\n\r\n  private acceptFraction(extend0: boolean, fraction: number, extend1: boolean) {\r\n    if (!extend0 && fraction < 0.0)\r\n      return false;\r\n    if (!extend1 && fraction > 1.0)\r\n      return false;\r\n    return true;\r\n  }\r\n\r\n  /** compute intersection of two line segments.\r\n   * filter by extension rules.\r\n   * reject if evaluated points do not match coordinates (e.g. close approach point)\r\n   * record with fraction mapping.\r\n   */\r\n  private recordPointWithLocalFractions(\r\n    localFractionA: number,\r\n    cpA: CurvePrimitive,\r\n    fractionA0: number,\r\n    fractionA1: number,\r\n    localFractionB: number,   // Computed intersection fraction\r\n    cpB: CurvePrimitive,\r\n    fractionB0: number,\r\n    fractionB1: number,\r\n    reversed: boolean,\r\n  ) {\r\n    const globalFractionA = Geometry.interpolate(fractionA0, localFractionA, fractionA1);\r\n    const globalFractionB = Geometry.interpolate(fractionB0, localFractionB, fractionB1);\r\n    // ignore duplicate of most recent point .  ..\r\n    const numPrevious = this._results.dataA.length;\r\n    if (numPrevious > 0) {\r\n      const topFractionA = this._results.dataA[numPrevious - 1].fraction;\r\n      const topFractionB = this._results.dataB[numPrevious - 1].fraction;\r\n      if (reversed) {\r\n        if (Geometry.isAlmostEqualNumber(topFractionA, globalFractionB) && Geometry.isAlmostEqualNumber(topFractionB, globalFractionA))\r\n          return;\r\n      } else {\r\n        if (Geometry.isAlmostEqualNumber(topFractionA, globalFractionA) && Geometry.isAlmostEqualNumber(topFractionB, globalFractionB))\r\n          return;\r\n      }\r\n    }\r\n    const pointA = cpA.fractionToPoint(globalFractionA);\r\n    const pointB = cpB.fractionToPoint(globalFractionB);\r\n    if (!pointA.isAlmostEqualMetric(pointB))\r\n      return;\r\n\r\n    const detailA = CurveLocationDetail.createCurveFractionPoint(cpA,\r\n      globalFractionA, pointA);\r\n    detailA.setIntervalRole(CurveIntervalRole.isolated);\r\n    const detailB = CurveLocationDetail.createCurveFractionPoint(cpB,\r\n      globalFractionB, pointB);\r\n    detailB.setIntervalRole(CurveIntervalRole.isolated);\r\n    if (reversed) {\r\n      this._results.dataA.push(detailB);\r\n      this._results.dataB.push(detailA);\r\n    } else {\r\n      this._results.dataA.push(detailA);\r\n      this._results.dataB.push(detailB);\r\n    }\r\n  }\r\n  /** compute intersection of two line segments.\r\n   * filter by extension rules.\r\n   * record with fraction mapping.\r\n   */\r\n  private computeSegmentSegment3D(\r\n    cpA: CurvePrimitive,\r\n    extendA0: boolean,\r\n    pointA0: Point3d,\r\n    fractionA0: number,\r\n    pointA1: Point3d,\r\n    fractionA1: number,\r\n    extendA1: boolean,\r\n    cpB: CurvePrimitive,\r\n    extendB0: boolean,\r\n    pointB0: Point3d,\r\n    fractionB0: number,\r\n    pointB1: Point3d,\r\n    fractionB1: number,\r\n    extendB1: boolean,\r\n    reversed: boolean,\r\n  ) {\r\n    const uv = CurveCurveIntersectXYZ._workVector2dA;\r\n    if (SmallSystem.lineSegment3dClosestApproachUnbounded(\r\n      pointA0, pointA1,\r\n      pointB0, pointB1, uv)\r\n      && this.acceptFraction(extendA0, uv.x, extendA1)\r\n      && this.acceptFraction(extendB0, uv.y, extendB1)\r\n    ) {\r\n      this.recordPointWithLocalFractions(uv.x, cpA, fractionA0, fractionA1, uv.y, cpB, fractionB0, fractionB1, reversed);\r\n    }\r\n  }\r\n  // Caller accesses data from a line segment and passes to here.\r\n  // (The line segment in question might be (a) a full line segment or (b) a fragment within a linestring.  The fraction and extend parameters\r\n  // allow all combinations to be passed in)\r\n  // This method applies transform.\r\n  private dispatchSegmentSegment(\r\n    cpA: CurvePrimitive,\r\n    extendA0: boolean,\r\n    pointA0: Point3d,\r\n    fractionA0: number,\r\n    pointA1: Point3d,\r\n    fractionA1: number,\r\n    extendA1: boolean,\r\n    cpB: CurvePrimitive,\r\n    extendB0: boolean,\r\n    pointB0: Point3d,\r\n    fractionB0: number,\r\n    pointB1: Point3d,\r\n    fractionB1: number,\r\n    extendB1: boolean,\r\n    reversed: boolean,\r\n  ) {\r\n    this.computeSegmentSegment3D(\r\n      cpA, extendA0, pointA0, fractionA0, pointA1, fractionA1, extendA1,\r\n      cpB, extendB0, pointB0, fractionB0, pointB1, fractionB1, extendB1,\r\n      reversed);\r\n  }\r\n  /**\r\n   * Create a plane whose normal is the \"better\" cross product: `vectorA.crossProduct(vectorB)` or `vectorA.crossProduct(vectorC)`\r\n   * * The heuristic for \"better\" is:\r\n   *   * first choice is cross product with `vectorB`, if `vectorA` and `vectorB` are sufficiently far from parallel (or anti-parallel).\r\n   *   * otherwise use vectorC\r\n   * @param origin plane origin\r\n   * @param vectorA vector which must be in the plane.\r\n   * @param cosineValue largest cosine of the angle theta between vectorA and vectorB to prefer their cross product,\r\n   *                    e.g. passing 0.94 ~ cos(20deg) will switch to using vectorC in the cross product if theta < ~20deg or theta > ~160deg.\r\n   * @param vectorB first candidate for additional in-plane vector\r\n   * @param vectorC second candidate for additional in-plane vector\r\n   */\r\n  public createPlaneWithPreferredPerpendicular(origin: Point3d, vectorA: Vector3d, cosineValue: number, vectorB: Vector3d, vectorC: Vector3d): Plane3dByOriginAndUnitNormal | undefined {\r\n    cosineValue = Geometry.restrictToInterval(Math.abs(cosineValue), 0.0, 1.0 - Geometry.smallFraction);\r\n    const dotAA = vectorA.magnitudeSquared();\r\n    const dotBB = vectorB.magnitudeSquared();\r\n    const dotAB = Math.abs(vectorA.dotProduct(vectorB));\r\n    const cross = vectorA.unitCrossProduct(dotAB * dotAB <= cosineValue * cosineValue * dotAA * dotBB ? vectorB : vectorC);\r\n    if (cross)\r\n      return Plane3dByOriginAndUnitNormal.create(origin, cross);\r\n    return undefined;\r\n  }\r\n  // Caller accesses data from a linestring or segment and passes it here.\r\n  // (The line segment in question might be (a) a full line segment or (b) a fragment within a linestring.  The fraction and extend parameters\r\n  // allow all combinations to be passed in)\r\n  private dispatchSegmentArc(\r\n    cpA: CurvePrimitive,\r\n    extendA0: boolean,\r\n    pointA0: Point3d,\r\n    fractionA0: number,\r\n    pointA1: Point3d,\r\n    fractionA1: number,\r\n    extendA1: boolean,\r\n    arc: Arc3d,\r\n    extendB0: boolean,\r\n    extendB1: boolean,\r\n    reversed: boolean,\r\n  ) {\r\n    const lineVector = Vector3d.createStartEnd(pointA0, pointA1);\r\n    const plane = this.createPlaneWithPreferredPerpendicular(pointA0, lineVector, 0.94, arc.perpendicularVector, arc.vector0);\r\n    if (plane !== undefined) {\r\n      const candidates: CurveLocationDetail[] = [];\r\n      arc.appendPlaneIntersectionPoints(plane, candidates);\r\n      let lineFraction;\r\n      let linePoint: Point3d | undefined;\r\n      for (const c of candidates) {\r\n        if (this.acceptFraction(extendB0, c.fraction, extendB1)) {\r\n          lineFraction = SmallSystem.lineSegment3dClosestPointUnbounded(pointA0, pointA1, c.point);\r\n          if (lineFraction !== undefined) {\r\n            linePoint = pointA0.interpolate(lineFraction, pointA1, linePoint);\r\n            if (linePoint.isAlmostEqualMetric(c.point)\r\n              && this.acceptFraction(extendA0, lineFraction, extendA1)) {\r\n              this.recordPointWithLocalFractions(lineFraction, cpA, fractionA0, fractionA1,\r\n                c.fraction, arc, 0, 1, reversed);\r\n\r\n            }\r\n          }\r\n        }\r\n      }\r\n    }\r\n  }\r\n\r\n  // Caller promises arcs are coplanar.\r\n  // Passes \"other\" as {center, vector0, vector90} in local xy space of cpA\r\n  // Solves the arc-arc equations for that local ellipse with unit circle.\r\n  // Solution fractions map directly to original arcs.\r\n  private dispatchArcArcInPlane(\r\n    cpA: Arc3d,\r\n    extendA: boolean,\r\n    cpB: Arc3d,\r\n    extendB: boolean,\r\n    reversed: boolean,\r\n  ) {\r\n    const otherVectors = cpA.otherArcAsLocalVectors(cpB);\r\n    if (otherVectors !== undefined) {\r\n      const ellipseRadians: number[] = [];\r\n      const circleRadians: number[] = [];\r\n      TrigPolynomial.solveUnitCircleHomogeneousEllipseIntersection(\r\n        otherVectors.center.x, otherVectors.center.y, 1.0,\r\n        otherVectors.vector0.x, otherVectors.vector0.y, 0.0,\r\n        otherVectors.vector90.x, otherVectors.vector90.y, 0.0,\r\n        ellipseRadians, circleRadians);\r\n      for (let i = 0; i < ellipseRadians.length; i++) {\r\n        const fractionA = cpA.sweep.radiansToSignedPeriodicFraction(circleRadians[i]);\r\n        const fractionB = cpA.sweep.radiansToSignedPeriodicFraction(ellipseRadians[i]);\r\n        // hm .. do we really need to check the fractions?  We know they are internal to the beziers\r\n        if (this.acceptFraction(extendA, fractionA, extendA) && this.acceptFraction(extendB, fractionB, extendB)) {\r\n          this.recordPointWithLocalFractions(fractionA, cpA, 0, 1,\r\n            fractionB, cpB, 0, 1, reversed);\r\n        }\r\n      }\r\n    }\r\n  }\r\n  // Caller accesses data from two arcs.\r\n  // Selects the best conditioned arc (in xy parts) as \"circle after inversion\"\r\n  // Solves the arc-arc equations\r\n  private dispatchArcArc(\r\n    cpA: Arc3d,\r\n    extendA: boolean,\r\n    cpB: Arc3d,\r\n    extendB: boolean,\r\n    reversed: boolean,\r\n  ) {\r\n    // If arcs are in different planes:\r\n    // 1) Intersect each plane with the other arc (quadratic)\r\n    // 2) accept points that appear in both intersection sets.\r\n    // If arcs are in parallel planes -- no intersections\r\n    // If arcs are in the same plane -- xy intersection in that plane.\r\n    const planeA = Plane3dByOriginAndUnitNormal.create(cpA.center, cpA.perpendicularVector);\r\n    const planeB = Plane3dByOriginAndUnitNormal.create(cpB.center, cpB.perpendicularVector);\r\n    if (planeA === undefined || planeB === undefined)\r\n      return;\r\n\r\n    if (planeA.getNormalRef().isParallelTo(planeB.getNormalRef())) {\r\n      if (planeA.isPointInPlane(planeB.getOriginRef()) && planeB.isPointInPlane(planeA.getOriginRef())) {\r\n        // coplanar !!!\r\n        this.dispatchArcArcInPlane(cpA, extendA, cpB, extendB, reversed);\r\n      }\r\n    } else {\r\n      const arcBPoints: CurveLocationDetail[] = [];\r\n      cpB.appendPlaneIntersectionPoints(planeA, arcBPoints);\r\n      const arcAPoints: CurveLocationDetail[] = [];\r\n      cpA.appendPlaneIntersectionPoints(planeB, arcAPoints);\r\n      for (const detailB of arcBPoints) {\r\n        for (const detailA of arcAPoints) {\r\n          if (detailA.point.isAlmostEqual(detailB.point)) {\r\n            if (this.acceptFraction(extendA, detailA.fraction, extendA)\r\n              && this.acceptFraction(extendB, detailB.fraction, extendB)) {\r\n              this.recordPointWithLocalFractions(detailA.fraction, cpA, 0, 1,\r\n                detailB.fraction, cpB, 0, 1, reversed);\r\n            }\r\n          }\r\n        }\r\n      }\r\n    }\r\n  }\r\n  // Caller accesses data from two arcs.\r\n  // Selects the best conditioned arc (in xy parts) as \"circle after inversion\"\r\n  // Solves the arc-arc equations\r\n  private dispatchArcBsplineCurve3d(\r\n    _: Arc3d,\r\n    _extendA: boolean,\r\n    _cpB: BSplineCurve3d,\r\n    _extendB: boolean,\r\n    _reversed: boolean,\r\n  ) {\r\n    /*\r\n    // Arc: X = C + cU + sV\r\n    // implicitize the arc as viewed.  This \"3d\" matrix is homogeneous \"XYW\" not \"xyz\"\r\n    let matrixA: Matrix3d;\r\n    if (this._worldToLocalPerspective) {\r\n      const dataA = cpA.toTransformedPoint4d(this._worldToLocalPerspective);\r\n      matrixA = Matrix3d.createColumnsXYW(dataA.vector0, dataA.vector0.w, dataA.vector90, dataA.vector90.w, dataA.center, dataA.center.w);\r\n    } else {\r\n      const dataA = cpA.toTransformedVectors(this._worldToLocalAffine);\r\n      matrixA = Matrix3d.createColumnsXYW(dataA.vector0, 0, dataA.vector90, 0, dataA.center, 1);\r\n    }\r\n    // The worldToLocal has moved the arc vectors into local space.\r\n    // matrixA captures the xyw parts (ignoring z)\r\n    // for any point in world space,\r\n    // THIS CODE ONLY WORKS FOR\r\n    const matrixAInverse = matrixA.inverse();\r\n    if (matrixAInverse) {\r\n      const orderF = cpB.order; // order of the beziers for simple coordinates\r\n      const orderG = 2 * orderF - 1;  // order of the (single) bezier for squared coordinates.\r\n      const coffF = new Float64Array(orderF);\r\n      const univariateBezierG = new UnivariateBezier(orderG);\r\n      const axx = matrixAInverse.at(0, 0); const axy = matrixAInverse.at(0, 1); const axz = 0.0; const axw = matrixAInverse.at(0, 2);\r\n      const ayx = matrixAInverse.at(1, 0); const ayy = matrixAInverse.at(1, 1); const ayz = 0.0; const ayw = matrixAInverse.at(1, 2);\r\n      const awx = matrixAInverse.at(2, 0); const awy = matrixAInverse.at(2, 1); const awz = 0.0; const aww = matrixAInverse.at(2, 2);\r\n\r\n      if (matrixAInverse) {\r\n        let bezier: BezierCurve3dH | undefined;\r\n        for (let spanIndex = 0; ; spanIndex++) {\r\n          bezier = cpB.getSaturatedBezierSpan3dH(spanIndex, bezier);\r\n          if (!bezier) break;\r\n          if (this._worldToLocalPerspective)\r\n            bezier.tryMultiplyMatrix4dInPlace(this._worldToLocalPerspective);\r\n          else if (this._worldToLocalAffine)\r\n            bezier.tryTransformInPlace(this._worldToLocalAffine);\r\n          univariateBezierG.zero();\r\n          bezier.poleProductsXYZW(coffF, axx, axy, axz, axw);\r\n          univariateBezierG.addSquaredSquaredBezier(coffF, 1.0);\r\n          bezier.poleProductsXYZW(coffF, ayx, ayy, ayz, ayw);\r\n          univariateBezierG.addSquaredSquaredBezier(coffF, 1.0);\r\n          bezier.poleProductsXYZW(coffF, awx, awy, awz, aww);\r\n          univariateBezierG.addSquaredSquaredBezier(coffF, -1.0);\r\n          const roots = univariateBezierG.roots(0.0, true);\r\n          if (roots) {\r\n            for (const root of roots) {\r\n              const fractionB = bezier.fractionToParentFraction(root);\r\n              // The univariate bezier (which has been transformed by the view transform) evaluates into xyw space\r\n              const bcurvePoint4d = bezier.fractionToPoint4d(root);\r\n              const c = bcurvePoint4d.dotProductXYZW(axx, axy, axz, axw);\r\n              const s = bcurvePoint4d.dotProductXYZW(ayx, ayy, ayz, ayw);\r\n              const arcFraction = cpA.sweep.radiansToSignedPeriodicFraction(Math.atan2(s, c));\r\n              if (this.acceptFraction(extendA, arcFraction, extendA) && this.acceptFraction(extendB, fractionB, extendB)) {\r\n                this.recordPointWithLocalFractions(arcFraction, cpA, 0, 1,\r\n                  fractionB, cpB, 0, 1, reversed);\r\n              }\r\n            }\r\n          }\r\n        }\r\n      }\r\n    }\r\n     */\r\n  }\r\n  /*\r\n  // apply the transformation to bezier curves. optionally construct ranges.\r\n  private transformBeziers(beziers: BezierCurve3dH[]) {\r\n    if (this._worldToLocalAffine) {\r\n      for (const bezier of beziers) bezier.tryTransformInPlace(this._worldToLocalAffine);\r\n    } else if (this._worldToLocalPerspective) {\r\n      for (const bezier of beziers) bezier.tryMultiplyMatrix4dInPlace(this._worldToLocalPerspective);\r\n    }\r\n  }\r\n   */\r\n  /*\r\n  private getRanges(beziers: BezierCurveBase[]): Range3d[] {\r\n    const ranges: Range3d[] = [];\r\n    ranges.length = 0;\r\n    for (const b of beziers) {\r\n      ranges.push(b.range());\r\n    }\r\n    return ranges;\r\n  }\r\n  private _xyzwA0?: Point4d;\r\n  private _xyzwA1?: Point4d;\r\n  private _xyzwPlane?: Point4d;\r\n  private _xyzwB?: Point4d;\r\n\r\n  private dispatchBezierBezierStrokeFirst(\r\n    bezierA: BezierCurve3dH,\r\n    bcurveA: BSplineCurve3dBase,\r\n    strokeCountA: number,\r\n    bezierB: BezierCurve3dH,\r\n    bcurveB: BSplineCurve3dBase,\r\n    _strokeCountB: number,\r\n    univariateBezierB: UnivariateBezier,  // caller-allocated for univariate coefficients.\r\n    reversed: boolean) {\r\n    if (!this._xyzwA0) this._xyzwA0 = Point4d.create();\r\n    if (!this._xyzwA1) this._xyzwA1 = Point4d.create();\r\n    if (!this._xyzwPlane) this._xyzwPlane = Point4d.create();\r\n    if (!this._xyzwB) this._xyzwB = Point4d.create();\r\n    /-*\r\n\r\n              const roots = univariateBezierG.roots(0.0, true);\r\n              if (roots) {\r\n                for (const root of roots) {\r\n                  const fractionB = bezier.fractionToParentFraction(root);\r\n                  // The univariate bezier (which has been transformed by the view transform) evaluates into xyw space\r\n                  const bcurvePoint4d = bezier.fractionToPoint4d(root);\r\n                  const c = bcurvePoint4d.dotProductXYZW(axx, axy, axz, axw);\r\n                  const s = bcurvePoint4d.dotProductXYZW(ayx, ayy, ayz, ayw);\r\n                  const arcFraction = cpA.sweep.radiansToSignedPeriodicFraction(Math.atan2(s, c));\r\n                  if (this.acceptFraction(extendA, arcFraction, extendA) && this.acceptFraction(extendB, fractionB, extendB)) {\r\n                    this.recordPointWithLocalFractions(arcFraction, cpA, 0, 1,\r\n                      fractionB, cpB, 0, 1, reversed);\r\n                  }\r\n                }\r\n   *-/\r\n    bezierA.fractionToPoint4d(0.0, this._xyzwA0);\r\n    let f0 = 0.0;\r\n    let f1 = 1.0;\r\n    const intervalTolerance = 1.0e-5;\r\n    const df = 1.0 / strokeCountA;\r\n    for (let i = 1; i <= strokeCountA; i++ , f0 = f1, this._xyzwA0.setFrom(this._xyzwA1)) {\r\n      f1 = i * df;\r\n      bezierA.fractionToPoint4d(f1, this._xyzwA1);\r\n      Point4d.createPlanePointPointZ(this._xyzwA0, this._xyzwA1, this._xyzwPlane);\r\n      bezierB.poleProductsXYZW(univariateBezierB.coffs, this._xyzwPlane.x, this._xyzwPlane.y, this._xyzwPlane.z, this._xyzwPlane.w);\r\n      let errors = 0;\r\n      const roots = univariateBezierB.roots(0.0, true);\r\n      if (roots)\r\n        for (const r of roots) {\r\n          const bezierBFraction = r;\r\n          bezierB.fractionToPoint4d(bezierBFraction, this._xyzwB);\r\n          const segmentAFraction = SmallSystem.lineSegment3dHXYClosestPointUnbounded(this._xyzwA0, this._xyzwA1, this._xyzwB);\r\n          if (segmentAFraction && Geometry.isIn01WithTolerance(segmentAFraction, intervalTolerance)) {\r\n            const bezierAFraction = Geometry.interpolate(f0, segmentAFraction, f1);\r\n            /*- TODO implement newton search\r\n            const xyMatchingFunction = new BezierBezierIntersectionXYRRToRRD(bezierA, bezierB);\r\n            const newtonSearcher = new Newton2dUnboundedWithDerivative(xyMatchingFunction);\r\n            newtonSearcher.setUV(bezierAFraction, bezierBFraction);\r\n            if (newtonSearcher.runIterations()) {\r\n              bezierAFraction = newtonSearcher.getU();\r\n              bezierBFraction = newtonSearcher.getV();\r\n            }\r\n             *-/\r\n              // We have a near intersection at fractions on the two beziers !!!\r\n              // Iterate on the curves for a true intersection ....\r\n              // NEEDS WORK -- just accept . . .\r\n            const bcurveAFraction = bezierA.fractionToParentFraction(bezierAFraction);\r\n            const bcurveBFraction = bezierB.fractionToParentFraction(bezierBFraction);\r\n            const xyzA0 = bezierA.fractionToPoint(bezierAFraction);\r\n            const xyzA1 = bcurveA.fractionToPoint(bcurveAFraction);\r\n            const xyzB0 = bezierB.fractionToPoint(bezierBFraction);\r\n            const xyzB1 = bcurveB.fractionToPoint(bcurveBFraction);\r\n            if (!xyzA0.isAlmostEqualXY(xyzA1))\r\n              errors++;\r\n            if (!xyzB0.isAlmostEqualXY(xyzB1))\r\n              errors++;\r\n            if (errors > 0 && !xyzA0.isAlmostEqual(xyzB0))\r\n              errors++;\r\n            if (errors > 0 && !xyzA1.isAlmostEqual(xyzB1))\r\n              errors++;\r\n            if (this.acceptFraction(false, bcurveAFraction, false) && this.acceptFraction(false, bcurveBFraction, false)) {\r\n              this.recordPointWithLocalFractions(bcurveAFraction, bcurveA, 0, 1,\r\n                bcurveBFraction, bcurveB, 0, 1, reversed);\r\n            }\r\n          }\r\n        }\r\n    }\r\n  }\r\n             */\r\n\r\n  // Caller accesses data from two arcs.\r\n  // Selects the best conditioned arc (in xy parts) as \"circle after inversion\"\r\n  // Solves the arc-arc equations\r\n  private dispatchBSplineCurve3dBSplineCurve3d(\r\n    _bcurveA: BSplineCurve3dBase,\r\n    _bcurveB: BSplineCurve3dBase,\r\n    _reversed: boolean) {\r\n    /*\r\n  const bezierSpanA = bcurveA.collectBezierSpans(true) as BezierCurve3dH[];\r\n  const bezierSpanB = bcurveB.collectBezierSpans(true) as BezierCurve3dH[];\r\n  const numA = bezierSpanA.length;\r\n  const numB = bezierSpanB.length;\r\n  this.transformBeziers(bezierSpanA);\r\n  this.transformBeziers(bezierSpanB);\r\n  const rangeA = this.getRanges(bezierSpanA);\r\n  const rangeB = this.getRanges(bezierSpanB);\r\n  const orderA = bcurveA.order;\r\n  const orderB = bcurveB.order;\r\n  const univariateCoffsA = new UnivariateBezier(orderA);\r\n  const univariateCoffsB = new UnivariateBezier(orderB);\r\n  for (let a = 0; a < numA; a++) {\r\n    for (let b = 0; b < numB; b++) {\r\n      if (rangeA[a].intersectsRangeXY(rangeB[b])) {\r\n        const strokeCountA = bezierSpanA[a].computeStrokeCountForOptions();\r\n        const strokeCountB = bezierSpanB[b].computeStrokeCountForOptions();\r\n        if (strokeCountA < strokeCountB)\r\n          this.dispatchBezierBezierStrokeFirst(bezierSpanA[a], bcurveA, strokeCountA, bezierSpanB[b], bcurveB, strokeCountB, univariateCoffsB, !_reversed);\r\n        else\r\n          this.dispatchBezierBezierStrokeFirst(bezierSpanB[b], bcurveB, strokeCountB, bezierSpanA[a], bcurveA, strokeCountA, univariateCoffsA, _reversed);\r\n      }\r\n    }\r\n  }\r\n               */\r\n  }\r\n\r\n  /**\r\n   * Apply the projection transform (if any) to (xyz, w)\r\n   * @param xyz xyz parts of input point.\r\n   * @param w   weight to use for homogeneous effects\r\n   */\r\n  /*\r\n  private projectPoint(xyz: XYAndZ, w: number = 1.0): Point4d {\r\n    if (this._worldToLocalPerspective)\r\n      return this._worldToLocalPerspective.multiplyPoint3d(xyz, w);\r\n    if (this._worldToLocalAffine)\r\n      return this._worldToLocalAffine.multiplyXYZW(xyz.x, xyz.y, xyz.z, w);\r\n    return Point4d.createFromPointAndWeight(xyz, w);\r\n  }\r\n  private mapNPCPlaneToWorld(npcPlane: Point4d, worldPlane: Point4d) {\r\n    // for NPC pointY, Y^ * H = 0 is \"on\" plane H.  (Hat is transpose)\r\n    // NPC Y is A*X for our transform A and worldPointX.\r\n    // hence (A X)^ * H = 0\r\n    // hence X^ * A^ * H = 0\r\n    // hence K = A^ * H\r\n    if (this._worldToLocalAffine) {\r\n      this._worldToLocalAffine.multiplyTransposeXYZW(npcPlane.x, npcPlane.y, npcPlane.z, npcPlane.w, worldPlane);\r\n    } else if (this._worldToLocalPerspective) {\r\n      this._worldToLocalPerspective.multiplyTransposePoint4d(npcPlane, worldPlane);\r\n    } else {\r\n      npcPlane.clone(worldPlane);\r\n    }\r\n  }\r\n   */\r\n  // Caller accesses data from segment and bsplineCurve\r\n  // Selects the best conditioned arc (in xy parts) as \"circle after inversion\"\r\n  // Solves the arc-arc equations\r\n  private dispatchSegmentBsplineCurve(\r\n    _cpA: CurvePrimitive,\r\n    _extendA0: boolean,\r\n    _pointA0: Point3d,\r\n    _fractionA0: number,\r\n    _pointA1: Point3d,\r\n    _fractionA1: number,\r\n    _extendA1: boolean,\r\n    _bcurve: BSplineCurve3d,\r\n    _extendB: boolean,\r\n    _reversed: boolean,\r\n  ) {\r\n    /*\r\n    const pointA0H = this.projectPoint(pointA0);\r\n    const pointA1H = this.projectPoint(pointA1);\r\n    const planeCoffs = Point4d.createPlanePointPointZ(pointA0H, pointA1H);\r\n    this.mapNPCPlaneToWorld(planeCoffs, planeCoffs);\r\n      // NOW .. we have a plane in world space.  Intersect it with the bspline:\r\n    const intersections: CurveLocationDetail[] = [];\r\n    bcurve.appendPlaneIntersectionPoints(planeCoffs, intersections);\r\n      // intersections has WORLD points with bspline fractions.   (The bspline fractions are all good 0..1 fractions within the spline.)\r\n      // accept those that are within the segment range.\r\n    for (const detail of intersections) {\r\n      const fractionB = detail.fraction;\r\n      const curvePoint = detail.point;\r\n      const curvePointH = this.projectPoint(curvePoint);\r\n      const lineFraction = SmallSystem.lineSegment3dHXYClosestPointUnbounded(pointA0H, pointA1H, curvePointH);\r\n      if (lineFraction !== undefined && this.acceptFraction(extendA0, lineFraction, extendA1) && this.acceptFraction(extendB, fractionB, extendB)) {\r\n        this.recordPointWithLocalFractions(lineFraction, cpA, fractionA0, fractionA1,\r\n          fractionB, bcurve, 0, 1, reversed);\r\n      }\r\n    }\r\n       */\r\n  }\r\n\r\n  private static _workPointAA0 = Point3d.create();\r\n  private static _workPointAA1 = Point3d.create();\r\n  private static _workPointBB0 = Point3d.create();\r\n  private static _workPointBB1 = Point3d.create();\r\n  /** low lever bspline curve -- STUB  .. */\r\n  public dispatchLineStringBSplineCurve(_lsA: LineString3d, _extendA: boolean, _curveB: BSplineCurve3d, _extendB: boolean, _reversed: boolean): any {\r\n    /*\r\n    const numA = lsA.numPoints();\r\n    if (numA > 1) {\r\n      const dfA = 1.0 / (numA - 1);\r\n      let fA0;\r\n      let fA1;\r\n      fA0 = 0.0;\r\n      const pointA0 = CurveCurveIntersectXYZ._workPointA0;\r\n      const pointA1 = CurveCurveIntersectXYZ._workPointA1;\r\n      lsA.pointAt(0, pointA0);\r\n      for (let iA = 1; iA < numA; iA++ , pointA0.setFrom(pointA1), fA0 = fA1) {\r\n        lsA.pointAt(iA, pointA1);\r\n        fA1 = iA * dfA;\r\n        this.dispatchSegmentBsplineCurve(\r\n          lsA, iA === 1 && extendA, pointA0, fA0, pointA1, fA1, (iA + 1) === numA && extendA,\r\n          curveB, extendB, reversed);\r\n      }\r\n    }\r\n    return undefined;\r\n     */\r\n  }\r\n\r\n  /** low lever segment intersect linestring .. */\r\n  public computeSegmentLineString(lsA: LineSegment3d, extendA: boolean, lsB: LineString3d, extendB: boolean, reversed: boolean): any {\r\n    const pointA0 = lsA.point0Ref;\r\n    const pointA1 = lsA.point1Ref;\r\n    const pointB0 = CurveCurveIntersectXYZ._workPointBB0;\r\n    const pointB1 = CurveCurveIntersectXYZ._workPointBB1;\r\n    const numB = lsB.numPoints();\r\n    if (numB > 1) {\r\n      const dfB = 1.0 / (numB - 1);\r\n      let fB0;\r\n      let fB1;\r\n      fB0 = 0.0;\r\n      lsB.pointAt(0, pointB0);\r\n      for (let ib = 1; ib < numB; ib++, pointB0.setFrom(pointB1), fB0 = fB1) {\r\n        lsB.pointAt(ib, pointB1);\r\n        fB1 = ib * dfB;\r\n        this.dispatchSegmentSegment(\r\n          lsA, extendA, pointA0, 0.0, pointA1, 1.0, extendA,\r\n          lsB, ib === 1 && extendB, pointB0, fB0, pointB1, fB1, (ib + 1) === numB && extendB,\r\n          reversed);\r\n      }\r\n    }\r\n    return undefined;\r\n  }\r\n\r\n  /** low lever arc intersect linestring .. */\r\n  public computeArcLineString(arcA: Arc3d, extendA: boolean, lsB: LineString3d, extendB: boolean, reversed: boolean): any {\r\n    const pointB0 = CurveCurveIntersectXYZ._workPointBB0;\r\n    const pointB1 = CurveCurveIntersectXYZ._workPointBB1;\r\n    const numB = lsB.numPoints();\r\n    if (numB > 1) {\r\n      const dfB = 1.0 / (numB - 1);\r\n      let fB0;\r\n      let fB1;\r\n      fB0 = 0.0;\r\n      lsB.pointAt(0, pointB0);\r\n      for (let ib = 1; ib < numB; ib++, pointB0.setFrom(pointB1), fB0 = fB1) {\r\n        lsB.pointAt(ib, pointB1);\r\n        fB1 = ib * dfB;\r\n        this.dispatchSegmentArc(\r\n          lsB, ib === 1 && extendB, pointB0, fB0, pointB1, fB1, (ib + 1) === numB && extendB,\r\n          arcA, extendA, extendA,\r\n          !reversed);\r\n      }\r\n    }\r\n    return undefined;\r\n  }\r\n\r\n  /** double dispatch handler for strongly typed segment.. */\r\n  public override handleLineSegment3d(segmentA: LineSegment3d): any {\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      const segmentB = this._geometryB;\r\n      this.dispatchSegmentSegment(\r\n        segmentA, this._extendA, segmentA.point0Ref, 0.0, segmentA.point1Ref, 1.0, this._extendA,\r\n        segmentB, this._extendB, segmentB.point0Ref, 0.0, segmentB.point1Ref, 1.0, this._extendB,\r\n        false);\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.computeSegmentLineString(segmentA, this._extendA, this._geometryB, this._extendB, false);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchSegmentArc(\r\n        segmentA, this._extendA, segmentA.point0Ref, 0.0, segmentA.point1Ref, 1.0, this._extendA,\r\n        this._geometryB, this._extendB, this._extendB, false);\r\n    } else if (this._geometryB instanceof BSplineCurve3d) {\r\n      this.dispatchSegmentBsplineCurve(\r\n        segmentA, this._extendA, segmentA.point0Ref, 0.0, segmentA.point1Ref, 1.0, this._extendA,\r\n        this._geometryB, this._extendB, false);\r\n    }\r\n  }\r\n  /** double dispatch handler for strongly typed linestring .. */\r\n  public override handleLineString3d(lsA: LineString3d): any {\r\n    if (this._geometryB instanceof LineString3d) {\r\n      const lsB = this._geometryB;\r\n      const pointA0 = CurveCurveIntersectXYZ._workPointAA0;\r\n      const pointA1 = CurveCurveIntersectXYZ._workPointAA1;\r\n      const pointB0 = CurveCurveIntersectXYZ._workPointBB0;\r\n      const pointB1 = CurveCurveIntersectXYZ._workPointBB1;\r\n      const numA = lsA.numPoints();\r\n      const numB = lsB.numPoints();\r\n      if (numA > 1 && numB > 1) {\r\n        lsA.pointAt(0, pointA0);\r\n        const dfA = 1.0 / (numA - 1);\r\n        const dfB = 1.0 / (numB - 1);\r\n        let fA0 = 0.0;\r\n        let fB0;\r\n        let fA1;\r\n        let fB1;\r\n        const extendA = this._extendA;\r\n        const extendB = this._extendB;\r\n        lsA.pointAt(0, pointA0);\r\n        for (let ia = 1; ia < numA; ia++, pointA0.setFrom(pointA1), fA0 = fA1) {\r\n          fA1 = ia * dfA;\r\n          fB0 = 0.0;\r\n          lsA.pointAt(ia, pointA1);\r\n          lsB.pointAt(0, pointB0);\r\n          for (let ib = 1; ib < numB; ib++, pointB0.setFrom(pointB1), fB0 = fB1) {\r\n            lsB.pointAt(ib, pointB1);\r\n            fB1 = ib * dfB;\r\n            this.dispatchSegmentSegment(\r\n              lsA, ia === 1 && extendA, pointA0, fA0, pointA1, fA1, (ia + 1) === numA && extendA,\r\n              lsB, ib === 1 && extendB, pointB0, fB0, pointB1, fB1, (ib + 1) === numB && extendB,\r\n              false);\r\n          }\r\n        }\r\n      }\r\n    } else if (this._geometryB instanceof LineSegment3d) {\r\n      this.computeSegmentLineString(this._geometryB, this._extendB, lsA, this._extendA, true);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.computeArcLineString(this._geometryB, this._extendB, lsA, this._extendA, true);\r\n    } else if (this._geometryB instanceof BSplineCurve3d) {\r\n      this.dispatchLineStringBSplineCurve(lsA, this._extendA, this._geometryB, this._extendB, false);\r\n    }\r\n    return undefined;\r\n  }\r\n  /** double dispatch handler for strongly typed arc .. */\r\n  public override handleArc3d(arc0: Arc3d): any {\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      this.dispatchSegmentArc(\r\n        this._geometryB, this._extendB, this._geometryB.point0Ref, 0.0, this._geometryB.point1Ref, 1.0, this._extendB,\r\n        arc0, this._extendA, this._extendA, true);\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.computeArcLineString(arc0, this._extendA, this._geometryB, this._extendB, false);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchArcArc(arc0, this._extendA, this._geometryB, this._extendB, false);\r\n    } else if (this._geometryB instanceof BSplineCurve3d) {\r\n      this.dispatchArcBsplineCurve3d(arc0, this._extendA, this._geometryB, this._extendB, false);\r\n    }\r\n    return undefined;\r\n  }\r\n  /** double dispatch handler for strongly typed bspline curve.. */\r\n  public override handleBSplineCurve3d(curve: BSplineCurve3d): any {\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      this.dispatchSegmentBsplineCurve(\r\n        this._geometryB, this._extendB, this._geometryB.point0Ref, 0.0, this._geometryB.point1Ref, 1.0, this._extendB,\r\n        curve, this._extendA, true);\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.dispatchLineStringBSplineCurve(this._geometryB, this._extendB, curve, this._extendA, true);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchArcBsplineCurve3d(this._geometryB, this._extendB, curve, this._extendA, true);\r\n    } else if (this._geometryB instanceof BSplineCurve3dBase) {\r\n      this.dispatchBSplineCurve3dBSplineCurve3d(curve, this._geometryB, false);\r\n    }\r\n    return undefined;\r\n  }\r\n  /** double dispatch handler for strongly typed homogeneous bspline curve. */\r\n  public override handleBSplineCurve3dH(_curve: BSplineCurve3dH): any {\r\n    /* NEEDS WORK -- make \"dispatch\" methods tolerant of both 3d and 3dH ...\"easy\" if both present BezierCurve3dH span loaders\r\n    if (this._geometryB instanceof LineSegment3d) {\r\n      this.dispatchSegmentBsplineCurve(\r\n        this._geometryB, this._extendB, this._geometryB.point0Ref, 0.0, this._geometryB.point1Ref, 1.0, this._extendB,\r\n        curve, this._extendA, true);\r\n    } else if (this._geometryB instanceof LineString3d) {\r\n      this.dispatchLineStringBSplineCurve(this._geometryB, this._extendB, curve, this._extendA, true);\r\n    } else if (this._geometryB instanceof Arc3d) {\r\n      this.dispatchArcBsplineCurve3d(this._geometryB, this._extendB, curve, this._extendA, true);\r\n    }\r\n     */\r\n    return undefined;\r\n  }\r\n}\r\n"]}

package/lib/esm/geometry3d/Matrix3d.d.ts

@@ -264,7 +264,7 @@ export declare class Matrix3d implements BeJSONFunctions {
     setIdentity(): void;
     /** Set the matrix to all zeros. */
     setZero(): void;
-    /** Copy contents from another matrix. */
+    /** Copy contents from the `other` matrix. If `other` is undefined, use identity matrix. */
     setFrom(other: Matrix3d | undefined): void;
     /**
      * Return a clone of this matrix.
@@ -348,6 +348,9 @@ export declare class Matrix3d implements BeJSONFunctions {
      * Construct a rigid matrix (orthogonal matrix with +1 determinant) using vectorA and its 2 perpendicular.
      * * If axisOrder is not passed then `AxisOrder = AxisOrder.ZXY` is used as default.
      * * This function internally uses createPerpendicularVectorFavorXYPlane and createRigidFromColumns.
+     * * If you want to rotate a given plane (which contains (0,0,0)) to the xy-plane, pass the normal vector of
+     * your plane into createRigidHeadsUp. The transpose of the returned Matrix3d can be used to rotate your plane
+     * to the xy-plane. If plane does not contain (0,0,0) then the plane is rotated to a plane parallel to the xy-plane.
      * * Visualization can be found at https://www.itwinjs.org/sandbox/SaeedTorabi/2PerpendicularVectorsTo1Vector
      */
     static createRigidHeadsUp(vectorA: Vector3d, axisOrder?: AxisOrder, result?: Matrix3d): Matrix3d;
@@ -370,6 +373,7 @@ export declare class Matrix3d implements BeJSONFunctions {
     /**
      * Replace current rows Ui and Uj with (c*Ui + s*Uj) and (c*Uj - s*Ui).
      * * There is no checking for i,j being 0,1,2.
+     * * The instance matrix A is multiplied in place on the left by a Givens rotation G, resulting in the matrix G*A.
      * @param i first row index. **must be 0,1,2** (unchecked)
      * @param j second row index. **must be 0,1,2** (unchecked)
      * @param c fist coefficient
@@ -379,6 +383,7 @@ export declare class Matrix3d implements BeJSONFunctions {
     /**
      * Replace current columns Ui and Uj with (c*Ui + s*Uj) and (c*Uj - s*Ui).
      * * There is no checking for i,j being 0,1,2.
+     * * The instance matrix A is multiplied in place on the right by a Givens rotation G, resulting in the matrix A*G.
      * * This is used in compute intensive inner loops
      * @param i first row index. **must be 0,1,2** (unchecked)
      * @param j second row index. **must be 0,1,2** (unchecked)
@@ -836,7 +841,7 @@ export declare class Matrix3d implements BeJSONFunctions {
      */
     multiplyInverseXYZAsVector3d(x: number, y: number, z: number, result?: Vector3d): Vector3d | undefined;
     /**
-     * Multiply `matrixInverse * [x,y,z]` and return result as `Point4d` with given weight.
+     * Multiply `matrixInverse * [x,y,z]` and return result as `Point4d` the with given weight as last element.
      * * Equivalent to solving `matrix * result = [x,y,z]` for an unknown `result`.
      * * Result is `undefined` if the matrix is singular (e.g. has parallel columns or a zero magnitude column)
      * @return result as a Point4d with the same weight.
@@ -1068,6 +1073,10 @@ export declare class Matrix3d implements BeJSONFunctions {
      * * column 1 is perpendicular to both. It is the "up" vector on the view plane.
      * * Multiplying the returned matrix times a local (view) vector gives the world vector.
      * * Multiplying transpose of the returned matrix times a world vector gives the local (view) vector.
+     * * If you want to rotate a given plane (which contains (0,0,0)) to the xy-plane, pass coordinates of the normal
+     * vector of your plane into createRigidViewAxesZTowardsEye. The transpose of the returned Matrix3d can be used
+     * to rotate your plane to the xy-plane. If plane does not contain (0,0,0) then the plane is rotated to a plane
+     * parallel to the xy-plane.
      * @param x eye x coordinate
      * @param y eye y coordinate
      * @param z eye z coordinate
@@ -1127,13 +1136,14 @@ export declare class Matrix3d implements BeJSONFunctions {
      * Test if all rows and columns are perpendicular to each other and have equal length.
      * If so, the length (or its negative) is the `scale` factor from a set of `orthonormal axes` to
      * the set of axes created by columns of `this` matrix. Otherwise, returns `undefined`.
+     * @param result optional pre-allocated object to populate and return
      * @returns returns `{ rigidAxes, scale }` where `rigidAxes` is a Matrix3d with its columns as the rigid axes
      * (with the scale factor removed) and `scale` is the scale factor.
      * * Note that determinant of a rigid matrix is +1.
      * * The context for this method is to determine if the matrix is the product a `rotation` matrix and a uniform
      * `scale` matrix (diagonal matrix with all diagonal entries the same nonzero number).
      */
-    factorRigidWithSignedScale(): {
+    factorRigidWithSignedScale(result?: Matrix3d): {
         rigidAxes: Matrix3d;
         scale: number;
     } | undefined;
@@ -1144,13 +1154,16 @@ export declare class Matrix3d implements BeJSONFunctions {
      * * columns are perpendicular and have unit length.
      * * transpose equals inverse.
      * * mirroring is removed.
+     * * This function internally uses `axisOrderCrossProductsInPlace` to make the matrix rigid.
      * @param axisOrder how to reorder the matrix columns
      * @return whether the adjusted matrix is `rigid` on return
      */
     makeRigid(axisOrder?: AxisOrder): boolean;
-    /** Create a new orthogonal matrix (perpendicular columns, unit length, transpose is inverse).
+    /**
+     * Create a new orthogonal matrix (perpendicular columns, unit length, transpose is inverse).
      * * Columns are taken from the source Matrix3d in order indicated by the axis order.
      * * Mirroring in the matrix is removed.
+     * * This function internally uses `axisOrderCrossProductsInPlace` to make the matrix rigid.
      */
     static createRigidFromMatrix3d(source: Matrix3d, axisOrder?: AxisOrder, result?: Matrix3d): Matrix3d | undefined;
     /**

package/lib/esm/geometry3d/Matrix3d.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"Matrix3d.d.ts","sourceRoot":"","sources":["../../../src/geometry3d/Matrix3d.ts"],"names":[],"mappings":"AAIA;;GAEG;AAEH,OAAO,EAAE,SAAS,EAAE,SAAS,EAAE,eAAe,EAAY,iBAAiB,EAAE,MAAM,aAAa,CAAC;AACjG,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC;AAChC,OAAO,EAAE,OAAO,EAAE,MAAM,mBAAmB,CAAC;AAC5C,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,GAAG,EAAE,MAAM,mBAAmB,CAAC;AAC3D,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AACxC,OAAO,EAAE,aAAa,EAAE,cAAc,EAAE,KAAK,EAAE,MAAM,EAAE,MAAM,YAAY,CAAC;AAI1E;;;;;;;;GAQG;AACH,qBAAa,iBAAiB;IAC5B;;;;;;;;;;;;OAYG;WACW,UAAU,CAAC,IAAI,EAAE,YAAY,EACzC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EACrC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EACrC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM;IAKvC;;;;OAIG;WACW,oBAAoB,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAczG;;;;OAIG;WACW,6BAA6B,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAclH;;;;OAIG;WACW,6BAA6B,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAclH,wCAAwC;WAC1B,gBAAgB,CAAC,CAAC,EAAE,YAAY;IAK9C;;;;OAIG;WACW,cAAc,CAAC,CAAC,EAAE,YAAY,EAAE,IAAI,CAAC,EAAE,YAAY,GAAG,YAAY;IAYhF,0CAA0C;WAC5B,IAAI,CAAC,CAAC,EAAE,YAAY,EAAE,IAAI,EAAE,YAAY,GAAG,YAAY;CAQtE;AAED;;;;;GAKG;AACH,oBAAY,kBAAkB;IAC5B;;;OAGG;IACH,OAAO,IAAA;IACP;;OAEG;IACH,aAAa,IAAA;IACb;;;OAGG;IACH,QAAQ,IAAA;CACT;AAED;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,qBAAa,QAAS,YAAW,eAAe;IAC9C,2EAA2E;IAC3E,OAAc,gBAAgB,UAAQ;IACtC,yEAAyE;IACzE,OAAc,WAAW,SAAK;IAC9B,2DAA2D;IAC3D,OAAc,eAAe,SAAK;IAClC;;;;;;;;OAQG;IACI,KAAK,EAAE,YAAY,CAAC;IAC3B;;;;;;;OAOG;IACI,YAAY,EAAE,YAAY,GAAG,SAAS,CAAC;IAC9C,oEAAoE;IAC7D,YAAY,EAAE,kBAAkB,CAAC;IACxC,0BAA0B;IAC1B,OAAO,CAAC,MAAM,CAAC,SAAS,CAAW;IACnC,+EAA+E;IAC/E,OAAO,CAAC,MAAM,CAAC,cAAc,CAAuB;IACpD,qEAAqE;IACrE,WAAkB,QAAQ,IAAI,QAAQ,CAOrC;IACD,4BAA4B;IACrB,MAAM,IAAI,QAAQ,CAAC,IAAI,CAAC;IAU/B;;;;OAIG;gBACgB,KAAK,CAAC,EAAE,YAAY;IAKvC;;;OAGG;IACI,MAAM,IAAI,aAAa;IAK9B;;;;;;;OAOG;IACI,WAAW,CAAC,IAAI,CAAC,EAAE,aAAa,GAAG,QAAQ,GAAG,IAAI;IA0CzD,4GAA4G;WAC9F,QAAQ,CAAC,IAAI,CAAC,EAAE,aAAa,GAAG,QAAQ;IAKtD;;;OAGG;IACI,aAAa,CAAC,KAAK,EAAE,QAAQ,EAAE,GAAG,CAAC,EAAE,MAAM,GAAG,OAAO;IAG5D;;;OAGG;IACI,mBAAmB,CAAC,WAAW,EAAE,SAAS,EAAE,KAAK,EAAE,QAAQ,EAAE,GAAG,CAAC,EAAE,MAAM,GAAG,OAAO;IAO1F;;;OAGG;IACI,sBAAsB,CAAC,WAAW,EAAE,SAAS,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,CAAC,EAAE,MAAM,GAAG,OAAO;IAOhH;;;;;OAKG;IACI,2BAA2B,CAAC,KAAK,EAAE,QAAQ,EAAE,GAAG,CAAC,EAAE,MAAM,GAAG,OAAO;IAyB1E,oDAAoD;IAC7C,YAAY,CAAC,KAAK,EAAE,QAAQ,GAAG,OAAO;IAG7C,iGAAiG;IACjG,IAAW,IAAI,IAAI,OAAO,CAMzB;IACD;;;OAGG;IACH,OAAO,CAAC,MAAM,CAAC,OAAO;IAGtB;;;;;;;;;;;;;;;;OAgBG;WACW,eAAe,CAC3B,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EACrC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EACrC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EACrC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAQ9B;;;;;;;OAOG;WACW,aAAa,CAAC,KAAK,EAAE,YAAY,EAAE,YAAY,CAAC,EAAE,YAAY,GAAG,QAAQ;IAUvF;;;;;;;;;;OAUG;WACW,wBAAwB,CAAC,SAAS,EAAE,SAAS,EAAE,OAAO,EAAE,QAAQ,GAAG,SAAS,EACxF,OAAO,EAAE,QAAQ,GAAG,SAAS,EAAE,OAAO,EAAE,QAAQ,GAAG,SAAS,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAiB5F;;;;OAIG;IACH,OAAO,CAAC,2BAA2B;IAMnC;;;OAGG;IACH,OAAO,CAAC,qBAAqB;IAS7B;;;;;;;;;;;OAWG;IACI,YAAY,CACjB,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EACrC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EACrC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,IAAI;IAM9C,qCAAqC;IAC9B,WAAW;IAIlB,mCAAmC;IAC5B,OAAO;IAId,yCAAyC;IAClC,OAAO,CAAC,KAAK,EAAE,QAAQ,GAAG,SAAS,GAAG,IAAI;IAoBjD;;;;OAIG;IACI,KAAK,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAKzC;;;;;;;;;OASG;WACW,UAAU,IAAI,QAAQ;IAKpC;;;;;;;;;;OAUG;WACW,cAAc,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAKzD;;;;;;OAMG;WACW,WAAW,CAAC,YAAY,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM,EACxF,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAoB9B;;;;;;OAMG;WACW,kBAAkB,CAAC,WAAW,EAAE,MAAM,GAAG,QAAQ;IAG/D;;;;;OAKG;WACW,qCAAqC,CAAC,OAAO,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAWnG;;;;;OAKG;WACW,8CAA8C,CAAC,OAAO,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAS5G;;;;;;OAMG;WACW,qBAAqB,CAAC,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EACzF,SAAS,EAAE,SAAS,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAOpD;;;;;;;;;OASG;WACW,sBAAsB,CAAC,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,SAAS,EAAE,SAAS,EAC7F,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAe1C;;;;;OAKG;WACW,kBAAkB,CAAC,OAAO,EAAE,QAAQ,EAAE,SAAS,GAAE,SAAyB,EACtF,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAS9B;;;;;;;OAOG;WACW,0BAA0B,CAAC,IAAI,EAAE,QAAQ,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAkB/G;;;;;OAKG;WACW,6BAA6B,CAAC,SAAS,EAAE,SAAS,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IA0B5G;;;;;;;OAOG;IACH,OAAO,CAAC,gBAAgB;IAWxB;;;;;;;;OAQG;IACI,mBAAmB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,IAAI;IAS5E;;;;;;OAMG;WACW,aAAa,CAAC,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IASjH;;;;;OAKG;WACW,gBAAgB,CAAC,OAAO,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,OAAO,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EACjF,OAAO,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IASzD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;OAkCG;WACW,gBAAgB,CAAC,WAAW,EAAE,QAAQ,EAAE,QAAQ,EAAE,QAAQ,EAAE,aAAa,GAAE,MAAU,EACjG,aAAa,GAAE,MAAU,GAAG,QAAQ,GAAG,SAAS;IA2BlD;;;;;;;;;;;;;;;OAeG;WACW,yBAAyB,CAAC,KAAK,EAAE,iBAAiB,EAAE,MAAM,GAAE,OAAe,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAkEvH;;;;;;;;;OASG;IACH,OAAO,CAAC,oBAAoB;IA4B5B;;;;;;;;;;;;;OAaG;IACI,oBAAoB,CAAC,gBAAgB,EAAE,QAAQ,EAAE,MAAM,EAAE,QAAQ,GAAG,OAAO;IAmBlF;;;;;;;;;;OAUG;IACH,OAAO,CAAC,wBAAwB;IAiChC;;;;;;;;;;;;OAYG;IACI,wBAAwB,CAAC,gBAAgB,EAAE,QAAQ,EAAE,MAAM,EAAE,QAAQ,GAAG,OAAO;IAgBtF;;;;OAIG;IACI,yBAAyB,IAAI;QAAE,IAAI,EAAE,QAAQ,CAAC;QAAC,KAAK,EAAE,KAAK,CAAC;QAAC,EAAE,EAAE,OAAO,CAAA;KAAE;IAgEjF;;;;;;;;OAQG;IACH,OAAO,CAAC,yBAAyB;IAmBjC;;;;;;;OAOG;IACI,0BAA0B,CAAC,QAAQ,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,GAAG,OAAO;IAejF;;;;;;;;;;OAUG;IACI,+BAA+B,CAAC,OAAO,EAAE,QAAQ,EAAE,KAAK,EAAE,OAAO,EAAE,OAAO,EAAE,QAAQ,GAAG,OAAO;IA+BrG;;;;;;OAMG;WACW,mCAAmC,CAAC,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,MAAM,EAAE,OAAO,EAAE,QAAQ,EACtG,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAoB1C,6DAA6D;WAC/C,4BAA4B,CAAC,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAGzH,0DAA0D;WAC5C,gCAAgC,CAAC,SAAS,EAAE,MAAM,GAAG,QAAQ;IAyB3E,sCAAsC;IAC/B,OAAO,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAG3C,sCAAsC;IAC/B,OAAO,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAG3C,sCAAsC;IAC/B,OAAO,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAG3C,4CAA4C;IACrC,uBAAuB,IAAI,MAAM;IAGxC,4CAA4C;IACrC,uBAAuB,IAAI,MAAM;IAGxC,4CAA4C;IACrC,uBAAuB,IAAI,MAAM;IAGxC,oCAAoC;IAC7B,gBAAgB,IAAI,MAAM;IAGjC,oCAAoC;IAC7B,gBAAgB,IAAI,MAAM;IAGjC,oCAAoC;IAC7B,gBAAgB,IAAI,MAAM;IAGjC,iDAAiD;IAC1C,6BAA6B,IAAI,MAAM;IAM9C,iCAAiC;IAC1B,aAAa,IAAI,MAAM;IAG9B,kCAAkC;IAC3B,aAAa,IAAI,MAAM;IAG9B,kCAAkC;IAC3B,aAAa,IAAI,MAAM;IAG9B,uDAAuD;IAChD,iBAAiB,IAAI,MAAM;IAKlC,uDAAuD;IAChD,iBAAiB,IAAI,MAAM;IAKlC,uDAAuD;IAChD,iBAAiB,IAAI,MAAM;IAKlC;;;;;;OAMG;IACI,YAAY,CAAC,WAAW,EAAE,SAAS,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM;IAG3E,mCAAmC;IAC5B,IAAI,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAGxC,mCAAmC;IAC5B,IAAI,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAGxC,mCAAmC;IAC5B,IAAI,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAGxC,wEAAwE;IACjE,UAAU,CAAC,MAAM,EAAE,GAAG,GAAG,MAAM;IAGtC,wEAAwE;IACjE,UAAU,CAAC,MAAM,EAAE,GAAG,GAAG,MAAM;IAGtC,wEAAwE;IACjE,UAAU,CAAC,MAAM,EAAE,GAAG,GAAG,MAAM;IAGtC,qEAAqE;IAC9D,OAAO,CAAC,MAAM,EAAE,GAAG,GAAG,MAAM;IAGnC,qEAAqE;IAC9D,OAAO,CAAC,MAAM,EAAE,GAAG,GAAG,MAAM;IAGnC,qEAAqE;IAC9D,OAAO,CAAC,MAAM,EAAE,GAAG,GAAG,MAAM;IAGnC,0DAA0D;IACnD,UAAU,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAG1D,0DAA0D;IACnD,UAAU,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAG1D,0DAA0D;IACnD,UAAU,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAG1D,0EAA0E;IACnE,kBAAkB,CAAC,MAAM,EAAE,GAAG,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAGnE,sEAAsE;IAC/D,oBAAoB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO;IAOhF;;;;OAIG;IACI,SAAS,CAAC,WAAW,EAAE,MAAM,EAAE,KAAK,EAAE,QAAQ,GAAG,SAAS;IAajE;;;;;;OAMG;IACI,UAAU,CAAC,OAAO,EAAE,QAAQ,GAAG,SAAS,EAAE,OAAO,EAAE,QAAQ,GAAG,SAAS,EAAE,OAAO,CAAC,EAAE,QAAQ;IAKlG;;;;OAIG;IACI,MAAM,CAAC,QAAQ,EAAE,MAAM,EAAE,KAAK,EAAE,QAAQ;IAO/C;;;;OAIG;IACI,SAAS,CAAC,WAAW,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IASlE;;;;OAIG;IACI,MAAM,CAAC,WAAW,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAS/D;;;;;;OAMG;WACW,UAAU,CAAC,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAS9G;;;;;;;;;;;;;;OAcG;WACW,sBAAsB,CAAC,SAAS,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAgBrG;;;;;;;OAOG;IACI,cAAc,CAAC,OAAO,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAWnE;;;;OAIG;IACI,0BAA0B,CAAC,IAAI,EAAE,GAAG,EAAE,GAAG,IAAI;IAOpD,yCAAyC;WAC3B,sBAAsB,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAWjH,iFAAiF;WACnE,mBAAmB,CAAC,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAS5G,6EAA6E;WAC/D,qBAAqB,CAAC,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAW7G,yFAAyF;WAC3E,4BAA4B,CAAC,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,cAAc,GAAG,IAAI;IAQvG,qGAAqG;WACvF,6BAA6B,CAAC,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,QAAQ,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EACxG,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAQ5B;;;;;;;;;;;;;;OAcG;WACW,qCAAqC,CAAC,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,QAAQ,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EACrG,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IASlD;;;;;;;;;;;;;;OAcG;WACW,mDAAmD,CAAC,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,QAAQ,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EACnH,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAS5D;;;;;;;;;;;;;OAaG;WACW,2CAA2C,CAAC,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,QAAQ,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAC3G,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAQjD;;;;;;;;;;;OAWG;IACI,uBAAuB,CAAC,MAAM,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAU7E;;;OAGG;IACI,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAOhF;;;;OAIG;IACI,gBAAgB,CAAC,GAAG,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG;IAS7C;;;OAGG;IACI,UAAU,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAOpE;;;OAGG;IACI,uBAAuB,CAAC,MAAM,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAQ5F;;;;OAIG;IACI,qBAAqB,CAAC,OAAO,EAAE,GAAG,GAAG,IAAI;IAQhD;;;;;OAKG;IACI,8BAA8B,CAAC,OAAO,EAAE,GAAG,GAAG,IAAI;IAQzD;;;;;;;;;;;;OAYG;IACI,oBAAoB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAOzF;;;;OAIG;IACI,eAAe,CAAC,MAAM,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAejF;;;;OAIG;IACI,wBAAwB,CAAC,MAAM,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAe1F;;;;;OAKG;IACI,4BAA4B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAY7G;;;;;OAKG;IACI,mBAAmB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAa7G;;;;OAIG;IACI,2BAA2B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAY1G;;;;;;OAMG;IACH,OAAO,CAAC,kBAAkB;IAgB1B;;;OAGG;IACI,oBAAoB,CAAC,KAAK,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAazE;;;OAGG;IACI,2BAA2B,CAAC,KAAK,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAY5F;;;OAGG;IACI,2BAA2B,CAAC,KAAK,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAY5F;;;;;;;OAOG;IACI,6BAA6B,CAAC,KAAK,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAWlF;;;;;;;OAOG;IACI,6BAA6B,CAAC,KAAK,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAWlF;;;;;;;;;;;;OAYG;IACI,uBAAuB,CAAC,KAAK,EAAE,SAAS,EAAE,MAAM,CAAC,EAAE,SAAS,GAAG,SAAS;IAW/E;;;;OAIG;IACI,SAAS,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAa7C;;OAEG;IACI,gBAAgB;IAKvB;;;;;;;;OAQG;IACI,OAAO,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAoBvD;;;OAGG;IACH,OAAO,CAAC,MAAM,CAAC,YAAY;IAK3B;;;;OAIG;IACH,OAAO,CAAC,MAAM,CAAC,sBAAsB;IAMrC;;;;OAIG;IACH,OAAO,CAAC,gCAAgC;IAMxC;;;;;;;;;;;;;;OAcG;IACI,6BAA6B,CAAC,SAAS,EAAE,SAAS,GAAG,IAAI;IAkChE;;;;;;OAMG;IACI,uBAAuB,CAAC,wBAAwB,CAAC,EAAE,QAAQ,GAAG,OAAO;IAW5E;;;;;;OAMG;IACI,oBAAoB,CAAC,qBAAqB,CAAC,EAAE,QAAQ,GAAG,OAAO;IAWtE;;OAEG;IACI,UAAU,IAAI,OAAO;IAG5B;;OAEG;IACI,YAAY,IAAI,IAAI;IAG3B;;;;;OAKG;IACI,oBAAoB,CAAC,mBAAmB,EAAE,OAAO,GAAG,OAAO;IAiClE;;;OAGG;WACW,WAAW,CAAC,GAAG,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,MAAM;IAG9D;;;OAGG;IACI,uBAAuB,CAAC,KAAK,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAIxF,kDAAkD;IAC3C,EAAE,CAAC,GAAG,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,MAAM;IAG9C,+CAA+C;IACxC,KAAK,CAAC,GAAG,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,GAAG,IAAI;IAI9D;;;;;OAKG;IACI,KAAK,CAAC,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAQxD;;;;;;OAMG;IACI,YAAY,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAQhG;;;;;OAKG;IACI,mBAAmB,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM;IAiBzE;;;;;;OAMG;IACI,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAQ7F;;;;;OAKG;IACI,gBAAgB,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM;IAiBtE;;;;OAIG;IACI,gBAAgB,CAAC,KAAK,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,GAAG,IAAI;IAK7D;;;;;;;;;;;;;;;OAeG;IACI,4BAA4B,CAAC,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,GAAG,IAAI;IAc9F;;;;;;;;;;;OAWG;WACW,8BAA8B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAyC1G,+CAA+C;IACxC,WAAW,IAAI,MAAM;IAQ5B;;;;OAIG;IACI,eAAe,IAAI,MAAM;IAQhC,+CAA+C;IACxC,UAAU,IAAI,MAAM;IAM3B,oDAAoD;IAC7C,kBAAkB,IAAI,MAAM;IAMnC,0DAA0D;IACnD,WAAW,IAAI,MAAM;IAG5B,4DAA4D;IACrD,MAAM,IAAI,MAAM;IAMvB,iGAAiG;IAC1F,OAAO,CAAC,KAAK,EAAE,QAAQ,GAAG,MAAM;IAMvC,uDAAuD;IACvD,IAAW,UAAU,IAAI,OAAO,CAE/B;IACD,2DAA2D;IAC3D,IAAW,UAAU,IAAI,OAAO,CAK/B;IACD,iHAAiH;IAC1G,cAAc,IAAI,MAAM;IAO/B,qDAAqD;IAC9C,WAAW,IAAI,OAAO;IAI7B,6DAA6D;IAC7D,IAAW,gBAAgB,IAAI,OAAO,CAErC;IACD,qEAAqE;IACrE,IAAW,iBAAiB,IAAI,OAAO,CAItC;IACD,qEAAqE;IACrE,IAAW,iBAAiB,IAAI,OAAO,CAItC;IACD;;;OAGG;IACI,iBAAiB,IAAI,MAAM,GAAG,SAAS;IAW9C;;;;OAIG;IACI,mCAAmC,IAAI,OAAO;IAIrD;;;;MAIE;IACK,OAAO,CAAC,WAAW,GAAE,OAAe,GAAG,OAAO;IAGrD;;;;;;;;;OASG;IACI,0BAA0B,IAAI;QAAE,SAAS,EAAE,QAAQ,CAAC;QAAC,KAAK,EAAE,MAAM,CAAA;KAAE,GAAG,SAAS;IAUvF,0FAA0F;IAC1F,IAAW,mBAAmB,IAAI,OAAO,CAkBxC;IACD;;;;;;;OAOG;IACI,SAAS,CAAC,SAAS,GAAE,SAAyB,GAAG,OAAO;IAS/D;;;OAGG;WACW,uBAAuB,CAAC,MAAM,EAAE,QAAQ,EAAE,SAAS,GAAE,SAAyB,EAC1F,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAM1C;;;;;;;OAOG;WACW,oBAAoB,CAAC,IAAI,EAAE,OAAO,GAAG,QAAQ;IA2B3D,wEAAwE;IACxE,OAAO,CAAC,MAAM,CAAC,eAAe;IAY9B;;;;;;;;OAQG;IACI,YAAY,IAAI,OAAO;CAkD/B"}
+{"version":3,"file":"Matrix3d.d.ts","sourceRoot":"","sources":["../../../src/geometry3d/Matrix3d.ts"],"names":[],"mappings":"AAIA;;GAEG;AAEH,OAAO,EAAE,SAAS,EAAE,SAAS,EAAE,eAAe,EAAY,iBAAiB,EAAE,MAAM,aAAa,CAAC;AACjG,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAChD,OAAO,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC;AAChC,OAAO,EAAE,OAAO,EAAE,MAAM,mBAAmB,CAAC;AAC5C,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,GAAG,EAAE,MAAM,mBAAmB,CAAC;AAC3D,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AACxC,OAAO,EAAE,aAAa,EAAE,cAAc,EAAE,KAAK,EAAE,MAAM,EAAE,MAAM,YAAY,CAAC;AAI1E;;;;;;;;GAQG;AACH,qBAAa,iBAAiB;IAC5B;;;;;;;;;;;;OAYG;WACW,UAAU,CAAC,IAAI,EAAE,YAAY,EACzC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EACrC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EACrC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM;IAKvC;;;;OAIG;WACW,oBAAoB,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAgBzG;;;;OAIG;WACW,6BAA6B,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAgBlH;;;;OAIG;WACW,6BAA6B,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAgBlH,wCAAwC;WAC1B,gBAAgB,CAAC,CAAC,EAAE,YAAY;IAK9C;;;;OAIG;WACW,cAAc,CAAC,CAAC,EAAE,YAAY,EAAE,IAAI,CAAC,EAAE,YAAY,GAAG,YAAY;IAYhF,0CAA0C;WAC5B,IAAI,CAAC,CAAC,EAAE,YAAY,EAAE,IAAI,EAAE,YAAY,GAAG,YAAY;CAQtE;AAED;;;;;GAKG;AACH,oBAAY,kBAAkB;IAC5B;;;OAGG;IACH,OAAO,IAAA;IACP;;OAEG;IACH,aAAa,IAAA;IACb;;;OAGG;IACH,QAAQ,IAAA;CACT;AAED;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,qBAAa,QAAS,YAAW,eAAe;IAC9C,2EAA2E;IAC3E,OAAc,gBAAgB,UAAQ;IACtC,yEAAyE;IACzE,OAAc,WAAW,SAAK;IAC9B,2DAA2D;IAC3D,OAAc,eAAe,SAAK;IAClC;;;;;;;;OAQG;IACI,KAAK,EAAE,YAAY,CAAC;IAC3B;;;;;;;OAOG;IACI,YAAY,EAAE,YAAY,GAAG,SAAS,CAAC;IAC9C,oEAAoE;IAC7D,YAAY,EAAE,kBAAkB,CAAC;IACxC,0BAA0B;IAC1B,OAAO,CAAC,MAAM,CAAC,SAAS,CAAW;IACnC,+EAA+E;IAC/E,OAAO,CAAC,MAAM,CAAC,cAAc,CAAuB;IACpD,qEAAqE;IACrE,WAAkB,QAAQ,IAAI,QAAQ,CAOrC;IACD,4BAA4B;IACrB,MAAM,IAAI,QAAQ,CAAC,IAAI,CAAC;IAU/B;;;;OAIG;gBACgB,KAAK,CAAC,EAAE,YAAY;IAKvC;;;OAGG;IACI,MAAM,IAAI,aAAa;IAK9B;;;;;;;OAOG;IACI,WAAW,CAAC,IAAI,CAAC,EAAE,aAAa,GAAG,QAAQ,GAAG,IAAI;IA0CzD,4GAA4G;WAC9F,QAAQ,CAAC,IAAI,CAAC,EAAE,aAAa,GAAG,QAAQ;IAKtD;;;OAGG;IACI,aAAa,CAAC,KAAK,EAAE,QAAQ,EAAE,GAAG,CAAC,EAAE,MAAM,GAAG,OAAO;IAG5D;;;OAGG;IACI,mBAAmB,CAAC,WAAW,EAAE,SAAS,EAAE,KAAK,EAAE,QAAQ,EAAE,GAAG,CAAC,EAAE,MAAM,GAAG,OAAO;IAO1F;;;OAGG;IACI,sBAAsB,CAAC,WAAW,EAAE,SAAS,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,CAAC,EAAE,MAAM,GAAG,OAAO;IAOhH;;;;;OAKG;IACI,2BAA2B,CAAC,KAAK,EAAE,QAAQ,EAAE,GAAG,CAAC,EAAE,MAAM,GAAG,OAAO;IAyB1E,oDAAoD;IAC7C,YAAY,CAAC,KAAK,EAAE,QAAQ,GAAG,OAAO;IAG7C,iGAAiG;IACjG,IAAW,IAAI,IAAI,OAAO,CAMzB;IACD;;;OAGG;IACH,OAAO,CAAC,MAAM,CAAC,OAAO;IAGtB;;;;;;;;;;;;;;;;OAgBG;WACW,eAAe,CAC3B,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EACrC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EACrC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EACrC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAQ9B;;;;;;;OAOG;WACW,aAAa,CAAC,KAAK,EAAE,YAAY,EAAE,YAAY,CAAC,EAAE,YAAY,GAAG,QAAQ;IAUvF;;;;;;;;;;OAUG;WACW,wBAAwB,CAAC,SAAS,EAAE,SAAS,EAAE,OAAO,EAAE,QAAQ,GAAG,SAAS,EACxF,OAAO,EAAE,QAAQ,GAAG,SAAS,EAAE,OAAO,EAAE,QAAQ,GAAG,SAAS,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAiB5F;;;;OAIG;IACH,OAAO,CAAC,2BAA2B;IAMnC;;;OAGG;IACH,OAAO,CAAC,qBAAqB;IAS7B;;;;;;;;;;;OAWG;IACI,YAAY,CACjB,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EACrC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EACrC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,IAAI;IAM9C,qCAAqC;IAC9B,WAAW;IAIlB,mCAAmC;IAC5B,OAAO;IAId,2FAA2F;IACpF,OAAO,CAAC,KAAK,EAAE,QAAQ,GAAG,SAAS,GAAG,IAAI;IAoBjD;;;;OAIG;IACI,KAAK,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAKzC;;;;;;;;;OASG;WACW,UAAU,IAAI,QAAQ;IAKpC;;;;;;;;;;OAUG;WACW,cAAc,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAKzD;;;;;;OAMG;WACW,WAAW,CAAC,YAAY,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM,EAAE,YAAY,EAAE,MAAM,EACxF,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAoB9B;;;;;;OAMG;WACW,kBAAkB,CAAC,WAAW,EAAE,MAAM,GAAG,QAAQ;IAG/D;;;;;OAKG;WACW,qCAAqC,CAAC,OAAO,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAWnG;;;;;OAKG;WACW,8CAA8C,CAAC,OAAO,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAS5G;;;;;;OAMG;WACW,qBAAqB,CAAC,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EACzF,SAAS,EAAE,SAAS,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAOpD;;;;;;;;;OASG;WACW,sBAAsB,CAAC,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,SAAS,EAAE,SAAS,EAC7F,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAe1C;;;;;;;;OAQG;WACW,kBAAkB,CAAC,OAAO,EAAE,QAAQ,EAAE,SAAS,GAAE,SAAyB,EACtF,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAS9B;;;;;;;OAOG;WACW,0BAA0B,CAAC,IAAI,EAAE,QAAQ,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAkB/G;;;;;OAKG;WACW,6BAA6B,CAAC,SAAS,EAAE,SAAS,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IA0B5G;;;;;;;;OAQG;IACH,OAAO,CAAC,gBAAgB;IAWxB;;;;;;;;;OASG;IACI,mBAAmB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,IAAI;IAS5E;;;;;;OAMG;WACW,aAAa,CAAC,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IASjH;;;;;OAKG;WACW,gBAAgB,CAAC,OAAO,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,OAAO,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EACjF,OAAO,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IASzD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;OAkCG;WACW,gBAAgB,CAAC,WAAW,EAAE,QAAQ,EAAE,QAAQ,EAAE,QAAQ,EAAE,aAAa,GAAE,MAAU,EACjG,aAAa,GAAE,MAAU,GAAG,QAAQ,GAAG,SAAS;IA2BlD;;;;;;;;;;;;;;;OAeG;WACW,yBAAyB,CAAC,KAAK,EAAE,iBAAiB,EAAE,MAAM,GAAE,OAAe,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAkEvH;;;;;;;;;OASG;IACH,OAAO,CAAC,oBAAoB;IA4B5B;;;;;;;;;;;;;OAaG;IACI,oBAAoB,CAAC,gBAAgB,EAAE,QAAQ,EAAE,MAAM,EAAE,QAAQ,GAAG,OAAO;IAmBlF;;;;;;;;;;OAUG;IACH,OAAO,CAAC,wBAAwB;IAiChC;;;;;;;;;;;;OAYG;IACI,wBAAwB,CAAC,gBAAgB,EAAE,QAAQ,EAAE,MAAM,EAAE,QAAQ,GAAG,OAAO;IAgBtF;;;;OAIG;IACI,yBAAyB,IAAI;QAAE,IAAI,EAAE,QAAQ,CAAC;QAAC,KAAK,EAAE,KAAK,CAAC;QAAC,EAAE,EAAE,OAAO,CAAA;KAAE;IAgEjF;;;;;;;;OAQG;IACH,OAAO,CAAC,yBAAyB;IAmBjC;;;;;;;OAOG;IACI,0BAA0B,CAAC,QAAQ,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,GAAG,OAAO;IAejF;;;;;;;;;;OAUG;IACI,+BAA+B,CAAC,OAAO,EAAE,QAAQ,EAAE,KAAK,EAAE,OAAO,EAAE,OAAO,EAAE,QAAQ,GAAG,OAAO;IA+BrG;;;;;;OAMG;WACW,mCAAmC,CAAC,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,MAAM,EAAE,OAAO,EAAE,QAAQ,EACtG,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAoB1C,6DAA6D;WAC/C,4BAA4B,CAAC,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAGzH,0DAA0D;WAC5C,gCAAgC,CAAC,SAAS,EAAE,MAAM,GAAG,QAAQ;IAyB3E,sCAAsC;IAC/B,OAAO,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAG3C,sCAAsC;IAC/B,OAAO,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAG3C,sCAAsC;IAC/B,OAAO,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAG3C,4CAA4C;IACrC,uBAAuB,IAAI,MAAM;IAGxC,4CAA4C;IACrC,uBAAuB,IAAI,MAAM;IAGxC,4CAA4C;IACrC,uBAAuB,IAAI,MAAM;IAGxC,oCAAoC;IAC7B,gBAAgB,IAAI,MAAM;IAGjC,oCAAoC;IAC7B,gBAAgB,IAAI,MAAM;IAGjC,oCAAoC;IAC7B,gBAAgB,IAAI,MAAM;IAGjC,iDAAiD;IAC1C,6BAA6B,IAAI,MAAM;IAM9C,iCAAiC;IAC1B,aAAa,IAAI,MAAM;IAG9B,kCAAkC;IAC3B,aAAa,IAAI,MAAM;IAG9B,kCAAkC;IAC3B,aAAa,IAAI,MAAM;IAG9B,uDAAuD;IAChD,iBAAiB,IAAI,MAAM;IAKlC,uDAAuD;IAChD,iBAAiB,IAAI,MAAM;IAKlC,uDAAuD;IAChD,iBAAiB,IAAI,MAAM;IAKlC;;;;;;OAMG;IACI,YAAY,CAAC,WAAW,EAAE,SAAS,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM;IAG3E,mCAAmC;IAC5B,IAAI,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAGxC,mCAAmC;IAC5B,IAAI,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAGxC,mCAAmC;IAC5B,IAAI,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAGxC,wEAAwE;IACjE,UAAU,CAAC,MAAM,EAAE,GAAG,GAAG,MAAM;IAGtC,wEAAwE;IACjE,UAAU,CAAC,MAAM,EAAE,GAAG,GAAG,MAAM;IAGtC,wEAAwE;IACjE,UAAU,CAAC,MAAM,EAAE,GAAG,GAAG,MAAM;IAGtC,qEAAqE;IAC9D,OAAO,CAAC,MAAM,EAAE,GAAG,GAAG,MAAM;IAGnC,qEAAqE;IAC9D,OAAO,CAAC,MAAM,EAAE,GAAG,GAAG,MAAM;IAGnC,qEAAqE;IAC9D,OAAO,CAAC,MAAM,EAAE,GAAG,GAAG,MAAM;IAGnC,0DAA0D;IACnD,UAAU,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAG1D,0DAA0D;IACnD,UAAU,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAG1D,0DAA0D;IACnD,UAAU,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAG1D,0EAA0E;IACnE,kBAAkB,CAAC,MAAM,EAAE,GAAG,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAGnE,sEAAsE;IAC/D,oBAAoB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO;IAOhF;;;;OAIG;IACI,SAAS,CAAC,WAAW,EAAE,MAAM,EAAE,KAAK,EAAE,QAAQ,GAAG,SAAS;IAajE;;;;;;OAMG;IACI,UAAU,CAAC,OAAO,EAAE,QAAQ,GAAG,SAAS,EAAE,OAAO,EAAE,QAAQ,GAAG,SAAS,EAAE,OAAO,CAAC,EAAE,QAAQ;IAKlG;;;;OAIG;IACI,MAAM,CAAC,QAAQ,EAAE,MAAM,EAAE,KAAK,EAAE,QAAQ;IAO/C;;;;OAIG;IACI,SAAS,CAAC,WAAW,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IASlE;;;;OAIG;IACI,MAAM,CAAC,WAAW,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAS/D;;;;;;OAMG;WACW,UAAU,CAAC,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAS9G;;;;;;;;;;;;;;OAcG;WACW,sBAAsB,CAAC,SAAS,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAgBrG;;;;;;;OAOG;IACI,cAAc,CAAC,OAAO,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAWnE;;;;OAIG;IACI,0BAA0B,CAAC,IAAI,EAAE,GAAG,EAAE,GAAG,IAAI;IAOpD,yCAAyC;WAC3B,sBAAsB,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAWjH,iFAAiF;WACnE,mBAAmB,CAAC,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAS5G,6EAA6E;WAC/D,qBAAqB,CAAC,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAW7G,yFAAyF;WAC3E,4BAA4B,CAAC,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,cAAc,GAAG,IAAI;IAQvG,qGAAqG;WACvF,6BAA6B,CAAC,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,QAAQ,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EACxG,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAQ5B;;;;;;;;;;;;;;OAcG;WACW,qCAAqC,CAAC,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,QAAQ,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EACrG,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IASlD;;;;;;;;;;;;;;OAcG;WACW,mDAAmD,CAAC,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,QAAQ,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EACnH,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAS5D;;;;;;;;;;;;;OAaG;WACW,2CAA2C,CAAC,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,QAAQ,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAC3G,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAQjD;;;;;;;;;;;OAWG;IACI,uBAAuB,CAAC,MAAM,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAU7E;;;OAGG;IACI,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAOhF;;;;OAIG;IACI,gBAAgB,CAAC,GAAG,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,GAAG,GAAG;IASnD;;;OAGG;IACI,UAAU,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAOpE;;;OAGG;IACI,uBAAuB,CAAC,MAAM,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAQ5F;;;;OAIG;IACI,qBAAqB,CAAC,OAAO,EAAE,GAAG,GAAG,IAAI;IAQhD;;;;;OAKG;IACI,8BAA8B,CAAC,OAAO,EAAE,GAAG,GAAG,IAAI;IAQzD;;;;;;;;;;;;OAYG;IACI,oBAAoB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAOzF;;;;OAIG;IACI,eAAe,CAAC,MAAM,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAejF;;;;OAIG;IACI,wBAAwB,CAAC,MAAM,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAe1F;;;;;OAKG;IACI,4BAA4B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAY7G;;;;;OAKG;IACI,mBAAmB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAa7G;;;;OAIG;IACI,2BAA2B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO,GAAG,SAAS;IAY1G;;;;;;OAMG;IACH,OAAO,CAAC,kBAAkB;IAgB1B;;;OAGG;IACI,oBAAoB,CAAC,KAAK,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAazE;;;OAGG;IACI,2BAA2B,CAAC,KAAK,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAY5F;;;OAGG;IACI,2BAA2B,CAAC,KAAK,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAY5F;;;;;;;OAOG;IACI,6BAA6B,CAAC,KAAK,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAWlF;;;;;;;OAOG;IACI,6BAA6B,CAAC,KAAK,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAWlF;;;;;;;;;;;;OAYG;IACI,uBAAuB,CAAC,KAAK,EAAE,SAAS,EAAE,MAAM,CAAC,EAAE,SAAS,GAAG,SAAS;IAW/E;;;;OAIG;IACI,SAAS,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAa7C;;OAEG;IACI,gBAAgB;IAKvB;;;;;;;;OAQG;IACI,OAAO,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAoBvD;;;OAGG;IACH,OAAO,CAAC,MAAM,CAAC,YAAY;IAK3B;;;;OAIG;IACH,OAAO,CAAC,MAAM,CAAC,sBAAsB;IAMrC;;;;OAIG;IACH,OAAO,CAAC,gCAAgC;IAMxC;;;;;;;;;;;;;;OAcG;IACI,6BAA6B,CAAC,SAAS,EAAE,SAAS,GAAG,IAAI;IAkChE;;;;;;OAMG;IACI,uBAAuB,CAAC,wBAAwB,CAAC,EAAE,QAAQ,GAAG,OAAO;IAW5E;;;;;;OAMG;IACI,oBAAoB,CAAC,qBAAqB,CAAC,EAAE,QAAQ,GAAG,OAAO;IAWtE;;OAEG;IACI,UAAU,IAAI,OAAO;IAG5B;;OAEG;IACI,YAAY,IAAI,IAAI;IAG3B;;;;;OAKG;IACI,oBAAoB,CAAC,mBAAmB,EAAE,OAAO,GAAG,OAAO;IAiClE;;;OAGG;WACW,WAAW,CAAC,GAAG,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,MAAM;IAG9D;;;OAGG;IACI,uBAAuB,CAAC,KAAK,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,OAAO,GAAG,OAAO;IAIxF,kDAAkD;IAC3C,EAAE,CAAC,GAAG,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,MAAM;IAG9C,+CAA+C;IACxC,KAAK,CAAC,GAAG,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,GAAG,IAAI;IAI9D;;;;;OAKG;IACI,KAAK,CAAC,KAAK,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAQxD;;;;;;OAMG;IACI,YAAY,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAQhG;;;;;OAKG;IACI,mBAAmB,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM;IAiBzE;;;;;;OAMG;IACI,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAQ7F;;;;;OAKG;IACI,gBAAgB,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM;IAiBtE;;;;OAIG;IACI,gBAAgB,CAAC,KAAK,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,GAAG,IAAI;IAK7D;;;;;;;;;;;;;;;OAeG;IACI,4BAA4B,CAAC,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,GAAG,IAAI;IAc9F;;;;;;;;;;;;;;;OAeG;WACW,8BAA8B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ;IAyC1G,+CAA+C;IACxC,WAAW,IAAI,MAAM;IAQ5B;;;;OAIG;IACI,eAAe,IAAI,MAAM;IAQhC,+CAA+C;IACxC,UAAU,IAAI,MAAM;IAM3B,oDAAoD;IAC7C,kBAAkB,IAAI,MAAM;IAMnC,0DAA0D;IACnD,WAAW,IAAI,MAAM;IAG5B,4DAA4D;IACrD,MAAM,IAAI,MAAM;IAMvB,iGAAiG;IAC1F,OAAO,CAAC,KAAK,EAAE,QAAQ,GAAG,MAAM;IAMvC,uDAAuD;IACvD,IAAW,UAAU,IAAI,OAAO,CAE/B;IACD,2DAA2D;IAC3D,IAAW,UAAU,IAAI,OAAO,CAK/B;IACD,iHAAiH;IAC1G,cAAc,IAAI,MAAM;IAO/B,qDAAqD;IAC9C,WAAW,IAAI,OAAO;IAI7B,6DAA6D;IAC7D,IAAW,gBAAgB,IAAI,OAAO,CAErC;IACD,qEAAqE;IACrE,IAAW,iBAAiB,IAAI,OAAO,CAItC;IACD,qEAAqE;IACrE,IAAW,iBAAiB,IAAI,OAAO,CAItC;IACD;;;OAGG;IACI,iBAAiB,IAAI,MAAM,GAAG,SAAS;IAW9C;;;;OAIG;IACI,mCAAmC,IAAI,OAAO;IAIrD;;;;MAIE;IACK,OAAO,CAAC,WAAW,GAAE,OAAe,GAAG,OAAO;IAGrD;;;;;;;;;;OAUG;IACI,0BAA0B,CAAC,MAAM,CAAC,EAAE,QAAQ,GAAG;QAAE,SAAS,EAAE,QAAQ,CAAC;QAAC,KAAK,EAAE,MAAM,CAAA;KAAE,GAAG,SAAS;IASxG,0FAA0F;IAC1F,IAAW,mBAAmB,IAAI,OAAO,CAkBxC;IACD;;;;;;;;OAQG;IACI,SAAS,CAAC,SAAS,GAAE,SAAyB,GAAG,OAAO;IAS/D;;;;;OAKG;WACW,uBAAuB,CAAC,MAAM,EAAE,QAAQ,EAAE,SAAS,GAAE,SAAyB,EAC1F,MAAM,CAAC,EAAE,QAAQ,GAAG,QAAQ,GAAG,SAAS;IAM1C;;;;;;;OAOG;WACW,oBAAoB,CAAC,IAAI,EAAE,OAAO,GAAG,QAAQ;IA2B3D,wEAAwE;IACxE,OAAO,CAAC,MAAM,CAAC,eAAe;IAY9B;;;;;;;;OAQG;IACI,YAAY,IAAI,OAAO;CAkD/B"}

package/lib/esm/geometry3d/Matrix3d.js

@@ -479,7 +479,7 @@ export class Matrix3d {
         this.setRowValues(0, 0, 0, 0, 0, 0, 0, 0, 0);
         this.inverseState = InverseMatrixState.singular;
     }
-    /** Copy contents from another matrix. */
+    /** Copy contents from the `other` matrix. If `other` is undefined, use identity matrix. */
     setFrom(other) {
         if (other === undefined) {
             this.setIdentity();
@@ -654,6 +654,9 @@ export class Matrix3d {
      * Construct a rigid matrix (orthogonal matrix with +1 determinant) using vectorA and its 2 perpendicular.
      * * If axisOrder is not passed then `AxisOrder = AxisOrder.ZXY` is used as default.
      * * This function internally uses createPerpendicularVectorFavorXYPlane and createRigidFromColumns.
+     * * If you want to rotate a given plane (which contains (0,0,0)) to the xy-plane, pass the normal vector of
+     * your plane into createRigidHeadsUp. The transpose of the returned Matrix3d can be used to rotate your plane
+     * to the xy-plane. If plane does not contain (0,0,0) then the plane is rotated to a plane parallel to the xy-plane.
      * * Visualization can be found at https://www.itwinjs.org/sandbox/SaeedTorabi/2PerpendicularVectorsTo1Vector
      */
     static createRigidHeadsUp(vectorA, axisOrder = AxisOrder.ZXY, result) {
@@ -711,6 +714,7 @@ export class Matrix3d {
     /**
      * Replace current rows Ui and Uj with (c*Ui + s*Uj) and (c*Uj - s*Ui).
      * * There is no checking for i,j being 0,1,2.
+     * * The instance matrix A is multiplied in place on the left by a Givens rotation G, resulting in the matrix G*A.
      * @param i first row index. **must be 0,1,2** (unchecked)
      * @param j second row index. **must be 0,1,2** (unchecked)
      * @param c fist coefficient
@@ -730,6 +734,7 @@ export class Matrix3d {
     /**
      * Replace current columns Ui and Uj with (c*Ui + s*Uj) and (c*Uj - s*Ui).
      * * There is no checking for i,j being 0,1,2.
+     * * The instance matrix A is multiplied in place on the right by a Givens rotation G, resulting in the matrix A*G.
      * * This is used in compute intensive inner loops
      * @param i first row index. **must be 0,1,2** (unchecked)
      * @param j second row index. **must be 0,1,2** (unchecked)
@@ -1755,7 +1760,7 @@ export class Matrix3d {
         return undefined;
     }
     /**
-     * Multiply `matrixInverse * [x,y,z]` and return result as `Point4d` with given weight.
+     * Multiply `matrixInverse * [x,y,z]` and return result as `Point4d` the with given weight as last element.
      * * Equivalent to solving `matrix * result = [x,y,z]` for an unknown `result`.
      * * Result is `undefined` if the matrix is singular (e.g. has parallel columns or a zero magnitude column)
      * @return result as a Point4d with the same weight.
@@ -2301,6 +2306,10 @@ export class Matrix3d {
      * * column 1 is perpendicular to both. It is the "up" vector on the view plane.
      * * Multiplying the returned matrix times a local (view) vector gives the world vector.
      * * Multiplying transpose of the returned matrix times a world vector gives the local (view) vector.
+     * * If you want to rotate a given plane (which contains (0,0,0)) to the xy-plane, pass coordinates of the normal
+     * vector of your plane into createRigidViewAxesZTowardsEye. The transpose of the returned Matrix3d can be used
+     * to rotate your plane to the xy-plane. If plane does not contain (0,0,0) then the plane is rotated to a plane
+     * parallel to the xy-plane.
      * @param x eye x coordinate
      * @param y eye y coordinate
      * @param z eye z coordinate
@@ -2469,21 +2478,21 @@ export class Matrix3d {
      * Test if all rows and columns are perpendicular to each other and have equal length.
      * If so, the length (or its negative) is the `scale` factor from a set of `orthonormal axes` to
      * the set of axes created by columns of `this` matrix. Otherwise, returns `undefined`.
+     * @param result optional pre-allocated object to populate and return
      * @returns returns `{ rigidAxes, scale }` where `rigidAxes` is a Matrix3d with its columns as the rigid axes
      * (with the scale factor removed) and `scale` is the scale factor.
      * * Note that determinant of a rigid matrix is +1.
      * * The context for this method is to determine if the matrix is the product a `rotation` matrix and a uniform
      * `scale` matrix (diagonal matrix with all diagonal entries the same nonzero number).
      */
-    factorRigidWithSignedScale() {
+    factorRigidWithSignedScale(result) {
         const product = this.multiplyMatrixMatrixTranspose(this);
         const scaleSquare = product.sameDiagonalScale();
         if (scaleSquare === undefined || scaleSquare <= 0.0)
             return undefined;
         const scale = this.determinant() > 0 ? Math.sqrt(scaleSquare) : -Math.sqrt(scaleSquare);
         const scaleInverse = 1.0 / scale;
-        const result = { rigidAxes: this.scaleColumns(scaleInverse, scaleInverse, scaleInverse), scale };
-        return result;
+        return { rigidAxes: this.scaleColumns(scaleInverse, scaleInverse, scaleInverse, result), scale };
     }
     /** Test if `this` matrix reorders and/or negates the columns of the `identity` matrix. */
     get isSignedPermutation() {
@@ -2512,6 +2521,7 @@ export class Matrix3d {
      * * columns are perpendicular and have unit length.
      * * transpose equals inverse.
      * * mirroring is removed.
+     * * This function internally uses `axisOrderCrossProductsInPlace` to make the matrix rigid.
      * @param axisOrder how to reorder the matrix columns
      * @return whether the adjusted matrix is `rigid` on return
      */
@@ -2524,9 +2534,11 @@ export class Matrix3d {
         this.axisOrderCrossProductsInPlace(axisOrder);
         return this.normalizeColumnsInPlace();
     }
-    /** Create a new orthogonal matrix (perpendicular columns, unit length, transpose is inverse).
+    /**
+     * Create a new orthogonal matrix (perpendicular columns, unit length, transpose is inverse).
      * * Columns are taken from the source Matrix3d in order indicated by the axis order.
      * * Mirroring in the matrix is removed.
+     * * This function internally uses `axisOrderCrossProductsInPlace` to make the matrix rigid.
      */
     static createRigidFromMatrix3d(source, axisOrder = AxisOrder.XYZ, result) {
         result = source.clone(result);