@itwin/core-geometry 4.3.0-dev.8 → 4.3.0

package/lib/esm/geometry3d/Angle.js.map

@@ -1 +1 @@
-{"version":3,"file":"Angle.js","sourceRoot":"","sources":["../../../src/geometry3d/Angle.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAEH,OAAO,EAA+B,QAAQ,EAAc,MAAM,aAAa,CAAC;AAEhF;;;;;;;;;GASG;AACH,MAAa,KAAK;IAiBhB,YAAoB,OAAO,GAAG,CAAC,EAAE,OAAgB;QAC/C,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;QACxB,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;IAC1B,CAAC;IACD,gDAAgD;IACzC,KAAK;QACV,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjD,CAAC;IACD,8CAA8C;IACvC,MAAM;QACX,OAAO,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;IAC7B,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,aAAa,CAAC,OAAe;QACzC,OAAO,IAAI,KAAK,CAAC,KAAK,CAAC,gBAAgB,CAAC,OAAO,CAAC,EAAE,OAAO,CAAC,CAAC;IAC7D,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,aAAa,CAAC,OAAe;QACzC,OAAO,IAAI,KAAK,CAAC,OAAO,CAAC,CAAC;IAC5B,CAAC;IACD;;OAEG;IACI,MAAM,CAAC,gBAAgB;QAC5B,OAAO,IAAI,KAAK,CAAC,QAAQ,CAAC,iBAAiB,CAAC,CAAC;IAC/C,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,iBAAiB,CAAC,MAAa,EAAE,QAAgB,EAAE,MAAa;QAC5E,OAAO,IAAI,KAAK,CAAC,QAAQ,CAAC,WAAW,CAAC,MAAM,CAAC,OAAO,EAAE,QAAQ,EAAE,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC;IACnF,CAAC;IACD;;;OAGG;IACI,WAAW,CAAC,KAAa;QAC9B,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,OAAO,GAAG,KAAK,CAAC,CAAC;IACzC,CAAC;IACD;;;OAGG;IACI,UAAU,CAAC,OAAe;QAC/B,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;QACxB,IAAI,CAAC,QAAQ,GAAG,SAAS,CAAC;IAC5B,CAAC;IACD;;;OAGG;IACI,UAAU,CAAC,OAAe;QAC/B,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC;QAChD,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;IAC1B,CAAC;IACD,yCAAyC;IAClC,MAAM,CAAC,SAAS;QACrB,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,EAAE,GAAG,GAAG,EAAE,KAAK,CAAC,CAAC;IACzC,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,WAAW,CAAC,SAAiB,EAAE,WAAmB;QAC9D,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,SAAS,EAAE,WAAW,CAAC,CAAC,CAAC;IACvD,CAAC;IACD;;;OAGG;IACI,OAAO,CAAC,KAAY;QACzB,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,QAAQ,CAAC;QAC/B,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,QAAQ,CAAC;IACjC,CAAC;IACD;;;;;;;OAOG;IACI,WAAW,CAAC,IAAiB,EAAE,iBAA0B;QAC9D,IAAI,CAAC,QAAQ,GAAG,iBAAiB,CAAC,CAAC,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1D,IAAI,CAAC,IAAI;YACP,OAAO;QACT,IAAI,OAAO,IAAI,KAAK,QAAQ,EAAE;YAC5B,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC;SACvB;aAAM,IAAI,OAAQ,IAAY,CAAC,OAAO,KAAK,QAAQ,EAAE;YACpD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,OAAO,CAAC,CAAC;SACxC;aAAM,IAAI,OAAQ,IAAY,CAAC,QAAQ,KAAK,QAAQ,EAAE;YACrD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,QAAQ,CAAC,CAAC;SACzC;aAAM,IAAI,OAAQ,IAAY,CAAC,OAAO,KAAK,QAAQ,EAAE;YACpD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,OAAO,CAAC,CAAC;SACxC;aAAM,IAAI,OAAQ,IAAY,CAAC,QAAQ,KAAK,QAAQ,EAAE;YACrD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,QAAQ,CAAC,CAAC;SACzC;IACH,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,QAAQ,CAAC,IAAiB,EAAE,iBAA0B;QAClE,MAAM,GAAG,GAAG,IAAI,KAAK,EAAE,CAAC;QACxB,GAAG,CAAC,WAAW,CAAC,IAAI,EAAE,iBAAiB,CAAC,CAAC;QACzC,OAAO,GAAG,CAAC;IACb,CAAC;IACD,+DAA+D;IACxD,MAAM;QACX,OAAO,IAAI,CAAC,OAAO,CAAC;IACtB,CAAC;IACD,yEAAyE;IAClE,aAAa;QAClB,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,CAAC;IACnC,CAAC;IACD,4CAA4C;IAC5C,IAAW,OAAO;QAChB,OAAO,IAAI,CAAC,QAAQ,CAAC;IACvB,CAAC;IACD,4CAA4C;IAC5C,IAAW,OAAO;QAChB,OAAO,IAAI,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,KAAK,CAAC,gBAAgB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAC7F,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,gBAAgB,CAAC,OAAe;QAC5C,OAAO,OAAO,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;IACjC,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,gBAAgB,CAAC,OAAe;QAC5C,IAAI,OAAO,GAAG,CAAC;YACb,OAAO,CAAC,KAAK,CAAC,gBAAgB,CAAC,CAAC,OAAO,CAAC,CAAC;QAC3C,8BAA8B;QAC9B,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QACnB,MAAM,MAAM,GAAG,KAAK,GAAG,EAAE,CAAC;QAC1B;;;;YAII;QACJ,IAAI,OAAO,IAAI,IAAI,GAAG,EAAE;YACtB,OAAO,MAAM,GAAG,OAAO,CAAC;QAC1B,IAAI,OAAO,GAAG,IAAI,GAAG,EAAE;YACrB,OAAO,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QAClD,IAAI,OAAO,IAAI,IAAI,GAAG,EAAE;YACtB,OAAO,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QAC7C,IAAI,OAAO,IAAI,IAAI,GAAG,EAAE;YACtB,OAAO,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QACnD,sDAAsD;QACtD,OAAO,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;IACnD,CAAC;IACD,qDAAqD;IAC9C,GAAG;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjC,CAAC;IACD,mDAAmD;IAC5C,GAAG;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjC,CAAC;IACD,sDAAsD;IAC/C,GAAG;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjC,CAAC;IACD,kEAAkE;IAC3D,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,OAAO,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,QAAQ,CAAC,gCAAgC,CAAC;IACxE,CAAC;IACD,iDAAiD;IAC1C,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,IAAI,QAAQ,CAAC,iBAAiB,CAAC;IAC7E,CAAC;IACD,yCAAyC;IACzC,IAAW,YAAY;QACrB,OAAO,KAAK,CAAC,mBAAmB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAClD,CAAC;IACD,+DAA+D;IAC/D,IAAW,YAAY;QACrB,OAAO,KAAK,CAAC,mBAAmB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAClD,CAAC;IACD,yDAAyD;IAClD,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,KAAK,CAAC;YACrB,IAAI,OAAO,GAAG,MAAM;gBAClB,OAAO,OAAO,CAAC;YACjB,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,MAAM,CAAC,CAAC;YAChD,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,MAAM,OAAO,GAAG,KAAK,CAAC,mBAAmB,CAAC,CAAC,OAAO,CAAC,CAAC;YACpD,OAAO,KAAK,GAAG,OAAO,CAAC;SACxB;QACD,sCAAsC;QACtC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,mDAAmD;IAC5C,MAAM,CAAC,sBAAsB,CAAC,OAAe;QAClD,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,KAAK;YAC5B,OAAO,OAAO,CAAC;QACjB,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,KAAK,CAAC;YACrB,MAAM,UAAU,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,OAAO,GAAG,KAAK,CAAC,GAAG,MAAM,CAAC,CAAC;YAC9D,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,OAAO,CAAC,KAAK,CAAC,sBAAsB,CAAC,CAAC,OAAO,CAAC,CAAC;SAChD;QACD,iCAAiC;QACjC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,yDAAyD;IAClD,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;YAC7B,IAAI,OAAO,GAAG,MAAM;gBAClB,OAAO,OAAO,CAAC;YACjB,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,MAAM,CAAC,CAAC;YAChD,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,OAAO,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,KAAK,CAAC,mBAAmB,CAAC,CAAC,OAAO,CAAC,CAAC;SAC5D;QACD,iCAAiC;QACjC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,0DAA0D;IACnD,MAAM,CAAC,0BAA0B,CAAC,OAAe;QACtD,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,IAAI,CAAC,EAAE;YAC9B,OAAO,OAAO,CAAC;QACjB,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;YAC7B,MAAM,UAAU,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,OAAO,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC;YAChE,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,OAAO,CAAC,KAAK,CAAC,0BAA0B,CAAC,CAAC,OAAO,CAAC,CAAC;SACpD;QACD,iCAAiC;QACjC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,mEAAmE;IAC5D,MAAM,CAAC,IAAI;QAChB,OAAO,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC;IACtB,CAAC;IACD,yCAAyC;IACzC,IAAW,WAAW;QACpB,OAAO,IAAI,CAAC,OAAO,KAAK,CAAC,CAAC;IAC5B,CAAC;IACD,uFAAuF;IACvF,IAAW,YAAY;QACrB,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,GAAG,QAAQ,CAAC,iBAAiB,CAAC;IAC7D,CAAC;IACD,wGAAwG;IACxG,IAAW,wBAAwB;QACjC,OAAO,KAAK,CAAC,mBAAmB,CAAC,IAAI,CAAC,OAAO,GAAG,GAAG,CAAC,CAAC;IACvD,CAAC;IACD,gEAAgE;IACzD,MAAM,CAAC,2BAA2B,CAAC,OAAe;QACvD,OAAO,KAAK,CAAC,aAAa,CAAC,KAAK,CAAC,mBAAmB,CAAC,OAAO,CAAC,CAAC,CAAC;IACjE,CAAC;IACD,mEAAmE;IAC5D,MAAM,CAAC,4BAA4B,CAAC,OAAe;QACxD,OAAO,KAAK,CAAC,aAAa,CAAC,KAAK,CAAC,sBAAsB,CAAC,OAAO,CAAC,CAAC,CAAC;IACpE,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,oCAAoC,CAAC,QAAgB,EAAE,QAAgB,EACnF,YAAoB,QAAQ,CAAC,iBAAiB;QAC9C,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,CAAC;QAC5C,IAAI,KAAK,IAAI,SAAS;YACpB,OAAO,IAAI,CAAC;QACd,MAAM,MAAM,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;QAC7B,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,IAAI,SAAS;YACvC,OAAO,IAAI,CAAC;QACd,MAAM,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC;QAC7C,MAAM,MAAM,GAAG,KAAK,GAAG,SAAS,GAAG,MAAM,CAAC;QAC1C,OAAO,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,SAAS,CAAC;IACvC,CAAC;IACD;;;KAGC;IACM,0BAA0B,CAAC,KAAY;QAC5C,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;IAC3D,CAAC;IAED;;;;OAIG;IACI,6BAA6B,CAAC,KAAY,EAAE,YAAoB,QAAQ,CAAC,iBAAiB;QAC/F,OAAO,KAAK,CAAC,oCAAoC,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;IAC9F,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,iCAAiC,CAAC,QAAgB,EAAE,QAAgB,EAChF,YAAoB,QAAQ,CAAC,iBAAiB;QAC9C,OAAO,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,SAAS,CAAC;IACnD,CAAC;IACD;;;;OAIG;IACI,0BAA0B,CAAC,KAAY,EAAE,YAAoB,QAAQ,CAAC,iBAAiB;QAC5F,OAAO,KAAK,CAAC,iCAAiC,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;IAC3F,CAAC;IACD;;;;;OAKG;IACI,aAAa,CAAC,KAAY,EAAE,YAAoB,QAAQ,CAAC,iBAAiB;QAC/E,OAAO,IAAI,CAAC,0BAA0B,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;IAC3D,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,qBAAqB,CAAC,KAAa,EAAE,KAAa,EAAE,KAAa;QAC7E,OAAO,KAAK,GAAG,QAAQ,CAAC,0BAA0B;eAC7C,KAAK,GAAG,QAAQ,CAAC,0BAA0B;eAC3C,KAAK,GAAG,KAAK,IAAI,QAAQ,CAAC,wBAAwB,GAAG,KAAK,GAAG,KAAK,CAAC;IAC1E,CAAC;IACD;;;;;;OAMG;IACI,MAAM,CAAC,+BAA+B,CAAC,MAAc,EAAE,MAAc;QAC1E,MAAM,CAAC,GAAG,QAAQ,CAAC,YAAY,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;QAChD,IAAI,CAAC,GAAG,QAAQ,CAAC,mBAAmB,EAAE;YACpC,OAAO,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,CAAC,YAAY;SACtD;aAAM;YACL;;;;eAIG;YACH,IAAI,IAAI,CAAC;YACT,IAAI,IAAI,GAAG,GAAG,CAAC;YACf,MAAM,KAAK,GAAG,MAAM,GAAG,CAAC,CAAC;YACzB,MAAM,KAAK,GAAG,MAAM,GAAG,CAAC,CAAC;YACzB,qEAAqE;YACrE,IAAI,KAAK,IAAI,GAAG,EAAE;gBAChB;;;;mBAIG;gBACH,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;gBACtC,uDAAuD;gBACvD,IAAI,GAAG,KAAK,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC;aAC7B;iBAAM;gBACL,2DAA2D;gBAC3D,IAAI,KAAK,GAAG,GAAG,EAAE;oBACf;;;;uBAIG;oBACH,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;oBACtC,2DAA2D;iBAC5D;qBAAM;oBACL,kDAAkD;oBAClD,IAAI,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;iBACxC;gBACD,uDAAuD;gBACvD,IAAI,GAAG,KAAK,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,kBAAkB;aAChD;YACD,OAAO,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,EAAE,OAAO,EAAE,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC;SAC9D;IACH,CAAC;IACD,8EAA8E;IACvE,MAAM,CAAC,gBAAgB,CAAC,KAAa,EAAE,YAAoB,QAAQ,CAAC,kBAAkB;QAC3F,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QACjC,IAAI,QAAQ,IAAI,SAAS;YACvB,OAAO,CAAC,CAAC;QACX,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,CAAC;QACjC,IAAI,CAAC,IAAI,SAAS;YAChB,OAAO,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;QAClC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,CAAC;QAC7B,IAAI,CAAC,IAAI,SAAS;YAChB,OAAO,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;QAClC,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;;;;;;;;;;;OAaG;IACI,MAAM,CAAC,gCAAgC,CAC5C,KAAa,EAAE,KAAa,EAAE,KAAa,EAAE,YAAqB,IAAI;QAGtE,MAAM,MAAM,GAAG,KAAK,GAAG,KAAK,CAAC;QAC7B,MAAM,MAAM,GAAG,GAAG,GAAG,KAAK,CAAC;QAC3B,IAAI,SAAS,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,QAAQ,CAAC,iBAAiB,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;YAClG,OAAO,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC;QAC1C,OAAO,KAAK,CAAC,+BAA+B,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC/D,CAAC;IACD;;;;;;;;;OASG;IACI,MAAM,CAAC,wBAAwB,CACpC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU;QAEtE,MAAM,KAAK,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC1C,OAAO,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,qBAAqB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACnF,CAAC;IACD;;;;;;;;;;;;;;;;;;OAkBG;IACI,MAAM,CAAC,gCAAgC,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EACnH,SAAiB,EAAE,SAAiB,EAAE,SAAiB,EACvD,mBAA4B,KAAK;QACjC,MAAM,KAAK,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC1C,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC7B,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC7B,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC7B,MAAM,MAAM,GAAG,SAAS,GAAG,EAAE,GAAG,SAAS,GAAG,EAAE,GAAG,SAAS,GAAG,EAAE,CAAC;QAChE,MAAM,cAAc,GAAG,QAAQ,CAAC,aAAa,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAC1D,IAAI,MAAM,GAAG,GAAG,EAAE;YAChB,IAAI,gBAAgB,EAAE;gBACpB,wGAAwG;gBACxG,4EAA4E;gBAC5E,OAAO,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,cAAc,EAAE,CAAC,KAAK,CAAC,CAAC;aACrD;iBAAM;gBACL,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,cAAc,EAAE,KAAK,CAAC,CAAC;aAC3C;SACF;aAAM;YACL,OAAO,IAAI,CAAC,KAAK,CAAC,cAAc,EAAE,KAAK,CAAC,CAAC;SAC1C;IACH,CAAC;IACD;;;OAGG;IACI,uBAAuB,CAAC,QAAgB;QAC7C,IAAI,IAAI,CAAC,QAAQ,KAAK,SAAS,EAAE;YAC/B,IAAI,CAAC,QAAQ,IAAI,QAAQ,GAAG,KAAK,CAAC;YAClC,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,gBAAgB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;SACvD;aAAM;YACL,IAAI,CAAC,QAAQ,IAAI,QAAQ,GAAG,KAAK,CAAC,UAAU,CAAC;SAC9C;IACH,CAAC;;AA7gBD,+DAA+D;AACxC,qBAAe,GAAG,mBAAmB,CAAC;AAC7D,8DAA8D;AACvC,oBAAc,GAAG,sBAAsB,CAAC;AAC/D,8DAA8D;AACvC,oBAAc,GAAG,uBAAuB,CAAC;AAChE,6DAA6D;AACtC,eAAS,GAAG,sBAAsB,CAAC;AAC1D,+DAA+D;AACxC,gBAAU,GAAG,sBAAsB,CAAC;AAC3D,qDAAqD;AAC9B,sBAAgB,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,cAAc,CAAC,CAAC;AACxE,qDAAqD;AAC9B,sBAAgB,GAAG,CAAC,KAAK,CAAC,cAAc,GAAG,IAAI,CAAC,CAAC;SAd7D,KAAK","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/** @packageDocumentation\r\n * @module CartesianGeometry\r\n */\r\n\r\nimport { AngleProps, BeJSONFunctions, Geometry, TrigValues } from \"../Geometry\";\r\n\r\n/**\r\n * An `Angle` carries the numeric value of an angle, with methods to allow (require!) callers to\r\n * be clear about whether their angle is degrees or radians.\r\n * * After the Angle object is created, the callers should not know or care whether it is stored in\r\n * `degrees` or `radians` because both are available if requested by caller.\r\n * * The various access method are named so that callers can specify whether untyped numbers passed in or\r\n * out are degrees or radians.\r\n * * Visualization can be found at https://www.itwinjs.org/sandbox/SaeedTorabi/AngleSweep\r\n * @public\r\n */\r\nexport class Angle implements BeJSONFunctions {\r\n  /** maximal accuracy value of pi/12 (15 degrees), in radians */\r\n  public static readonly piOver12Radians = 0.26179938779914946;\r\n  /** maximal accuracy value of pi/4 (45 degrees), in radians */\r\n  public static readonly piOver4Radians = 7.853981633974483e-001;\r\n  /** maximal accuracy value of pi/2 (90 degrees), in radians */\r\n  public static readonly piOver2Radians = 1.5707963267948966e+000;\r\n  /** maximal accuracy value of pi (180 degrees), in radians */\r\n  public static readonly piRadians = 3.141592653589793e+000;\r\n  /** maximal accuracy value of 2*pi (360 degrees), in radians */\r\n  public static readonly pi2Radians = 6.283185307179586e+000;\r\n  /** scale factor for converting radians to degrees */\r\n  public static readonly degreesPerRadian = (45.0 / Angle.piOver4Radians);\r\n  /** scale factor for converting degrees to radians */\r\n  public static readonly radiansPerDegree = (Angle.piOver4Radians / 45.0);\r\n  private _radians: number;\r\n  private _degrees?: number;\r\n  private constructor(radians = 0, degrees?: number) {\r\n    this._radians = radians;\r\n    this._degrees = degrees;\r\n  }\r\n  /** Return a new angle with the same content. */\r\n  public clone(): Angle {\r\n    return new Angle(this._radians, this._degrees);\r\n  }\r\n  /** Freeze this instance so it is read-only */\r\n  public freeze(): Readonly<this> {\r\n    return Object.freeze(this);\r\n  }\r\n  /**\r\n   * Return a new Angle object for angle given in degrees.\r\n   * @param degrees angle in degrees\r\n   */\r\n  public static createDegrees(degrees: number): Angle {\r\n    return new Angle(Angle.degreesToRadians(degrees), degrees);\r\n  }\r\n  /**\r\n   * Return a (new) Angle object for a value given in radians.\r\n   * @param radians angle in radians\r\n   */\r\n  public static createRadians(radians: number): Angle {\r\n    return new Angle(radians);\r\n  }\r\n  /**\r\n   * Return a new `Angle` object with the default \"small\" angle measurement specified by [[Geometry.smallAngleRadians]].\r\n   */\r\n  public static createSmallAngle(): Angle {\r\n    return new Angle(Geometry.smallAngleRadians);\r\n  }\r\n  /**\r\n   * Return a (new) Angle object that is interpolated between two inputs (based on a fraction)\r\n   * @param angle0 first angle in radians\r\n   * @param fraction the interpolation fraction\r\n   * @param angle1 second angle in radians\r\n   */\r\n  public static createInterpolate(angle0: Angle, fraction: number, angle1: Angle): Angle {\r\n    return new Angle(Geometry.interpolate(angle0.radians, fraction, angle1.radians));\r\n  }\r\n  /**\r\n   * Return a (new) Angle object, with angle scaled from existing angle.\r\n   * @param scale scale factor to apply to angle.\r\n   */\r\n  public cloneScaled(scale: number): Angle {\r\n    return new Angle(this.radians * scale);\r\n  }\r\n  /**\r\n   * Set this angle to a value given in radians.\r\n   * @param radians angle given in radians\r\n   */\r\n  public setRadians(radians: number) {\r\n    this._radians = radians;\r\n    this._degrees = undefined;\r\n  }\r\n  /**\r\n   * Set this angle to a value given in degrees.\r\n   * @param degrees angle given in degrees.\r\n   */\r\n  public setDegrees(degrees: number) {\r\n    this._radians = Angle.degreesToRadians(degrees);\r\n    this._degrees = degrees;\r\n  }\r\n  /** Create an angle for a full circle. */\r\n  public static create360(): Angle {\r\n    return new Angle(Math.PI * 2.0, 360.0);\r\n  }\r\n  /**\r\n   * Create a (strongly typed) Angle whose tangent is `numerator/denominator`, using the signs of both in\r\n   * determining the (otherwise ambiguous) quadrant.\r\n   * @param numerator numerator for tangent\r\n   * @param denominator denominator for tangent\r\n   */\r\n  public static createAtan2(numerator: number, denominator: number): Angle {\r\n    return new Angle(Math.atan2(numerator, denominator));\r\n  }\r\n  /**\r\n   * Copy all contents of `other` to this Angle.\r\n   * @param other source data\r\n   */\r\n  public setFrom(other: Angle) {\r\n    this._radians = other._radians;\r\n    this._degrees = other._degrees;\r\n  }\r\n  /**\r\n   * Set an Angle from a JSON object\r\n   * * A simple number is considered as degrees.\r\n   * * specified `json.degrees` or `json._degrees` is degree value.\r\n   * * specified `son.radians` or `json._radians` is radians value.\r\n   * @param json object from JSON.parse. If a number, value is in *DEGREES*\r\n   * @param defaultValRadians if json is undefined, default value in radians.\r\n   */\r\n  public setFromJSON(json?: AngleProps, defaultValRadians?: number) {\r\n    this._radians = defaultValRadians ? defaultValRadians : 0;\r\n    if (!json)\r\n      return;\r\n    if (typeof json === \"number\") {\r\n      this.setDegrees(json);\r\n    } else if (typeof (json as any).degrees === \"number\") {\r\n      this.setDegrees((json as any).degrees);\r\n    } else if (typeof (json as any)._degrees === \"number\") {\r\n      this.setDegrees((json as any)._degrees);\r\n    } else if (typeof (json as any).radians === \"number\") {\r\n      this.setRadians((json as any).radians);\r\n    } else if (typeof (json as any)._radians === \"number\") {\r\n      this.setRadians((json as any)._radians);\r\n    }\r\n  }\r\n  /**\r\n   * Create an Angle from a JSON object\r\n   * @param json object from JSON.parse. If a number, value is in *DEGREES*\r\n   * @param defaultValRadians if json is undefined, default value in radians.\r\n   * @return a new Angle\r\n   */\r\n  public static fromJSON(json?: AngleProps, defaultValRadians?: number): Angle {\r\n    const val = new Angle();\r\n    val.setFromJSON(json, defaultValRadians);\r\n    return val;\r\n  }\r\n  /** Convert an Angle to a JSON object as a number in degrees */\r\n  public toJSON(): AngleProps {\r\n    return this.degrees;\r\n  }\r\n  /** Return a json object with radians keyword, e.g. `{ radians: 0.10}` */\r\n  public toJSONRadians(): AngleProps {\r\n    return { radians: this.radians };\r\n  }\r\n  /** Return the angle measured in radians. */\r\n  public get radians(): number {\r\n    return this._radians;\r\n  }\r\n  /** Return the angle measured in degrees. */\r\n  public get degrees(): number {\r\n    return this._degrees !== undefined ? this._degrees : Angle.radiansToDegrees(this._radians);\r\n  }\r\n  /**\r\n   * Convert an angle in degrees to radians.\r\n   * @param degrees angle in degrees\r\n   */\r\n  public static degreesToRadians(degrees: number): number {\r\n    return degrees * Math.PI / 180;\r\n  }\r\n  /**\r\n   * Convert an angle in radians to degrees.\r\n   * @param degrees angle in radians\r\n   */\r\n  public static radiansToDegrees(radians: number): number {\r\n    if (radians < 0)\r\n      return -Angle.radiansToDegrees(-radians);\r\n    // Now radians is positive ...\r\n    const pi = Math.PI;\r\n    const factor = 180.0 / pi;\r\n    /* the following if statements are for round-off reasons. The problem is that no IEEE number is\r\n      * an exact hit for any primary multiple of pi (90, 180, etc). The following is supposed to have\r\n      * a better chance that if the input was computed by direct assignment from 90, 180, etc degrees\r\n      * it will return exactly 90,180 etc.\r\n      */\r\n    if (radians <= 0.25 * pi)\r\n      return factor * radians;\r\n    if (radians < 0.75 * pi)\r\n      return 90.0 + 180 * ((radians - 0.5 * pi) / pi);\r\n    if (radians <= 1.25 * pi)\r\n      return 180.0 + 180 * ((radians - pi) / pi);\r\n    if (radians <= 1.75 * pi)\r\n      return 270.0 + 180 * ((radians - 1.5 * pi) / pi);\r\n    // all larger radians reference from 360 degrees (2PI)\r\n    return 360.0 + 180 * ((radians - 2.0 * pi) / pi);\r\n  }\r\n  /** Return the cosine of this Angle object's angle */\r\n  public cos(): number {\r\n    return Math.cos(this._radians);\r\n  }\r\n  /** Return the sine of this Angle object's angle */\r\n  public sin(): number {\r\n    return Math.sin(this._radians);\r\n  }\r\n  /** Return the tangent of this Angle object's angle */\r\n  public tan(): number {\r\n    return Math.tan(this._radians);\r\n  }\r\n  /** Test if a radians (absolute) value is nearly 2PI or larger! */\r\n  public static isFullCircleRadians(radians: number): boolean {\r\n    return Math.abs(radians) >= Geometry.fullCircleRadiansMinusSmallAngle;\r\n  }\r\n  /** Test if the radians value is a half circle */\r\n  public static isHalfCircleRadians(radians: number): boolean {\r\n    return Math.abs(Math.abs(radians) - Math.PI) <= Geometry.smallAngleRadians;\r\n  }\r\n  /** Test if the angle is a full circle */\r\n  public get isFullCircle(): boolean {\r\n    return Angle.isFullCircleRadians(this._radians);\r\n  }\r\n  /** Test if the angle is a half circle (in either direction) */\r\n  public get isHalfCircle(): boolean {\r\n    return Angle.isHalfCircleRadians(this._radians);\r\n  }\r\n  /** Adjust a radians value so it is positive in 0..360 */\r\n  public static adjustDegrees0To360(degrees: number): number {\r\n    if (degrees >= 0) {\r\n      const period = 360.0;\r\n      if (degrees < period)\r\n        return degrees;\r\n      const numPeriods = Math.floor(degrees / period);\r\n      return degrees - numPeriods * period;\r\n    } else if (degrees < 0) {\r\n      // negative angle ...\r\n      const radians = Angle.adjustDegrees0To360(-degrees);\r\n      return 360.0 - radians;\r\n    }\r\n    // fall through for Nan (disaster) !!!\r\n    return 0;\r\n  }\r\n  /** Adjust a radians value so it is in -180..180 */\r\n  public static adjustDegreesSigned180(degrees: number): number {\r\n    if (Math.abs(degrees) <= 180.0)\r\n      return degrees;\r\n    if (degrees >= 0) {\r\n      const period = 360.0;\r\n      const numPeriods = 1 + Math.floor((degrees - 180.0) / period);\r\n      return degrees - numPeriods * period;\r\n    } else if (degrees < 0) {\r\n      // negative angle ...\r\n      return -Angle.adjustDegreesSigned180(-degrees);\r\n    }\r\n    // fall through for NaN disaster.\r\n    return 0;\r\n  }\r\n  /** Adjust a radians value so it is positive in 0..2Pi */\r\n  public static adjustRadians0To2Pi(radians: number): number {\r\n    if (radians >= 0) {\r\n      const period = Math.PI * 2.0;\r\n      if (radians < period)\r\n        return radians;\r\n      const numPeriods = Math.floor(radians / period);\r\n      return radians - numPeriods * period;\r\n    } else if (radians < 0) {\r\n      // negative angle ...\r\n      return Math.PI * 2.0 - Angle.adjustRadians0To2Pi(-radians);\r\n    }\r\n    // fall through for NaN disaster.\r\n    return 0;\r\n  }\r\n  /** Adjust a radians value so it is positive in -PI..PI */\r\n  public static adjustRadiansMinusPiPlusPi(radians: number): number {\r\n    if (Math.abs(radians) <= Math.PI)\r\n      return radians;\r\n    if (radians >= 0) {\r\n      const period = Math.PI * 2.0;\r\n      const numPeriods = 1 + Math.floor((radians - Math.PI) / period);\r\n      return radians - numPeriods * period;\r\n    } else if (radians < 0) {\r\n      // negative angle ...\r\n      return -Angle.adjustRadiansMinusPiPlusPi(-radians);\r\n    }\r\n    // fall through for NaN disaster.\r\n    return 0;\r\n  }\r\n  /** Return a (newly allocated) Angle object with value 0 radians */\r\n  public static zero(): Angle {\r\n    return new Angle(0);\r\n  }\r\n  /** Test if the angle is exactly zero. */\r\n  public get isExactZero(): boolean {\r\n    return this.radians === 0;\r\n  }\r\n  /** Test if the angle is almost zero (within tolerance `Geometry.smallAngleRadians`) */\r\n  public get isAlmostZero(): boolean {\r\n    return Math.abs(this.radians) < Geometry.smallAngleRadians;\r\n  }\r\n  /** Test if the angle is almost a north or south pole (within tolerance `Geometry.smallAngleRadians`) */\r\n  public get isAlmostNorthOrSouthPole(): boolean {\r\n    return Angle.isHalfCircleRadians(this.radians * 2.0);\r\n  }\r\n  /** Create an angle object with degrees adjusted into 0..360. */\r\n  public static createDegreesAdjustPositive(degrees: number): Angle {\r\n    return Angle.createDegrees(Angle.adjustDegrees0To360(degrees));\r\n  }\r\n  /** Create an angle object with degrees adjusted into -180..180. */\r\n  public static createDegreesAdjustSigned180(degrees: number): Angle {\r\n    return Angle.createDegrees(Angle.adjustDegreesSigned180(degrees));\r\n  }\r\n  /**\r\n   * Test if two radians values are equivalent, allowing shift by full circle (i.e. by a multiple of `2*PI`)\r\n   * @param radiansA first radians value\r\n   * @param radiansB second radians value\r\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\r\n   */\r\n  public static isAlmostEqualRadiansAllowPeriodShift(radiansA: number, radiansB: number,\r\n    radianTol: number = Geometry.smallAngleRadians): boolean {\r\n    const delta = Math.abs(radiansA - radiansB);\r\n    if (delta <= radianTol)\r\n      return true;\r\n    const period = Math.PI * 2.0;\r\n    if (Math.abs(delta - period) <= radianTol)\r\n      return true;\r\n    const numPeriod = Math.round(delta / period);\r\n    const delta1 = delta - numPeriod * period;\r\n    return Math.abs(delta1) <= radianTol;\r\n  }\r\n  /**\r\n * Test if this angle has magnitude no greater than that of `other`.\r\n * @param other the other angle\r\n */\r\n  public isMagnitudeLessThanOrEqual(other: Angle): boolean {\r\n    return Math.abs(this.radians) <= Math.abs(other.radians);\r\n  }\r\n\r\n  /**\r\n   * Test if this angle and other are equivalent, allowing shift by full circle (i.e., multiples of `2 * PI`).\r\n   * @param other the other angle\r\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\r\n   */\r\n  public isAlmostEqualAllowPeriodShift(other: Angle, radianTol: number = Geometry.smallAngleRadians): boolean {\r\n    return Angle.isAlmostEqualRadiansAllowPeriodShift(this._radians, other._radians, radianTol);\r\n  }\r\n  /**\r\n   * Test if two angle (in radians)  almost equal, NOT allowing shift by full circle (i.e., multiples of `2 * PI`).\r\n   * @param radiansA first radians value\r\n   * @param radiansB second radians value\r\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\r\n   */\r\n  public static isAlmostEqualRadiansNoPeriodShift(radiansA: number, radiansB: number,\r\n    radianTol: number = Geometry.smallAngleRadians): boolean {\r\n    return Math.abs(radiansA - radiansB) < radianTol;\r\n  }\r\n  /**\r\n   * Test if two this angle and other are almost equal, NOT allowing shift by full circle (i.e., multiples of `2 * PI`).\r\n   * @param other the other angle\r\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\r\n   */\r\n  public isAlmostEqualNoPeriodShift(other: Angle, radianTol: number = Geometry.smallAngleRadians): boolean {\r\n    return Angle.isAlmostEqualRadiansNoPeriodShift(this._radians, other._radians, radianTol);\r\n  }\r\n  /**\r\n   * Test if two this angle and other are almost equal, NOT allowing shift by full circle (i.e., multiples of `2 * PI`).\r\n   * * This function is same as isAlmostEqualRadiansNoPeriodShift. Please use isAlmostEqualRadiansNoPeriodShift.\r\n   * @param other the other angle\r\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\r\n   */\r\n  public isAlmostEqual(other: Angle, radianTol: number = Geometry.smallAngleRadians): boolean {\r\n    return this.isAlmostEqualNoPeriodShift(other, radianTol);\r\n  }\r\n  /**\r\n   * Test if dot product values indicate non-zero length perpendicular vectors.\r\n   * @param dotUU dot product of vectorU with itself\r\n   * @param dotVV dot product of vectorV with itself\r\n   * @param dotUV dot product of vectorU with vectorV\r\n   */\r\n  public static isPerpendicularDotSet(dotUU: number, dotVV: number, dotUV: number): boolean {\r\n    return dotUU > Geometry.smallMetricDistanceSquared\r\n      && dotVV > Geometry.smallMetricDistanceSquared\r\n      && dotUV * dotUV <= Geometry.smallAngleRadiansSquared * dotUU * dotVV;\r\n  }\r\n  /**\r\n   * Return cosine, sine, and radians for the half angle of a \"cosine,sine\" pair.\r\n   * * This function assumes the input arguments are related to an angle between -PI and PI\r\n   * * This function returns an angle between -PI and PI\r\n   * @param rCos2A cosine value (scaled by radius) for initial angle.\r\n   * @param rSin2A sine value (scaled by radius) for final angle.\r\n   */\r\n  public static trigValuesToHalfAngleTrigValues(rCos2A: number, rSin2A: number): TrigValues {\r\n    const r = Geometry.hypotenuseXY(rCos2A, rSin2A);\r\n    if (r < Geometry.smallMetricDistance) {\r\n      return { c: 1.0, s: 0.0, radians: 0.0 }; // angle = 0\r\n    } else {\r\n      /* If the caller really gave you sine and cosine values, r should be 1.  However,\r\n       * to allow scaled values -- e.g. the x and y components of any vector -- we normalize\r\n       * right here. This adds an extra sqrt and 2 divides to the whole process, but improves\r\n       * both the usefulness and robustness of the computation.\r\n       */\r\n      let cosA;\r\n      let sinA = 0.0;\r\n      const cos2A = rCos2A / r;\r\n      const sin2A = rSin2A / r;\r\n      // Original angle in NE and SE quadrants. Half angle in same quadrant\r\n      if (cos2A >= 0.0) {\r\n        /*\r\n         * We know cos2A = (cosA)^2 - (sinA)^2 and 1 = (cosA)^2 + (sinA)^2\r\n         * so 1 + cos2A = 2(cosA)^2 and therefore, cosA = sqrt((1+cos2A)/2)\r\n         * cosine is positive in NE and SE quadrants so we use +sqrt\r\n         */\r\n        cosA = Math.sqrt(0.5 * (1.0 + cos2A));\r\n        // We know sin2A = 2 sinA cosA so sinA = sin2A/(2*cosA)\r\n        sinA = sin2A / (2.0 * cosA);\r\n      } else {\r\n        // Original angle in NW quadrant. Half angle in NE quadrant\r\n        if (sin2A > 0.0) {\r\n          /*\r\n           * We know cos2A = (cosA)^2 - (sinA)^2 and 1 = (cosA)^2 + (sinA)^2\r\n           * so 1 - cos2A = 2(sinA)^2 and therefore, sinA = sqrt((1-cos2A)/2)\r\n           * sine is positive in NE quadrant so we use +sqrt\r\n           */\r\n          sinA = Math.sqrt(0.5 * (1.0 - cos2A));\r\n          // Original angle in SW quadrant. Half angle in SE quadrant\r\n        } else {\r\n          // sine is negative in SE quadrant so we use -sqrt\r\n          sinA = -Math.sqrt(0.5 * (1.0 - cos2A));\r\n        }\r\n        // We know sin2A = 2 sinA cosA so cosA = sin2A/(2*sinA)\r\n        cosA = sin2A / (2.0 * sinA); // always positive\r\n      }\r\n      return { c: cosA, s: sinA, radians: Math.atan2(sinA, cosA) };\r\n    }\r\n  }\r\n  /** If value is close to -1, -0.5, 0, 0.5, 1, adjust it to the exact value. */\r\n  public static cleanupTrigValue(value: number, tolerance: number = Geometry.smallFloatingPoint): number {\r\n    const absValue = Math.abs(value);\r\n    if (absValue <= tolerance)\r\n      return 0;\r\n    let a = Math.abs(absValue - 0.5);\r\n    if (a <= tolerance)\r\n      return value < 0.0 ? -0.5 : 0.5;\r\n    a = Math.abs(absValue - 1.0);\r\n    if (a <= tolerance)\r\n      return value < 0.0 ? -1.0 : 1.0;\r\n    return value;\r\n  }\r\n  /**\r\n   * Return the half angle cosine, sine, and radians for given dot products between vectors. The vectors define\r\n   * an ellipse using x(t) = c + U cos(t) + V sin(t) so U and V are at angle t=0 degree and t=90 degree. The\r\n   * half angle t0 is an angle such that x(t0) is one of the ellipse semi-axis.\r\n   * * This construction arises e.g. in `Arc3d.toScaledMatrix3d`.\r\n   * * Given ellipse x(t) = c + U cos(t) + V sin(t), find t0 such that radial vector W(t0) = x(t0) - c is\r\n   * perpendicular to the ellipse.\r\n   * * Then 0 = W(t0).x'(t0) = (U cos(t0) + V sin(t0)).(V cos(t0) - U sin(t0)) = U.V cos(2t0) + 0.5 (V.V - U.U) sin(2t0)\r\n   * implies sin(2t0) / cos(2t0) = 2 U.V / (U.U - V.V), i.e., t0 can be computed given the three dot products on the RHS.\r\n   * math details can be found at docs/learning/geometry/Angle.md\r\n   * @param dotUU dot product of vectorU with itself\r\n   * @param dotVV dot product of vectorV with itself\r\n   * @param dotUV dot product of vectorU with vectorV\r\n   */\r\n  public static dotProductsToHalfAngleTrigValues(\r\n    dotUU: number, dotVV: number, dotUV: number, favorZero: boolean = true,\r\n  ): TrigValues {\r\n\r\n    const cos2t0 = dotUU - dotVV;\r\n    const sin2t0 = 2.0 * dotUV;\r\n    if (favorZero && Math.abs(sin2t0) < Geometry.smallAngleRadians * (Math.abs(dotUU) + Math.abs(dotVV)))\r\n      return { c: 1.0, s: 0.0, radians: 0.0 };\r\n    return Angle.trigValuesToHalfAngleTrigValues(cos2t0, sin2t0);\r\n  }\r\n  /**\r\n   * Returns the angle between two vectors, with the vectors given as xyz components\r\n   * * The returned angle is between 0 and PI\r\n   * @param ux x component of vector u\r\n   * @param uy y component of vector u\r\n   * @param uz z component of vector u\r\n   * @param vx x component of vector v\r\n   * @param vy y component of vector v\r\n   * @param vz z component of vector v\r\n   */\r\n  public static radiansBetweenVectorsXYZ(\r\n    ux: number, uy: number, uz: number, vx: number, vy: number, vz: number,\r\n  ): number {\r\n    const uDotV = ux * vx + uy * vy + uz * vz;\r\n    return Math.atan2(Geometry.crossProductMagnitude(ux, uy, uz, vx, vy, vz), uDotV);\r\n  }\r\n  /**\r\n   * Returns the angle between two vectors, with the vectors given as xyz components, and an up vector to resolve\r\n   * angle to a full 2PI range.\r\n   * * The returned angle is (-PI < radians <= PI) or (0 <= radians < 2 * PI)\r\n   * * The angle is in the plane of the U and V vectors.\r\n   * * The upVector determines a positive side of the plane but need not be strictly perpendicular to the plane.\r\n   *\r\n   * @param ux x component of vector u\r\n   * @param uy y component of vector u\r\n   * @param uz z component of vector u\r\n   * @param vx x component of vector v\r\n   * @param vy y component of vector v\r\n   * @param vz z component of vector v\r\n   * @param upVectorX x component of vector to positive side of plane.\r\n   * @param upVectorY y component of vector to positive side of plane.\r\n   * @param upVectorZ z component of vector to positive side of plane.\r\n   * @param adjustToAllPositive if true, return strictly non-negative sweep (0 <= radians < 2*PI). If false, return\r\n   * signed (-PI < radians <= PI)\r\n   */\r\n  public static orientedRadiansBetweenVectorsXYZ(ux: number, uy: number, uz: number, vx: number, vy: number, vz: number,\r\n    upVectorX: number, upVectorY: number, upVectorZ: number,\r\n    adjustToPositive: boolean = false): number {\r\n    const uDotV = ux * vx + uy * vy + uz * vz;\r\n    const wx = uy * vz - uz * vy;\r\n    const wy = uz * vx - ux * vz;\r\n    const wz = ux * vy - uy * vx;\r\n    const upDotW = upVectorX * wx + upVectorY * wy + upVectorZ * wz;\r\n    const crossMagnitude = Geometry.hypotenuseXYZ(wx, wy, wz);\r\n    if (upDotW < 0.0) {\r\n      if (adjustToPositive) {\r\n        // The turn is greater than 180 degrees.  Take a peculiarly oriented atan2 to get the excess-180 part as\r\n        // addition to PI. This gives the smoothest numerical transition passing PI.\r\n        return Math.PI + Math.atan2(crossMagnitude, -uDotV);\r\n      } else {\r\n        return -Math.atan2(crossMagnitude, uDotV);\r\n      }\r\n    } else {\r\n      return Math.atan2(crossMagnitude, uDotV);\r\n    }\r\n  }\r\n  /**\r\n   * Add a multiple of a full circle angle (360 degrees, 2PI) in place.\r\n   * @param multiple multiplier factor\r\n   */\r\n  public addMultipleOf2PiInPlace(multiple: number) {\r\n    if (this._degrees !== undefined) {\r\n      this._degrees += multiple * 360.0;\r\n      this._radians = Angle.degreesToRadians(this._degrees);\r\n    } else {\r\n      this._radians += multiple * Angle.pi2Radians;\r\n    }\r\n  }\r\n}\r\n"]}
+{"version":3,"file":"Angle.js","sourceRoot":"","sources":["../../../src/geometry3d/Angle.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAEH,OAAO,EAA+B,QAAQ,EAAc,MAAM,aAAa,CAAC;AAEhF;;;;;;;;;GASG;AACH,MAAa,KAAK;IAiBhB,YAAoB,OAAO,GAAG,CAAC,EAAE,OAAgB;QAC/C,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;QACxB,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;IAC1B,CAAC;IACD,gDAAgD;IACzC,KAAK;QACV,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjD,CAAC;IACD,8CAA8C;IACvC,MAAM;QACX,OAAO,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;IAC7B,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,aAAa,CAAC,OAAe;QACzC,OAAO,IAAI,KAAK,CAAC,KAAK,CAAC,gBAAgB,CAAC,OAAO,CAAC,EAAE,OAAO,CAAC,CAAC;IAC7D,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,aAAa,CAAC,OAAe;QACzC,OAAO,IAAI,KAAK,CAAC,OAAO,CAAC,CAAC;IAC5B,CAAC;IACD;;OAEG;IACI,MAAM,CAAC,gBAAgB;QAC5B,OAAO,IAAI,KAAK,CAAC,QAAQ,CAAC,iBAAiB,CAAC,CAAC;IAC/C,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,iBAAiB,CAAC,MAAa,EAAE,QAAgB,EAAE,MAAa;QAC5E,OAAO,IAAI,KAAK,CAAC,QAAQ,CAAC,WAAW,CAAC,MAAM,CAAC,OAAO,EAAE,QAAQ,EAAE,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC;IACnF,CAAC;IACD;;;OAGG;IACI,WAAW,CAAC,KAAa;QAC9B,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,OAAO,GAAG,KAAK,CAAC,CAAC;IACzC,CAAC;IACD;;;OAGG;IACI,UAAU,CAAC,OAAe;QAC/B,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;QACxB,IAAI,CAAC,QAAQ,GAAG,SAAS,CAAC;IAC5B,CAAC;IACD;;;OAGG;IACI,UAAU,CAAC,OAAe;QAC/B,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC;QAChD,IAAI,CAAC,QAAQ,GAAG,OAAO,CAAC;IAC1B,CAAC;IACD,yCAAyC;IAClC,MAAM,CAAC,SAAS;QACrB,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,EAAE,GAAG,GAAG,EAAE,KAAK,CAAC,CAAC;IACzC,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,WAAW,CAAC,SAAiB,EAAE,WAAmB;QAC9D,OAAO,IAAI,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,SAAS,EAAE,WAAW,CAAC,CAAC,CAAC;IACvD,CAAC;IACD;;;OAGG;IACI,OAAO,CAAC,KAAY;QACzB,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,QAAQ,CAAC;QAC/B,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,QAAQ,CAAC;IACjC,CAAC;IACD;;;;;;;OAOG;IACI,WAAW,CAAC,IAAiB,EAAE,iBAA0B;QAC9D,IAAI,CAAC,QAAQ,GAAG,iBAAiB,CAAC,CAAC,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1D,IAAI,CAAC,IAAI;YACP,OAAO;QACT,IAAI,OAAO,IAAI,KAAK,QAAQ,EAAE;YAC5B,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC;SACvB;aAAM,IAAI,OAAQ,IAAY,CAAC,OAAO,KAAK,QAAQ,EAAE;YACpD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,OAAO,CAAC,CAAC;SACxC;aAAM,IAAI,OAAQ,IAAY,CAAC,QAAQ,KAAK,QAAQ,EAAE;YACrD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,QAAQ,CAAC,CAAC;SACzC;aAAM,IAAI,OAAQ,IAAY,CAAC,OAAO,KAAK,QAAQ,EAAE;YACpD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,OAAO,CAAC,CAAC;SACxC;aAAM,IAAI,OAAQ,IAAY,CAAC,QAAQ,KAAK,QAAQ,EAAE;YACrD,IAAI,CAAC,UAAU,CAAE,IAAY,CAAC,QAAQ,CAAC,CAAC;SACzC;IACH,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,QAAQ,CAAC,IAAiB,EAAE,iBAA0B;QAClE,MAAM,GAAG,GAAG,IAAI,KAAK,EAAE,CAAC;QACxB,GAAG,CAAC,WAAW,CAAC,IAAI,EAAE,iBAAiB,CAAC,CAAC;QACzC,OAAO,GAAG,CAAC;IACb,CAAC;IACD,+DAA+D;IACxD,MAAM;QACX,OAAO,IAAI,CAAC,OAAO,CAAC;IACtB,CAAC;IACD,yEAAyE;IAClE,aAAa;QAClB,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,OAAO,EAAE,CAAC;IACnC,CAAC;IACD,4CAA4C;IAC5C,IAAW,OAAO;QAChB,OAAO,IAAI,CAAC,QAAQ,CAAC;IACvB,CAAC;IACD,4CAA4C;IAC5C,IAAW,OAAO;QAChB,OAAO,IAAI,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,KAAK,CAAC,gBAAgB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAC7F,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,gBAAgB,CAAC,OAAe;QAC5C,OAAO,OAAO,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;IACjC,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,gBAAgB,CAAC,OAAe;QAC5C,IAAI,OAAO,GAAG,CAAC;YACb,OAAO,CAAC,KAAK,CAAC,gBAAgB,CAAC,CAAC,OAAO,CAAC,CAAC;QAC3C,8BAA8B;QAC9B,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QACnB,MAAM,MAAM,GAAG,KAAK,GAAG,EAAE,CAAC;QAC1B;;;;YAII;QACJ,IAAI,OAAO,IAAI,IAAI,GAAG,EAAE;YACtB,OAAO,MAAM,GAAG,OAAO,CAAC;QAC1B,IAAI,OAAO,GAAG,IAAI,GAAG,EAAE;YACrB,OAAO,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QAClD,IAAI,OAAO,IAAI,IAAI,GAAG,EAAE;YACtB,OAAO,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QAC7C,IAAI,OAAO,IAAI,IAAI,GAAG,EAAE;YACtB,OAAO,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QACnD,sDAAsD;QACtD,OAAO,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC,OAAO,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;IACnD,CAAC;IACD,qDAAqD;IAC9C,GAAG;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjC,CAAC;IACD,mDAAmD;IAC5C,GAAG;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjC,CAAC;IACD,sDAAsD;IAC/C,GAAG;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IACjC,CAAC;IACD,kEAAkE;IAC3D,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,OAAO,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,QAAQ,CAAC,gCAAgC,CAAC;IACxE,CAAC;IACD,iDAAiD;IAC1C,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,IAAI,QAAQ,CAAC,iBAAiB,CAAC;IAC7E,CAAC;IACD,yCAAyC;IACzC,IAAW,YAAY;QACrB,OAAO,KAAK,CAAC,mBAAmB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAClD,CAAC;IACD,+DAA+D;IAC/D,IAAW,YAAY;QACrB,OAAO,KAAK,CAAC,mBAAmB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAClD,CAAC;IACD,yDAAyD;IAClD,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,KAAK,CAAC;YACrB,IAAI,OAAO,GAAG,MAAM;gBAClB,OAAO,OAAO,CAAC;YACjB,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,MAAM,CAAC,CAAC;YAChD,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,MAAM,OAAO,GAAG,KAAK,CAAC,mBAAmB,CAAC,CAAC,OAAO,CAAC,CAAC;YACpD,OAAO,KAAK,GAAG,OAAO,CAAC;SACxB;QACD,sCAAsC;QACtC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,mDAAmD;IAC5C,MAAM,CAAC,sBAAsB,CAAC,OAAe;QAClD,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,KAAK;YAC5B,OAAO,OAAO,CAAC;QACjB,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,KAAK,CAAC;YACrB,MAAM,UAAU,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,OAAO,GAAG,KAAK,CAAC,GAAG,MAAM,CAAC,CAAC;YAC9D,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,OAAO,CAAC,KAAK,CAAC,sBAAsB,CAAC,CAAC,OAAO,CAAC,CAAC;SAChD;QACD,iCAAiC;QACjC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,yDAAyD;IAClD,MAAM,CAAC,mBAAmB,CAAC,OAAe;QAC/C,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;YAC7B,IAAI,OAAO,GAAG,MAAM;gBAClB,OAAO,OAAO,CAAC;YACjB,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,OAAO,GAAG,MAAM,CAAC,CAAC;YAChD,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,OAAO,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,KAAK,CAAC,mBAAmB,CAAC,CAAC,OAAO,CAAC,CAAC;SAC5D;QACD,iCAAiC;QACjC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,0DAA0D;IACnD,MAAM,CAAC,0BAA0B,CAAC,OAAe;QACtD,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,IAAI,CAAC,EAAE;YAC9B,OAAO,OAAO,CAAC;QACjB,IAAI,OAAO,IAAI,CAAC,EAAE;YAChB,MAAM,MAAM,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;YAC7B,MAAM,UAAU,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,OAAO,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,CAAC;YAChE,OAAO,OAAO,GAAG,UAAU,GAAG,MAAM,CAAC;SACtC;aAAM,IAAI,OAAO,GAAG,CAAC,EAAE;YACtB,qBAAqB;YACrB,OAAO,CAAC,KAAK,CAAC,0BAA0B,CAAC,CAAC,OAAO,CAAC,CAAC;SACpD;QACD,iCAAiC;QACjC,OAAO,CAAC,CAAC;IACX,CAAC;IACD,mEAAmE;IAC5D,MAAM,CAAC,IAAI;QAChB,OAAO,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC;IACtB,CAAC;IACD,yCAAyC;IACzC,IAAW,WAAW;QACpB,OAAO,IAAI,CAAC,OAAO,KAAK,CAAC,CAAC;IAC5B,CAAC;IACD,uFAAuF;IACvF,IAAW,YAAY;QACrB,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,GAAG,QAAQ,CAAC,iBAAiB,CAAC;IAC7D,CAAC;IACD,wGAAwG;IACxG,IAAW,wBAAwB;QACjC,OAAO,KAAK,CAAC,mBAAmB,CAAC,IAAI,CAAC,OAAO,GAAG,GAAG,CAAC,CAAC;IACvD,CAAC;IACD,gEAAgE;IACzD,MAAM,CAAC,2BAA2B,CAAC,OAAe;QACvD,OAAO,KAAK,CAAC,aAAa,CAAC,KAAK,CAAC,mBAAmB,CAAC,OAAO,CAAC,CAAC,CAAC;IACjE,CAAC;IACD,mEAAmE;IAC5D,MAAM,CAAC,4BAA4B,CAAC,OAAe;QACxD,OAAO,KAAK,CAAC,aAAa,CAAC,KAAK,CAAC,sBAAsB,CAAC,OAAO,CAAC,CAAC,CAAC;IACpE,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,oCAAoC,CAAC,QAAgB,EAAE,QAAgB,EACnF,YAAoB,QAAQ,CAAC,iBAAiB;QAC9C,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,CAAC;QAC5C,IAAI,KAAK,IAAI,SAAS;YACpB,OAAO,IAAI,CAAC;QACd,MAAM,MAAM,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC;QAC7B,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,IAAI,SAAS;YACvC,OAAO,IAAI,CAAC;QACd,MAAM,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC;QAC7C,MAAM,MAAM,GAAG,KAAK,GAAG,SAAS,GAAG,MAAM,CAAC;QAC1C,OAAO,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,SAAS,CAAC;IACvC,CAAC;IACD;;;KAGC;IACM,0BAA0B,CAAC,KAAY;QAC5C,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;IAC3D,CAAC;IAED;;;;OAIG;IACI,6BAA6B,CAAC,KAAY,EAAE,YAAoB,QAAQ,CAAC,iBAAiB;QAC/F,OAAO,KAAK,CAAC,oCAAoC,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;IAC9F,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,iCAAiC,CAAC,QAAgB,EAAE,QAAgB,EAChF,YAAoB,QAAQ,CAAC,iBAAiB;QAC9C,OAAO,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,SAAS,CAAC;IACnD,CAAC;IACD;;;;OAIG;IACI,0BAA0B,CAAC,KAAY,EAAE,YAAoB,QAAQ,CAAC,iBAAiB;QAC5F,OAAO,KAAK,CAAC,iCAAiC,CAAC,IAAI,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;IAC3F,CAAC;IACD;;;;;OAKG;IACI,aAAa,CAAC,KAAY,EAAE,YAAoB,QAAQ,CAAC,iBAAiB;QAC/E,OAAO,IAAI,CAAC,0BAA0B,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;IAC3D,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,qBAAqB,CAAC,KAAa,EAAE,KAAa,EAAE,KAAa;QAC7E,OAAO,KAAK,GAAG,QAAQ,CAAC,0BAA0B;eAC7C,KAAK,GAAG,QAAQ,CAAC,0BAA0B;eAC3C,KAAK,GAAG,KAAK,IAAI,QAAQ,CAAC,wBAAwB,GAAG,KAAK,GAAG,KAAK,CAAC;IAC1E,CAAC;IACD;;;;;;OAMG;IACI,MAAM,CAAC,+BAA+B,CAAC,MAAc,EAAE,MAAc;QAC1E,MAAM,CAAC,GAAG,QAAQ,CAAC,YAAY,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;QAChD,IAAI,CAAC,GAAG,QAAQ,CAAC,mBAAmB,EAAE;YACpC,OAAO,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC,CAAC,YAAY;SACtD;aAAM;YACL;;;;eAIG;YACH,IAAI,IAAI,CAAC;YACT,IAAI,IAAI,GAAG,GAAG,CAAC;YACf,MAAM,KAAK,GAAG,MAAM,GAAG,CAAC,CAAC;YACzB,MAAM,KAAK,GAAG,MAAM,GAAG,CAAC,CAAC;YACzB,qEAAqE;YACrE,IAAI,KAAK,IAAI,GAAG,EAAE;gBAChB;;;;mBAIG;gBACH,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;gBACtC,uDAAuD;gBACvD,IAAI,GAAG,KAAK,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC;aAC7B;iBAAM;gBACL,2DAA2D;gBAC3D,IAAI,KAAK,GAAG,GAAG,EAAE;oBACf;;;;uBAIG;oBACH,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;oBACtC,2DAA2D;iBAC5D;qBAAM;oBACL,kDAAkD;oBAClD,IAAI,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;iBACxC;gBACD,uDAAuD;gBACvD,IAAI,GAAG,KAAK,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,kBAAkB;aAChD;YACD,OAAO,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,EAAE,OAAO,EAAE,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC;SAC9D;IACH,CAAC;IACD,8EAA8E;IACvE,MAAM,CAAC,gBAAgB,CAAC,KAAa,EAAE,YAAoB,QAAQ,CAAC,kBAAkB;QAC3F,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QACjC,IAAI,QAAQ,IAAI,SAAS;YACvB,OAAO,CAAC,CAAC;QACX,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,CAAC;QACjC,IAAI,CAAC,IAAI,SAAS;YAChB,OAAO,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;QAClC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,CAAC;QAC7B,IAAI,CAAC,IAAI,SAAS;YAChB,OAAO,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;QAClC,OAAO,KAAK,CAAC;IACf,CAAC;IACD;;;;;;;;;;;;;OAaG;IACI,MAAM,CAAC,gCAAgC,CAC5C,KAAa,EAAE,KAAa,EAAE,KAAa,EAAE,YAAqB,IAAI;QAGtE,MAAM,MAAM,GAAG,KAAK,GAAG,KAAK,CAAC;QAC7B,MAAM,MAAM,GAAG,GAAG,GAAG,KAAK,CAAC;QAC3B,IAAI,SAAS,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,QAAQ,CAAC,iBAAiB,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;YAClG,OAAO,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,EAAE,CAAC;QAC1C,OAAO,KAAK,CAAC,+BAA+B,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC/D,CAAC;IACD;;;;;;;;;OASG;IACI,MAAM,CAAC,wBAAwB,CACpC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU;QAEtE,MAAM,KAAK,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC1C,OAAO,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,qBAAqB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACnF,CAAC;IACD;;;;;;;;;;;;;;;;;OAiBG;IACI,MAAM,CAAC,gCAAgC,CAC5C,EAAU,EAAE,EAAU,EAAE,EAAU,EAClC,EAAU,EAAE,EAAU,EAAE,EAAU,EAClC,SAAiB,EAAE,SAAiB,EAAE,SAAiB,EACvD,mBAA4B,KAAK;QAEjC,MAAM,KAAK,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC1C,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC7B,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC7B,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC7B,MAAM,MAAM,GAAG,SAAS,GAAG,EAAE,GAAG,SAAS,GAAG,EAAE,GAAG,SAAS,GAAG,EAAE,CAAC;QAChE,MAAM,cAAc,GAAG,QAAQ,CAAC,aAAa,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAC1D,IAAI,MAAM,GAAG,GAAG,EAAE;YAChB,IAAI,gBAAgB,EAAE;gBACpB,wGAAwG;gBACxG,4EAA4E;gBAC5E,OAAO,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,cAAc,EAAE,CAAC,KAAK,CAAC,CAAC;aACrD;iBAAM;gBACL,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,cAAc,EAAE,KAAK,CAAC,CAAC;aAC3C;SACF;aAAM;YACL,OAAO,IAAI,CAAC,KAAK,CAAC,cAAc,EAAE,KAAK,CAAC,CAAC;SAC1C;IACH,CAAC;IACD;;;OAGG;IACI,uBAAuB,CAAC,QAAgB;QAC7C,IAAI,IAAI,CAAC,QAAQ,KAAK,SAAS,EAAE;YAC/B,IAAI,CAAC,QAAQ,IAAI,QAAQ,GAAG,KAAK,CAAC;YAClC,IAAI,CAAC,QAAQ,GAAG,KAAK,CAAC,gBAAgB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;SACvD;aAAM;YACL,IAAI,CAAC,QAAQ,IAAI,QAAQ,GAAG,KAAK,CAAC,UAAU,CAAC;SAC9C;IACH,CAAC;;AA/gBD,+DAA+D;AACxC,qBAAe,GAAG,mBAAmB,CAAC;AAC7D,8DAA8D;AACvC,oBAAc,GAAG,sBAAsB,CAAC;AAC/D,8DAA8D;AACvC,oBAAc,GAAG,uBAAuB,CAAC;AAChE,6DAA6D;AACtC,eAAS,GAAG,sBAAsB,CAAC;AAC1D,+DAA+D;AACxC,gBAAU,GAAG,sBAAsB,CAAC;AAC3D,qDAAqD;AAC9B,sBAAgB,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,cAAc,CAAC,CAAC;AACxE,qDAAqD;AAC9B,sBAAgB,GAAG,CAAC,KAAK,CAAC,cAAc,GAAG,IAAI,CAAC,CAAC;SAd7D,KAAK","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/** @packageDocumentation\r\n * @module CartesianGeometry\r\n */\r\n\r\nimport { AngleProps, BeJSONFunctions, Geometry, TrigValues } from \"../Geometry\";\r\n\r\n/**\r\n * An `Angle` carries the numeric value of an angle, with methods to allow (require!) callers to\r\n * be clear about whether their angle is degrees or radians.\r\n * * After the Angle object is created, the callers should not know or care whether it is stored in\r\n * `degrees` or `radians` because both are available if requested by caller.\r\n * * The various access method are named so that callers can specify whether untyped numbers passed in or\r\n * out are degrees or radians.\r\n * * Visualization can be found at https://www.itwinjs.org/sandbox/SaeedTorabi/AngleSweep\r\n * @public\r\n */\r\nexport class Angle implements BeJSONFunctions {\r\n  /** maximal accuracy value of pi/12 (15 degrees), in radians */\r\n  public static readonly piOver12Radians = 0.26179938779914946;\r\n  /** maximal accuracy value of pi/4 (45 degrees), in radians */\r\n  public static readonly piOver4Radians = 7.853981633974483e-001;\r\n  /** maximal accuracy value of pi/2 (90 degrees), in radians */\r\n  public static readonly piOver2Radians = 1.5707963267948966e+000;\r\n  /** maximal accuracy value of pi (180 degrees), in radians */\r\n  public static readonly piRadians = 3.141592653589793e+000;\r\n  /** maximal accuracy value of 2*pi (360 degrees), in radians */\r\n  public static readonly pi2Radians = 6.283185307179586e+000;\r\n  /** scale factor for converting radians to degrees */\r\n  public static readonly degreesPerRadian = (45.0 / Angle.piOver4Radians);\r\n  /** scale factor for converting degrees to radians */\r\n  public static readonly radiansPerDegree = (Angle.piOver4Radians / 45.0);\r\n  private _radians: number;\r\n  private _degrees?: number;\r\n  private constructor(radians = 0, degrees?: number) {\r\n    this._radians = radians;\r\n    this._degrees = degrees;\r\n  }\r\n  /** Return a new angle with the same content. */\r\n  public clone(): Angle {\r\n    return new Angle(this._radians, this._degrees);\r\n  }\r\n  /** Freeze this instance so it is read-only */\r\n  public freeze(): Readonly<this> {\r\n    return Object.freeze(this);\r\n  }\r\n  /**\r\n   * Return a new Angle object for angle given in degrees.\r\n   * @param degrees angle in degrees\r\n   */\r\n  public static createDegrees(degrees: number): Angle {\r\n    return new Angle(Angle.degreesToRadians(degrees), degrees);\r\n  }\r\n  /**\r\n   * Return a (new) Angle object for a value given in radians.\r\n   * @param radians angle in radians\r\n   */\r\n  public static createRadians(radians: number): Angle {\r\n    return new Angle(radians);\r\n  }\r\n  /**\r\n   * Return a new `Angle` object with the default \"small\" angle measurement specified by [[Geometry.smallAngleRadians]].\r\n   */\r\n  public static createSmallAngle(): Angle {\r\n    return new Angle(Geometry.smallAngleRadians);\r\n  }\r\n  /**\r\n   * Return a (new) Angle object that is interpolated between two inputs (based on a fraction)\r\n   * @param angle0 first angle in radians\r\n   * @param fraction the interpolation fraction\r\n   * @param angle1 second angle in radians\r\n   */\r\n  public static createInterpolate(angle0: Angle, fraction: number, angle1: Angle): Angle {\r\n    return new Angle(Geometry.interpolate(angle0.radians, fraction, angle1.radians));\r\n  }\r\n  /**\r\n   * Return a (new) Angle object, with angle scaled from existing angle.\r\n   * @param scale scale factor to apply to angle.\r\n   */\r\n  public cloneScaled(scale: number): Angle {\r\n    return new Angle(this.radians * scale);\r\n  }\r\n  /**\r\n   * Set this angle to a value given in radians.\r\n   * @param radians angle given in radians\r\n   */\r\n  public setRadians(radians: number) {\r\n    this._radians = radians;\r\n    this._degrees = undefined;\r\n  }\r\n  /**\r\n   * Set this angle to a value given in degrees.\r\n   * @param degrees angle given in degrees.\r\n   */\r\n  public setDegrees(degrees: number) {\r\n    this._radians = Angle.degreesToRadians(degrees);\r\n    this._degrees = degrees;\r\n  }\r\n  /** Create an angle for a full circle. */\r\n  public static create360(): Angle {\r\n    return new Angle(Math.PI * 2.0, 360.0);\r\n  }\r\n  /**\r\n   * Create a (strongly typed) Angle whose tangent is `numerator/denominator`, using the signs of both in\r\n   * determining the (otherwise ambiguous) quadrant.\r\n   * @param numerator numerator for tangent\r\n   * @param denominator denominator for tangent\r\n   */\r\n  public static createAtan2(numerator: number, denominator: number): Angle {\r\n    return new Angle(Math.atan2(numerator, denominator));\r\n  }\r\n  /**\r\n   * Copy all contents of `other` to this Angle.\r\n   * @param other source data\r\n   */\r\n  public setFrom(other: Angle) {\r\n    this._radians = other._radians;\r\n    this._degrees = other._degrees;\r\n  }\r\n  /**\r\n   * Set an Angle from a JSON object\r\n   * * A simple number is considered as degrees.\r\n   * * specified `json.degrees` or `json._degrees` is degree value.\r\n   * * specified `son.radians` or `json._radians` is radians value.\r\n   * @param json object from JSON.parse. If a number, value is in *DEGREES*\r\n   * @param defaultValRadians if json is undefined, default value in radians.\r\n   */\r\n  public setFromJSON(json?: AngleProps, defaultValRadians?: number) {\r\n    this._radians = defaultValRadians ? defaultValRadians : 0;\r\n    if (!json)\r\n      return;\r\n    if (typeof json === \"number\") {\r\n      this.setDegrees(json);\r\n    } else if (typeof (json as any).degrees === \"number\") {\r\n      this.setDegrees((json as any).degrees);\r\n    } else if (typeof (json as any)._degrees === \"number\") {\r\n      this.setDegrees((json as any)._degrees);\r\n    } else if (typeof (json as any).radians === \"number\") {\r\n      this.setRadians((json as any).radians);\r\n    } else if (typeof (json as any)._radians === \"number\") {\r\n      this.setRadians((json as any)._radians);\r\n    }\r\n  }\r\n  /**\r\n   * Create an Angle from a JSON object\r\n   * @param json object from JSON.parse. If a number, value is in *DEGREES*\r\n   * @param defaultValRadians if json is undefined, default value in radians.\r\n   * @return a new Angle\r\n   */\r\n  public static fromJSON(json?: AngleProps, defaultValRadians?: number): Angle {\r\n    const val = new Angle();\r\n    val.setFromJSON(json, defaultValRadians);\r\n    return val;\r\n  }\r\n  /** Convert an Angle to a JSON object as a number in degrees */\r\n  public toJSON(): AngleProps {\r\n    return this.degrees;\r\n  }\r\n  /** Return a json object with radians keyword, e.g. `{ radians: 0.10}` */\r\n  public toJSONRadians(): AngleProps {\r\n    return { radians: this.radians };\r\n  }\r\n  /** Return the angle measured in radians. */\r\n  public get radians(): number {\r\n    return this._radians;\r\n  }\r\n  /** Return the angle measured in degrees. */\r\n  public get degrees(): number {\r\n    return this._degrees !== undefined ? this._degrees : Angle.radiansToDegrees(this._radians);\r\n  }\r\n  /**\r\n   * Convert an angle in degrees to radians.\r\n   * @param degrees angle in degrees\r\n   */\r\n  public static degreesToRadians(degrees: number): number {\r\n    return degrees * Math.PI / 180;\r\n  }\r\n  /**\r\n   * Convert an angle in radians to degrees.\r\n   * @param degrees angle in radians\r\n   */\r\n  public static radiansToDegrees(radians: number): number {\r\n    if (radians < 0)\r\n      return -Angle.radiansToDegrees(-radians);\r\n    // Now radians is positive ...\r\n    const pi = Math.PI;\r\n    const factor = 180.0 / pi;\r\n    /* the following if statements are for round-off reasons. The problem is that no IEEE number is\r\n      * an exact hit for any primary multiple of pi (90, 180, etc). The following is supposed to have\r\n      * a better chance that if the input was computed by direct assignment from 90, 180, etc degrees\r\n      * it will return exactly 90,180 etc.\r\n      */\r\n    if (radians <= 0.25 * pi)\r\n      return factor * radians;\r\n    if (radians < 0.75 * pi)\r\n      return 90.0 + 180 * ((radians - 0.5 * pi) / pi);\r\n    if (radians <= 1.25 * pi)\r\n      return 180.0 + 180 * ((radians - pi) / pi);\r\n    if (radians <= 1.75 * pi)\r\n      return 270.0 + 180 * ((radians - 1.5 * pi) / pi);\r\n    // all larger radians reference from 360 degrees (2PI)\r\n    return 360.0 + 180 * ((radians - 2.0 * pi) / pi);\r\n  }\r\n  /** Return the cosine of this Angle object's angle */\r\n  public cos(): number {\r\n    return Math.cos(this._radians);\r\n  }\r\n  /** Return the sine of this Angle object's angle */\r\n  public sin(): number {\r\n    return Math.sin(this._radians);\r\n  }\r\n  /** Return the tangent of this Angle object's angle */\r\n  public tan(): number {\r\n    return Math.tan(this._radians);\r\n  }\r\n  /** Test if a radians (absolute) value is nearly 2PI or larger! */\r\n  public static isFullCircleRadians(radians: number): boolean {\r\n    return Math.abs(radians) >= Geometry.fullCircleRadiansMinusSmallAngle;\r\n  }\r\n  /** Test if the radians value is a half circle */\r\n  public static isHalfCircleRadians(radians: number): boolean {\r\n    return Math.abs(Math.abs(radians) - Math.PI) <= Geometry.smallAngleRadians;\r\n  }\r\n  /** Test if the angle is a full circle */\r\n  public get isFullCircle(): boolean {\r\n    return Angle.isFullCircleRadians(this._radians);\r\n  }\r\n  /** Test if the angle is a half circle (in either direction) */\r\n  public get isHalfCircle(): boolean {\r\n    return Angle.isHalfCircleRadians(this._radians);\r\n  }\r\n  /** Adjust a radians value so it is positive in 0..360 */\r\n  public static adjustDegrees0To360(degrees: number): number {\r\n    if (degrees >= 0) {\r\n      const period = 360.0;\r\n      if (degrees < period)\r\n        return degrees;\r\n      const numPeriods = Math.floor(degrees / period);\r\n      return degrees - numPeriods * period;\r\n    } else if (degrees < 0) {\r\n      // negative angle ...\r\n      const radians = Angle.adjustDegrees0To360(-degrees);\r\n      return 360.0 - radians;\r\n    }\r\n    // fall through for Nan (disaster) !!!\r\n    return 0;\r\n  }\r\n  /** Adjust a radians value so it is in -180..180 */\r\n  public static adjustDegreesSigned180(degrees: number): number {\r\n    if (Math.abs(degrees) <= 180.0)\r\n      return degrees;\r\n    if (degrees >= 0) {\r\n      const period = 360.0;\r\n      const numPeriods = 1 + Math.floor((degrees - 180.0) / period);\r\n      return degrees - numPeriods * period;\r\n    } else if (degrees < 0) {\r\n      // negative angle ...\r\n      return -Angle.adjustDegreesSigned180(-degrees);\r\n    }\r\n    // fall through for NaN disaster.\r\n    return 0;\r\n  }\r\n  /** Adjust a radians value so it is positive in 0..2Pi */\r\n  public static adjustRadians0To2Pi(radians: number): number {\r\n    if (radians >= 0) {\r\n      const period = Math.PI * 2.0;\r\n      if (radians < period)\r\n        return radians;\r\n      const numPeriods = Math.floor(radians / period);\r\n      return radians - numPeriods * period;\r\n    } else if (radians < 0) {\r\n      // negative angle ...\r\n      return Math.PI * 2.0 - Angle.adjustRadians0To2Pi(-radians);\r\n    }\r\n    // fall through for NaN disaster.\r\n    return 0;\r\n  }\r\n  /** Adjust a radians value so it is positive in -PI..PI */\r\n  public static adjustRadiansMinusPiPlusPi(radians: number): number {\r\n    if (Math.abs(radians) <= Math.PI)\r\n      return radians;\r\n    if (radians >= 0) {\r\n      const period = Math.PI * 2.0;\r\n      const numPeriods = 1 + Math.floor((radians - Math.PI) / period);\r\n      return radians - numPeriods * period;\r\n    } else if (radians < 0) {\r\n      // negative angle ...\r\n      return -Angle.adjustRadiansMinusPiPlusPi(-radians);\r\n    }\r\n    // fall through for NaN disaster.\r\n    return 0;\r\n  }\r\n  /** Return a (newly allocated) Angle object with value 0 radians */\r\n  public static zero(): Angle {\r\n    return new Angle(0);\r\n  }\r\n  /** Test if the angle is exactly zero. */\r\n  public get isExactZero(): boolean {\r\n    return this.radians === 0;\r\n  }\r\n  /** Test if the angle is almost zero (within tolerance `Geometry.smallAngleRadians`) */\r\n  public get isAlmostZero(): boolean {\r\n    return Math.abs(this.radians) < Geometry.smallAngleRadians;\r\n  }\r\n  /** Test if the angle is almost a north or south pole (within tolerance `Geometry.smallAngleRadians`) */\r\n  public get isAlmostNorthOrSouthPole(): boolean {\r\n    return Angle.isHalfCircleRadians(this.radians * 2.0);\r\n  }\r\n  /** Create an angle object with degrees adjusted into 0..360. */\r\n  public static createDegreesAdjustPositive(degrees: number): Angle {\r\n    return Angle.createDegrees(Angle.adjustDegrees0To360(degrees));\r\n  }\r\n  /** Create an angle object with degrees adjusted into -180..180. */\r\n  public static createDegreesAdjustSigned180(degrees: number): Angle {\r\n    return Angle.createDegrees(Angle.adjustDegreesSigned180(degrees));\r\n  }\r\n  /**\r\n   * Test if two radians values are equivalent, allowing shift by full circle (i.e. by a multiple of `2*PI`)\r\n   * @param radiansA first radians value\r\n   * @param radiansB second radians value\r\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\r\n   */\r\n  public static isAlmostEqualRadiansAllowPeriodShift(radiansA: number, radiansB: number,\r\n    radianTol: number = Geometry.smallAngleRadians): boolean {\r\n    const delta = Math.abs(radiansA - radiansB);\r\n    if (delta <= radianTol)\r\n      return true;\r\n    const period = Math.PI * 2.0;\r\n    if (Math.abs(delta - period) <= radianTol)\r\n      return true;\r\n    const numPeriod = Math.round(delta / period);\r\n    const delta1 = delta - numPeriod * period;\r\n    return Math.abs(delta1) <= radianTol;\r\n  }\r\n  /**\r\n * Test if this angle has magnitude no greater than that of `other`.\r\n * @param other the other angle\r\n */\r\n  public isMagnitudeLessThanOrEqual(other: Angle): boolean {\r\n    return Math.abs(this.radians) <= Math.abs(other.radians);\r\n  }\r\n\r\n  /**\r\n   * Test if this angle and other are equivalent, allowing shift by full circle (i.e., multiples of `2 * PI`).\r\n   * @param other the other angle\r\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\r\n   */\r\n  public isAlmostEqualAllowPeriodShift(other: Angle, radianTol: number = Geometry.smallAngleRadians): boolean {\r\n    return Angle.isAlmostEqualRadiansAllowPeriodShift(this._radians, other._radians, radianTol);\r\n  }\r\n  /**\r\n   * Test if two angle (in radians)  almost equal, NOT allowing shift by full circle (i.e., multiples of `2 * PI`).\r\n   * @param radiansA first radians value\r\n   * @param radiansB second radians value\r\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\r\n   */\r\n  public static isAlmostEqualRadiansNoPeriodShift(radiansA: number, radiansB: number,\r\n    radianTol: number = Geometry.smallAngleRadians): boolean {\r\n    return Math.abs(radiansA - radiansB) < radianTol;\r\n  }\r\n  /**\r\n   * Test if two this angle and other are almost equal, NOT allowing shift by full circle (i.e., multiples of `2 * PI`).\r\n   * @param other the other angle\r\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\r\n   */\r\n  public isAlmostEqualNoPeriodShift(other: Angle, radianTol: number = Geometry.smallAngleRadians): boolean {\r\n    return Angle.isAlmostEqualRadiansNoPeriodShift(this._radians, other._radians, radianTol);\r\n  }\r\n  /**\r\n   * Test if two this angle and other are almost equal, NOT allowing shift by full circle (i.e., multiples of `2 * PI`).\r\n   * * This function is same as isAlmostEqualRadiansNoPeriodShift. Please use isAlmostEqualRadiansNoPeriodShift.\r\n   * @param other the other angle\r\n   * @param radianTol radian tolerance with default value of Geometry.smallAngleRadians\r\n   */\r\n  public isAlmostEqual(other: Angle, radianTol: number = Geometry.smallAngleRadians): boolean {\r\n    return this.isAlmostEqualNoPeriodShift(other, radianTol);\r\n  }\r\n  /**\r\n   * Test if dot product values indicate non-zero length perpendicular vectors.\r\n   * @param dotUU dot product of vectorU with itself\r\n   * @param dotVV dot product of vectorV with itself\r\n   * @param dotUV dot product of vectorU with vectorV\r\n   */\r\n  public static isPerpendicularDotSet(dotUU: number, dotVV: number, dotUV: number): boolean {\r\n    return dotUU > Geometry.smallMetricDistanceSquared\r\n      && dotVV > Geometry.smallMetricDistanceSquared\r\n      && dotUV * dotUV <= Geometry.smallAngleRadiansSquared * dotUU * dotVV;\r\n  }\r\n  /**\r\n   * Return cosine, sine, and radians for the half angle of a \"cosine,sine\" pair.\r\n   * * This function assumes the input arguments are related to an angle between -PI and PI\r\n   * * This function returns an angle between -PI and PI\r\n   * @param rCos2A cosine value (scaled by radius) for initial angle.\r\n   * @param rSin2A sine value (scaled by radius) for final angle.\r\n   */\r\n  public static trigValuesToHalfAngleTrigValues(rCos2A: number, rSin2A: number): TrigValues {\r\n    const r = Geometry.hypotenuseXY(rCos2A, rSin2A);\r\n    if (r < Geometry.smallMetricDistance) {\r\n      return { c: 1.0, s: 0.0, radians: 0.0 }; // angle = 0\r\n    } else {\r\n      /* If the caller really gave you sine and cosine values, r should be 1.  However,\r\n       * to allow scaled values -- e.g. the x and y components of any vector -- we normalize\r\n       * right here. This adds an extra sqrt and 2 divides to the whole process, but improves\r\n       * both the usefulness and robustness of the computation.\r\n       */\r\n      let cosA;\r\n      let sinA = 0.0;\r\n      const cos2A = rCos2A / r;\r\n      const sin2A = rSin2A / r;\r\n      // Original angle in NE and SE quadrants. Half angle in same quadrant\r\n      if (cos2A >= 0.0) {\r\n        /*\r\n         * We know cos2A = (cosA)^2 - (sinA)^2 and 1 = (cosA)^2 + (sinA)^2\r\n         * so 1 + cos2A = 2(cosA)^2 and therefore, cosA = sqrt((1+cos2A)/2)\r\n         * cosine is positive in NE and SE quadrants so we use +sqrt\r\n         */\r\n        cosA = Math.sqrt(0.5 * (1.0 + cos2A));\r\n        // We know sin2A = 2 sinA cosA so sinA = sin2A/(2*cosA)\r\n        sinA = sin2A / (2.0 * cosA);\r\n      } else {\r\n        // Original angle in NW quadrant. Half angle in NE quadrant\r\n        if (sin2A > 0.0) {\r\n          /*\r\n           * We know cos2A = (cosA)^2 - (sinA)^2 and 1 = (cosA)^2 + (sinA)^2\r\n           * so 1 - cos2A = 2(sinA)^2 and therefore, sinA = sqrt((1-cos2A)/2)\r\n           * sine is positive in NE quadrant so we use +sqrt\r\n           */\r\n          sinA = Math.sqrt(0.5 * (1.0 - cos2A));\r\n          // Original angle in SW quadrant. Half angle in SE quadrant\r\n        } else {\r\n          // sine is negative in SE quadrant so we use -sqrt\r\n          sinA = -Math.sqrt(0.5 * (1.0 - cos2A));\r\n        }\r\n        // We know sin2A = 2 sinA cosA so cosA = sin2A/(2*sinA)\r\n        cosA = sin2A / (2.0 * sinA); // always positive\r\n      }\r\n      return { c: cosA, s: sinA, radians: Math.atan2(sinA, cosA) };\r\n    }\r\n  }\r\n  /** If value is close to -1, -0.5, 0, 0.5, 1, adjust it to the exact value. */\r\n  public static cleanupTrigValue(value: number, tolerance: number = Geometry.smallFloatingPoint): number {\r\n    const absValue = Math.abs(value);\r\n    if (absValue <= tolerance)\r\n      return 0;\r\n    let a = Math.abs(absValue - 0.5);\r\n    if (a <= tolerance)\r\n      return value < 0.0 ? -0.5 : 0.5;\r\n    a = Math.abs(absValue - 1.0);\r\n    if (a <= tolerance)\r\n      return value < 0.0 ? -1.0 : 1.0;\r\n    return value;\r\n  }\r\n  /**\r\n   * Return the half angle cosine, sine, and radians for given dot products between vectors. The vectors define\r\n   * an ellipse using x(t) = c + U cos(t) + V sin(t) so U and V are at angle t=0 degree and t=90 degree. The\r\n   * half angle t0 is an angle such that x(t0) is one of the ellipse semi-axis.\r\n   * * This construction arises e.g. in `Arc3d.toScaledMatrix3d`.\r\n   * * Given ellipse x(t) = c + U cos(t) + V sin(t), find t0 such that radial vector W(t0) = x(t0) - c is\r\n   * perpendicular to the ellipse.\r\n   * * Then 0 = W(t0).x'(t0) = (U cos(t0) + V sin(t0)).(V cos(t0) - U sin(t0)) = U.V cos(2t0) + 0.5 (V.V - U.U) sin(2t0)\r\n   * implies sin(2t0) / cos(2t0) = 2 U.V / (U.U - V.V), i.e., t0 can be computed given the three dot products on the RHS.\r\n   * math details can be found at docs/learning/geometry/Angle.md\r\n   * @param dotUU dot product of vectorU with itself\r\n   * @param dotVV dot product of vectorV with itself\r\n   * @param dotUV dot product of vectorU with vectorV\r\n   */\r\n  public static dotProductsToHalfAngleTrigValues(\r\n    dotUU: number, dotVV: number, dotUV: number, favorZero: boolean = true,\r\n  ): TrigValues {\r\n\r\n    const cos2t0 = dotUU - dotVV;\r\n    const sin2t0 = 2.0 * dotUV;\r\n    if (favorZero && Math.abs(sin2t0) < Geometry.smallAngleRadians * (Math.abs(dotUU) + Math.abs(dotVV)))\r\n      return { c: 1.0, s: 0.0, radians: 0.0 };\r\n    return Angle.trigValuesToHalfAngleTrigValues(cos2t0, sin2t0);\r\n  }\r\n  /**\r\n   * Returns the angle between two vectors, with the vectors given as xyz components\r\n   * * The returned angle is between 0 and PI\r\n   * @param ux x component of vector u\r\n   * @param uy y component of vector u\r\n   * @param uz z component of vector u\r\n   * @param vx x component of vector v\r\n   * @param vy y component of vector v\r\n   * @param vz z component of vector v\r\n   */\r\n  public static radiansBetweenVectorsXYZ(\r\n    ux: number, uy: number, uz: number, vx: number, vy: number, vz: number,\r\n  ): number {\r\n    const uDotV = ux * vx + uy * vy + uz * vz;\r\n    return Math.atan2(Geometry.crossProductMagnitude(ux, uy, uz, vx, vy, vz), uDotV);\r\n  }\r\n  /**\r\n   * Returns the angle between two vectors, with the vectors given as xyz components, and an up vector to resolve\r\n   * angle to a full 2PI range.\r\n   * * The returned angle is (-PI < radians <= PI) or (0 <= radians < 2 * PI)\r\n   * * The angle is in the plane of the U and V vectors.\r\n   * * The upVector determines a positive side of the plane but need not be strictly perpendicular to the plane.\r\n   * @param ux x component of vector u\r\n   * @param uy y component of vector u\r\n   * @param uz z component of vector u\r\n   * @param vx x component of vector v\r\n   * @param vy y component of vector v\r\n   * @param vz z component of vector v\r\n   * @param upVectorX x component of vector to positive side of plane.\r\n   * @param upVectorY y component of vector to positive side of plane.\r\n   * @param upVectorZ z component of vector to positive side of plane.\r\n   * @param adjustToAllPositive if true, return strictly non-negative sweep (0 <= radians < 2*PI). If false, return\r\n   * signed (-PI < radians <= PI)\r\n   */\r\n  public static orientedRadiansBetweenVectorsXYZ(\r\n    ux: number, uy: number, uz: number,\r\n    vx: number, vy: number, vz: number,\r\n    upVectorX: number, upVectorY: number, upVectorZ: number,\r\n    adjustToPositive: boolean = false,\r\n  ): number {\r\n    const uDotV = ux * vx + uy * vy + uz * vz;\r\n    const wx = uy * vz - uz * vy;\r\n    const wy = uz * vx - ux * vz;\r\n    const wz = ux * vy - uy * vx;\r\n    const upDotW = upVectorX * wx + upVectorY * wy + upVectorZ * wz;\r\n    const crossMagnitude = Geometry.hypotenuseXYZ(wx, wy, wz);\r\n    if (upDotW < 0.0) {\r\n      if (adjustToPositive) {\r\n        // The turn is greater than 180 degrees.  Take a peculiarly oriented atan2 to get the excess-180 part as\r\n        // addition to PI. This gives the smoothest numerical transition passing PI.\r\n        return Math.PI + Math.atan2(crossMagnitude, -uDotV);\r\n      } else {\r\n        return -Math.atan2(crossMagnitude, uDotV);\r\n      }\r\n    } else {\r\n      return Math.atan2(crossMagnitude, uDotV);\r\n    }\r\n  }\r\n  /**\r\n   * Add a multiple of a full circle angle (360 degrees, 2PI) in place.\r\n   * @param multiple multiplier factor\r\n   */\r\n  public addMultipleOf2PiInPlace(multiple: number) {\r\n    if (this._degrees !== undefined) {\r\n      this._degrees += multiple * 360.0;\r\n      this._radians = Angle.degreesToRadians(this._degrees);\r\n    } else {\r\n      this._radians += multiple * Angle.pi2Radians;\r\n    }\r\n  }\r\n}\r\n"]}

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

@@ -21,17 +21,33 @@ export declare class CoincidentGeometryQuery {
      * * reevaluate the points as simple interpolation between given points.
      */
     static assignDetailInterpolatedFractionsAndPoints(detail: CurveLocationDetail, f0: number, f1: number, pointA: Point3d, pointB: Point3d, swap?: boolean): void;
-    /** Return a curve location detail with projection of a `spacePoint` to the line segment with `pointA` and `pointB`
-     *
-     */
+    /** Return a curve location detail with projection of a `spacePoint` to the line segment with `pointA` and `pointB` */
     projectPointToSegmentXY(spacePoint: Point3d, pointA: Point3d, pointB: Point3d): CurveLocationDetail;
     /**
-     * * project `pointA0` and `pointA1` onto the segment with `pointB0` and `pointB1`
-     * * In the returned detail pair, the `detailA` has fractions along segmentA, and `detailB` has fractions along segment B
+     * Given a detail pair representing the projection of each of two colinear line segments onto the other,
+     * clamp the details (in place) to the line segments' endpoints according to the given flags.
+     * @param overlap segment overlap as returned by [[coincidentSegmentRangeXY]], modified on return
+     * @param pointA0 start point of segment A
+     * @param pointA1 end point of segment A
+     * @param pointB0 start point of segment B
+     * @param pointB1 end point of segment B
+     * @param extendA0 whether to extend segment A beyond its start
+     * @param extendA1 whether to extend segment A beyond its end
+     * @param extendB0 whether to extend segment B beyond its start
+     * @param extendB1 whether to extend segment B beyond its end
+     * @return reference to the input clamped in place, or undefined (leaving interval untouched) if clamping would result in empty interval.
+     */
+    clampCoincidentOverlapToSegmentBounds(overlap: CurveLocationDetailPair, pointA0: Point3d, pointA1: Point3d, pointB0: Point3d, pointB1: Point3d, extendA0?: boolean, extendA1?: boolean, extendB0?: boolean, extendB1?: boolean): CurveLocationDetailPair | undefined;
+    /**
+     * Compute whether two line segments have a coincident overlap in xy.
+     * * Project `pointA0` and `pointA1` onto the line formed by `pointB0` and `pointB1` and vice versa
+     * * If all projection distances are sufficiently small, return a detail pair recording the coincident interval, optionally clipped to segment bounds.
      * @param pointA0 start point of segment A
      * @param pointA1 end point of segment A
      * @param pointB0 start point of segment B
      * @param pointB1 end point of segment B
+     * @param restrictToBounds whether to clip the coincident segment details to the segment bounds
+     * @return detail pair for the coincident interval (`detailA` has fractions along segment A, and `detailB` has fractions along segment B), or undefined if no coincidence
      */
     coincidentSegmentRangeXY(pointA0: Point3d, pointA1: Point3d, pointB0: Point3d, pointB1: Point3d, restrictToBounds?: boolean): CurveLocationDetailPair | undefined;
     /**

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

@@ -1 +1 @@
-{"version":3,"file":"CoincidentGeometryOps.d.ts","sourceRoot":"","sources":["../../../src/geometry3d/CoincidentGeometryOps.ts"],"names":[],"mappings":"AASA,OAAO,EAAE,KAAK,EAAE,MAAM,gBAAgB,CAAC;AACvC,OAAO,EAAE,mBAAmB,EAAE,uBAAuB,EAAE,MAAM,8BAA8B,CAAC;AAI5F,OAAO,EAAE,OAAO,EAAY,MAAM,mBAAmB,CAAC;AAGtD;;;;;GAKG;AACH,qBAAa,uBAAuB;IAClC,OAAO,CAAC,QAAQ,CAAC,CAAW;IAC5B,OAAO,CAAC,QAAQ,CAAC,CAAW;IAC5B,OAAO,CAAC,OAAO,CAAC,CAAU;IAC1B,OAAO,CAAC,OAAO,CAAC,CAAU;IAC1B,OAAO,CAAC,UAAU,CAAS;IAC3B,IAAW,SAAS,IAAI,MAAM,CAE7B;IACD,OAAO;WAGO,MAAM,CAAC,SAAS,GAAE,MAAqC,GAAG,uBAAuB;IAG/F;;;OAGG;WACW,0CAA0C,CAAC,MAAM,EAAE,mBAAmB,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAC1G,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,IAAI,GAAE,OAAe;IAazD;;OAEG;IACI,uBAAuB,CAAC,UAAU,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,GAAG,mBAAmB;IAS1G;;;;;;;OAOG;IACI,wBAAwB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,gBAAgB,GAAE,OAAc,GAAG,uBAAuB,GAAG,SAAS;IA2D9K;;;;;;;;OAQG;IACH,OAAO,CAAC,gBAAgB;IAaxB,OAAO,CAAC,gBAAgB;IAQxB;;;;;;OAMG;IACI,2BAA2B,CAAC,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,iBAAiB,GAAE,OAAc,GAAG,uBAAuB,EAAE,GAAG,SAAS;CAgEvI"}
+{"version":3,"file":"CoincidentGeometryOps.d.ts","sourceRoot":"","sources":["../../../src/geometry3d/CoincidentGeometryOps.ts"],"names":[],"mappings":"AASA,OAAO,EAAE,KAAK,EAAE,MAAM,gBAAgB,CAAC;AACvC,OAAO,EAAE,mBAAmB,EAAE,uBAAuB,EAAE,MAAM,8BAA8B,CAAC;AAI5F,OAAO,EAAE,OAAO,EAAY,MAAM,mBAAmB,CAAC;AAGtD;;;;;GAKG;AACH,qBAAa,uBAAuB;IAClC,OAAO,CAAC,QAAQ,CAAC,CAAW;IAC5B,OAAO,CAAC,QAAQ,CAAC,CAAW;IAC5B,OAAO,CAAC,OAAO,CAAC,CAAU;IAC1B,OAAO,CAAC,OAAO,CAAC,CAAU;IAC1B,OAAO,CAAC,UAAU,CAAS;IAC3B,IAAW,SAAS,IAAI,MAAM,CAE7B;IACD,OAAO;WAGO,MAAM,CAAC,SAAS,GAAE,MAAqC,GAAG,uBAAuB;IAG/F;;;OAGG;WACW,0CAA0C,CAAC,MAAM,EAAE,mBAAmB,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAC1G,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,IAAI,GAAE,OAAe;IAazD,sHAAsH;IAC/G,uBAAuB,CAAC,UAAU,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,GAAG,mBAAmB;IAS1G;;;;;;;;;;;;;OAaG;IACI,qCAAqC,CAAC,OAAO,EAAE,uBAAuB,EAC3E,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EACtE,QAAQ,GAAE,OAAe,EAAE,QAAQ,GAAE,OAAe,EAAE,QAAQ,GAAE,OAAe,EAAE,QAAQ,GAAE,OAAe,GACzG,uBAAuB,GAAG,SAAS;IAoDtC;;;;;;;;;;OAUG;IACI,wBAAwB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,gBAAgB,GAAE,OAAc,GAAG,uBAAuB,GAAG,SAAS;IAuB9K;;;;;;;;OAQG;IACH,OAAO,CAAC,gBAAgB;IAaxB,OAAO,CAAC,gBAAgB;IAQxB;;;;;;OAMG;IACI,2BAA2B,CAAC,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,iBAAiB,GAAE,OAAc,GAAG,uBAAuB,EAAE,GAAG,SAAS;CAgEvI"}

package/lib/esm/geometry3d/CoincidentGeometryOps.js

@@ -43,9 +43,7 @@ export class CoincidentGeometryQuery {
         detail.point = pointA.interpolate(detail.fraction, pointB, detail.point);
         detail.point1 = pointA.interpolate(detail.fraction1, pointB, detail.point1);
     }
-    /** Return a curve location detail with projection of a `spacePoint` to the line segment with `pointA` and `pointB`
-     *
-     */
+    /** Return a curve location detail with projection of a `spacePoint` to the line segment with `pointA` and `pointB` */
     projectPointToSegmentXY(spacePoint, pointA, pointB) {
         this._vectorU = Vector3d.createStartEnd(pointA, pointB, this._vectorU);
         this._vectorV = Vector3d.createStartEnd(pointA, spacePoint, this._vectorV);
@@ -55,12 +53,78 @@ export class CoincidentGeometryQuery {
         return CurveLocationDetail.createCurveFractionPoint(undefined, fraction, pointA.interpolate(fraction, pointB));
     }
     /**
-     * * project `pointA0` and `pointA1` onto the segment with `pointB0` and `pointB1`
-     * * In the returned detail pair, the `detailA` has fractions along segmentA, and `detailB` has fractions along segment B
+     * Given a detail pair representing the projection of each of two colinear line segments onto the other,
+     * clamp the details (in place) to the line segments' endpoints according to the given flags.
+     * @param overlap segment overlap as returned by [[coincidentSegmentRangeXY]], modified on return
+     * @param pointA0 start point of segment A
+     * @param pointA1 end point of segment A
+     * @param pointB0 start point of segment B
+     * @param pointB1 end point of segment B
+     * @param extendA0 whether to extend segment A beyond its start
+     * @param extendA1 whether to extend segment A beyond its end
+     * @param extendB0 whether to extend segment B beyond its start
+     * @param extendB1 whether to extend segment B beyond its end
+     * @return reference to the input clamped in place, or undefined (leaving interval untouched) if clamping would result in empty interval.
+     */
+    clampCoincidentOverlapToSegmentBounds(overlap, pointA0, pointA1, pointB0, pointB1, extendA0 = false, extendA1 = false, extendB0 = false, extendB1 = false) {
+        const rangeA = Segment1d.create(overlap.detailA.fraction, overlap.detailA.hasFraction1 ? overlap.detailA.fraction1 : overlap.detailA.fraction);
+        const rangeB = Segment1d.create(overlap.detailB.fraction, overlap.detailB.hasFraction1 ? overlap.detailB.fraction1 : overlap.detailB.fraction);
+        const reversed = rangeA.signedDelta() < 0.0;
+        const updateIntervalFromRangesAndInterpolatedPoints = () => {
+            const a0 = rangeA.x0;
+            const a1 = rangeA.x1;
+            const b0 = rangeB.x0;
+            const b1 = rangeB.x1;
+            CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(overlap.detailA, a0, a1, pointA0, pointA1, a0 > a1);
+            CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(overlap.detailB, b0, b1, pointB0, pointB1, b0 > b1);
+            return overlap;
+        };
+        const haveIntervalA = rangeA.clampDirectedTo01(!extendA0, !extendA1, false);
+        const haveIntervalB = rangeB.clampDirectedTo01(!extendB0, !extendB1, false);
+        if (haveIntervalA && haveIntervalB) {
+            if (Geometry.isAlmostEqualNumber(rangeA.absoluteDelta(), rangeB.absoluteDelta(), Geometry.smallFraction))
+                return updateIntervalFromRangesAndInterpolatedPoints(); // intersection of partially clamped ranges
+            else if (rangeA.clampDirectedTo01(true, true, false) && rangeB.clampDirectedTo01(true, true, false))
+                return updateIntervalFromRangesAndInterpolatedPoints(); // intersection of fully clamped ranges
+        }
+        const collapseToSingleton = (pointA, pointB, atStartA, atStartB) => {
+            pointA.clone(overlap.detailA.point);
+            pointB.clone(overlap.detailB.point);
+            overlap.detailA.fraction = atStartA ? 0.0 : 1.0;
+            overlap.detailB.fraction = atStartB ? 0.0 : 1.0;
+            overlap.detailA.collapseToStart();
+            overlap.detailB.collapseToStart();
+            return overlap;
+        };
+        const haveSingletonA = rangeA.clampDirectedTo01(true, true, true);
+        const haveSingletonB = rangeB.clampDirectedTo01(true, true, true);
+        if (haveSingletonA && haveSingletonB) { // intersection is a single point
+            const point1 = overlap.detailA.point1 ?? overlap.detailA.point;
+            if (reversed) {
+                if (overlap.detailA.point.isAlmostEqual(pointA0, this.tolerance))
+                    return collapseToSingleton(pointA0, pointB0, true, true);
+                else if (point1.isAlmostEqual(pointA1, this.tolerance))
+                    return collapseToSingleton(pointA1, pointB1, false, false);
+            }
+            else {
+                if (point1.isAlmostEqual(pointA0, this.tolerance))
+                    return collapseToSingleton(pointA0, pointB1, true, false);
+                else if (overlap.detailA.point.isAlmostEqual(pointA1, this.tolerance))
+                    return collapseToSingleton(pointA1, pointB0, false, true);
+            }
+        }
+        return undefined; // no intersection
+    }
+    /**
+     * Compute whether two line segments have a coincident overlap in xy.
+     * * Project `pointA0` and `pointA1` onto the line formed by `pointB0` and `pointB1` and vice versa
+     * * If all projection distances are sufficiently small, return a detail pair recording the coincident interval, optionally clipped to segment bounds.
      * @param pointA0 start point of segment A
      * @param pointA1 end point of segment A
      * @param pointB0 start point of segment B
      * @param pointB1 end point of segment B
+     * @param restrictToBounds whether to clip the coincident segment details to the segment bounds
+     * @return detail pair for the coincident interval (`detailA` has fractions along segment A, and `detailB` has fractions along segment B), or undefined if no coincidence
      */
     coincidentSegmentRangeXY(pointA0, pointA1, pointB0, pointB1, restrictToBounds = true) {
         const detailA0OnB = this.projectPointToSegmentXY(pointA0, pointB0, pointB1);
@@ -79,46 +143,8 @@ export class CoincidentGeometryQuery {
         detailA0OnB.point1 = detailA1OnB.point; // capture -- detailA1OnB is not reused.
         detailB0OnA.fraction1 = detailB1OnA.fraction;
         detailB0OnA.point1 = detailB1OnA.point;
-        if (!restrictToBounds)
-            return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);
-        const segment = Segment1d.create(detailB0OnA.fraction, detailB0OnA.fraction1);
-        if (segment.clampDirectedTo01()) {
-            const f0 = segment.x0;
-            const f1 = segment.x1;
-            const h0 = detailB0OnA.inverseInterpolateFraction(f0);
-            const h1 = detailB0OnA.inverseInterpolateFraction(f1);
-            // recompute fractions and points..
-            CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(detailB0OnA, f0, f1, pointA0, pointA1, f0 > f1);
-            CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(detailA0OnB, h0, h1, pointB0, pointB1, h0 > h1);
-            return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);
-        }
-        else {
-            if (segment.signedDelta() < 0.0) {
-                if (detailB0OnA.point.isAlmostEqual(pointA0, this.tolerance)) {
-                    detailB0OnA.collapseToStart();
-                    detailA0OnB.collapseToStart();
-                    return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);
-                }
-                if (detailB0OnA.point1.isAlmostEqual(pointA1, this.tolerance)) {
-                    detailB0OnA.collapseToEnd();
-                    detailA0OnB.collapseToEnd();
-                    return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);
-                }
-            }
-            else {
-                if (detailB0OnA.point.isAlmostEqual(pointA1, this.tolerance)) {
-                    detailB0OnA.collapseToStart();
-                    detailA0OnB.collapseToEnd();
-                    return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);
-                }
-                if (detailB0OnA.point1.isAlmostEqual(pointA0, this.tolerance)) {
-                    detailB0OnA.collapseToEnd();
-                    detailA0OnB.collapseToStart();
-                    return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);
-                }
-            }
-        }
-        return undefined;
+        const interval = CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);
+        return restrictToBounds ? this.clampCoincidentOverlapToSegmentBounds(interval, pointA0, pointA1, pointB0, pointB1) : interval;
     }
     /**
      * Create a CurveLocationDetailPair for a coincident interval of two overlapping curves

package/lib/esm/geometry3d/CoincidentGeometryOps.js.map

@@ -1 +1 @@
-{"version":3,"file":"CoincidentGeometryOps.js","sourceRoot":"","sources":["../../../src/geometry3d/CoincidentGeometryOps.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AACH,OAAO,EAAE,MAAM,EAAE,MAAM,qBAAqB,CAAC;AAE7C,OAAO,EAAE,mBAAmB,EAAE,uBAAuB,EAAE,MAAM,8BAA8B,CAAC;AAE5F,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AACvC,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAC1C,OAAO,EAAW,QAAQ,EAAE,MAAM,mBAAmB,CAAC;AACtD,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AAExC;;;;;GAKG;AACH,MAAM,OAAO,uBAAuB;IAMlC,IAAW,SAAS;QAClB,OAAO,IAAI,CAAC,UAAU,CAAC;IACzB,CAAC;IACD,YAAoB,YAAoB,QAAQ,CAAC,mBAAmB;QAClE,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;IAC9B,CAAC;IACM,MAAM,CAAC,MAAM,CAAC,YAAoB,QAAQ,CAAC,mBAAmB;QACnE,OAAO,IAAI,uBAAuB,CAAC,SAAS,CAAC,CAAC;IAChD,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,0CAA0C,CAAC,MAA2B,EAAE,EAAU,EAAE,EAAU,EAC1G,MAAe,EAAE,MAAe,EAAE,OAAgB,KAAK;QACvD,IAAI,IAAI,EAAE;YACR,MAAM,CAAC,QAAQ,GAAG,EAAE,CAAC;YACrB,MAAM,CAAC,SAAS,GAAG,EAAE,CAAC;SACvB;aAAM;YACL,MAAM,CAAC,QAAQ,GAAG,EAAE,CAAC;YACrB,MAAM,CAAC,SAAS,GAAG,EAAE,CAAC;SACvB;QACD,MAAM,CAAC,KAAK,GAAG,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,QAAQ,EAAE,MAAM,EAAE,MAAM,CAAC,KAAK,CAAC,CAAC;QACzE,MAAM,CAAC,MAAM,GAAG,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,SAAS,EAAE,MAAM,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;IAE9E,CAAC;IAED;;OAEG;IACI,uBAAuB,CAAC,UAAmB,EAAE,MAAe,EAAE,MAAe;QAClF,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC,cAAc,CAAC,MAAM,EAAE,MAAM,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QACvE,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC,cAAc,CAAC,MAAM,EAAE,UAAU,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QAC3E,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,YAAY,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QACxD,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,YAAY,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QACxD,MAAM,QAAQ,GAAG,QAAQ,CAAC,kBAAkB,CAAC,KAAK,EAAE,KAAK,EAAE,GAAG,CAAC,CAAC;QAChE,OAAO,mBAAmB,CAAC,wBAAwB,CAAC,SAAS,EAAE,QAAQ,EACrE,MAAM,CAAC,WAAW,CAAC,QAAQ,EAAE,MAAM,CAAC,CAAC,CAAC;IAC1C,CAAC;IACD;;;;;;;OAOG;IACI,wBAAwB,CAAC,OAAgB,EAAE,OAAgB,EAAE,OAAgB,EAAE,OAAgB,EAAE,mBAA4B,IAAI;QACtI,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QACnB,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QAEnB,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QACnB,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QAEnB,WAAW,CAAC,SAAS,GAAG,WAAW,CAAC,QAAQ,CAAC;QAC7C,WAAW,CAAC,MAAM,GAAG,WAAW,CAAC,KAAK,CAAC,CAAE,wCAAwC;QACjF,WAAW,CAAC,SAAS,GAAG,WAAW,CAAC,QAAQ,CAAC;QAC7C,WAAW,CAAC,MAAM,GAAG,WAAW,CAAC,KAAK,CAAC;QACvC,IAAI,CAAC,gBAAgB;YACnB,OAAO,uBAAuB,CAAC,aAAa,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;QAEzE,MAAM,OAAO,GAAG,SAAS,CAAC,MAAM,CAAC,WAAW,CAAC,QAAQ,EAAE,WAAW,CAAC,SAAS,CAAC,CAAC;QAC9E,IAAI,OAAO,CAAC,iBAAiB,EAAE,EAAE;YAC/B,MAAM,EAAE,GAAG,OAAO,CAAC,EAAE,CAAC;YACtB,MAAM,EAAE,GAAG,OAAO,CAAC,EAAE,CAAC;YACtB,MAAM,EAAE,GAAG,WAAW,CAAC,0BAA0B,CAAC,EAAE,CAAC,CAAC;YACtD,MAAM,EAAE,GAAG,WAAW,CAAC,0BAA0B,CAAC,EAAE,CAAC,CAAC;YACtD,mCAAmC;YACnC,uBAAuB,CAAC,0CAA0C,CAAC,WAAW,EAAE,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;YACnH,uBAAuB,CAAC,0CAA0C,CAAC,WAAW,EAAE,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;YACnH,OAAO,uBAAuB,CAAC,aAAa,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;SACxE;aAAM;YACL,IAAI,OAAO,CAAC,WAAW,EAAE,GAAG,GAAG,EAAE;gBAC/B,IAAI,WAAW,CAAC,KAAK,CAAC,aAAa,CAAC,OAAO,EAAE,IAAI,CAAC,SAAS,CAAC,EAAE;oBAC5D,WAAW,CAAC,eAAe,EAAE,CAAC;oBAC9B,WAAW,CAAC,eAAe,EAAE,CAAC;oBAC9B,OAAO,uBAAuB,CAAC,aAAa,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;iBACxE;gBACD,IAAI,WAAW,CAAC,MAAM,CAAC,aAAa,CAAC,OAAO,EAAE,IAAI,CAAC,SAAS,CAAC,EAAE;oBAC7D,WAAW,CAAC,aAAa,EAAE,CAAC;oBAC5B,WAAW,CAAC,aAAa,EAAE,CAAC;oBAC5B,OAAO,uBAAuB,CAAC,aAAa,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;iBACxE;aACF;iBAAM;gBACL,IAAI,WAAW,CAAC,KAAK,CAAC,aAAa,CAAC,OAAO,EAAE,IAAI,CAAC,SAAS,CAAC,EAAE;oBAC5D,WAAW,CAAC,eAAe,EAAE,CAAC;oBAC9B,WAAW,CAAC,aAAa,EAAE,CAAC;oBAC5B,OAAO,uBAAuB,CAAC,aAAa,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;iBACxE;gBACD,IAAI,WAAW,CAAC,MAAM,CAAC,aAAa,CAAC,OAAO,EAAE,IAAI,CAAC,SAAS,CAAC,EAAE;oBAC7D,WAAW,CAAC,aAAa,EAAE,CAAC;oBAC5B,WAAW,CAAC,eAAe,EAAE,CAAC;oBAC9B,OAAO,uBAAuB,CAAC,aAAa,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;iBACxE;aACF;SACF;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD;;;;;;;;OAQG;IACK,gBAAgB,CAAC,GAAmB,EAAE,GAAmB,EAAE,YAAuB,EAAE,UAAkB,EAAE,UAAkB,EAAE,OAAgB;QAClJ,MAAM,MAAM,GAAG,UAAU,GAAG,UAAU,CAAC;QACvC,MAAM,EAAE,GAAG,QAAQ,CAAC,yBAAyB,CAAC,YAAY,CAAC,EAAE,GAAG,UAAU,EAAE,MAAM,CAAC,CAAC;QACpF,MAAM,EAAE,GAAG,QAAQ,CAAC,yBAAyB,CAAC,YAAY,CAAC,EAAE,GAAG,UAAU,EAAE,MAAM,CAAC,CAAC;QACpF,IAAI,EAAE,KAAK,SAAS,IAAI,EAAE,KAAK,SAAS,EAAE;YACxC,MAAM,OAAO,GAAG,mBAAmB,CAAC,oCAAoC,CAAC,GAAG,EAAE,YAAY,CAAC,EAAE,EAAE,YAAY,CAAC,EAAE,CAAC,CAAC;YAChH,MAAM,OAAO,GAAG,mBAAmB,CAAC,oCAAoC,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;YACtF,IAAI,OAAO;gBACT,OAAO,CAAC,sBAAsB,EAAE,CAAC;YACnC,OAAO,uBAAuB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;SAChE;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACO,gBAAgB,CAAC,MAA6C,EAAE,IAAyC;QAC/G,IAAI,IAAI,KAAK,SAAS;YACpB,OAAO,MAAM,CAAC;QAChB,IAAI,MAAM,KAAK,SAAS;YACtB,OAAO,CAAC,IAAI,CAAC,CAAC;QAChB,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAClB,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;;;;OAMG;IACI,2BAA2B,CAAC,IAAW,EAAE,IAAW,EAAE,oBAA6B,IAAI;QAC5F,IAAI,MAA6C,CAAC;QAClD,IAAI,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,SAAS,CAAC,EAAE;YAC1D,MAAM,UAAU,GAAG,IAAI,CAAC,SAAS,CAAC,2BAA2B,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;YAC9E,IAAI,UAAU,EAAE;gBACd,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAAC,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAC/D,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAAC,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAC/D,MAAM,EAAE,GAAG,QAAQ,CAAC,YAAY,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;gBACzC,MAAM,EAAE,GAAG,QAAQ,CAAC,YAAY,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;gBACzC,MAAM,GAAG,GAAG,QAAQ,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBACpD,MAAM,KAAK,GAAG,QAAQ,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBACxD,IAAI,QAAQ,CAAC,mBAAmB,CAAC,EAAE,EAAE,GAAG,CAAC;uBACpC,QAAQ,CAAC,mBAAmB,CAAC,EAAE,EAAE,GAAG,CAAC;uBACrC,QAAQ,CAAC,mBAAmB,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE;oBACzC,MAAM,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAQ,iDAAiD;oBACnG,MAAM,cAAc,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAI,6EAA6E;oBAC/H,MAAM,gBAAgB,GAAG,cAAc,GAAG,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,CAAC,CAAG,qCAAqC;oBAC3H,MAAM,cAAc,GAAG,cAAc,GAAG,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,CAAC,CAAO,mCAAmC;oBACzH,MAAM,sBAAsB,GAAG,CAAC,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC;oBACxG,MAAM,MAAM,GAAG,UAAU,CAAC,qBAAqB,CAAC,gBAAgB,EAAE,cAAc,CAAC,CAAC;oBAClF,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC;oBAC1B,MAAM,eAAe,GAAG,MAAM,CAAC,cAAc,EAAE,CAAC;oBAChD,MAAM,UAAU,GAAG,MAAM,CAAC,iCAAiC,CAAC,MAAM,CAAC,YAAY,CAAC,CAAC,CAAG,oCAAoC;oBACxH,MAAM,CAAC,UAAU,IAAI,GAAG,CAAC,CAAC;oBAC1B,MAAM,aAAa,GAAG,MAAM,CAAC,YAAY,GAAG,MAAM,CAAC,YAAY,CAAC,CAAoB,oCAAoC;oBACxH,MAAM,UAAU,GAAG,UAAU,GAAG,aAAa,CAAC,CAAsC,kCAAkC;oBACtH,MAAM,cAAc,GAAG,SAAS,CAAC,MAAM,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;oBAEhE;;;;uBAIG;oBACH,MAAM,4BAA4B,GAAG,CAAC,uBAAkC,EAAE,eAAwB,EAAW,EAAE;wBAC7G,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;wBACxC,MAAM,SAAS,GAAG,uBAAuB,CAAC,EAAE,CAAC;wBAC7C,MAAM,OAAO,GAAG,uBAAuB,CAAC,EAAE,CAAC;wBAC3C,IAAI,uBAAuB,CAAC,iBAAiB,EAAE,IAAI,CAAC,QAAQ,CAAC,eAAe,CAAC,uBAAuB,CAAC,aAAa,EAAE,CAAC,EAAE;4BACrH,MAAM,GAAG,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,IAAI,CAAC,gBAAgB,CAAC,IAAI,EAAE,IAAI,EAAE,uBAAuB,EAAE,SAAS,EAAE,OAAO,EAAE,sBAAsB,CAAC,CAAC,CAAC;yBAChJ;6BAAM,EAAG,6BAA6B;4BACrC,MAAM,eAAe,GAAG,sBAAsB,CAAC,CAAC,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC,eAAe,CAAC;4BACpF,MAAM,MAAM,GAAG,IAAI,CAAC,OAAO,GAAG,eAAe,CAAC,CAAC,CAAC,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAA,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;4BAC3G,MAAM,MAAM,GAAG,IAAI,CAAC,OAAO,GAAG,eAAe,CAAC,CAAC,CAAC,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAA,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;4BAC3G,IAAI,MAAM,CAAC,aAAa,CAAC,MAAM,EAAE,IAAI,CAAC,SAAS,CAAC,EAAE;gCAChD,MAAM,OAAO,GAAG,mBAAmB,CAAC,wBAAwB,CAAC,IAAI,EAAE,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;gCACpG,MAAM,OAAO,GAAG,mBAAmB,CAAC,wBAAwB,CAAC,IAAI,EAAE,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;gCACpG,MAAM,GAAG,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,uBAAuB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;6BACjG;yBACF;wBACD,OAAO,MAAM,KAAK,SAAS,IAAI,MAAM,CAAC,MAAM,GAAG,IAAI,CAAC;oBACtD,CAAC,CAAC;oBAEF,4BAA4B,CAAC,cAAc,EAAE,KAAK,CAAC,CAAC,CAAE,wDAAwD;oBAE9G,yEAAyE;oBACzE,IAAI,UAAU,IAAI,eAAe;wBAC/B,4BAA4B,CAAC,SAAS,CAAC,MAAM,CAAC,UAAU,GAAG,eAAe,EAAE,UAAU,GAAG,eAAe,CAAC,EAAE,IAAI,CAAC,CAAC;yBAC9G,IAAI,UAAU,KAAK,GAAG;wBACzB,4BAA4B,CAAC,SAAS,CAAC,MAAM,CAAC,UAAU,GAAG,eAAe,EAAE,UAAU,GAAG,eAAe,CAAC,EAAE,IAAI,CAAC,CAAC;iBACpH;aACF;SACF;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;CACF","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n\r\n/** @packageDocumentation\r\n * @module CartesianGeometry\r\n */\r\nimport { assert } from \"@itwin/core-bentley\";\r\nimport { Arc3d } from \"../curve/Arc3d\";\r\nimport { CurveLocationDetail, CurveLocationDetailPair } from \"../curve/CurveLocationDetail\";\r\nimport { CurvePrimitive } from \"../curve/CurvePrimitive\";\r\nimport { Geometry } from \"../Geometry\";\r\nimport { AngleSweep } from \"./AngleSweep\";\r\nimport { Point3d, Vector3d } from \"./Point3dVector3d\";\r\nimport { Segment1d } from \"./Segment1d\";\r\n\r\n/**\r\n * `CoincidentGeometryQuery` has methods useful in testing for overlapping geometry.\r\n * * Each instance carries tolerance information that can be reused over extended call sequences.\r\n * * These methods are expected to be called internally by CurveCurve intersection methods.\r\n * @internal\r\n */\r\nexport class CoincidentGeometryQuery {\r\n  private _vectorU?: Vector3d;\r\n  private _vectorV?: Vector3d;\r\n  private _point0?: Point3d;\r\n  private _point1?: Point3d;\r\n  private _tolerance: number;\r\n  public get tolerance(): number {\r\n    return this._tolerance;\r\n  }\r\n  private constructor(tolerance: number = Geometry.smallMetricDistance) {\r\n    this._tolerance = tolerance;\r\n  }\r\n  public static create(tolerance: number = Geometry.smallMetricDistance): CoincidentGeometryQuery {\r\n    return new CoincidentGeometryQuery(tolerance);\r\n  }\r\n  /**\r\n   * * Assign both the fraction and fraction1 values in the detail, possibly swapped.\r\n   * * reevaluate the points as simple interpolation between given points.\r\n   */\r\n  public static assignDetailInterpolatedFractionsAndPoints(detail: CurveLocationDetail, f0: number, f1: number,\r\n    pointA: Point3d, pointB: Point3d, swap: boolean = false) {\r\n    if (swap) {\r\n      detail.fraction = f1;\r\n      detail.fraction1 = f0;\r\n    } else {\r\n      detail.fraction = f0;\r\n      detail.fraction1 = f1;\r\n    }\r\n    detail.point = pointA.interpolate(detail.fraction, pointB, detail.point);\r\n    detail.point1 = pointA.interpolate(detail.fraction1, pointB, detail.point1);\r\n\r\n  }\r\n\r\n  /** Return a curve location detail with projection of a `spacePoint` to the line segment with `pointA` and `pointB`\r\n   *\r\n   */\r\n  public projectPointToSegmentXY(spacePoint: Point3d, pointA: Point3d, pointB: Point3d): CurveLocationDetail {\r\n    this._vectorU = Vector3d.createStartEnd(pointA, pointB, this._vectorU);\r\n    this._vectorV = Vector3d.createStartEnd(pointA, spacePoint, this._vectorV);\r\n    const uDotU = this._vectorU.dotProductXY(this._vectorU);\r\n    const uDotV = this._vectorU.dotProductXY(this._vectorV);\r\n    const fraction = Geometry.safeDivideFraction(uDotV, uDotU, 0.0);\r\n    return CurveLocationDetail.createCurveFractionPoint(undefined, fraction,\r\n      pointA.interpolate(fraction, pointB));\r\n  }\r\n  /**\r\n   * * project `pointA0` and `pointA1` onto the segment with `pointB0` and `pointB1`\r\n   * * In the returned detail pair, the `detailA` has fractions along segmentA, and `detailB` has fractions along segment B\r\n   * @param pointA0 start point of segment A\r\n   * @param pointA1 end point of segment A\r\n   * @param pointB0 start point of segment B\r\n   * @param pointB1 end point of segment B\r\n   */\r\n  public coincidentSegmentRangeXY(pointA0: Point3d, pointA1: Point3d, pointB0: Point3d, pointB1: Point3d, restrictToBounds: boolean = true): CurveLocationDetailPair | undefined {\r\n    const detailA0OnB = this.projectPointToSegmentXY(pointA0, pointB0, pointB1);\r\n    if (pointA0.distanceXY(detailA0OnB.point) > this._tolerance)\r\n      return undefined;\r\n    const detailA1OnB = this.projectPointToSegmentXY(pointA1, pointB0, pointB1);\r\n    if (pointA1.distanceXY(detailA1OnB.point) > this._tolerance)\r\n      return undefined;\r\n\r\n    const detailB0OnA = this.projectPointToSegmentXY(pointB0, pointA0, pointA1);\r\n    if (pointB0.distanceXY(detailB0OnA.point) > this._tolerance)\r\n      return undefined;\r\n    const detailB1OnA = this.projectPointToSegmentXY(pointB1, pointA0, pointA1);\r\n    if (pointB1.distanceXY(detailB1OnA.point) > this._tolerance)\r\n      return undefined;\r\n\r\n    detailA0OnB.fraction1 = detailA1OnB.fraction;\r\n    detailA0OnB.point1 = detailA1OnB.point;  // capture -- detailA1OnB is not reused.\r\n    detailB0OnA.fraction1 = detailB1OnA.fraction;\r\n    detailB0OnA.point1 = detailB1OnA.point;\r\n    if (!restrictToBounds)\r\n      return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);\r\n\r\n    const segment = Segment1d.create(detailB0OnA.fraction, detailB0OnA.fraction1);\r\n    if (segment.clampDirectedTo01()) {\r\n      const f0 = segment.x0;\r\n      const f1 = segment.x1;\r\n      const h0 = detailB0OnA.inverseInterpolateFraction(f0);\r\n      const h1 = detailB0OnA.inverseInterpolateFraction(f1);\r\n      // recompute fractions and points..\r\n      CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(detailB0OnA, f0, f1, pointA0, pointA1, f0 > f1);\r\n      CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(detailA0OnB, h0, h1, pointB0, pointB1, h0 > h1);\r\n      return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);\r\n    } else {\r\n      if (segment.signedDelta() < 0.0) {\r\n        if (detailB0OnA.point.isAlmostEqual(pointA0, this.tolerance)) {\r\n          detailB0OnA.collapseToStart();\r\n          detailA0OnB.collapseToStart();\r\n          return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);\r\n        }\r\n        if (detailB0OnA.point1.isAlmostEqual(pointA1, this.tolerance)) {\r\n          detailB0OnA.collapseToEnd();\r\n          detailA0OnB.collapseToEnd();\r\n          return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);\r\n        }\r\n      } else {\r\n        if (detailB0OnA.point.isAlmostEqual(pointA1, this.tolerance)) {\r\n          detailB0OnA.collapseToStart();\r\n          detailA0OnB.collapseToEnd();\r\n          return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);\r\n        }\r\n        if (detailB0OnA.point1.isAlmostEqual(pointA0, this.tolerance)) {\r\n          detailB0OnA.collapseToEnd();\r\n          detailA0OnB.collapseToStart();\r\n          return CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);\r\n        }\r\n      }\r\n    }\r\n    return undefined;\r\n  }\r\n  /**\r\n   * Create a CurveLocationDetailPair for a coincident interval of two overlapping curves\r\n   * @param cpA curveA\r\n   * @param cpB curveB\r\n   * @param fractionsOnA coincident interval of curveB in fraction space of curveA\r\n   * @param fractionB0 curveB start in fraction space of curveA\r\n   * @param fractionB1 curveB end in fraction space of curveA\r\n   * @param reverse whether curveB and curveA have opposite direction\r\n   */\r\n  private createDetailPair(cpA: CurvePrimitive, cpB: CurvePrimitive, fractionsOnA: Segment1d, fractionB0: number, fractionB1: number, reverse: boolean): CurveLocationDetailPair | undefined {\r\n    const deltaB = fractionB1 - fractionB0;\r\n    const g0 = Geometry.conditionalDivideFraction(fractionsOnA.x0 - fractionB0, deltaB);\r\n    const g1 = Geometry.conditionalDivideFraction(fractionsOnA.x1 - fractionB0, deltaB);\r\n    if (g0 !== undefined && g1 !== undefined) {\r\n      const detailA = CurveLocationDetail.createCurveEvaluatedFractionFraction(cpA, fractionsOnA.x0, fractionsOnA.x1);\r\n      const detailB = CurveLocationDetail.createCurveEvaluatedFractionFraction(cpB, g0, g1);\r\n      if (reverse)\r\n        detailA.swapFractionsAndPoints();\r\n      return CurveLocationDetailPair.createCapture(detailA, detailB);\r\n    }\r\n    return undefined;\r\n  }\r\n  private appendDetailPair(result: CurveLocationDetailPair[] | undefined, pair: CurveLocationDetailPair | undefined): CurveLocationDetailPair[] | undefined {\r\n    if (pair === undefined)\r\n      return result;\r\n    if (result === undefined)\r\n      return [pair];\r\n    result.push(pair);\r\n    return result;\r\n  }\r\n  /**\r\n   * Test if 2 arcs have coinciding portions.\r\n   * @param arcA\r\n   * @param arcB\r\n   * @param _restrictToBounds\r\n   * @return 0, 1, or 2 overlap points/intervals\r\n   */\r\n  public coincidentArcIntersectionXY(arcA: Arc3d, arcB: Arc3d, _restrictToBounds: boolean = true): CurveLocationDetailPair[] | undefined {\r\n    let result: CurveLocationDetailPair[] | undefined;\r\n    if (arcA.center.isAlmostEqual(arcB.center, this.tolerance)) {\r\n      const matrixBToA = arcA.matrixRef.multiplyMatrixInverseMatrix(arcB.matrixRef);\r\n      if (matrixBToA) {\r\n        const ux = matrixBToA.at(0, 0); const uy = matrixBToA.at(1, 0);\r\n        const vx = matrixBToA.at(0, 1); const vy = matrixBToA.at(1, 1);\r\n        const ru = Geometry.hypotenuseXY(ux, uy);\r\n        const rv = Geometry.hypotenuseXY(vx, vy);\r\n        const dot = Geometry.dotProductXYXY(ux, uy, vx, vy);\r\n        const cross = Geometry.crossProductXYXY(ux, uy, vx, vy);\r\n        if (Geometry.isAlmostEqualNumber(ru, 1.0)\r\n          && Geometry.isAlmostEqualNumber(rv, 1.0)\r\n          && Geometry.isAlmostEqualNumber(0, dot)) {\r\n          const alphaB0Radians = Math.atan2(uy, ux);        // angular position of arcB 0 point in arcA sweep\r\n          const sweepDirection = cross > 0 ? 1.0 : -1.0;    // 1 if arcB parameter space sweeps in same direction as arcA, -1 if opposite\r\n          const betaStartRadians = alphaB0Radians + sweepDirection * arcB.sweep.startRadians;   // arcB start in arcA parameter space\r\n          const betaEndRadians = alphaB0Radians + sweepDirection * arcB.sweep.endRadians;       // arcB end in arcA parameter space\r\n          const fractionSpacesReversed = (sweepDirection * arcA.sweep.sweepRadians * arcB.sweep.sweepRadians) < 0;\r\n          const sweepB = AngleSweep.createStartEndRadians(betaStartRadians, betaEndRadians);\r\n          const sweepA = arcA.sweep;\r\n          const fractionPeriodA = sweepA.fractionPeriod();\r\n          const fractionB0 = sweepA.radiansToPositivePeriodicFraction(sweepB.startRadians);   // arcB start in arcA fraction space\r\n          assert(fractionB0 >= 0.0);\r\n          const fractionSweep = sweepB.sweepRadians / sweepA.sweepRadians;                    // arcB sweep in arcA fraction space\r\n          const fractionB1 = fractionB0 + fractionSweep;                                      // arcB end in arcA fraction space\r\n          const fractionSweepB = Segment1d.create(fractionB0, fractionB1);\r\n\r\n          /** lambda to add a coincident interval or isolated intersection, given inputs in arcA fraction space\r\n           * @param arcBInArcAFractionSpace span of arcB in arcA fraction space. On return, clamped to [0,1] if nontrivial.\r\n           * @param testStartOfArcA if no nontrivial coincident interval was found, look for an isolated intersection at the start (true) or end (false) of arcA\r\n           * @returns whether a detail pair was appended to result\r\n           */\r\n          const appendCoincidentIntersection = (arcBInArcAFractionSpace: Segment1d, testStartOfArcA: boolean): boolean => {\r\n            const size = result ? result.length : 0;\r\n            const arcBStart = arcBInArcAFractionSpace.x0;\r\n            const arcBEnd = arcBInArcAFractionSpace.x1;\r\n            if (arcBInArcAFractionSpace.clampDirectedTo01() && !Geometry.isSmallRelative(arcBInArcAFractionSpace.absoluteDelta())) {\r\n              result = this.appendDetailPair(result, this.createDetailPair(arcA, arcB, arcBInArcAFractionSpace, arcBStart, arcBEnd, fractionSpacesReversed));\r\n            } else {  // test isolated intersection\r\n              const testStartOfArcB = fractionSpacesReversed ? testStartOfArcA : !testStartOfArcA;\r\n              const arcAPt = this._point0 = testStartOfArcA ? arcA.startPoint(this._point0): arcA.endPoint(this._point0);\r\n              const arcBPt = this._point1 = testStartOfArcB ? arcB.startPoint(this._point1): arcB.endPoint(this._point1);\r\n              if (arcAPt.isAlmostEqual(arcBPt, this.tolerance)) {\r\n                const detailA = CurveLocationDetail.createCurveFractionPoint(arcA, testStartOfArcA ? 0 : 1, arcAPt);\r\n                const detailB = CurveLocationDetail.createCurveFractionPoint(arcB, testStartOfArcB ? 0 : 1, arcBPt);\r\n                result = this.appendDetailPair(result, CurveLocationDetailPair.createCapture(detailA, detailB));\r\n              }\r\n            }\r\n            return result !== undefined && result.length > size;\r\n          };\r\n\r\n          appendCoincidentIntersection(fractionSweepB, false);  // compute overlap in strict interior, or at end of arcA\r\n\r\n          // check overlap at start of arcA with a periodic shift of fractionSweepB\r\n          if (fractionB1 >= fractionPeriodA)\r\n            appendCoincidentIntersection(Segment1d.create(fractionB0 - fractionPeriodA, fractionB1 - fractionPeriodA), true);\r\n          else if (fractionB0 === 0.0)\r\n            appendCoincidentIntersection(Segment1d.create(fractionB0 + fractionPeriodA, fractionB1 + fractionPeriodA), true);\r\n        }\r\n      }\r\n    }\r\n    return result;\r\n  }\r\n}\r\n"]}
+{"version":3,"file":"CoincidentGeometryOps.js","sourceRoot":"","sources":["../../../src/geometry3d/CoincidentGeometryOps.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AACH,OAAO,EAAE,MAAM,EAAE,MAAM,qBAAqB,CAAC;AAE7C,OAAO,EAAE,mBAAmB,EAAE,uBAAuB,EAAE,MAAM,8BAA8B,CAAC;AAE5F,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AACvC,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAC1C,OAAO,EAAW,QAAQ,EAAE,MAAM,mBAAmB,CAAC;AACtD,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AAExC;;;;;GAKG;AACH,MAAM,OAAO,uBAAuB;IAMlC,IAAW,SAAS;QAClB,OAAO,IAAI,CAAC,UAAU,CAAC;IACzB,CAAC;IACD,YAAoB,YAAoB,QAAQ,CAAC,mBAAmB;QAClE,IAAI,CAAC,UAAU,GAAG,SAAS,CAAC;IAC9B,CAAC;IACM,MAAM,CAAC,MAAM,CAAC,YAAoB,QAAQ,CAAC,mBAAmB;QACnE,OAAO,IAAI,uBAAuB,CAAC,SAAS,CAAC,CAAC;IAChD,CAAC;IACD;;;OAGG;IACI,MAAM,CAAC,0CAA0C,CAAC,MAA2B,EAAE,EAAU,EAAE,EAAU,EAC1G,MAAe,EAAE,MAAe,EAAE,OAAgB,KAAK;QACvD,IAAI,IAAI,EAAE;YACR,MAAM,CAAC,QAAQ,GAAG,EAAE,CAAC;YACrB,MAAM,CAAC,SAAS,GAAG,EAAE,CAAC;SACvB;aAAM;YACL,MAAM,CAAC,QAAQ,GAAG,EAAE,CAAC;YACrB,MAAM,CAAC,SAAS,GAAG,EAAE,CAAC;SACvB;QACD,MAAM,CAAC,KAAK,GAAG,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,QAAQ,EAAE,MAAM,EAAE,MAAM,CAAC,KAAK,CAAC,CAAC;QACzE,MAAM,CAAC,MAAM,GAAG,MAAM,CAAC,WAAW,CAAC,MAAM,CAAC,SAAS,EAAE,MAAM,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;IAE9E,CAAC;IAED,sHAAsH;IAC/G,uBAAuB,CAAC,UAAmB,EAAE,MAAe,EAAE,MAAe;QAClF,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC,cAAc,CAAC,MAAM,EAAE,MAAM,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QACvE,IAAI,CAAC,QAAQ,GAAG,QAAQ,CAAC,cAAc,CAAC,MAAM,EAAE,UAAU,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QAC3E,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,YAAY,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QACxD,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,YAAY,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QACxD,MAAM,QAAQ,GAAG,QAAQ,CAAC,kBAAkB,CAAC,KAAK,EAAE,KAAK,EAAE,GAAG,CAAC,CAAC;QAChE,OAAO,mBAAmB,CAAC,wBAAwB,CAAC,SAAS,EAAE,QAAQ,EACrE,MAAM,CAAC,WAAW,CAAC,QAAQ,EAAE,MAAM,CAAC,CAAC,CAAC;IAC1C,CAAC;IACD;;;;;;;;;;;;;OAaG;IACI,qCAAqC,CAAC,OAAgC,EAC3E,OAAgB,EAAE,OAAgB,EAAE,OAAgB,EAAE,OAAgB,EACtE,WAAoB,KAAK,EAAE,WAAoB,KAAK,EAAE,WAAoB,KAAK,EAAE,WAAoB,KAAK;QAE1G,MAAM,MAAM,GAAG,SAAS,CAAC,MAAM,CAAC,OAAO,CAAC,OAAO,CAAC,QAAQ,EAAE,OAAO,CAAC,OAAO,CAAC,YAAY,CAAC,CAAC,CAAC,OAAO,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,OAAO,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC;QAC/I,MAAM,MAAM,GAAG,SAAS,CAAC,MAAM,CAAC,OAAO,CAAC,OAAO,CAAC,QAAQ,EAAE,OAAO,CAAC,OAAO,CAAC,YAAY,CAAC,CAAC,CAAC,OAAO,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,OAAO,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC;QAC/I,MAAM,QAAQ,GAAG,MAAM,CAAC,WAAW,EAAE,GAAG,GAAG,CAAC;QAE5C,MAAM,6CAA6C,GAAG,GAA4B,EAAE;YAClF,MAAM,EAAE,GAAG,MAAM,CAAC,EAAE,CAAC;YACrB,MAAM,EAAE,GAAG,MAAM,CAAC,EAAE,CAAC;YACrB,MAAM,EAAE,GAAG,MAAM,CAAC,EAAE,CAAC;YACrB,MAAM,EAAE,GAAG,MAAM,CAAC,EAAE,CAAC;YACrB,uBAAuB,CAAC,0CAA0C,CAAC,OAAO,CAAC,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;YACvH,uBAAuB,CAAC,0CAA0C,CAAC,OAAO,CAAC,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,OAAO,EAAE,OAAO,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;YACvH,OAAO,OAAO,CAAC;QACjB,CAAC,CAAC;QAEF,MAAM,aAAa,GAAG,MAAM,CAAC,iBAAiB,CAAC,CAAC,QAAQ,EAAE,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;QAC5E,MAAM,aAAa,GAAG,MAAM,CAAC,iBAAiB,CAAC,CAAC,QAAQ,EAAE,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC;QAC5E,IAAI,aAAa,IAAI,aAAa,EAAE;YAClC,IAAI,QAAQ,CAAC,mBAAmB,CAAC,MAAM,CAAC,aAAa,EAAE,EAAE,MAAM,CAAC,aAAa,EAAE,EAAE,QAAQ,CAAC,aAAa,CAAC;gBACtG,OAAO,6CAA6C,EAAE,CAAC,CAAE,2CAA2C;iBACjG,IAAI,MAAM,CAAC,iBAAiB,CAAC,IAAI,EAAE,IAAI,EAAE,KAAK,CAAC,IAAI,MAAM,CAAC,iBAAiB,CAAC,IAAI,EAAE,IAAI,EAAE,KAAK,CAAC;gBACjG,OAAO,6CAA6C,EAAE,CAAC,CAAE,uCAAuC;SACnG;QAED,MAAM,mBAAmB,GAAG,CAAC,MAAe,EAAE,MAAe,EAAE,QAAiB,EAAE,QAAiB,EAA2B,EAAE;YAC9H,MAAM,CAAC,KAAK,CAAC,OAAO,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;YACpC,MAAM,CAAC,KAAK,CAAC,OAAO,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;YACpC,OAAO,CAAC,OAAO,CAAC,QAAQ,GAAG,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;YAChD,OAAO,CAAC,OAAO,CAAC,QAAQ,GAAG,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;YAChD,OAAO,CAAC,OAAO,CAAC,eAAe,EAAE,CAAC;YAClC,OAAO,CAAC,OAAO,CAAC,eAAe,EAAE,CAAC;YAClC,OAAO,OAAO,CAAC;QACjB,CAAC,CAAC;QAEF,MAAM,cAAc,GAAG,MAAM,CAAC,iBAAiB,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAClE,MAAM,cAAc,GAAG,MAAM,CAAC,iBAAiB,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAClE,IAAI,cAAc,IAAI,cAAc,EAAE,EAAE,iCAAiC;YACvE,MAAM,MAAM,GAAG,OAAO,CAAC,OAAO,CAAC,MAAM,IAAI,OAAO,CAAC,OAAO,CAAC,KAAK,CAAC;YAC/D,IAAI,QAAQ,EAAE;gBACZ,IAAI,OAAO,CAAC,OAAO,CAAC,KAAK,CAAC,aAAa,CAAC,OAAO,EAAE,IAAI,CAAC,SAAS,CAAC;oBAC9D,OAAO,mBAAmB,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;qBACtD,IAAI,MAAM,CAAC,aAAa,CAAC,OAAO,EAAE,IAAI,CAAC,SAAS,CAAC;oBACpD,OAAO,mBAAmB,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC;aAC9D;iBAAM;gBACL,IAAI,MAAM,CAAC,aAAa,CAAC,OAAO,EAAE,IAAI,CAAC,SAAS,CAAC;oBAC/C,OAAO,mBAAmB,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC;qBACvD,IAAI,OAAO,CAAC,OAAO,CAAC,KAAK,CAAC,aAAa,CAAC,OAAO,EAAE,IAAI,CAAC,SAAS,CAAC;oBACnE,OAAO,mBAAmB,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,EAAE,IAAI,CAAC,CAAC;aAC7D;SACF;QACD,OAAO,SAAS,CAAC,CAAC,kBAAkB;IACtC,CAAC;IACD;;;;;;;;;;OAUG;IACI,wBAAwB,CAAC,OAAgB,EAAE,OAAgB,EAAE,OAAgB,EAAE,OAAgB,EAAE,mBAA4B,IAAI;QACtI,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QACnB,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QAEnB,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QACnB,MAAM,WAAW,GAAG,IAAI,CAAC,uBAAuB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;QAC5E,IAAI,OAAO,CAAC,UAAU,CAAC,WAAW,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,UAAU;YACzD,OAAO,SAAS,CAAC;QAEnB,WAAW,CAAC,SAAS,GAAG,WAAW,CAAC,QAAQ,CAAC;QAC7C,WAAW,CAAC,MAAM,GAAG,WAAW,CAAC,KAAK,CAAC,CAAE,wCAAwC;QACjF,WAAW,CAAC,SAAS,GAAG,WAAW,CAAC,QAAQ,CAAC;QAC7C,WAAW,CAAC,MAAM,GAAG,WAAW,CAAC,KAAK,CAAC;QACvC,MAAM,QAAQ,GAAG,uBAAuB,CAAC,aAAa,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;QAEjF,OAAO,gBAAgB,CAAC,CAAC,CAAC,IAAI,CAAC,qCAAqC,CAAC,QAAQ,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC;IAChI,CAAC;IACD;;;;;;;;OAQG;IACK,gBAAgB,CAAC,GAAmB,EAAE,GAAmB,EAAE,YAAuB,EAAE,UAAkB,EAAE,UAAkB,EAAE,OAAgB;QAClJ,MAAM,MAAM,GAAG,UAAU,GAAG,UAAU,CAAC;QACvC,MAAM,EAAE,GAAG,QAAQ,CAAC,yBAAyB,CAAC,YAAY,CAAC,EAAE,GAAG,UAAU,EAAE,MAAM,CAAC,CAAC;QACpF,MAAM,EAAE,GAAG,QAAQ,CAAC,yBAAyB,CAAC,YAAY,CAAC,EAAE,GAAG,UAAU,EAAE,MAAM,CAAC,CAAC;QACpF,IAAI,EAAE,KAAK,SAAS,IAAI,EAAE,KAAK,SAAS,EAAE;YACxC,MAAM,OAAO,GAAG,mBAAmB,CAAC,oCAAoC,CAAC,GAAG,EAAE,YAAY,CAAC,EAAE,EAAE,YAAY,CAAC,EAAE,CAAC,CAAC;YAChH,MAAM,OAAO,GAAG,mBAAmB,CAAC,oCAAoC,CAAC,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;YACtF,IAAI,OAAO;gBACT,OAAO,CAAC,sBAAsB,EAAE,CAAC;YACnC,OAAO,uBAAuB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;SAChE;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACO,gBAAgB,CAAC,MAA6C,EAAE,IAAyC;QAC/G,IAAI,IAAI,KAAK,SAAS;YACpB,OAAO,MAAM,CAAC;QAChB,IAAI,MAAM,KAAK,SAAS;YACtB,OAAO,CAAC,IAAI,CAAC,CAAC;QAChB,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAClB,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;;;;OAMG;IACI,2BAA2B,CAAC,IAAW,EAAE,IAAW,EAAE,oBAA6B,IAAI;QAC5F,IAAI,MAA6C,CAAC;QAClD,IAAI,IAAI,CAAC,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,SAAS,CAAC,EAAE;YAC1D,MAAM,UAAU,GAAG,IAAI,CAAC,SAAS,CAAC,2BAA2B,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;YAC9E,IAAI,UAAU,EAAE;gBACd,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAAC,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAC/D,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAAC,MAAM,EAAE,GAAG,UAAU,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBAC/D,MAAM,EAAE,GAAG,QAAQ,CAAC,YAAY,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;gBACzC,MAAM,EAAE,GAAG,QAAQ,CAAC,YAAY,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;gBACzC,MAAM,GAAG,GAAG,QAAQ,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBACpD,MAAM,KAAK,GAAG,QAAQ,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBACxD,IAAI,QAAQ,CAAC,mBAAmB,CAAC,EAAE,EAAE,GAAG,CAAC;uBACpC,QAAQ,CAAC,mBAAmB,CAAC,EAAE,EAAE,GAAG,CAAC;uBACrC,QAAQ,CAAC,mBAAmB,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE;oBACzC,MAAM,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAQ,iDAAiD;oBACnG,MAAM,cAAc,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAI,6EAA6E;oBAC/H,MAAM,gBAAgB,GAAG,cAAc,GAAG,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,CAAC,CAAG,qCAAqC;oBAC3H,MAAM,cAAc,GAAG,cAAc,GAAG,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,CAAC,CAAO,mCAAmC;oBACzH,MAAM,sBAAsB,GAAG,CAAC,cAAc,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,GAAG,IAAI,CAAC,KAAK,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC;oBACxG,MAAM,MAAM,GAAG,UAAU,CAAC,qBAAqB,CAAC,gBAAgB,EAAE,cAAc,CAAC,CAAC;oBAClF,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC;oBAC1B,MAAM,eAAe,GAAG,MAAM,CAAC,cAAc,EAAE,CAAC;oBAChD,MAAM,UAAU,GAAG,MAAM,CAAC,iCAAiC,CAAC,MAAM,CAAC,YAAY,CAAC,CAAC,CAAG,oCAAoC;oBACxH,MAAM,CAAC,UAAU,IAAI,GAAG,CAAC,CAAC;oBAC1B,MAAM,aAAa,GAAG,MAAM,CAAC,YAAY,GAAG,MAAM,CAAC,YAAY,CAAC,CAAoB,oCAAoC;oBACxH,MAAM,UAAU,GAAG,UAAU,GAAG,aAAa,CAAC,CAAsC,kCAAkC;oBACtH,MAAM,cAAc,GAAG,SAAS,CAAC,MAAM,CAAC,UAAU,EAAE,UAAU,CAAC,CAAC;oBAEhE;;;;uBAIG;oBACH,MAAM,4BAA4B,GAAG,CAAC,uBAAkC,EAAE,eAAwB,EAAW,EAAE;wBAC7G,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;wBACxC,MAAM,SAAS,GAAG,uBAAuB,CAAC,EAAE,CAAC;wBAC7C,MAAM,OAAO,GAAG,uBAAuB,CAAC,EAAE,CAAC;wBAC3C,IAAI,uBAAuB,CAAC,iBAAiB,EAAE,IAAI,CAAC,QAAQ,CAAC,eAAe,CAAC,uBAAuB,CAAC,aAAa,EAAE,CAAC,EAAE;4BACrH,MAAM,GAAG,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,IAAI,CAAC,gBAAgB,CAAC,IAAI,EAAE,IAAI,EAAE,uBAAuB,EAAE,SAAS,EAAE,OAAO,EAAE,sBAAsB,CAAC,CAAC,CAAC;yBAChJ;6BAAM,EAAG,6BAA6B;4BACrC,MAAM,eAAe,GAAG,sBAAsB,CAAC,CAAC,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC,eAAe,CAAC;4BACpF,MAAM,MAAM,GAAG,IAAI,CAAC,OAAO,GAAG,eAAe,CAAC,CAAC,CAAC,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;4BAC5G,MAAM,MAAM,GAAG,IAAI,CAAC,OAAO,GAAG,eAAe,CAAC,CAAC,CAAC,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;4BAC5G,IAAI,MAAM,CAAC,aAAa,CAAC,MAAM,EAAE,IAAI,CAAC,SAAS,CAAC,EAAE;gCAChD,MAAM,OAAO,GAAG,mBAAmB,CAAC,wBAAwB,CAAC,IAAI,EAAE,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;gCACpG,MAAM,OAAO,GAAG,mBAAmB,CAAC,wBAAwB,CAAC,IAAI,EAAE,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;gCACpG,MAAM,GAAG,IAAI,CAAC,gBAAgB,CAAC,MAAM,EAAE,uBAAuB,CAAC,aAAa,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,CAAC;6BACjG;yBACF;wBACD,OAAO,MAAM,KAAK,SAAS,IAAI,MAAM,CAAC,MAAM,GAAG,IAAI,CAAC;oBACtD,CAAC,CAAC;oBAEF,4BAA4B,CAAC,cAAc,EAAE,KAAK,CAAC,CAAC,CAAE,wDAAwD;oBAE9G,yEAAyE;oBACzE,IAAI,UAAU,IAAI,eAAe;wBAC/B,4BAA4B,CAAC,SAAS,CAAC,MAAM,CAAC,UAAU,GAAG,eAAe,EAAE,UAAU,GAAG,eAAe,CAAC,EAAE,IAAI,CAAC,CAAC;yBAC9G,IAAI,UAAU,KAAK,GAAG;wBACzB,4BAA4B,CAAC,SAAS,CAAC,MAAM,CAAC,UAAU,GAAG,eAAe,EAAE,UAAU,GAAG,eAAe,CAAC,EAAE,IAAI,CAAC,CAAC;iBACpH;aACF;SACF;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;CACF","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n\r\n/** @packageDocumentation\r\n * @module CartesianGeometry\r\n */\r\nimport { assert } from \"@itwin/core-bentley\";\r\nimport { Arc3d } from \"../curve/Arc3d\";\r\nimport { CurveLocationDetail, CurveLocationDetailPair } from \"../curve/CurveLocationDetail\";\r\nimport { CurvePrimitive } from \"../curve/CurvePrimitive\";\r\nimport { Geometry } from \"../Geometry\";\r\nimport { AngleSweep } from \"./AngleSweep\";\r\nimport { Point3d, Vector3d } from \"./Point3dVector3d\";\r\nimport { Segment1d } from \"./Segment1d\";\r\n\r\n/**\r\n * `CoincidentGeometryQuery` has methods useful in testing for overlapping geometry.\r\n * * Each instance carries tolerance information that can be reused over extended call sequences.\r\n * * These methods are expected to be called internally by CurveCurve intersection methods.\r\n * @internal\r\n */\r\nexport class CoincidentGeometryQuery {\r\n  private _vectorU?: Vector3d;\r\n  private _vectorV?: Vector3d;\r\n  private _point0?: Point3d;\r\n  private _point1?: Point3d;\r\n  private _tolerance: number;\r\n  public get tolerance(): number {\r\n    return this._tolerance;\r\n  }\r\n  private constructor(tolerance: number = Geometry.smallMetricDistance) {\r\n    this._tolerance = tolerance;\r\n  }\r\n  public static create(tolerance: number = Geometry.smallMetricDistance): CoincidentGeometryQuery {\r\n    return new CoincidentGeometryQuery(tolerance);\r\n  }\r\n  /**\r\n   * * Assign both the fraction and fraction1 values in the detail, possibly swapped.\r\n   * * reevaluate the points as simple interpolation between given points.\r\n   */\r\n  public static assignDetailInterpolatedFractionsAndPoints(detail: CurveLocationDetail, f0: number, f1: number,\r\n    pointA: Point3d, pointB: Point3d, swap: boolean = false) {\r\n    if (swap) {\r\n      detail.fraction = f1;\r\n      detail.fraction1 = f0;\r\n    } else {\r\n      detail.fraction = f0;\r\n      detail.fraction1 = f1;\r\n    }\r\n    detail.point = pointA.interpolate(detail.fraction, pointB, detail.point);\r\n    detail.point1 = pointA.interpolate(detail.fraction1, pointB, detail.point1);\r\n\r\n  }\r\n\r\n  /** Return a curve location detail with projection of a `spacePoint` to the line segment with `pointA` and `pointB` */\r\n  public projectPointToSegmentXY(spacePoint: Point3d, pointA: Point3d, pointB: Point3d): CurveLocationDetail {\r\n    this._vectorU = Vector3d.createStartEnd(pointA, pointB, this._vectorU);\r\n    this._vectorV = Vector3d.createStartEnd(pointA, spacePoint, this._vectorV);\r\n    const uDotU = this._vectorU.dotProductXY(this._vectorU);\r\n    const uDotV = this._vectorU.dotProductXY(this._vectorV);\r\n    const fraction = Geometry.safeDivideFraction(uDotV, uDotU, 0.0);\r\n    return CurveLocationDetail.createCurveFractionPoint(undefined, fraction,\r\n      pointA.interpolate(fraction, pointB));\r\n  }\r\n  /**\r\n   * Given a detail pair representing the projection of each of two colinear line segments onto the other,\r\n   * clamp the details (in place) to the line segments' endpoints according to the given flags.\r\n   * @param overlap segment overlap as returned by [[coincidentSegmentRangeXY]], modified on return\r\n   * @param pointA0 start point of segment A\r\n   * @param pointA1 end point of segment A\r\n   * @param pointB0 start point of segment B\r\n   * @param pointB1 end point of segment B\r\n   * @param extendA0 whether to extend segment A beyond its start\r\n   * @param extendA1 whether to extend segment A beyond its end\r\n   * @param extendB0 whether to extend segment B beyond its start\r\n   * @param extendB1 whether to extend segment B beyond its end\r\n   * @return reference to the input clamped in place, or undefined (leaving interval untouched) if clamping would result in empty interval.\r\n   */\r\n  public clampCoincidentOverlapToSegmentBounds(overlap: CurveLocationDetailPair,\r\n    pointA0: Point3d, pointA1: Point3d, pointB0: Point3d, pointB1: Point3d,\r\n    extendA0: boolean = false, extendA1: boolean = false, extendB0: boolean = false, extendB1: boolean = false,\r\n  ): CurveLocationDetailPair | undefined {\r\n    const rangeA = Segment1d.create(overlap.detailA.fraction, overlap.detailA.hasFraction1 ? overlap.detailA.fraction1 : overlap.detailA.fraction);\r\n    const rangeB = Segment1d.create(overlap.detailB.fraction, overlap.detailB.hasFraction1 ? overlap.detailB.fraction1 : overlap.detailB.fraction);\r\n    const reversed = rangeA.signedDelta() < 0.0;\r\n\r\n    const updateIntervalFromRangesAndInterpolatedPoints = (): CurveLocationDetailPair => {\r\n      const a0 = rangeA.x0;\r\n      const a1 = rangeA.x1;\r\n      const b0 = rangeB.x0;\r\n      const b1 = rangeB.x1;\r\n      CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(overlap.detailA, a0, a1, pointA0, pointA1, a0 > a1);\r\n      CoincidentGeometryQuery.assignDetailInterpolatedFractionsAndPoints(overlap.detailB, b0, b1, pointB0, pointB1, b0 > b1);\r\n      return overlap;\r\n    };\r\n\r\n    const haveIntervalA = rangeA.clampDirectedTo01(!extendA0, !extendA1, false);\r\n    const haveIntervalB = rangeB.clampDirectedTo01(!extendB0, !extendB1, false);\r\n    if (haveIntervalA && haveIntervalB) {\r\n      if (Geometry.isAlmostEqualNumber(rangeA.absoluteDelta(), rangeB.absoluteDelta(), Geometry.smallFraction))\r\n        return updateIntervalFromRangesAndInterpolatedPoints();  // intersection of partially clamped ranges\r\n      else if (rangeA.clampDirectedTo01(true, true, false) && rangeB.clampDirectedTo01(true, true, false))\r\n        return updateIntervalFromRangesAndInterpolatedPoints();  // intersection of fully clamped ranges\r\n    }\r\n\r\n    const collapseToSingleton = (pointA: Point3d, pointB: Point3d, atStartA: boolean, atStartB: boolean): CurveLocationDetailPair => {\r\n      pointA.clone(overlap.detailA.point);\r\n      pointB.clone(overlap.detailB.point);\r\n      overlap.detailA.fraction = atStartA ? 0.0 : 1.0;\r\n      overlap.detailB.fraction = atStartB ? 0.0 : 1.0;\r\n      overlap.detailA.collapseToStart();\r\n      overlap.detailB.collapseToStart();\r\n      return overlap;\r\n    };\r\n\r\n    const haveSingletonA = rangeA.clampDirectedTo01(true, true, true);\r\n    const haveSingletonB = rangeB.clampDirectedTo01(true, true, true);\r\n    if (haveSingletonA && haveSingletonB) { // intersection is a single point\r\n      const point1 = overlap.detailA.point1 ?? overlap.detailA.point;\r\n      if (reversed) {\r\n        if (overlap.detailA.point.isAlmostEqual(pointA0, this.tolerance))\r\n          return collapseToSingleton(pointA0, pointB0, true, true);\r\n        else if (point1.isAlmostEqual(pointA1, this.tolerance))\r\n          return collapseToSingleton(pointA1, pointB1, false, false);\r\n      } else {\r\n        if (point1.isAlmostEqual(pointA0, this.tolerance))\r\n          return collapseToSingleton(pointA0, pointB1, true, false);\r\n        else if (overlap.detailA.point.isAlmostEqual(pointA1, this.tolerance))\r\n          return collapseToSingleton(pointA1, pointB0, false, true);\r\n      }\r\n    }\r\n    return undefined; // no intersection\r\n  }\r\n  /**\r\n   * Compute whether two line segments have a coincident overlap in xy.\r\n   * * Project `pointA0` and `pointA1` onto the line formed by `pointB0` and `pointB1` and vice versa\r\n   * * If all projection distances are sufficiently small, return a detail pair recording the coincident interval, optionally clipped to segment bounds.\r\n   * @param pointA0 start point of segment A\r\n   * @param pointA1 end point of segment A\r\n   * @param pointB0 start point of segment B\r\n   * @param pointB1 end point of segment B\r\n   * @param restrictToBounds whether to clip the coincident segment details to the segment bounds\r\n   * @return detail pair for the coincident interval (`detailA` has fractions along segment A, and `detailB` has fractions along segment B), or undefined if no coincidence\r\n   */\r\n  public coincidentSegmentRangeXY(pointA0: Point3d, pointA1: Point3d, pointB0: Point3d, pointB1: Point3d, restrictToBounds: boolean = true): CurveLocationDetailPair | undefined {\r\n    const detailA0OnB = this.projectPointToSegmentXY(pointA0, pointB0, pointB1);\r\n    if (pointA0.distanceXY(detailA0OnB.point) > this._tolerance)\r\n      return undefined;\r\n    const detailA1OnB = this.projectPointToSegmentXY(pointA1, pointB0, pointB1);\r\n    if (pointA1.distanceXY(detailA1OnB.point) > this._tolerance)\r\n      return undefined;\r\n\r\n    const detailB0OnA = this.projectPointToSegmentXY(pointB0, pointA0, pointA1);\r\n    if (pointB0.distanceXY(detailB0OnA.point) > this._tolerance)\r\n      return undefined;\r\n    const detailB1OnA = this.projectPointToSegmentXY(pointB1, pointA0, pointA1);\r\n    if (pointB1.distanceXY(detailB1OnA.point) > this._tolerance)\r\n      return undefined;\r\n\r\n    detailA0OnB.fraction1 = detailA1OnB.fraction;\r\n    detailA0OnB.point1 = detailA1OnB.point;  // capture -- detailA1OnB is not reused.\r\n    detailB0OnA.fraction1 = detailB1OnA.fraction;\r\n    detailB0OnA.point1 = detailB1OnA.point;\r\n    const interval = CurveLocationDetailPair.createCapture(detailB0OnA, detailA0OnB);\r\n\r\n    return restrictToBounds ? this.clampCoincidentOverlapToSegmentBounds(interval, pointA0, pointA1, pointB0, pointB1) : interval;\r\n  }\r\n  /**\r\n   * Create a CurveLocationDetailPair for a coincident interval of two overlapping curves\r\n   * @param cpA curveA\r\n   * @param cpB curveB\r\n   * @param fractionsOnA coincident interval of curveB in fraction space of curveA\r\n   * @param fractionB0 curveB start in fraction space of curveA\r\n   * @param fractionB1 curveB end in fraction space of curveA\r\n   * @param reverse whether curveB and curveA have opposite direction\r\n   */\r\n  private createDetailPair(cpA: CurvePrimitive, cpB: CurvePrimitive, fractionsOnA: Segment1d, fractionB0: number, fractionB1: number, reverse: boolean): CurveLocationDetailPair | undefined {\r\n    const deltaB = fractionB1 - fractionB0;\r\n    const g0 = Geometry.conditionalDivideFraction(fractionsOnA.x0 - fractionB0, deltaB);\r\n    const g1 = Geometry.conditionalDivideFraction(fractionsOnA.x1 - fractionB0, deltaB);\r\n    if (g0 !== undefined && g1 !== undefined) {\r\n      const detailA = CurveLocationDetail.createCurveEvaluatedFractionFraction(cpA, fractionsOnA.x0, fractionsOnA.x1);\r\n      const detailB = CurveLocationDetail.createCurveEvaluatedFractionFraction(cpB, g0, g1);\r\n      if (reverse)\r\n        detailA.swapFractionsAndPoints();\r\n      return CurveLocationDetailPair.createCapture(detailA, detailB);\r\n    }\r\n    return undefined;\r\n  }\r\n  private appendDetailPair(result: CurveLocationDetailPair[] | undefined, pair: CurveLocationDetailPair | undefined): CurveLocationDetailPair[] | undefined {\r\n    if (pair === undefined)\r\n      return result;\r\n    if (result === undefined)\r\n      return [pair];\r\n    result.push(pair);\r\n    return result;\r\n  }\r\n  /**\r\n   * Test if 2 arcs have coinciding portions.\r\n   * @param arcA\r\n   * @param arcB\r\n   * @param _restrictToBounds\r\n   * @return 0, 1, or 2 overlap points/intervals\r\n   */\r\n  public coincidentArcIntersectionXY(arcA: Arc3d, arcB: Arc3d, _restrictToBounds: boolean = true): CurveLocationDetailPair[] | undefined {\r\n    let result: CurveLocationDetailPair[] | undefined;\r\n    if (arcA.center.isAlmostEqual(arcB.center, this.tolerance)) {\r\n      const matrixBToA = arcA.matrixRef.multiplyMatrixInverseMatrix(arcB.matrixRef);\r\n      if (matrixBToA) {\r\n        const ux = matrixBToA.at(0, 0); const uy = matrixBToA.at(1, 0);\r\n        const vx = matrixBToA.at(0, 1); const vy = matrixBToA.at(1, 1);\r\n        const ru = Geometry.hypotenuseXY(ux, uy);\r\n        const rv = Geometry.hypotenuseXY(vx, vy);\r\n        const dot = Geometry.dotProductXYXY(ux, uy, vx, vy);\r\n        const cross = Geometry.crossProductXYXY(ux, uy, vx, vy);\r\n        if (Geometry.isAlmostEqualNumber(ru, 1.0)\r\n          && Geometry.isAlmostEqualNumber(rv, 1.0)\r\n          && Geometry.isAlmostEqualNumber(0, dot)) {\r\n          const alphaB0Radians = Math.atan2(uy, ux);        // angular position of arcB 0 point in arcA sweep\r\n          const sweepDirection = cross > 0 ? 1.0 : -1.0;    // 1 if arcB parameter space sweeps in same direction as arcA, -1 if opposite\r\n          const betaStartRadians = alphaB0Radians + sweepDirection * arcB.sweep.startRadians;   // arcB start in arcA parameter space\r\n          const betaEndRadians = alphaB0Radians + sweepDirection * arcB.sweep.endRadians;       // arcB end in arcA parameter space\r\n          const fractionSpacesReversed = (sweepDirection * arcA.sweep.sweepRadians * arcB.sweep.sweepRadians) < 0;\r\n          const sweepB = AngleSweep.createStartEndRadians(betaStartRadians, betaEndRadians);\r\n          const sweepA = arcA.sweep;\r\n          const fractionPeriodA = sweepA.fractionPeriod();\r\n          const fractionB0 = sweepA.radiansToPositivePeriodicFraction(sweepB.startRadians);   // arcB start in arcA fraction space\r\n          assert(fractionB0 >= 0.0);\r\n          const fractionSweep = sweepB.sweepRadians / sweepA.sweepRadians;                    // arcB sweep in arcA fraction space\r\n          const fractionB1 = fractionB0 + fractionSweep;                                      // arcB end in arcA fraction space\r\n          const fractionSweepB = Segment1d.create(fractionB0, fractionB1);\r\n\r\n          /** lambda to add a coincident interval or isolated intersection, given inputs in arcA fraction space\r\n           * @param arcBInArcAFractionSpace span of arcB in arcA fraction space. On return, clamped to [0,1] if nontrivial.\r\n           * @param testStartOfArcA if no nontrivial coincident interval was found, look for an isolated intersection at the start (true) or end (false) of arcA\r\n           * @returns whether a detail pair was appended to result\r\n           */\r\n          const appendCoincidentIntersection = (arcBInArcAFractionSpace: Segment1d, testStartOfArcA: boolean): boolean => {\r\n            const size = result ? result.length : 0;\r\n            const arcBStart = arcBInArcAFractionSpace.x0;\r\n            const arcBEnd = arcBInArcAFractionSpace.x1;\r\n            if (arcBInArcAFractionSpace.clampDirectedTo01() && !Geometry.isSmallRelative(arcBInArcAFractionSpace.absoluteDelta())) {\r\n              result = this.appendDetailPair(result, this.createDetailPair(arcA, arcB, arcBInArcAFractionSpace, arcBStart, arcBEnd, fractionSpacesReversed));\r\n            } else {  // test isolated intersection\r\n              const testStartOfArcB = fractionSpacesReversed ? testStartOfArcA : !testStartOfArcA;\r\n              const arcAPt = this._point0 = testStartOfArcA ? arcA.startPoint(this._point0) : arcA.endPoint(this._point0);\r\n              const arcBPt = this._point1 = testStartOfArcB ? arcB.startPoint(this._point1) : arcB.endPoint(this._point1);\r\n              if (arcAPt.isAlmostEqual(arcBPt, this.tolerance)) {\r\n                const detailA = CurveLocationDetail.createCurveFractionPoint(arcA, testStartOfArcA ? 0 : 1, arcAPt);\r\n                const detailB = CurveLocationDetail.createCurveFractionPoint(arcB, testStartOfArcB ? 0 : 1, arcBPt);\r\n                result = this.appendDetailPair(result, CurveLocationDetailPair.createCapture(detailA, detailB));\r\n              }\r\n            }\r\n            return result !== undefined && result.length > size;\r\n          };\r\n\r\n          appendCoincidentIntersection(fractionSweepB, false);  // compute overlap in strict interior, or at end of arcA\r\n\r\n          // check overlap at start of arcA with a periodic shift of fractionSweepB\r\n          if (fractionB1 >= fractionPeriodA)\r\n            appendCoincidentIntersection(Segment1d.create(fractionB0 - fractionPeriodA, fractionB1 - fractionPeriodA), true);\r\n          else if (fractionB0 === 0.0)\r\n            appendCoincidentIntersection(Segment1d.create(fractionB0 + fractionPeriodA, fractionB1 + fractionPeriodA), true);\r\n        }\r\n      }\r\n    }\r\n    return result;\r\n  }\r\n}\r\n"]}

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

@@ -34,28 +34,35 @@ export declare class PackedMatrix3dOps {
     static loadMatrix(dest: Float64Array, a00: number, a01: number, a02: number, a10: number, a11: number, a12: number, a20: number, a21: number, a22: number): void;
     /**
      * Multiply 3x3 matrix `a*b`, store in `result`.
-     * * All params assumed length 9, allocated by caller.
-     * * c may alias either input.
+     * @param a left matrix in product. ASSUMED length 9.
+     * @param b right matrix in product. ASSUMED length 9.
+     * @param result optional destination array of length 9. If insufficient length, a new array is returned. May refer to same array as `a` or `b`.
+     * @return matrix product `a*b`
      */
     static multiplyMatrixMatrix(a: Float64Array, b: Float64Array, result?: Float64Array): Float64Array;
     /**
      * Multiply 3x3 matrix `a*bTranspose`, store in `result`.
-     * * All params assumed length 9, allocated by caller.
-     * * c may alias either input.
+     * @param a left matrix in product. ASSUMED length 9.
+     * @param b transpose of right matrix in product. ASSUMED length 9.
+     * @param result optional destination array of length 9. If insufficient length, a new array is returned. May refer to same array as `a` or `b`.
+     * @return matrix product `a*b^T`
      */
     static multiplyMatrixMatrixTranspose(a: Float64Array, b: Float64Array, result?: Float64Array): Float64Array;
     /**
      * Multiply 3x3 matrix `aTranspose*b`, store in `result`.
-     * * All params assumed length 9, allocated by caller.
-     * * c may alias either input.
+     * @param a transpose of left matrix in product. ASSUMED length 9.
+     * @param b right matrix in product. ASSUMED length 9.
+     * @param result optional destination array of length 9. If insufficient length, a new array is returned. May refer to same array as `a` or `b`.
+     * @return matrix product `a^T*b`
      */
     static multiplyMatrixTransposeMatrix(a: Float64Array, b: Float64Array, result?: Float64Array): Float64Array;
     /** Transpose 3x3 matrix `a` in place */
     static transposeInPlace(a: Float64Array): void;
     /**
-     * Returns the transpose of 3x3 matrix `a`
-     * * If `dest` is passed as argument, then the function copies the transpose of 3x3 matrix `a` into `dest`
-     * * `a` is not changed unless also passed as the dest, i.e., copyTransposed(a,a) transposes `a` in place
+     * Compute transpose of 3x3 matrix `a`, store in `dest`.
+     * @param a source matrix. ASSUMED length 9. Note that `a` is not changed unless also passed as `dest`, i.e., `copyTransposed(a,a)` transposes `a` in place.
+     * @param dest optional destination array of length 9. If insufficient length, a new array is returned. May refer to same array as `a`.
+     * @return matrix `a^T`
      */
     static copyTransposed(a: Float64Array, dest?: Float64Array): Float64Array;
     /** Copy matrix `a` entries into `dest` */
@@ -146,8 +153,7 @@ export declare class Matrix3d implements BeJSONFunctions {
     freeze(): Readonly<this>;
     /**
      * Constructor
-     * @param coffs optional coefficient array.
-     * * **WARNING:** coffs is captured (i.e., is now owned by the Matrix3d object and can be modified by it).
+     * @param coffs optional coefficient array of length at least 9 (CAPTURED). If undefined or insufficient length, a new array is created.
      */
     constructor(coffs?: Float64Array);
     /**
@@ -219,6 +225,7 @@ export declare class Matrix3d implements BeJSONFunctions {
      * Create a Matrix3d with caller-supplied coefficients and optional inverse coefficients.
      * * The inputs are captured into (i.e., owned by) the new Matrix3d.
      * * The caller is responsible for validity of the inverse coefficients.
+     * * If either array is insufficiently sized, it is ignored.
      * @param coffs (required) array of 9 coefficients.
      * @param inverseCoffs (optional) array of 9 coefficients.
      * @returns a Matrix3d populated by a coffs array.
@@ -769,7 +776,7 @@ export declare class Matrix3d implements BeJSONFunctions {
      * @param y y part of multiplied point
      * @param z z part of multiplied point
      * @param w w part of multiplied point
-     * @param result optional preallocated result.
+     * @param result optional preallocated result. If length < 4, a new array is returned.
      */
     static xyzPlusMatrixTimesWeightedCoordinatesToFloat64Array(origin: XYZ, matrix: Matrix3d, x: number, y: number, z: number, w: number, result?: Float64Array): Float64Array;
     /**
@@ -784,7 +791,7 @@ export declare class Matrix3d implements BeJSONFunctions {
      * @param x x part of multiplied point
      * @param y y part of multiplied point
      * @param z z part of multiplied point
-     * @param result optional preallocated result.
+     * @param result optional preallocated result. If length < 3, a new array is returned.
      */
     static xyzPlusMatrixTimesCoordinatesToFloat64Array(origin: XYZ, matrix: Matrix3d, x: number, y: number, z: number, result?: Float64Array): Float64Array;
     /**