@itwin/core-geometry 3.1.0-dev.8 → 3.2.0-dev.3

package/lib/esm/bspline/BSplineCurveOps.js.map

@@ -1 +1 @@
-{"version":3,"file":"BSplineCurveOps.js","sourceRoot":"","sources":["../../../src/bspline/BSplineCurveOps.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAEH,8EAA8E;AAE9E,OAAO,EAAE,gBAAgB,EAAE,MAAM,gCAAgC,CAAC;AAClE,OAAO,EAAE,oBAAoB,EAAE,MAAM,oCAAoC,CAAC;AAC1E,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AACvC,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AAClE,OAAO,EAAE,YAAY,EAAE,MAAM,0BAA0B,CAAC;AACxD,OAAO,EAAE,cAAc,EAAE,MAAM,gBAAgB,CAAC;AAChD,OAAO,EAAE,eAAe,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAE3D,OAAO,EAAE,YAAY,EAAE,MAAM,4BAA4B,CAAC;AAE1D;;;GAGG;AACH,MAAM,OAAO,eAAe;IAC1B;;;;OAIG;IACI,MAAM,CAAC,mBAAmB,CAAC,MAAwC,EAAE,KAAa;QACvF,MAAM,SAAS,GAAG,MAAM,CAAC,MAAM,CAAC;QAChC,IAAI,KAAK,GAAG,SAAS,IAAI,KAAK,GAAG,CAAC;YAChC,OAAO,SAAS,CAAC;QACnB,MAAM,MAAM,GAAG,KAAK,GAAG,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,CAAI,oDAAoD;QAClF,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,SAAS,CAAC,CAAC;QAChD,MAAM,cAAc,GAAG,IAAI,YAAY,CAAC,KAAK,CAAC,CAAC;QAC/C,MAAM,GAAG,GAAG,IAAI,gBAAgB,EAAE,CAAC;QACnC,MAAM,KAAK,GAAG,UAAU,CAAC,oBAAoB,CAAC,SAAS,EAAE,KAAK,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;QAC9E,MAAM,GAAG,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;QAC7B,KAAK,IAAI,cAAc,GAAG,CAAC,EAAE,cAAc,GAAG,SAAS,EAAE,cAAc,EAAE,EAAE;YACzE,MAAM,CAAC,GAAG,KAAK,CAAC,YAAY,CAAC,cAAc,CAAC,CAAC;YAC7C,MAAM,SAAS,GAAG,KAAK,CAAC,mBAAmB,CAAC,CAAC,CAAC,CAAC;YAC/C,KAAK,CAAC,sBAAsB,CAAC,SAAS,EAAE,CAAC,EAAE,cAAc,CAAC,CAAC;YAC3D,oEAAoE;YACpE,IAAI,QAAQ,GAAG,CAAC,CAAC;YACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE;gBAC5B,IAAI,cAAc,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,QAAQ,CAAC;oBAC9C,QAAQ,GAAG,CAAC,CAAC;YACjB,MAAM,2BAA2B,GAAG,MAAM,GAAG,QAAQ,CAAC;YACtD,MAAM,QAAQ,GAAG,cAAc,GAAG,EAAE,CAAC;YACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE;gBAC9B,MAAM,UAAU,GAAG,cAAc,GAAG,MAAM,GAAG,2BAA2B,GAAG,CAAC,CAAC;gBAC7E,IAAI,QAAQ,GAAG,UAAU,IAAI,CAAC,IAAI,UAAU,GAAG,SAAS;oBACtD,MAAM,CAAC,QAAQ,GAAG,2BAA2B,GAAG,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC,CAAC;aAC1E;YACD,IAAI,MAAM,YAAY,oBAAoB,EAAE;gBAC1C,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,+BAA+B,CAAC,cAAc,EAAE,GAAG,CAAC,CAAC,CAAC;aACvE;iBAAM;gBACL,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,cAAc,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC;aAC1C;SACF;QACD,MAAM,KAAK,GAAG,YAAY,CAAC,4BAA4B,CAAC,SAAS,EAAE,EAAE,EAAE,MAAM,EAAE,CAAC,EAAE,GAAG,CAAC,WAAW,EAAE,CAAC,CAAC;QACrG,IAAI,KAAK,EAAE;YACT,OAAO,cAAc,CAAC,MAAM,CAAC,KAAK,EAAE,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;SACzD;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IAED;;;OAGG;IACI,MAAM,CAAC,0BAA0B,CAAC,OAAoC;QAC3E,+GAA+G;QAC/G,wEAAwE;QACxE,MAAM,gBAAgB,GAAG,OAAO,CAAC,KAAK,EAAE,CAAC;QACzC,IAAI,CAAC,IAAI,CAAC,UAAU,CAAC,eAAe,CAAC,gBAAgB,CAAC;YACpD,OAAO,SAAS,CAAC;QAEnB,MAAM,KAAK,GAAG,IAAI,CAAC,UAAU,CAAC,cAAc,CAAC,gBAAgB,CAAC,CAAC;QAC/D,IAAI,SAAS,KAAK,KAAK;YACrB,OAAO,SAAS,CAAC;QAEnB,MAAM,SAAS,GAAG,IAAI,CAAC,UAAU,CAAC,iCAAiC,CAAC,gBAAgB,CAAC,KAAK,EAAE,gBAAgB,CAAC,MAAM,CAAC,CAAC;QACrH,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,SAAS,CAAC;QAEnB,MAAM,WAAW,GAAG,cAAc,CAAC,MAAM,CAAC,KAAK,EAAE,SAAS,EAAE,gBAAgB,CAAC,KAAK,CAAC,CAAC;QAEpF,IAAI,gBAAgB,CAAC,MAAM;YACzB,WAAW,aAAX,WAAW,uBAAX,WAAW,CAAE,YAAY,CAAC,eAAe,CAAC,yBAAyB,CAAC,CAAC;QAEvE,OAAO,WAAW,CAAC;IACrB,CAAC;CACF;AAED;;;GAGG;AACH,WAAiB,eAAe;IAC9B;;;;;OAKG;IACH,MAAa,UAAU;QACrB;;WAEG;QACK,MAAM,CAAC,cAAc,CAAC,KAA2B;YACvD,IAAI,SAAS,KAAK,KAAK,IAAI,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE;gBAC3C,KAAK,GAAG,SAAS,CAAC;gBAClB,OAAO,KAAK,CAAC;aACd;YACD,MAAM,OAAO,GAAG,UAAU,CAAC,MAAM,CAAC,KAAK,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;YACnD,IAAI,OAAO,CAAC,YAAY,GAAG,UAAU,CAAC,aAAa,EAAE;gBACnD,KAAK,GAAG,SAAS,CAAC;gBAClB,OAAO,KAAK,CAAC;aACd;YACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,QAAQ,EAAE,EAAE,CAAC,EAAE;gBACzC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,sBAAsB,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;aACnD;YACD,KAAK,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;YAC9B,OAAO,IAAI,CAAC;QACd,CAAC;QAED,qEAAqE;QAC7D,MAAM,CAAC,8BAA8B,CAAC,SAAoB;YAChE,IAAI,SAAS,CAAC,MAAM,GAAG,CAAC;gBACtB,OAAO,SAAS,CAAC;YACnB,MAAM,MAAM,GAAa,CAAC,GAAG,CAAC,CAAC;YAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,CAAC;gBACvC,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YACtE,IAAI,CAAC,IAAI,CAAC,cAAc,CAAC,MAAM,CAAC;gBAC9B,OAAO,SAAS,CAAC;YACnB,OAAO,MAAM,CAAC;QAChB,CAAC;QAED,gEAAgE;QACxD,MAAM,CAAC,0BAA0B,CAAC,SAAiB;YACzD,IAAI,SAAS,GAAG,CAAC;gBACf,OAAO,SAAS,CAAC;YACnB,MAAM,KAAK,GAAG,UAAU,CAAC,oBAAoB,CAAC,SAAS,GAAG,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;YAC1E,MAAM,MAAM,GAAa,EAAE,CAAC;YAC5B,KAAK,IAAI,CAAC,GAAG,KAAK,CAAC,aAAa,EAAE,CAAC,IAAI,KAAK,CAAC,cAAc,EAAE,EAAE,CAAC;gBAC9D,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YAC9B,OAAO,MAAM,CAAC;QAChB,CAAC;QAED,qEAAqE;QAC7D,MAAM,CAAC,wBAAwB,CAAC,OAAoC;YAC1E,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK,IAAI,OAAO,CAAC,KAAK,CAAC,MAAM,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM;gBAClF,OAAO,CAAC,KAAK,GAAG,SAAS,CAAC;YAE5B,gDAAgD;YAChD,MAAM,MAAM,GAAG,gBAAgB,CAAC,MAAM,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;YAC1D,MAAM,OAAO,GAAG,MAAM,CAAC,qBAAqB,EAAE,CAAC;YAC/C,MAAM,CAAC,KAAK,EAAE,CAAC;YAEf,8BAA8B;YAC9B,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,MAAM,GAAG,CAAC,EAAE,KAAK,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,EAAE,EAAE,KAAK,EAAE;gBACzE,IAAI,KAAK,KAAK,OAAO,CAAC,MAAM,CAAC;oBAC3B,EAAE,MAAM,CAAC,CAAC,qBAAqB;;oBAE/B,MAAM,CAAC,IAAI,CAAC,OAAO,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC;aACjD;YACD,OAAO,CAAC,SAAS,GAAG,MAAM,CAAC,eAAe,EAAE,CAAC;YAE7C,oDAAoD;YACpD,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK,EAAE;gBAC/B,MAAM,QAAQ,GAAa,EAAE,CAAC;gBAC9B,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,MAAM,GAAG,CAAC,EAAE,KAAK,GAAG,OAAO,CAAC,KAAK,CAAC,MAAM,EAAE,EAAE,KAAK,EAAE;oBACrE,IAAI,KAAK,KAAK,OAAO,CAAC,MAAM,CAAC;wBAC3B,EAAE,MAAM,CAAC,CAAC,OAAO;;wBAEjB,QAAQ,CAAC,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC;iBACvC;gBACD,OAAO,CAAC,KAAK,GAAG,QAAQ,CAAC,KAAK,EAAE,CAAC;aAClC;YACD,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;;UAKE;QACM,MAAM,CAAC,gCAAgC,CAAC,SAAoB,EAAE,aAAiC,EAAE,MAA2B;YAClI,IAAI,MAA4B,CAAC;YACjC,IAAI,aAAa,IAAI,CAAC,MAAM;gBAC1B,MAAM,GAAG,IAAI,CAAC,8BAA8B,CAAC,SAAS,CAAC,CAAC;YAC1D,IAAI,SAAS,KAAK,MAAM;gBACtB,MAAM,GAAG,IAAI,CAAC,0BAA0B,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC;YAC7D,OAAO,MAAM,CAAC;QAChB,CAAC;QAED;;;WAGG;QACK,MAAM,CAAC,sBAAsB,CAAC,OAAoC;;YACxE,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK;gBAC7B,OAAO,CAAC,KAAK,GAAG,IAAI,CAAC,gCAAgC,CAAC,OAAO,CAAC,SAAS,EAAE,OAAO,CAAC,eAAe,EAAE,OAAO,CAAC,MAAM,CAAC,CAAC;YACpH,OAAO,CAAA,MAAA,OAAO,CAAC,KAAK,0CAAE,MAAM,MAAK,OAAO,CAAC,SAAS,CAAC,MAAM,CAAC;QAC5D,CAAC;QAED;;;;;WAKG;QACK,MAAM,CAAC,qBAAqB,CAAC,KAAe,EAAE,IAAc,EAAE,KAAe,EAAE,KAAa,EAAE,WAAmB,EAAE,MAAc,EAAE,YAAoB,EAAE,YAAoB;YACnL,IAAI,eAAe,GAAG,GAAG,GAAG,CAAC,YAAY,GAAG,YAAY,GAAG,MAAM,CAAC,CAAC;YACnE,KAAK,CAAC,KAAK,CAAC,GAAG,MAAM,GAAG,MAAM,GAAG,eAAe,CAAC;YACjD,IAAI,CAAC,KAAK,CAAC,GAAG,MAAM,GAAG,CAAC,YAAY,GAAG,YAAY,CAAC,GAAG,eAAe,CAAC;YAEvE,eAAe,GAAG,GAAG,GAAG,CAAC,YAAY,GAAG,MAAM,GAAG,WAAW,CAAC,CAAC;YAC9D,IAAI,CAAC,KAAK,CAAC,IAAI,YAAY,GAAG,CAAC,MAAM,GAAG,WAAW,CAAC,GAAG,eAAe,CAAC;YACvE,KAAK,CAAC,KAAK,CAAC,GAAG,YAAY,GAAG,YAAY,GAAG,eAAe,CAAC;YAE7D,eAAe,GAAG,GAAG,GAAG,CAAC,YAAY,GAAG,MAAM,CAAC,CAAC;YAChD,KAAK,CAAC,KAAK,CAAC,IAAI,eAAe,CAAC;YAChC,IAAI,CAAC,KAAK,CAAC,IAAI,eAAe,CAAC;YAC/B,KAAK,CAAC,KAAK,CAAC,IAAI,eAAe,CAAC;QAClC,CAAC;QAED;;;;WAIG;QACK,MAAM,CAAC,kBAAkB,CAAC,KAAe,EAAE,IAAc,EAAE,KAAe;YAChF,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC;gBAC/D,OAAO,KAAK,CAAC;YAEf,kBAAkB;YAClB,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;YAChD,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;YACxB,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;;;;WAOG;QACM,MAAM,CAAC,kBAAkB,CAAC,KAAe,EAAE,IAAc,EAAE,KAAe,EAAE,OAAoC,EAAE,sBAA+B,EAAE,oBAA6B;YACvL,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK;gBAC7B,OAAO,KAAK,CAAC;YACf,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC;gBAC/D,OAAO,KAAK,CAAC;YACf,IAAI,OAAO,CAAC,KAAK,CAAC,MAAM,KAAK,CAAC,IAAI,OAAO,CAAC,SAAS,CAAC,MAAM,KAAK,CAAC;gBAC9D,OAAO,KAAK,CAAC;YAEf,IAAI,WAAW,GAAG,CAAC,EAAE,MAAM,GAAG,CAAC,EAAE,YAAY,GAAG,CAAC,EAAE,YAAY,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,aAAa,GAAG,CAAC,CAAC;YAEhG,YAAY;YACZ,IAAI,sBAAsB,EAAE;gBAC1B,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;gBACf,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAC7C,WAAW,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAClD,GAAG,GAAG,MAAM,GAAG,WAAW,CAAC;gBAC3B,aAAa,GAAG,GAAG,GAAG,GAAG,CAAC;gBAC1B,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,GAAG,aAAa,CAAC;gBACzC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,aAAa,CAAC;aACpC;iBAAM;gBACL,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;gBAC1B,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;aACf;YAED,aAAa;YACb,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YACnD,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YAC7C,aAAa,GAAG,GAAG,GAAG,CAAC,YAAY,GAAG,MAAM,CAAC,CAAC;YAC9C,aAAa,IAAI,aAAa,CAAC;YAC/B,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,MAAM,GAAG,aAAa,CAAC;YAC3C,IAAI,CAAC,CAAC,CAAC,GAAI,GAAG,GAAG,CAAC,MAAM,GAAG,YAAY,CAAC,GAAG,aAAa,CAAC;YACzD,KAAK,CAAC,CAAC,CAAC,GAAG,YAAY,GAAG,YAAY,GAAG,aAAa,CAAC;YAEvD,WAAW;YACX,IAAI,oBAAoB,EAAE;gBACxB,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBACnD,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBACnD,GAAG,GAAG,YAAY,GAAG,YAAY,CAAC;gBAClC,aAAa,GAAG,GAAG,GAAG,GAAG,CAAC;gBAC1B,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,YAAY,GAAG,aAAa,CAAC;gBACzC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,YAAY,GAAG,GAAG,CAAC,GAAG,aAAa,CAAC;gBAC/C,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;aAChB;iBAAM;gBACL,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;gBAC1B,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;aACf;YACD,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;;;;WAOG;QACK,MAAM,CAAC,wBAAwB,CAAC,KAAe,EAAE,IAAc,EAAE,KAAe,EAAE,OAAoC,EAAE,sBAA+B,EAAE,oBAA6B;YAC5L,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK;gBAC7B,OAAO,KAAK,CAAC;YACf,IAAI,KAAK,CAAC,MAAM,KAAK,IAAI,CAAC,MAAM,IAAI,KAAK,CAAC,MAAM,KAAK,KAAK,CAAC,MAAM,IAAI,KAAK,CAAC,MAAM,KAAK,OAAO,CAAC,KAAK,CAAC,MAAM;gBACxG,OAAO,KAAK,CAAC;YACf,IAAI,OAAO,CAAC,KAAK,CAAC,MAAM,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM;gBACnD,OAAO,KAAK,CAAC;YAEf,MAAM,YAAY,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;YAClD,MAAM,kBAAkB,GAAG,YAAY,GAAG,CAAC,CAAC;YAC5C,IAAI,WAAW,GAAG,CAAC,EAAE,MAAM,GAAG,CAAC,EAAE,YAAY,GAAG,CAAC,EAAE,YAAY,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,aAAa,GAAG,CAAC,CAAC;YAEhG,IAAI,OAAO,CAAC,MAAM,EAAE;gBAClB,YAAY;gBACZ,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAC7C,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,CAAC;gBACzF,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,CAAC,CAAC;gBACzF,WAAW,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAClD,IAAI,CAAC,qBAAqB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,WAAW,EAAE,MAAM,EAAE,YAAY,EAAE,YAAY,CAAC,CAAC;gBAEnG,aAAa;gBACb,YAAY,GAAG,YAAY,CAAC;gBAC5B,YAAY,GAAG,MAAM,CAAC;gBACtB,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAC7C,WAAW,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAClD,IAAI,CAAC,qBAAqB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,WAAW,EAAE,MAAM,EAAE,YAAY,EAAE,YAAY,CAAC,CAAC;gBAEnG,qDAAqD;gBACrD,WAAW,GAAG,YAAY,CAAC;gBAC3B,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,CAAC,CAAC;gBACnF,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,CAAC;gBAC7F,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,CAAC;gBACzF,IAAI,CAAC,qBAAqB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,kBAAkB,EAAE,WAAW,EAAE,MAAM,EAAE,YAAY,EAAE,YAAY,CAAC,CAAC;aACrH;iBAAM,EAAE,OAAO;gBACd,YAAY;gBACZ,IAAI,sBAAsB,EAAE;oBAC1B,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;oBACf,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;oBAC7C,WAAW,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;oBAClD,GAAG,GAAG,MAAM,GAAG,WAAW,CAAC;oBAC3B,aAAa,GAAG,GAAG,GAAG,GAAG,CAAC;oBAC1B,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,GAAG,aAAa,CAAC;oBACzC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,aAAa,CAAC;iBACpC;qBAAM;oBACL,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;oBAC1B,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;iBACf;gBAED,aAAa;gBACb,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAC7C,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBACnD,YAAY,GAAG,GAAG,CAAC;gBACnB,WAAW,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAClD,IAAI,CAAC,qBAAqB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,WAAW,EAAE,MAAM,EAAE,YAAY,EAAE,YAAY,CAAC,CAAC;gBAEnG,kBAAkB;gBAClB,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,CAAC,CAAC;gBACnF,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,CAAC;gBACzF,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,CAAC;gBAC7F,WAAW,GAAG,GAAG,CAAC;gBAClB,IAAI,CAAC,qBAAqB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,kBAAkB,EAAE,WAAW,EAAE,MAAM,EAAE,YAAY,EAAE,YAAY,CAAC,CAAC;gBAEpH,WAAW;gBACX,IAAI,oBAAoB,EAAE;oBACxB,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,YAAY,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC;oBAC7E,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC;oBACjF,GAAG,GAAG,YAAY,GAAG,YAAY,CAAC;oBAClC,aAAa,GAAG,GAAG,GAAG,GAAG,CAAC;oBAC1B,KAAK,CAAC,YAAY,CAAC,GAAG,CAAC,YAAY,GAAG,aAAa,CAAC;oBACpD,IAAI,CAAC,YAAY,CAAC,GAAG,CAAC,YAAY,GAAG,GAAG,CAAC,GAAG,aAAa,CAAC;oBAC1D,KAAK,CAAC,YAAY,CAAC,GAAG,GAAG,CAAC;iBAC3B;qBAAM;oBACL,KAAK,CAAC,YAAY,CAAC,GAAG,KAAK,CAAC,YAAY,CAAC,GAAG,GAAG,CAAC;oBAChD,IAAI,CAAC,YAAY,CAAC,GAAG,GAAG,CAAC;iBAC1B;aACF;YAED,cAAc;YACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,kBAAkB,EAAE,EAAE,CAAC,EAAE;gBAC3C,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBACjD,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC3D,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBACvD,WAAW,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC1D,IAAI,CAAC,qBAAqB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,WAAW,EAAE,MAAM,EAAE,YAAY,EAAE,YAAY,CAAC,CAAC;aACpG;YACD,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;;WAKG;QACK,MAAM,CAAC,WAAW,CAAC,KAAe,EAAE,IAAc,EAAE,KAAe,EAAE,OAAoC;YAC/G,IAAI,sBAAsB,GAAG,KAAK,CAAC;YACnC,IAAI,oBAAoB,GAAG,KAAK,CAAC;YACjC,IAAI,OAAO,CAAC,iBAAiB,IAAI,CAAC,OAAO,CAAC,MAAM,EAAE;gBAChD,sBAAsB,GAAG,CAAC,SAAS,KAAK,OAAO,CAAC,YAAY,CAAC,CAAC;gBAC9D,oBAAoB,GAAG,CAAC,SAAS,KAAK,OAAO,CAAC,UAAU,CAAC,CAAC;aAC3D;YAED,IAAI,SAAS,GAAG,KAAK,CAAC;YACtB,IAAI,CAAC,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM;gBAChC,SAAS,GAAG,IAAI,CAAC,kBAAkB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC;iBACrD,IAAI,CAAC,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM;gBACrC,SAAS,GAAG,IAAI,CAAC,kBAAkB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,OAAO,EAAE,sBAAsB,EAAE,oBAAoB,CAAC,CAAC;iBAC5G,IAAI,CAAC,IAAI,OAAO,CAAC,SAAS,CAAC,MAAM;gBACpC,SAAS,GAAG,IAAI,CAAC,wBAAwB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,OAAO,EAAE,sBAAsB,EAAE,oBAAoB,CAAC,CAAC;YAEvH,OAAO,SAAS,CAAC;QACnB,CAAC;QAED;;;;;WAKG;QACK,MAAM,CAAC,qBAAqB,CAAC,OAAkB,EAAE,OAAoC,EAAE,OAAgB;YAC7G,IAAI,OAAO,CAAC,MAAM,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC;gBACjD,OAAO,KAAK,CAAC;YACf,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK;gBAC7B,OAAO,KAAK,CAAC;YAEf,MAAM,KAAK,GAAG,GAAG,GAAC,GAAG,CAAC;YACtB,MAAM,YAAY,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;YAElD,IAAI,CAAC,KAAK,YAAY,EAAE,EAAE,gBAAgB;gBACxC,IAAI,OAAO;oBACT,OAAO,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;;oBAEtD,OAAO,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;gBACxD,OAAO,IAAI,CAAC;aACb;YAED,IAAI,CAAC,KAAK,YAAY,EAAE;gBACtB,MAAM,KAAK,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC5F,MAAM,IAAI,GAAG,GAAG,GAAG,KAAK,CAAC;gBACzB,MAAM,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,KAAK,GAAG,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;gBAC1F,IAAI,OAAO;oBACT,OAAO,CAAC,gBAAgB,CAAC,IAAI,EAAE,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,WAAW,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;;oBAElH,OAAO,CAAC,gBAAgB,CAAC,IAAI,EAAE,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,WAAW,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;gBAClH,OAAO,IAAI,CAAC;aACb;YAED,mBAAmB;YACnB,IAAI,OAAO,EAAE;gBACT,MAAM,KAAK,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC5F,MAAM,IAAI,GAAG,GAAG,GAAG,KAAK,CAAC;gBACzB,MAAM,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,KAAK,GAAG,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;gBAC1F,OAAO,CAAC,gBAAgB,CAAC,IAAI,EAAE,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,WAAW,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;aACrH;iBAAM;gBACH,MAAM,KAAK,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,YAAY,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,YAAY,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC;gBAChJ,MAAM,IAAI,GAAG,GAAG,GAAG,KAAK,CAAC;gBACzB,MAAM,IAAI,GAAG,OAAO,CAAC,YAAY,CAAC,CAAC,WAAW,CAAC,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,EAAE,CAAC,KAAK,GAAG,KAAK,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;gBACnI,OAAO,CAAC,gBAAgB,CAAC,IAAI,EAAE,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,WAAW,CAAC,KAAK,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC;aAChK;YACD,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;;;WAMG;QACK,MAAM,CAAC,sBAAsB,CAAC,OAAkB,EAAE,OAAoC,EAAE,OAAgB;YAC9G,IAAI,OAAO,CAAC,MAAM,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC;gBACjD,OAAO,KAAK,CAAC;YAEf,MAAM,YAAY,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;YAClD,IAAI,CAAC,KAAK,YAAY;gBACpB,OAAO,IAAI,CAAC,qBAAqB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;YAE/D,IAAI,OAAO;gBACT,OAAO,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;;gBAExB,OAAO,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;YAC5D,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;;;WAMG;QACK,MAAM,CAAC,gCAAgC,CAAC,OAAkB,EAAE,OAAoC,EAAE,OAAgB;YACxH,IAAI,OAAO,CAAC,MAAM,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC;gBACjD,OAAO,KAAK,CAAC;YAEf,MAAM,OAAO,GAAG,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC,YAAY,CAAC,CAAC,CAAC,OAAO,CAAC,UAAU,CAAC;YACpE,IAAI,SAAS,KAAK,OAAO;gBACvB,OAAO,KAAK,CAAC;YAEf,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,+BAA+B;YAC7C,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,2FAA2F;YACzG,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,oCAAoC;YAElD,MAAM,YAAY,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;YAClD,IAAI,CAAC,KAAK,YAAY,EAAE,EAAE,yBAAyB;gBACjD,IAAI,OAAO,EAAE;oBACX,IAAI,GAAG,CAAC,CAAC;oBACT,IAAI,GAAG,CAAC,CAAC;oBACT,IAAI,GAAG,CAAC,CAAC;iBACV;qBAAM;oBACL,IAAI,GAAG,CAAC,CAAC;oBACT,IAAI,GAAG,CAAC,CAAC;oBACT,IAAI,GAAG,CAAC,CAAC;iBACV;aACF;iBAAM;gBACL,IAAI,OAAO,EAAE;oBACX,IAAI,GAAG,CAAC,CAAC;oBACT,IAAI,GAAG,CAAC,CAAC;oBACT,IAAI,GAAG,CAAC,CAAC;iBACV;qBAAM;oBACL,IAAI,GAAG,YAAY,GAAG,CAAC,CAAC;oBACxB,IAAI,GAAG,YAAY,GAAG,CAAC,CAAC;oBACxB,IAAI,GAAG,YAAY,CAAC;iBACrB;aACF;YAED,kCAAkC;YAClC,MAAM,WAAW,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC;YAC1D,OAAO,CAAC,IAAI,CAAC,CAAC,UAAU,CAAC,OAAO,EAAE,WAAW,GAAG,GAAG,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC;YACpE,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;;WAKG;QACK,MAAM,CAAC,iCAAiC,CAAC,OAAkB,EAAE,OAAoC,EAAE,OAAgB;YACzH,IAAI,OAAO,CAAC,MAAM,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC;gBACjD,OAAO,KAAK,CAAC;YAEf,MAAM,OAAO,GAAG,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC,YAAY,CAAC,CAAC,CAAC,OAAO,CAAC,UAAU,CAAC;YACpE,IAAI,SAAS,KAAK,OAAO;gBACvB,OAAO,KAAK,CAAC;YAEf,uCAAuC;YACvC,IAAI,CAAC,IAAI,CAAC,qBAAqB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC;gBACxD,OAAO,KAAK,CAAC;YAEf,MAAM,YAAY,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;YAClD,MAAM,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,+BAA+B;YAC5E,MAAM,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,2FAA2F;YAExI,kFAAkF;YAClF,OAAO,CAAC,IAAI,CAAC,CAAC,UAAU,CAAC,OAAO,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC;YACxF,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;WAIG;QACK,MAAM,CAAC,+BAA+B,CAAC,OAAkB,EAAE,OAAoC;YACrG,MAAM,YAAY,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;YAClD,IAAI,CAAC,OAAO,CAAC,kBAAkB;mBACxB,YAAY,IAAI,CAAC;mBACjB,CAAC,SAAS,KAAK,OAAO,CAAC,YAAY,IAAI,SAAS,KAAK,OAAO,CAAC,UAAU,CAAC;mBACxE,OAAO,CAAC,iBAAiB;mBACzB,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,EAAE;gBAC3D,OAAO,IAAI,CAAC;aACX;YACH,qFAAqF;YACrF,IAAI,SAAS,KAAK,OAAO,CAAC,YAAY,EAAE,EAAE,4DAA4D;gBACpG,MAAM,mBAAmB,GAAG,QAAQ,CAAC,cAAc,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;gBACxF,IAAI,SAAS,KAAK,mBAAmB,EAAE;oBACrC,MAAM,aAAa,GAAG,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC;oBACpF,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,UAAU,CAAC,mBAAmB,EAAE,aAAa,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC;iBACrG;aACF;iBAAM,IAAI,SAAS,KAAK,OAAO,CAAC,UAAU,EAAE,EAAG,4DAA4D;gBAC1G,MAAM,iBAAiB,GAAG,QAAQ,CAAC,cAAc,CAAC,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;gBACpH,IAAI,SAAS,KAAK,iBAAiB,EAAE;oBACnC,MAAM,eAAe,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;oBACxD,OAAO,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,iBAAiB,EAAE,eAAe,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;iBACvE;aACF;iBAAM,EAAG,8DAA8D;gBACtE,MAAM,aAAa,GAAG,QAAQ,CAAC,cAAc,CAAC,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;gBACjG,IAAI,SAAS,KAAK,aAAa,EAAE;oBAC/B,MAAM,eAAe,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;oBACxD,OAAO,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,aAAa,EAAE,eAAe,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;oBAClE,MAAM,aAAa,GAAG,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC;oBACpF,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,UAAU,CAAC,aAAa,EAAE,CAAC,aAAa,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC;iBAChG;aACF;YACD,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;WAIG;QACK,MAAM,CAAC,gBAAgB,CAAC,OAAkB,EAAE,OAAoC;YACtF,IAAI,OAAO,CAAC,MAAM,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM;gBAC7C,OAAO,KAAK,CAAC;YAEf,2CAA2C;YAC3C,MAAM,MAAM,GAAG,OAAO,CAAC,UAAU,EAAE,CAAC;YACpC,MAAM,MAAM,GAAG,OAAO,CAAC,UAAU,EAAE,CAAC;YACpC,OAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;YAC7B,OAAO,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;YAE9C,IAAI,SAAS,GAAG,KAAK,CAAC;YACtB,IAAI,SAAS,KAAK,OAAO,CAAC,YAAY,EAAE;gBACtC,IAAI,OAAO,CAAC,iBAAiB;oBAC3B,SAAS,GAAG,IAAI,CAAC,sBAAsB,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;;oBAEhE,SAAS,GAAG,IAAI,CAAC,qBAAqB,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;aAClE;iBAAM,EAAE,qBAAqB;gBAC5B,IAAI,OAAO,CAAC,kBAAkB;oBAC5B,SAAS,GAAG,IAAI,CAAC,gCAAgC,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;;oBAE1E,SAAS,GAAG,IAAI,CAAC,iCAAiC,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;aAC9E;YAED,IAAI,SAAS,KAAK,OAAO,CAAC,UAAU,EAAE;gBACpC,IAAI,OAAO,CAAC,iBAAiB;oBAC3B,SAAS,GAAG,IAAI,CAAC,sBAAsB,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;;oBAEjE,SAAS,GAAG,IAAI,CAAC,qBAAqB,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;aACnE;iBAAM,EAAE,mBAAmB;gBAC1B,IAAI,OAAO,CAAC,kBAAkB;oBAC5B,SAAS,GAAG,IAAI,CAAC,gCAAgC,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;;oBAE3E,SAAS,GAAG,IAAI,CAAC,iCAAiC,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;aAC/E;YAED,IAAI,SAAS;gBACX,SAAS,GAAG,IAAI,CAAC,+BAA+B,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;YAErE,OAAO,SAAS,CAAC;QACnB,CAAC;QAED;;WAEG;QACK,MAAM,CAAC,oBAAoB,CAAC,MAAiB,EAAE,KAAe,EAAE,IAAc,EAAE,KAAe;YACrG,IAAI,KAAK,CAAC,MAAM,KAAK,IAAI,CAAC,MAAM,IAAI,KAAK,CAAC,MAAM,KAAK,KAAK,CAAC,MAAM,IAAI,KAAK,CAAC,MAAM,KAAK,MAAM,CAAC,MAAM;gBACjG,OAAO,SAAS,CAAC;YACnB,MAAM,KAAK,GAAc,EAAE,CAAC;YAC5B,MAAM,YAAY,GAAG,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC;YACvC,MAAM,OAAO,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAE,wBAAwB;YAC9D,IAAI,GAAG,GAAuB,GAAG,CAAC;YAElC,6BAA6B;YAC7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,YAAY,EAAE,EAAE,CAAC,EAAE;gBACrC,IAAI,SAAS,KAAK,CAAC,GAAG,GAAG,QAAQ,CAAC,yBAAyB,CAAC,CAAE,KAAK,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;oBACnF,OAAO,SAAS,CAAC;gBACnB,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC9B,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC9B,OAAO,CAAC,CAAC,CAAC,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;aAClD;YAED,8BAA8B;YAC9B,IAAI,SAAS,KAAK,CAAC,GAAG,GAAG,QAAQ,CAAC,yBAAyB,CAAC,GAAG,EAAE,IAAI,CAAC,YAAY,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC;gBACjH,OAAO,SAAS,CAAC;YACnB,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC;YAC/B,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC;YAC5C,KAAK,IAAI,CAAC,GAAG,YAAY,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,EAAE;gBAC1C,IAAI,SAAS,KAAK,CAAC,GAAG,GAAG,QAAQ,CAAC,yBAAyB,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;oBACxE,OAAO,SAAS,CAAC;gBACnB,OAAO,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAE,KAAK,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,EAAE,CAAE,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;gBAChI,KAAK,CAAC,CAAC,CAAC,GAAG,CAAE,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;aACnF;YAED,8BAA8B;YAC9B,IAAI,SAAS,KAAK,CAAC,GAAG,GAAG,QAAQ,CAAC,yBAAyB,CAAC,GAAG,EAAE,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC/E,OAAO,SAAS,CAAC;YACnB,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;YACjD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,YAAY,EAAE,EAAE,CAAC,EAAE;gBACrC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;aACzD;YACD,OAAO,KAAK,CAAC;QACf,CAAC;QAED;;WAEG;QACI,MAAM,CAAC,eAAe,CAAC,OAAoC;YAChE,OAAO,CAAC,KAAK,GAAG,CAAC,CAAC;YAElB,yGAAyG;YACzG,OAAO,CAAC,KAAK,GAAG,IAAI,CAAC,iCAAiC,CAAC,OAAO,CAAC,KAAK,EAAE,OAAO,CAAC,SAAS,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;YAEpG,oEAAoE;YACtE,IAAI,CAAC,IAAI,CAAC,wBAAwB,CAAC,OAAO,CAAC;gBACzC,OAAO,KAAK,CAAC;YAEf,uDAAuD;YACvD,IAAI,eAAe,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC,OAAO,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;YAC1G,IAAI,CAAC,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM,IAAI,eAAe,EAAE;gBACrD,OAAO,CAAC,SAAS,CAAC,GAAG,EAAE,CAAC;gBACxB,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK;oBAC7B,OAAO,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC;gBACtB,eAAe,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC,OAAO,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;aACvG;YACD,IAAI,OAAO,CAAC,SAAS,CAAC,MAAM,IAAI,CAAC,EAAE;gBACjC,IAAI,eAAe;oBACjB,OAAO,KAAK,CAAC;gBACf,OAAO,CAAC,MAAM,GAAG,KAAK,CAAC;aACxB;YAED,kCAAkC;YAClC,IAAI,OAAO,CAAC,MAAM,EAAE;gBAClB,IAAI,CAAC,eAAe,EAAE;oBACpB,OAAO,CAAC,SAAS,CAAC,IAAI,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC;oBACrD,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK,EAAE,EAAG,4BAA4B;wBAC9D,OAAO,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;qBACzJ;iBACF;gBACD,IAAI,OAAO,CAAC,SAAS,CAAC,MAAM,IAAI,CAAC;oBAC/B,OAAO,CAAC,MAAM,GAAG,KAAK,CAAC,CAAC,kDAAkD;aAC7E;YAED,IAAI,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC;gBAC9B,OAAO,KAAK,CAAC;YAEf,oCAAoC;YACpC,IAAI,SAAS,KAAK,OAAO,CAAC,YAAY,EAAE;gBACtC,IAAI,OAAO,CAAC,YAAY,CAAC,YAAY;oBACnC,OAAO,CAAC,YAAY,GAAG,SAAS,CAAC;;oBAEjC,OAAO,CAAC,YAAY,CAAC,gBAAgB,EAAE,CAAC;aAC3C;YACD,IAAI,SAAS,KAAK,OAAO,CAAC,UAAU,EAAE;gBACpC,IAAI,OAAO,CAAC,UAAU,CAAC,YAAY;oBACjC,OAAO,CAAC,UAAU,GAAG,SAAS,CAAC;;oBAE/B,OAAO,CAAC,UAAU,CAAC,gBAAgB,EAAE,CAAC;aACzC;YAED,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;YAII;QACG,MAAM,CAAC,iCAAiC,CAAC,KAA2B,EAAE,YAAoB,EAAE,SAAmB;YACpH,IAAI,MAAM,GAAG,KAAK,aAAL,KAAK,uBAAL,KAAK,CAAE,KAAK,EAAE,CAAC;YAC5B,IAAI,SAAS,KAAK,MAAM,EAAE;gBACxB,MAAM,aAAa,GAAG,MAAM,CAAC,MAAM,GAAG,YAAY,CAAC;gBACnD,QAAQ,aAAa,EAAE;oBACrB,KAAK,CAAC,CAAC,CAAC,EAAI,yCAAyC;wBACnD,MAAM;qBACP;oBACD,KAAK,CAAC,CAAC,CAAK,0BAA0B;oBACtC,KAAK,CAAC,CAAC,CAAC,EAAI,0BAA0B;wBACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,aAAa,GAAG,CAAC,EAAE,EAAE,CAAC,EAAE;4BAC1C,MAAM,CAAC,GAAG,EAAE,CAAC;4BACb,MAAM,CAAC,KAAK,EAAE,CAAC;yBAChB;wBACD,MAAM;qBACP;oBACD,OAAO,CAAC,CAAC,EAAG,yCAAyC;wBACnD,MAAM,GAAG,SAAS,CAAC;wBACnB,MAAM;qBACP;iBACF;gBACD,IAAI,SAAS,IAAI,CAAC,IAAI,CAAC,cAAc,CAAC,MAAM,CAAC;oBAC3C,MAAM,GAAG,SAAS,CAAC;aACtB;YACD,OAAO,MAAM,CAAC;QAChB,CAAC;QAED;;;YAGI;QACG,MAAM,CAAC,iCAAiC,CAAC,MAA4B,EAAE,MAAgB,EAAE,MAAgB;YAC9G,MAAM,KAAK,GAAG,MAAM,aAAN,MAAM,uBAAN,MAAM,CAAE,KAAK,EAAE,CAAC;YAC9B,IAAI,SAAS,KAAK,KAAK,EAAE;gBACvB,MAAM,aAAa,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACrC,IAAI,MAAM,EAAE;oBACV,MAAM,KAAK,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC;oBAC/B,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,IAAI,aAAa,EAAE,KAAK,IAAI,CAAC,EAAE;wBACtD,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,CAAC,CAAE,oBAAoB;wBACxD,KAAK,CAAC,IAAI,CAAC,GAAG,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAK,0BAA0B;qBAC/D;iBACF;qBAAM;oBACL,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,aAAa,GAAG,CAAC,EAAE,EAAE,CAAC,EAAE;wBAC1C,KAAK,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;wBACnB,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;qBACjB;iBACF;aACF;YACD,OAAO,KAAK,CAAC;QACf,CAAC;QAED,sDAAsD;QAC/C,MAAM,CAAC,kBAAkB,CAAC,KAAgC;YAC/D,IAAI,SAAS,KAAK,KAAK,CAAC,KAAK,EAAE;gBAC7B,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,iCAAiC,CAAC,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,SAAS,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;gBACjG,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,iCAAiC,CAAC,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;aACvF;iBAAM;gBACL,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,gCAAgC,CAAC,YAAY,CAAC,iBAAiB,CAAC,KAAK,CAAC,SAAS,CAAC,EAAE,KAAK,CAAC,eAAe,EAAE,KAAK,CAAC,MAAM,CAAC,CAAC;gBAC1I,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,iCAAiC,CAAC,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;aACvF;QACH,CAAC;QAED;;WAEG;QACI,MAAM,CAAC,cAAc,CAAC,OAAoC;YAC/D,IAAI,CAAC,IAAI,CAAC,sBAAsB,CAAC,OAAO,CAAC,IAAI,CAAC,SAAS,KAAK,OAAO,CAAC,KAAK,CAAC;gBACxE,OAAO,SAAS,CAAC;YAEnB,MAAM,MAAM,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,CAAC;YACxC,MAAM,KAAK,GAAa,KAAK,CAAC,MAAM,CAAC,CAAC;YACtC,MAAM,IAAI,GAAa,KAAK,CAAC,MAAM,CAAC,CAAC;YACrC,MAAM,KAAK,GAAa,KAAK,CAAC,MAAM,CAAC,CAAC;YACtC,IAAI,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,OAAO,CAAC;gBAChD,OAAO,SAAS,CAAC;YAEnB,IAAI,KAAK,GAAyC,EAAE,CAAC;YACrD,IAAI,CAAC,OAAO,CAAC,MAAM,EAAE;gBACnB,MAAM,OAAO,GAAG,OAAO,CAAC,SAAS,CAAC,KAAK,EAAE,CAAC;gBAC1C,IAAI,CAAC,IAAI,CAAC,gBAAgB,CAAC,OAAO,EAAE,OAAO,CAAC;oBAC1C,OAAO,SAAS,CAAC;gBACnB,IAAI,OAAO,CAAC,MAAM,KAAK,MAAM,GAAG,CAAC;oBAC/B,OAAO,SAAS,CAAC,CAAC,+EAA+E;gBAEnG,iDAAiD;gBACjD,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;gBAC5C,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;gBACzC,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,WAAW,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,EAAE,IAAI,GAAG,MAAM,EAAE,EAAE,IAAI,EAAE;oBACvE,MAAM,CAAC,WAAW,EAAE,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC;oBACpC,MAAM,CAAC,WAAW,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC;oBACnC,MAAM,CAAC,WAAW,EAAE,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC;oBACpC,GAAG,CAAC,QAAQ,EAAE,CAAC,GAAG,OAAO,CAAC,IAAI,GAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACpC,GAAG,CAAC,QAAQ,EAAE,CAAC,GAAG,OAAO,CAAC,IAAI,GAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACpC,GAAG,CAAC,QAAQ,EAAE,CAAC,GAAG,OAAO,CAAC,IAAI,GAAC,CAAC,CAAC,CAAC,CAAC,CAAC;iBACrC;gBAED,MAAM,QAAQ,GAAG,YAAY,CAAC,4BAA4B,CAAC,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;gBACtF,IAAI,SAAS,KAAK,QAAQ;oBACxB,OAAO,SAAS,CAAC;gBAEnB,wCAAwC;gBACxC,KAAK,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,QAAQ,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;gBAClD,IAAI,MAAM,GAAG,CAAC,CAAC;gBACf,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACzC,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACzC,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACzC,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,QAAQ,CAAC,MAAM,GAAI;oBAC7C,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,QAAQ,CAAC,KAAK,EAAE,CAAC,CAAC;iBACrC;gBACD,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBACpE,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBACpE,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;aACrE;iBAAM,EAAE,SAAS;gBAChB,IAAI,SAAS,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,oBAAoB,CAAC,OAAO,CAAC,SAAS,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC,EAAE;oBAC5F,IAAI,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE;wBACpB,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,GAAG,EAAG,CAAC,CAAC,CAAE,8CAA8C;wBAC5E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,KAAK,GAAG,CAAC,EAAE,EAAE,CAAC;4BACxC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,uCAAuC;qBACxE;iBACF;aACF;YACD,OAAO,KAAK,CAAC;QACf,CAAC;KACF;IAxwBY,0BAAU,aAwwBtB,CAAA;AACH,CAAC,EAhxBgB,eAAe,KAAf,eAAe,QAgxB/B","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 Bspline\r\n */\r\n\r\n/* eslint-disable @typescript-eslint/naming-convention, no-empty, no-console*/\r\n\r\nimport { GrowableXYZArray } from \"../geometry3d/GrowableXYZArray\";\r\nimport { IndexedXYZCollection } from \"../geometry3d/IndexedXYZCollection\";\r\nimport { Geometry } from \"../Geometry\";\r\nimport { Point3d, Vector3d } from \"../geometry3d/Point3dVector3d\";\r\nimport { BandedSystem } from \"../numerics/BandedSystem\";\r\nimport { BSplineCurve3d } from \"./BSplineCurve\";\r\nimport { BSplineWrapMode, KnotVector } from \"./KnotVector\";\r\nimport { InterpolationCurve3dOptions, InterpolationCurve3dProps } from \"./InterpolationCurve3d\";\r\nimport { Point3dArray } from \"../geometry3d/PointHelpers\";\r\n\r\n/**\r\n * A class with static methods for creating B-spline curves.\r\n * @public\r\n */\r\nexport class BSplineCurveOps {\r\n  /**\r\n   * Greville knot algorithm: alternative to traditional c2 cubic algorithm allowing for any order, but no end conditions or periodicity.\r\n   * @param points pass-through points.\r\n   * @param order bspline order (1 more than degree)\r\n   */\r\n  public static createThroughPoints(points: IndexedXYZCollection | Point3d[], order: number): BSplineCurve3d | undefined {\r\n    const numPoints = points.length;\r\n    if (order > numPoints || order < 2)\r\n      return undefined;\r\n    const degree = order - 1;\r\n    const bw = 1 + 2 * degree;    // probably less than that . . just zeros at fringe.\r\n    const matrix = new Float64Array(bw * numPoints);\r\n    const basisFunctions = new Float64Array(order);\r\n    const rhs = new GrowableXYZArray();\r\n    const knots = KnotVector.createUniformClamped(numPoints, order - 1, 0.0, 1.0);\r\n    const xyz = Point3d.create();\r\n    for (let basePointIndex = 0; basePointIndex < numPoints; basePointIndex++) {\r\n      const u = knots.grevilleKnot(basePointIndex);\r\n      const spanIndex = knots.knotToLeftKnotIndex(u);\r\n      knots.evaluateBasisFunctions(spanIndex, u, basisFunctions);\r\n      // puzzlement .. how do the max points shift within the order spots?\r\n      let maxIndex = 0;\r\n      for (let i = 1; i < order; i++)\r\n        if (basisFunctions[i] > basisFunctions[maxIndex])\r\n          maxIndex = i;\r\n      const basisFunctionStartWithinRow = degree - maxIndex;\r\n      const rowStart = basePointIndex * bw;\r\n      for (let i = 0; i < order; i++) {\r\n        const realColumn = basePointIndex - degree + basisFunctionStartWithinRow + i;\r\n        if (rowStart + realColumn >= 0 && realColumn < numPoints)\r\n          matrix[rowStart + basisFunctionStartWithinRow + i] = basisFunctions[i];\r\n      }\r\n      if (points instanceof IndexedXYZCollection) {\r\n        rhs.push(points.getPoint3dAtUncheckedPointIndex(basePointIndex, xyz));\r\n      } else {\r\n        rhs.push(points[basePointIndex].clone());\r\n      }\r\n    }\r\n    const poles = BandedSystem.solveBandedSystemMultipleRHS(numPoints, bw, matrix, 3, rhs.float64Data());\r\n    if (poles) {\r\n      return BSplineCurve3d.create(poles, knots.knots, order);\r\n    }\r\n    return undefined;\r\n  }\r\n\r\n  /**\r\n   * Construct BSplineCurve3d that fit points using the C2 cubic algorithm.\r\n   * @param options curve definition, unmodified\r\n   */\r\n  public static createThroughPointsC2Cubic(options: InterpolationCurve3dOptions): BSplineCurve3d | undefined {\r\n    // Work on a copy rather than installing computed knots/tangents that could become stale when fit points change\r\n    // Knots/tangents that come in, however, are used without recomputation.\r\n    const validatedOptions = options.clone();\r\n    if (!this.C2CubicFit.validateOptions(validatedOptions))\r\n      return undefined;\r\n\r\n    const poles = this.C2CubicFit.constructPoles(validatedOptions);\r\n    if (undefined === poles)\r\n      return undefined;\r\n\r\n    const fullKnots = this.C2CubicFit.convertFitParamsToCubicKnotVector(validatedOptions.knots, validatedOptions.closed);\r\n    if (undefined === fullKnots)\r\n      return undefined;\r\n\r\n    const interpolant = BSplineCurve3d.create(poles, fullKnots, validatedOptions.order);\r\n\r\n    if (validatedOptions.closed)\r\n      interpolant?.setWrappable(BSplineWrapMode.OpenByAddingControlPoints);\r\n\r\n    return interpolant;\r\n  }\r\n}\r\n\r\n/**\r\n * Namespace for collecting curve fit API\r\n * @public\r\n */\r\nexport namespace BSplineCurveOps {\r\n  /**\r\n   * A helper class for creating C2 cubic fit curves.\r\n   * Knots herein are understood to be *interior* knots (including the start/end knot),\r\n   * so that there is one knot per fit point. In other words, the knots are fit parameters.\r\n   * @private\r\n   */\r\n  export class C2CubicFit {\r\n    /** Transform fit parameters to span [0,1]\r\n     * @param knots fit parameters, normalized in place\r\n     */\r\n    private static normalizeKnots(knots: number[] | undefined): boolean {\r\n      if (undefined === knots || knots.length < 2) {\r\n        knots = undefined;\r\n        return false;\r\n      }\r\n      const myKnots = KnotVector.create(knots, 1, false);\r\n      if (myKnots.knotLength01 < KnotVector.knotTolerance) {\r\n        knots = undefined;\r\n        return false;\r\n      }\r\n      for (let i = 0; i < myKnots.numSpans; ++i) {\r\n        knots[i] = myKnots.spanFractionToFraction(i, 0.0);\r\n      }\r\n      knots[knots.length - 1] = 1.0;\r\n      return true;\r\n    }\r\n\r\n    /** Compute chord-length fit parameters for C2 cubic fit algorithm */\r\n    private static constructChordLengthParameters(fitPoints: Point3d[]): number[] | undefined {\r\n      if (fitPoints.length < 2)\r\n        return undefined;\r\n      const params: number[] = [0.0];\r\n      for (let i = 1; i < fitPoints.length; ++i)\r\n        params[i] = params[i - 1] + fitPoints[i].distance(fitPoints[i - 1]);\r\n      if (!this.normalizeKnots(params))\r\n        return undefined;\r\n      return params;\r\n    }\r\n\r\n    /** Compute uniform fit parameters for C2 cubic fit algorithm */\r\n    private static constructUniformParameters(numParams: number): number[] | undefined {\r\n      if (numParams < 2)\r\n        return undefined;\r\n      const knots = KnotVector.createUniformClamped(numParams + 2, 3, 0.0, 1.0);\r\n      const params: number[] = [];\r\n      for (let i = knots.leftKnotIndex; i <= knots.rightKnotIndex; ++i)\r\n        params.push(knots.knots[i]);\r\n      return params;\r\n    }\r\n\r\n    /** Remove duplicate fit points, and their given knots in parallel */\r\n    private static removeDuplicateFitPoints(options: InterpolationCurve3dOptions): boolean {\r\n      if (undefined !== options.knots && options.knots.length !== options.fitPoints.length)\r\n        options.knots = undefined;\r\n\r\n      // get indices of duplicate points to be removed\r\n      const newPts = GrowableXYZArray.create(options.fitPoints);\r\n      const indices = newPts.findOrderedDuplicates();\r\n      newPts.clear();\r\n\r\n      // remove duplicate fit points\r\n      for (let iRead = 0, iIndex = 0; iRead < options.fitPoints.length; ++iRead) {\r\n        if (iRead === indices[iIndex])\r\n          ++iIndex; // skip the duplicate\r\n        else\r\n          newPts.push(options.fitPoints[iRead].clone());\r\n      }\r\n      options.fitPoints = newPts.getPoint3dArray();\r\n\r\n      // remove params corresponding to removed fit points\r\n      if (undefined !== options.knots) {\r\n        const newKnots: number[] = [];\r\n        for (let iRead = 0, iIndex = 0; iRead < options.knots.length; ++iRead) {\r\n          if (iRead === indices[iIndex])\r\n            ++iIndex; // skip\r\n          else\r\n            newKnots.push(options.knots[iRead]);\r\n        }\r\n        options.knots = newKnots.slice();\r\n      }\r\n      return true;\r\n    }\r\n\r\n    /** Construct fit parameters for the c2 cubic fit algorithm.\r\n     * @param fitPoints validated fit points (should not contain duplicates)\r\n     * @param isChordLength whether knots are computed using distances between successive fit points\r\n     * @param closed whether curve is periodically defined\r\n     * @return fit parameters, one per fit point\r\n    */\r\n     public static constructFitParametersFromPoints(fitPoints: Point3d[], isChordLength: number | undefined, closed: boolean | undefined): number[] | undefined {\r\n      let params: number[] | undefined;\r\n      if (isChordLength || !closed)\r\n        params = this.constructChordLengthParameters(fitPoints);\r\n      if (undefined === params)\r\n        params = this.constructUniformParameters(fitPoints.length);\r\n      return params;\r\n    }\r\n\r\n    /** Construct fit parameters for the c2 cubic fit algorithm, if they are missing.\r\n     * @param options validated as per validateOptions, possibly modified\r\n     * @return whether fit parameters are valid\r\n     */\r\n     public static constructFitParameters(options: InterpolationCurve3dOptions): boolean {\r\n      if (undefined === options.knots)\r\n        options.knots = this.constructFitParametersFromPoints(options.fitPoints, options.isChordLenKnots, options.closed);\r\n      return options.knots?.length === options.fitPoints.length;\r\n    }\r\n\r\n    /** Compute a row of the tridiagonal system matrix from Farin.\r\n     * @param alpha sub-diagonal, length = # fit points\r\n     * @param beta diagonal, length = # fit points\r\n     * @param gamma super-diagonal, length = # fit points\r\n     * @param index 0-based row index to set\r\n     */\r\n    private static computeAlphaBetaGamma(alpha: number[], beta: number[], gamma: number[], index: number, deltaIPlus1: number, deltaI: number, deltaIMinus1: number, deltaIMinus2: number) {\r\n      let denomReciprocal = 1.0 / (deltaIMinus2 + deltaIMinus1 + deltaI);\r\n      alpha[index] = deltaI * deltaI * denomReciprocal;\r\n      beta[index] = deltaI * (deltaIMinus2 + deltaIMinus1) * denomReciprocal;\r\n\r\n      denomReciprocal = 1.0 / (deltaIMinus1 + deltaI + deltaIPlus1);\r\n      beta[index] += deltaIMinus1 * (deltaI + deltaIPlus1) * denomReciprocal;\r\n      gamma[index] = deltaIMinus1 * deltaIMinus1 * denomReciprocal;\r\n\r\n      denomReciprocal = 1.0 / (deltaIMinus1 + deltaI);\r\n      alpha[index] *= denomReciprocal;\r\n      beta[index] *= denomReciprocal;\r\n      gamma[index] *= denomReciprocal;\r\n    }\r\n\r\n    /** Setup tridiagonal system for 2 fit points\r\n     * @param alpha sub-diagonal, length = 2\r\n     * @param beta diagonal, length = 2\r\n     * @param gamma super-diagonal, length = 2\r\n     */\r\n    private static setUpSystem2Points(alpha: number[], beta: number[], gamma: number[]): boolean {\r\n      if (alpha.length !== 2 || beta.length !== 2 || gamma.length !== 2)\r\n        return false;\r\n\r\n      // identity matrix\r\n      alpha[0] = alpha[1] = gamma[0] = gamma[1] = 0.0;\r\n      beta[0] = beta[1] = 1.0;\r\n      return true;\r\n    }\r\n\r\n    /** Setup tridiagonal system for 3 fit points\r\n     * @param alpha sub-diagonal, length = 3\r\n     * @param beta diagonal, length = 3\r\n     * @param gamma super-diagonal, length = 3\r\n     * @param options validated as per validateOptions, unmodified\r\n     * @param useNaturalStartTangent whether to bake the natural end condition into the first row\r\n     * @param useNaturalEndTangent whether to bake the natural end condition into the last row\r\n     */\r\n     private static setUpSystem3Points(alpha: number[], beta: number[], gamma: number[], options: InterpolationCurve3dOptions, useNaturalStartTangent: boolean, useNaturalEndTangent: boolean): boolean {\r\n      if (undefined === options.knots)\r\n        return false;\r\n      if (alpha.length !== 3 || beta.length !== 3 || gamma.length !== 3)\r\n        return false;\r\n      if (options.knots.length !== 3 || options.fitPoints.length !== 3)\r\n        return false;\r\n\r\n      let deltaIPlus1 = 0, deltaI = 0, deltaIMinus1 = 0, deltaIMinus2 = 0, sum = 0, sumReciprocal = 0;\r\n\r\n      // first row\r\n      if (useNaturalStartTangent) {\r\n        alpha[0] = 0.0;\r\n        deltaI = options.knots[1] - options.knots[0];\r\n        deltaIPlus1 = options.knots[2] - options.knots[1];\r\n        sum = deltaI + deltaIPlus1;\r\n        sumReciprocal = 1.0 / sum;\r\n        beta[0] = (deltaI + sum) * sumReciprocal;\r\n        gamma[0] = -deltaI * sumReciprocal;\r\n      } else {\r\n        alpha[0] = gamma[0] = 0.0;\r\n        beta[0] = 1.0;\r\n      }\r\n\r\n      // middle row\r\n      deltaIMinus1 = options.knots[1] - options.knots[0];\r\n      deltaI = options.knots[2] - options.knots[1];\r\n      sumReciprocal = 1.0 / (deltaIMinus1 + deltaI);\r\n      sumReciprocal *= sumReciprocal;\r\n      alpha[1] = deltaI * deltaI * sumReciprocal;\r\n      beta[1]  = 2.0 * (deltaI * deltaIMinus1) * sumReciprocal;\r\n      gamma[1] = deltaIMinus1 * deltaIMinus1 * sumReciprocal;\r\n\r\n      // last row\r\n      if (useNaturalEndTangent) {\r\n        deltaIMinus1 = options.knots[2] - options.knots[1];\r\n        deltaIMinus2 = options.knots[1] - options.knots[0];\r\n        sum = deltaIMinus2 + deltaIMinus1;\r\n        sumReciprocal = 1.0 / sum;\r\n        alpha[2] = -deltaIMinus1 * sumReciprocal;\r\n        beta[2] = (deltaIMinus1 + sum) * sumReciprocal;\r\n        gamma[2] = 0.0;\r\n      } else {\r\n        alpha[2] = gamma[2] = 0.0;\r\n        beta[2] = 1.0;\r\n      }\r\n      return true;\r\n    }\r\n\r\n    /** Setup tridiagonal system for 4 fit points or more\r\n     * @param alpha sub-diagonal, length = # fit points\r\n     * @param beta diagonal, length = # fit points\r\n     * @param gamma super-diagonal, length = # fit points\r\n     * @param options validated as per validateOptions, unmodified\r\n     * @param useNaturalStartTangent whether to bake the natural end condition into the first row\r\n     * @param useNaturalEndTangent whether to bake the natural end condition into the last row\r\n     */\r\n    private static setUpSystem4PointsOrMore(alpha: number[], beta: number[], gamma: number[], options: InterpolationCurve3dOptions, useNaturalStartTangent: boolean, useNaturalEndTangent: boolean): boolean {\r\n      if (undefined === options.knots)\r\n        return false;\r\n      if (alpha.length !== beta.length || alpha.length !== gamma.length || alpha.length !== options.knots.length)\r\n        return false;\r\n      if (options.knots.length !== options.fitPoints.length)\r\n        return false;\r\n\r\n      const numIntervals = options.fitPoints.length - 1;\r\n      const numIntervalsMinus1 = numIntervals - 1;\r\n      let deltaIPlus1 = 0, deltaI = 0, deltaIMinus1 = 0, deltaIMinus2 = 0, sum = 0, sumReciprocal = 0;\r\n\r\n      if (options.closed) {\r\n        // first row\r\n        deltaI = options.knots[1] - options.knots[0];\r\n        deltaIMinus2 = options.knots[numIntervalsMinus1] - options.knots[numIntervalsMinus1 - 1];\r\n        deltaIMinus1 = options.knots[numIntervalsMinus1 + 1] - options.knots[numIntervalsMinus1];\r\n        deltaIPlus1 = options.knots[2] - options.knots[1];\r\n        this.computeAlphaBetaGamma(alpha, beta, gamma, 0, deltaIPlus1, deltaI, deltaIMinus1, deltaIMinus2);\r\n\r\n        // second row\r\n        deltaIMinus2 = deltaIMinus1;\r\n        deltaIMinus1 = deltaI;\r\n        deltaI = options.knots[2] - options.knots[1];\r\n        deltaIPlus1 = options.knots[3] - options.knots[2];\r\n        this.computeAlphaBetaGamma(alpha, beta, gamma, 1, deltaIPlus1, deltaI, deltaIMinus1, deltaIMinus2);\r\n\r\n        // last row; there's one less equation than open case\r\n        deltaIPlus1 = deltaIMinus1;\r\n        deltaI = options.knots[numIntervalsMinus1 + 1] - options.knots[numIntervalsMinus1];\r\n        deltaIMinus2 = options.knots[numIntervalsMinus1 - 1] - options.knots[numIntervalsMinus1 - 2];\r\n        deltaIMinus1 = options.knots[numIntervalsMinus1] - options.knots[numIntervalsMinus1 - 1];\r\n        this.computeAlphaBetaGamma(alpha, beta, gamma, numIntervalsMinus1, deltaIPlus1, deltaI, deltaIMinus1, deltaIMinus2);\r\n      } else { // open\r\n        // first row\r\n        if (useNaturalStartTangent) {\r\n          alpha[0] = 0.0;\r\n          deltaI = options.knots[1] - options.knots[0];\r\n          deltaIPlus1 = options.knots[2] - options.knots[1];\r\n          sum = deltaI + deltaIPlus1;\r\n          sumReciprocal = 1.0 / sum;\r\n          beta[0] = (deltaI + sum) * sumReciprocal;\r\n          gamma[0] = -deltaI * sumReciprocal;\r\n        } else {\r\n          alpha[0] = gamma[0] = 0.0;\r\n          beta[0] = 1.0;\r\n        }\r\n\r\n        // second row\r\n        deltaI = options.knots[2] - options.knots[1];\r\n        deltaIMinus1 = options.knots[1] - options.knots[0];\r\n        deltaIMinus2 = 0.0;\r\n        deltaIPlus1 = options.knots[3] - options.knots[2];\r\n        this.computeAlphaBetaGamma(alpha, beta, gamma, 1, deltaIPlus1, deltaI, deltaIMinus1, deltaIMinus2);\r\n\r\n        // penultimate row\r\n        deltaI = options.knots[numIntervalsMinus1 + 1] - options.knots[numIntervalsMinus1];\r\n        deltaIMinus1 = options.knots[numIntervalsMinus1] - options.knots[numIntervalsMinus1 - 1];\r\n        deltaIMinus2 = options.knots[numIntervalsMinus1 - 1] - options.knots[numIntervalsMinus1 - 2];\r\n        deltaIPlus1 = 0.0;\r\n        this.computeAlphaBetaGamma(alpha, beta, gamma, numIntervalsMinus1, deltaIPlus1, deltaI, deltaIMinus1, deltaIMinus2);\r\n\r\n        // last row\r\n        if (useNaturalEndTangent) {\r\n          deltaIMinus1 = options.knots[numIntervals] - options.knots[numIntervals - 1];\r\n          deltaIMinus2 = options.knots[numIntervals - 1] - options.knots[numIntervals - 2];\r\n          sum = deltaIMinus2 + deltaIMinus1;\r\n          sumReciprocal = 1.0 / sum;\r\n          alpha[numIntervals] = -deltaIMinus1 * sumReciprocal;\r\n          beta[numIntervals] = (deltaIMinus1 + sum) * sumReciprocal;\r\n          gamma[numIntervals] = 0.0;\r\n        } else {\r\n          alpha[numIntervals] = gamma[numIntervals] = 0.0;\r\n          beta[numIntervals] = 1.0;\r\n        }\r\n      }\r\n\r\n      // middle rows\r\n      for (let i = 2; i < numIntervalsMinus1; ++i) {\r\n        deltaI = options.knots[i + 1] - options.knots[i];\r\n        deltaIMinus2 = options.knots[i - 1] - options.knots[i - 2];\r\n        deltaIMinus1 = options.knots[i] - options.knots[i - 1];\r\n        deltaIPlus1 = options.knots[i + 2] - options.knots[i + 1];\r\n        this.computeAlphaBetaGamma(alpha, beta, gamma, i, deltaIPlus1, deltaI, deltaIMinus1, deltaIMinus2);\r\n      }\r\n      return true;\r\n    }\r\n\r\n    /** Setup tridiagonal system\r\n     * @param alpha sub-diagonal, length = # fitPoints\r\n     * @param beta diagonal, length = # fitPoints\r\n     * @param gamma super-diagonal, length = # fitPoints\r\n     * @param options validated as per validateOptions, unmodified\r\n     */\r\n    private static setUpSystem(alpha: number[], beta: number[], gamma: number[], options: InterpolationCurve3dOptions): boolean {\r\n      let useNaturalStartTangent = false;\r\n      let useNaturalEndTangent = false;\r\n      if (options.isNaturalTangents && !options.closed) {\r\n        useNaturalStartTangent = (undefined === options.startTangent);\r\n        useNaturalEndTangent = (undefined === options.endTangent);\r\n      }\r\n\r\n      let succeeded = false;\r\n      if (2 === options.fitPoints.length)\r\n        succeeded = this.setUpSystem2Points(alpha, beta, gamma);\r\n      else if (3 === options.fitPoints.length)\r\n        succeeded = this.setUpSystem3Points(alpha, beta, gamma, options, useNaturalStartTangent, useNaturalEndTangent);\r\n      else if (4 <= options.fitPoints.length)\r\n        succeeded = this.setUpSystem4PointsOrMore(alpha, beta, gamma, options, useNaturalStartTangent, useNaturalEndTangent);\r\n\r\n      return succeeded;\r\n    }\r\n\r\n    /** Set the Bessel end condition for the linear system.\r\n     * @param dataPts array whose middle is the system rhs (augmented with first/last fitPoint at beginning/end);\r\n     *                2nd or penultimate point is set by this function.\r\n     * @param options validated as per validateOptions, unmodified\r\n     * @param atStart whether end condition is for start of curve (false: end of curve)\r\n     */\r\n    private static setBesselEndCondition(dataPts: Point3d[], options: InterpolationCurve3dOptions, atStart: boolean): boolean {\r\n      if (dataPts.length !== options.fitPoints.length + 2)\r\n        return false;\r\n      if (undefined === options.knots)\r\n        return false;\r\n\r\n      const scale = 1.0/3.0;\r\n      const numIntervals = options.fitPoints.length - 1;\r\n\r\n      if (1 === numIntervals) { // linear Bezier\r\n        if (atStart)\r\n          dataPts[0].interpolate(scale, dataPts[3], dataPts[1]);\r\n        else\r\n          dataPts[3].interpolate(scale, dataPts[0], dataPts[2]);\r\n        return true;\r\n      }\r\n\r\n      if (2 === numIntervals) {\r\n        const alpha = (options.knots[2] - options.knots[1]) / (options.knots[2] - options.knots[0]);\r\n        const beta = 1.0 - alpha;\r\n        const temp = dataPts[2].plus2Scaled(dataPts[0], -alpha * alpha, dataPts[4], -beta * beta);\r\n        if (atStart)\r\n          Point3d.createAdd2Scaled(temp, 1.0 / (2.0 * alpha), dataPts[0], alpha).interpolate(scale, dataPts[0], dataPts[1]);\r\n        else\r\n          Point3d.createAdd2Scaled(temp, 1.0 / (2.0 * beta), dataPts[4], beta).interpolate(scale, dataPts[4], dataPts[3]);\r\n        return true;\r\n      }\r\n\r\n      // numIntervals > 2\r\n      if (atStart) {\r\n          const alpha = (options.knots[2] - options.knots[1]) / (options.knots[2] - options.knots[0]);\r\n          const beta = 1.0 - alpha;\r\n          const temp = dataPts[2].plus2Scaled(dataPts[0], -alpha * alpha, dataPts[3], -beta * beta);\r\n          Point3d.createAdd2Scaled(temp, 1.0 / (2.0 * alpha), dataPts[0], alpha).interpolate(scale, dataPts[0], dataPts[1]);\r\n      } else {\r\n          const alpha = (options.knots[numIntervals] - options.knots[numIntervals - 1]) / (options.knots[numIntervals] - options.knots[numIntervals - 2]);\r\n          const beta = 1.0 - alpha;\r\n          const temp = dataPts[numIntervals].plus2Scaled(dataPts[numIntervals - 1], -alpha * alpha, dataPts[numIntervals + 2], -beta * beta);\r\n          Point3d.createAdd2Scaled(temp, 1.0 / (2.0 * beta), dataPts[numIntervals + 2], beta).interpolate(scale, dataPts[numIntervals + 2], dataPts[numIntervals + 1]);\r\n      }\r\n      return true;\r\n    }\r\n\r\n    /** Set the natural end condition for the linear system.\r\n     *  This is the end condition used by ADSK for fit-splines with a given zero tangent.\r\n     * @param dataPts array whose middle is the system rhs (augmented with first/last fitPoint at beginning/end);\r\n     *                2nd or penultimate point is set by this function.\r\n     * @param options validated as per validateOptions, unmodified\r\n     * @param atStart whether end condition is for start of curve (false: end of curve)\r\n     */\r\n    private static setNaturalEndCondition(dataPts: Point3d[], options: InterpolationCurve3dOptions, atStart: boolean): boolean {\r\n      if (dataPts.length !== options.fitPoints.length + 2)\r\n        return false;\r\n\r\n      const numIntervals = options.fitPoints.length - 1;\r\n      if (1 === numIntervals)\r\n        return this.setBesselEndCondition(dataPts, options, atStart);\r\n\r\n      if (atStart)\r\n        dataPts[1] = dataPts[0];\r\n      else\r\n        dataPts[dataPts.length - 2] = dataPts[dataPts.length - 1];\r\n      return true;\r\n    }\r\n\r\n    /** Set the end condition for the linear system to the given tangent, scaled by chord length.\r\n     *  This is the end condition used by ADSK for fit-splines with a given nonzero tangent.\r\n     * @param dataPts array whose middle is the system rhs (augmented with first/last fitPoint at beginning/end);\r\n     *                2nd or penultimate point is set by this function.\r\n     * @param options validated as per validateOptions, unmodified\r\n     * @param atStart whether end condition is for start of curve (false: end of curve)\r\n     */\r\n    private static setChordLengthScaledEndCondition(dataPts: Point3d[], options: InterpolationCurve3dOptions, atStart: boolean): boolean {\r\n      if (dataPts.length !== options.fitPoints.length + 2)\r\n        return false;\r\n\r\n      const tangent = atStart ? options.startTangent : options.endTangent;\r\n      if (undefined === tangent)\r\n        return false;\r\n\r\n      let iExt = 0; // index of first/last fitPoint\r\n      let iSet = 0; // index of 2nd/penultimate Bezier point to set (determines start/end tangent of the curve)\r\n      let iInt = 0; // index of 2nd/penultimate fitPoint\r\n\r\n      const numIntervals = options.fitPoints.length - 1;\r\n      if (1 === numIntervals) { // no interior fit points\r\n        if (atStart) {\r\n          iExt = 0;\r\n          iSet = 1;\r\n          iInt = 3;\r\n        } else {\r\n          iExt = 3;\r\n          iSet = 2;\r\n          iInt = 0;\r\n        }\r\n      } else {\r\n        if (atStart) {\r\n          iExt = 0;\r\n          iSet = 1;\r\n          iInt = 2;\r\n        } else {\r\n          iExt = numIntervals + 2;\r\n          iSet = numIntervals + 1;\r\n          iInt = numIntervals;\r\n        }\r\n      }\r\n\r\n      // NOTE: tangent points INTO curve\r\n      const chordLength = dataPts[iInt].distance(dataPts[iExt]);\r\n      dataPts[iExt].plusScaled(tangent, chordLength / 3.0, dataPts[iSet]);\r\n      return true;\r\n    }\r\n\r\n    /** Set the end condition for the linear system to the given tangent, scaled by bessel length.\r\n     * @param dataPts array whose middle is the system rhs (augmented with first/last fitPoint at beginning/end);\r\n     *                2nd or penultimate point is set by this function.\r\n     * @param options validated as per validateOptions, unmodified\r\n     * @param atStart whether end condition is for start of curve (false: end of curve)\r\n     */\r\n    private static setBesselLengthScaledEndCondition(dataPts: Point3d[], options: InterpolationCurve3dOptions, atStart: boolean): boolean {\r\n      if (dataPts.length !== options.fitPoints.length + 2)\r\n        return false;\r\n\r\n      const tangent = atStart ? options.startTangent : options.endTangent;\r\n      if (undefined === tangent)\r\n        return false;\r\n\r\n      // temporarily set bessel end condition\r\n      if (!this.setBesselEndCondition(dataPts, options, atStart))\r\n        return false;\r\n\r\n      const numIntervals = options.fitPoints.length - 1;\r\n      const iExt = atStart ? 0 : numIntervals + 2; // index of first/last fitPoint\r\n      const iSet = atStart ? 1 : numIntervals + 1; // index of 2nd/penultimate Bezier point to set (determines start/end tangent of the curve)\r\n\r\n      // reset end condition with our tangent, but scaled to the bessel tangent's length\r\n      dataPts[iExt].plusScaled(tangent, dataPts[iExt].distance(dataPts[iSet]), dataPts[iSet]);\r\n      return true;\r\n    }\r\n\r\n    /** Set the end condition for a physically closed (non-periodic) interpolant.\r\n     * @param dataPts array whose middle is the system rhs (augmented with first/last fitPoint at beginning/end);\r\n     *                2nd or penultimate point is set by this function.\r\n     * @param options validated as per validateOptions, unmodified\r\n     */\r\n    private static setPhysicallyClosedEndCondition(dataPts: Point3d[], options: InterpolationCurve3dOptions): boolean {\r\n      const numIntervals = options.fitPoints.length - 1;\r\n      if (!options.isColinearTangents\r\n          || numIntervals <= 2\r\n          || (undefined !== options.startTangent && undefined !== options.endTangent)\r\n          || options.isNaturalTangents\r\n          || !dataPts[0].isAlmostEqual(dataPts[numIntervals + 2])) {\r\n        return true;\r\n        }\r\n      // force parallel start/end tangents, using chord length scale for undefined tangents\r\n      if (undefined !== options.startTangent) { // start tangent is supplied; compute a parallel end tangent\r\n        const outwardStartTangent = Vector3d.createStartEnd(dataPts[1], dataPts[0]).normalize();\r\n        if (undefined !== outwardStartTangent) {\r\n          const endTangentMag = dataPts[numIntervals + 2].distance(dataPts[numIntervals + 1]);\r\n          dataPts[numIntervals + 2].plusScaled(outwardStartTangent, endTangentMag, dataPts[numIntervals + 1]);\r\n        }\r\n      } else if (undefined !== options.endTangent) {  // end tangent is supplied; compute a parallel start tangent\r\n        const outwardEndTangent = Vector3d.createStartEnd(dataPts[numIntervals + 1], dataPts[numIntervals + 2]).normalize();\r\n        if (undefined !== outwardEndTangent) {\r\n          const startTangentMag = dataPts[0].distance(dataPts[1]);\r\n          dataPts[0].plusScaled(outwardEndTangent, startTangentMag, dataPts[1]);\r\n        }\r\n      } else {  // neither tangent is supplied, compute both along same vector\r\n        const commonTangent = Vector3d.createStartEnd(dataPts[numIntervals + 1], dataPts[1]).normalize();\r\n        if (undefined !== commonTangent) {\r\n          const startTangentMag = dataPts[0].distance(dataPts[1]);\r\n          dataPts[0].plusScaled(commonTangent, startTangentMag, dataPts[1]);\r\n          const endTangentMag = dataPts[numIntervals + 2].distance(dataPts[numIntervals + 1]);\r\n          dataPts[numIntervals + 2].plusScaled(commonTangent, -endTangentMag, dataPts[numIntervals + 1]);\r\n        }\r\n      }\r\n      return true;\r\n    }\r\n\r\n    /** Set end conditions for the linear system to solve for the poles of the open interpolant, as per Farin 3e/4e.\r\n     * @param dataPts array whose interior is the system rhs and whose first/last entries are the first/last fitPoints;\r\n     *                points are inserted to become the 2nd and penultimate dataPts, the first/last rows of the system rhs.\r\n     * @param options validated as per validateOptions, unmodified\r\n     */\r\n    private static setEndConditions(dataPts: Point3d[], options: InterpolationCurve3dOptions): boolean {\r\n      if (dataPts.length !== options.fitPoints.length)\r\n        return false;\r\n\r\n      // insert dummy points to be computed below\r\n      const dummy0 = Point3d.createZero();\r\n      const dummy1 = Point3d.createZero();\r\n      dataPts.splice(1, 0, dummy0);\r\n      dataPts.splice(dataPts.length - 1, 0, dummy1);\r\n\r\n      let succeeded = false;\r\n      if (undefined === options.startTangent) {\r\n        if (options.isNaturalTangents)\r\n          succeeded = this.setNaturalEndCondition(dataPts, options, true);\r\n        else\r\n          succeeded = this.setBesselEndCondition(dataPts, options, true);\r\n      } else { // scale startTangent\r\n        if (options.isChordLenTangents)\r\n          succeeded = this.setChordLengthScaledEndCondition(dataPts, options, true);\r\n        else\r\n          succeeded = this.setBesselLengthScaledEndCondition(dataPts, options, true);\r\n      }\r\n\r\n      if (undefined === options.endTangent) {\r\n        if (options.isNaturalTangents)\r\n          succeeded = this.setNaturalEndCondition(dataPts, options, false);\r\n        else\r\n          succeeded = this.setBesselEndCondition(dataPts, options, false);\r\n      } else { // scale endTangent\r\n        if (options.isChordLenTangents)\r\n          succeeded = this.setChordLengthScaledEndCondition(dataPts, options, false);\r\n        else\r\n          succeeded = this.setBesselLengthScaledEndCondition(dataPts, options, false);\r\n      }\r\n\r\n      if (succeeded)\r\n        succeeded = this.setPhysicallyClosedEndCondition(dataPts, options);\r\n\r\n      return succeeded;\r\n    }\r\n\r\n    /** Solve the near tridiagonal system for a periodic C2 cubic interpolant.\r\n     * Alpha, beta, gamma are computed by setUpSystem, have same length as fitPts, and are overwritten.\r\n     */\r\n    private static solveNearTridiagonal(fitPts: Point3d[], alpha: number[], beta: number[], gamma: number[]): Point3d[] | undefined {\r\n      if (alpha.length !== beta.length || alpha.length !== gamma.length || alpha.length !== fitPts.length)\r\n        return undefined;\r\n      const poles: Point3d[] = [];\r\n      const numIntervals = fitPts.length - 1;\r\n      const leftPts = fitPts.slice(0, -1);  // last fitPt is ignored\r\n      let tmp: number | undefined = 0.0;\r\n\r\n      // first forward substitution\r\n      for (let i = 1; i < numIntervals; ++i) {\r\n        if (undefined === (tmp = Geometry.conditionalDivideFraction(- alpha[i], beta[i - 1])))\r\n          return undefined;\r\n        beta[i] += tmp * gamma[i - 1];\r\n        alpha[i] = tmp * alpha[i - 1];\r\n        leftPts[i].addScaledInPlace(leftPts[i - 1], tmp);\r\n      }\r\n\r\n      // first backward substitution\r\n      if (undefined === (tmp = Geometry.conditionalDivideFraction(1.0, beta[numIntervals - 1] + alpha[numIntervals - 1])))\r\n        return undefined;\r\n      gamma[numIntervals - 1] *= tmp;\r\n      leftPts[numIntervals - 1].scaleInPlace(tmp);\r\n      for (let i = numIntervals - 2; i >= 0; --i) {\r\n        if (undefined === (tmp = Geometry.conditionalDivideFraction(1.0, beta[i])))\r\n          return undefined;\r\n        Point3d.createScale(leftPts[i].plus2Scaled(leftPts[i + 1], - gamma[i], leftPts[numIntervals - 1], - alpha[i]), tmp, leftPts[i]);\r\n        gamma[i] = - (gamma[i] * gamma[i + 1] + alpha[i] * gamma[numIntervals - 1]) * tmp;\r\n      }\r\n\r\n      // second forward substitution\r\n      if (undefined === (tmp = Geometry.conditionalDivideFraction(1.0, 1.0 + gamma[0])))\r\n        return undefined;\r\n      poles.push(Point3d.createScale(leftPts[0], tmp));\r\n      for (let i = 1; i < numIntervals; ++i) {\r\n        poles.push(leftPts[i].plusScaled(poles[0], - gamma[i]));\r\n      }\r\n      return poles;\r\n    }\r\n\r\n    /** Adjust options by correcting invalid combinations\r\n     * @param options curve definition, possibly modified\r\n     */\r\n    public static validateOptions(options: InterpolationCurve3dOptions): boolean {\r\n      options.order = 4;\r\n\r\n      // remove relevant exterior knots so knots and fit points align *before* we start compressing fit points.\r\n      options.knots = this.convertCubicKnotVectorToFitParams(options.knots, options.fitPoints.length, true);\r\n\r\n        // compress out duplicate fit points (and their corresponding knots)\r\n      if (!this.removeDuplicateFitPoints(options))\r\n        return false;\r\n\r\n      // if only 2 unique points, then must create open curve\r\n      let hasClosurePoint = options.fitPoints[0].isAlmostEqual(options.fitPoints[options.fitPoints.length - 1]);\r\n      if (3 === options.fitPoints.length && hasClosurePoint) {\r\n        options.fitPoints.pop();\r\n        if (undefined !== options.knots)\r\n          options.knots.pop();\r\n        hasClosurePoint = options.fitPoints[0].isAlmostEqual(options.fitPoints[options.fitPoints.length - 1]);\r\n      }\r\n      if (options.fitPoints.length <= 2) {\r\n        if (hasClosurePoint)\r\n          return false;\r\n        options.closed = false;\r\n      }\r\n\r\n      // append closure point if missing\r\n      if (options.closed) {\r\n        if (!hasClosurePoint) {\r\n          options.fitPoints.push(options.fitPoints[0].clone());\r\n          if (undefined !== options.knots) {  // best guess: uniform knots\r\n            options.knots.push(options.knots[options.knots.length - 1] + (options.knots[options.knots.length - 1] - options.knots[0]) / (options.knots.length - 1));\r\n          }\r\n        }\r\n        if (options.fitPoints.length <= 4)\r\n          options.closed = false; // can't fit cubic closed curve to 3 unique points\r\n      }\r\n\r\n      if (options.fitPoints.length < 2)\r\n        return false;\r\n\r\n      // ASSUME: tangents point INTO curve\r\n      if (undefined !== options.startTangent) {\r\n        if (options.startTangent.isAlmostZero)\r\n          options.startTangent = undefined;\r\n        else\r\n          options.startTangent.normalizeInPlace();\r\n      }\r\n      if (undefined !== options.endTangent) {\r\n        if (options.endTangent.isAlmostZero)\r\n          options.endTangent = undefined;\r\n        else\r\n          options.endTangent.normalizeInPlace();\r\n      }\r\n\r\n      return true;\r\n    }\r\n\r\n    /** Converts a full cubic knot vector of expected length into fit parameters, by removing extraneous exterior knots.\r\n     * @param knots cubic knot vector, unmodified\r\n     * @param numFitPoints number of fit points\r\n     * @return fit parameters, or undefined if unexpected input\r\n     **/\r\n    public static convertCubicKnotVectorToFitParams(knots: number[] | undefined, numFitPoints: number, normalize?: boolean): number[] | undefined {\r\n      let params = knots?.slice();\r\n      if (undefined !== params) {\r\n        const numExtraKnots = params.length - numFitPoints;\r\n        switch (numExtraKnots) {\r\n          case 0: {   // ASSUME caller passed in interior knots\r\n            break;\r\n          }\r\n          case 4:     // modern full cubic knots\r\n          case 6: {   // legacy full cubic knots\r\n            for (let i = 0; i < numExtraKnots / 2; ++i) {\r\n              params.pop();\r\n              params.shift();\r\n            }\r\n            break;\r\n          }\r\n          default: {  // other knot configurations are unusable\r\n            params = undefined;\r\n            break;\r\n          }\r\n        }\r\n        if (normalize && !this.normalizeKnots(params))\r\n          params = undefined;\r\n      }\r\n      return params;\r\n    }\r\n\r\n    /** Return fit parameters augmented to a full cubic knot vector.\r\n     * @param params fit parameters, unmodified\r\n     * @param legacy whether to create a legacy (DGN) full knot vector, or modern vector with two less knots\r\n     **/\r\n    public static convertFitParamsToCubicKnotVector(params: number[] | undefined, closed?: boolean, legacy?: boolean): number[] | undefined {\r\n      const knots = params?.slice();\r\n      if (undefined !== knots) {\r\n        const numExtraKnots = legacy ? 6 : 4;\r\n        if (closed) {\r\n          const iTail = knots.length - 2;\r\n          for (let iHead = 2; iHead <= numExtraKnots; iHead += 2) {\r\n            knots.unshift(knots[iTail] - 1.0);  // index is constant\r\n            knots.push(1.0 + knots[iHead]);     // index increments by two\r\n          }\r\n        } else {\r\n          for (let i = 0; i < numExtraKnots / 2; ++i) {\r\n            knots.unshift(0.0);\r\n            knots.push(1.0);\r\n          }\r\n        }\r\n      }\r\n      return knots;\r\n    }\r\n\r\n    /** Ensure full legacy knot vector for JSON export **/\r\n    public static convertToJsonKnots(props: InterpolationCurve3dProps) {\r\n      if (undefined !== props.knots) {\r\n        props.knots = this.convertCubicKnotVectorToFitParams(props.knots, props.fitPoints.length, false);\r\n        props.knots = this.convertFitParamsToCubicKnotVector(props.knots, props.closed, true);\r\n      } else {\r\n        props.knots = this.constructFitParametersFromPoints(Point3dArray.clonePoint3dArray(props.fitPoints), props.isChordLenKnots, props.closed);\r\n        props.knots = this.convertFitParamsToCubicKnotVector(props.knots, props.closed, true);\r\n      }\r\n    }\r\n\r\n    /** Construct the control points for the c2 cubic fit algorithm\r\n     * @param options validated as per validateOptions, possibly modified\r\n     */\r\n    public static constructPoles(options: InterpolationCurve3dOptions): Point3d[] | Float64Array | undefined {\r\n      if (!this.constructFitParameters(options) || (undefined === options.knots))\r\n        return undefined;\r\n\r\n      const numRow = options.fitPoints.length;\r\n      const alpha: number[] = Array(numRow);\r\n      const beta: number[] = Array(numRow);\r\n      const gamma: number[] = Array(numRow);\r\n      if (!this.setUpSystem(alpha, beta, gamma, options))\r\n        return undefined;\r\n\r\n      let poles: Point3d[] | Float64Array | undefined = [];\r\n      if (!options.closed) {\r\n        const dataPts = options.fitPoints.slice();\r\n        if (!this.setEndConditions(dataPts, options))\r\n          return undefined;\r\n        if (dataPts.length !== numRow + 2)\r\n          return undefined; // sanity check: we added 2nd/penultimate points, rhs is middle numRow entries.\r\n\r\n        // construct tridiagonal banded system components\r\n        const matrix = new Float64Array(numRow * 3);\r\n        const rhs = new Float64Array(numRow * 3);\r\n        for (let iRow = 0, iMatrixRead = 0, iRhsRead = 0; iRow < numRow; ++iRow) {\r\n          matrix[iMatrixRead++] = alpha[iRow];\r\n          matrix[iMatrixRead++] = beta[iRow];\r\n          matrix[iMatrixRead++] = gamma[iRow];\r\n          rhs[iRhsRead++] = dataPts[iRow+1].x;\r\n          rhs[iRhsRead++] = dataPts[iRow+1].y;\r\n          rhs[iRhsRead++] = dataPts[iRow+1].z;\r\n        }\r\n\r\n        const solution = BandedSystem.solveBandedSystemMultipleRHS(numRow, 3, matrix, 3, rhs);\r\n        if (undefined === solution)\r\n          return undefined;\r\n\r\n        // pre/append first/last poles/fitPoints\r\n        poles = new Float64Array(3 + solution.length + 3);\r\n        let iWrite = 0;\r\n        poles[iWrite++] = options.fitPoints[0].x;\r\n        poles[iWrite++] = options.fitPoints[0].y;\r\n        poles[iWrite++] = options.fitPoints[0].z;\r\n        for (let iRead = 0; iRead < solution.length; ) {\r\n          poles[iWrite++] = solution[iRead++];\r\n        }\r\n        poles[iWrite++] = options.fitPoints[options.fitPoints.length - 1].x;\r\n        poles[iWrite++] = options.fitPoints[options.fitPoints.length - 1].y;\r\n        poles[iWrite++] = options.fitPoints[options.fitPoints.length - 1].z;\r\n      } else { // closed\r\n        if (undefined !== (poles = this.solveNearTridiagonal(options.fitPoints, alpha, beta, gamma))) {\r\n          if (poles.length > 2) {\r\n            poles.unshift(poles.pop()!);  // shift poles right to line up with the knots\r\n            for (let i = 0; i < options.order - 1; ++i)\r\n              poles.push(poles[i].clone()); // periodically extend (the modern way)\r\n          }\r\n        }\r\n      }\r\n      return poles;\r\n    }\r\n  }\r\n}\r\n"]}
+{"version":3,"file":"BSplineCurveOps.js","sourceRoot":"","sources":["../../../src/bspline/BSplineCurveOps.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAEH,8EAA8E;AAE9E,OAAO,EAAE,gBAAgB,EAAE,MAAM,gCAAgC,CAAC;AAClE,OAAO,EAAE,oBAAoB,EAAE,MAAM,oCAAoC,CAAC;AAC1E,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AACvC,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AAClE,OAAO,EAAE,YAAY,EAAE,MAAM,0BAA0B,CAAC;AACxD,OAAO,EAAE,cAAc,EAAE,MAAM,gBAAgB,CAAC;AAChD,OAAO,EAAE,eAAe,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAE3D,OAAO,EAAE,YAAY,EAAE,MAAM,4BAA4B,CAAC;AAE1D;;;GAGG;AACH,MAAM,OAAO,eAAe;IAC1B;;;;OAIG;IACI,MAAM,CAAC,mBAAmB,CAAC,MAAwC,EAAE,KAAa;QACvF,MAAM,SAAS,GAAG,MAAM,CAAC,MAAM,CAAC;QAChC,IAAI,KAAK,GAAG,SAAS,IAAI,KAAK,GAAG,CAAC;YAChC,OAAO,SAAS,CAAC;QACnB,MAAM,MAAM,GAAG,KAAK,GAAG,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,CAAI,oDAAoD;QAClF,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,SAAS,CAAC,CAAC;QAChD,MAAM,cAAc,GAAG,IAAI,YAAY,CAAC,KAAK,CAAC,CAAC;QAC/C,MAAM,GAAG,GAAG,IAAI,gBAAgB,EAAE,CAAC;QACnC,MAAM,KAAK,GAAG,UAAU,CAAC,oBAAoB,CAAC,SAAS,EAAE,KAAK,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;QAC9E,MAAM,GAAG,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC;QAC7B,KAAK,IAAI,cAAc,GAAG,CAAC,EAAE,cAAc,GAAG,SAAS,EAAE,cAAc,EAAE,EAAE;YACzE,MAAM,CAAC,GAAG,KAAK,CAAC,YAAY,CAAC,cAAc,CAAC,CAAC;YAC7C,MAAM,SAAS,GAAG,KAAK,CAAC,mBAAmB,CAAC,CAAC,CAAC,CAAC;YAC/C,KAAK,CAAC,sBAAsB,CAAC,SAAS,EAAE,CAAC,EAAE,cAAc,CAAC,CAAC;YAC3D,oEAAoE;YACpE,IAAI,QAAQ,GAAG,CAAC,CAAC;YACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE;gBAC5B,IAAI,cAAc,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,QAAQ,CAAC;oBAC9C,QAAQ,GAAG,CAAC,CAAC;YACjB,MAAM,2BAA2B,GAAG,MAAM,GAAG,QAAQ,CAAC;YACtD,MAAM,QAAQ,GAAG,cAAc,GAAG,EAAE,CAAC;YACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE;gBAC9B,MAAM,UAAU,GAAG,cAAc,GAAG,MAAM,GAAG,2BAA2B,GAAG,CAAC,CAAC;gBAC7E,IAAI,QAAQ,GAAG,UAAU,IAAI,CAAC,IAAI,UAAU,GAAG,SAAS;oBACtD,MAAM,CAAC,QAAQ,GAAG,2BAA2B,GAAG,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC,CAAC;aAC1E;YACD,IAAI,MAAM,YAAY,oBAAoB,EAAE;gBAC1C,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,+BAA+B,CAAC,cAAc,EAAE,GAAG,CAAC,CAAC,CAAC;aACvE;iBAAM;gBACL,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,cAAc,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC;aAC1C;SACF;QACD,MAAM,KAAK,GAAG,YAAY,CAAC,4BAA4B,CAAC,SAAS,EAAE,EAAE,EAAE,MAAM,EAAE,CAAC,EAAE,GAAG,CAAC,WAAW,EAAE,CAAC,CAAC;QACrG,IAAI,KAAK,EAAE;YACT,OAAO,cAAc,CAAC,MAAM,CAAC,KAAK,EAAE,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;SACzD;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IAED;;;OAGG;IACI,MAAM,CAAC,0BAA0B,CAAC,OAAoC;QAC3E,+GAA+G;QAC/G,wEAAwE;QACxE,MAAM,gBAAgB,GAAG,OAAO,CAAC,KAAK,EAAE,CAAC;QACzC,IAAI,CAAC,IAAI,CAAC,UAAU,CAAC,eAAe,CAAC,gBAAgB,CAAC;YACpD,OAAO,SAAS,CAAC;QAEnB,MAAM,KAAK,GAAG,IAAI,CAAC,UAAU,CAAC,cAAc,CAAC,gBAAgB,CAAC,CAAC;QAC/D,IAAI,SAAS,KAAK,KAAK;YACrB,OAAO,SAAS,CAAC;QAEnB,MAAM,SAAS,GAAG,IAAI,CAAC,UAAU,CAAC,iCAAiC,CAAC,gBAAgB,CAAC,KAAK,EAAE,gBAAgB,CAAC,MAAM,CAAC,CAAC;QACrH,IAAI,SAAS,KAAK,SAAS;YACzB,OAAO,SAAS,CAAC;QAEnB,MAAM,WAAW,GAAG,cAAc,CAAC,MAAM,CAAC,KAAK,EAAE,SAAS,EAAE,gBAAgB,CAAC,KAAK,CAAC,CAAC;QAEpF,IAAI,gBAAgB,CAAC,MAAM;YACzB,WAAW,aAAX,WAAW,uBAAX,WAAW,CAAE,YAAY,CAAC,eAAe,CAAC,yBAAyB,CAAC,CAAC;QAEvE,OAAO,WAAW,CAAC;IACrB,CAAC;CACF;AAED;;;GAGG;AACH,WAAiB,eAAe;IAC9B;;;;;OAKG;IACH,MAAa,UAAU;QACrB;;WAEG;QACK,MAAM,CAAC,cAAc,CAAC,KAA2B;YACvD,IAAI,SAAS,KAAK,KAAK,IAAI,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE;gBAC3C,KAAK,GAAG,SAAS,CAAC;gBAClB,OAAO,KAAK,CAAC;aACd;YACD,MAAM,OAAO,GAAG,UAAU,CAAC,MAAM,CAAC,KAAK,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,8BAA8B;YAClF,IAAI,CAAC,OAAO,CAAC,SAAS,EAAE,EAAE;gBACxB,KAAK,GAAG,SAAS,CAAC;gBAClB,OAAO,KAAK,CAAC;aACd;YACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,CAAC,MAAM,EAAE,EAAE,CAAC;gBACnC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YAC9B,OAAO,IAAI,CAAC;QACd,CAAC;QAED,qEAAqE;QAC7D,MAAM,CAAC,8BAA8B,CAAC,SAAoB;YAChE,IAAI,SAAS,CAAC,MAAM,GAAG,CAAC;gBACtB,OAAO,SAAS,CAAC;YACnB,MAAM,MAAM,GAAa,CAAC,GAAG,CAAC,CAAC;YAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,EAAE,CAAC;gBACvC,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YACtE,IAAI,CAAC,IAAI,CAAC,cAAc,CAAC,MAAM,CAAC;gBAC9B,OAAO,SAAS,CAAC;YACnB,OAAO,MAAM,CAAC;QAChB,CAAC;QAED,gEAAgE;QACxD,MAAM,CAAC,0BAA0B,CAAC,SAAiB;YACzD,IAAI,SAAS,GAAG,CAAC;gBACf,OAAO,SAAS,CAAC;YACnB,MAAM,KAAK,GAAG,UAAU,CAAC,oBAAoB,CAAC,SAAS,GAAG,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;YAC1E,MAAM,MAAM,GAAa,EAAE,CAAC;YAC5B,KAAK,IAAI,CAAC,GAAG,KAAK,CAAC,aAAa,EAAE,CAAC,IAAI,KAAK,CAAC,cAAc,EAAE,EAAE,CAAC;gBAC9D,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YAC9B,OAAO,MAAM,CAAC;QAChB,CAAC;QAED,qEAAqE;QAC7D,MAAM,CAAC,wBAAwB,CAAC,OAAoC;YAC1E,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK,IAAI,OAAO,CAAC,KAAK,CAAC,MAAM,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM;gBAClF,OAAO,CAAC,KAAK,GAAG,SAAS,CAAC;YAE5B,gDAAgD;YAChD,MAAM,MAAM,GAAG,gBAAgB,CAAC,MAAM,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;YAC1D,MAAM,OAAO,GAAG,MAAM,CAAC,qBAAqB,EAAE,CAAC;YAC/C,MAAM,CAAC,KAAK,EAAE,CAAC;YAEf,8BAA8B;YAC9B,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,MAAM,GAAG,CAAC,EAAE,KAAK,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,EAAE,EAAE,KAAK,EAAE;gBACzE,IAAI,KAAK,KAAK,OAAO,CAAC,MAAM,CAAC;oBAC3B,EAAE,MAAM,CAAC,CAAC,qBAAqB;;oBAE/B,MAAM,CAAC,IAAI,CAAC,OAAO,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC;aACjD;YACD,OAAO,CAAC,SAAS,GAAG,MAAM,CAAC,eAAe,EAAE,CAAC;YAE7C,oDAAoD;YACpD,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK,EAAE;gBAC/B,MAAM,QAAQ,GAAa,EAAE,CAAC;gBAC9B,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,MAAM,GAAG,CAAC,EAAE,KAAK,GAAG,OAAO,CAAC,KAAK,CAAC,MAAM,EAAE,EAAE,KAAK,EAAE;oBACrE,IAAI,KAAK,KAAK,OAAO,CAAC,MAAM,CAAC;wBAC3B,EAAE,MAAM,CAAC,CAAC,OAAO;;wBAEjB,QAAQ,CAAC,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC;iBACvC;gBACD,OAAO,CAAC,KAAK,GAAG,QAAQ,CAAC,KAAK,EAAE,CAAC;aAClC;YACD,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;;UAKE;QACM,MAAM,CAAC,gCAAgC,CAAC,SAAoB,EAAE,aAAiC,EAAE,MAA2B;YAClI,IAAI,MAA4B,CAAC;YACjC,IAAI,aAAa,IAAI,CAAC,MAAM;gBAC1B,MAAM,GAAG,IAAI,CAAC,8BAA8B,CAAC,SAAS,CAAC,CAAC;YAC1D,IAAI,SAAS,KAAK,MAAM;gBACtB,MAAM,GAAG,IAAI,CAAC,0BAA0B,CAAC,SAAS,CAAC,MAAM,CAAC,CAAC;YAC7D,OAAO,MAAM,CAAC;QAChB,CAAC;QAED;;;WAGG;QACK,MAAM,CAAC,sBAAsB,CAAC,OAAoC;;YACxE,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK;gBAC7B,OAAO,CAAC,KAAK,GAAG,IAAI,CAAC,gCAAgC,CAAC,OAAO,CAAC,SAAS,EAAE,OAAO,CAAC,eAAe,EAAE,OAAO,CAAC,MAAM,CAAC,CAAC;YACpH,OAAO,CAAA,MAAA,OAAO,CAAC,KAAK,0CAAE,MAAM,MAAK,OAAO,CAAC,SAAS,CAAC,MAAM,CAAC;QAC5D,CAAC;QAED;;;;;WAKG;QACK,MAAM,CAAC,qBAAqB,CAAC,KAAe,EAAE,IAAc,EAAE,KAAe,EAAE,KAAa,EAAE,WAAmB,EAAE,MAAc,EAAE,YAAoB,EAAE,YAAoB;YACnL,IAAI,eAAe,GAAG,GAAG,GAAG,CAAC,YAAY,GAAG,YAAY,GAAG,MAAM,CAAC,CAAC;YACnE,KAAK,CAAC,KAAK,CAAC,GAAG,MAAM,GAAG,MAAM,GAAG,eAAe,CAAC;YACjD,IAAI,CAAC,KAAK,CAAC,GAAG,MAAM,GAAG,CAAC,YAAY,GAAG,YAAY,CAAC,GAAG,eAAe,CAAC;YAEvE,eAAe,GAAG,GAAG,GAAG,CAAC,YAAY,GAAG,MAAM,GAAG,WAAW,CAAC,CAAC;YAC9D,IAAI,CAAC,KAAK,CAAC,IAAI,YAAY,GAAG,CAAC,MAAM,GAAG,WAAW,CAAC,GAAG,eAAe,CAAC;YACvE,KAAK,CAAC,KAAK,CAAC,GAAG,YAAY,GAAG,YAAY,GAAG,eAAe,CAAC;YAE7D,eAAe,GAAG,GAAG,GAAG,CAAC,YAAY,GAAG,MAAM,CAAC,CAAC;YAChD,KAAK,CAAC,KAAK,CAAC,IAAI,eAAe,CAAC;YAChC,IAAI,CAAC,KAAK,CAAC,IAAI,eAAe,CAAC;YAC/B,KAAK,CAAC,KAAK,CAAC,IAAI,eAAe,CAAC;QAClC,CAAC;QAED;;;;WAIG;QACK,MAAM,CAAC,kBAAkB,CAAC,KAAe,EAAE,IAAc,EAAE,KAAe;YAChF,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC;gBAC/D,OAAO,KAAK,CAAC;YAEf,kBAAkB;YAClB,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;YAChD,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;YACxB,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;;;;WAOG;QACM,MAAM,CAAC,kBAAkB,CAAC,KAAe,EAAE,IAAc,EAAE,KAAe,EAAE,OAAoC,EAAE,sBAA+B,EAAE,oBAA6B;YACvL,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK;gBAC7B,OAAO,KAAK,CAAC;YACf,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC;gBAC/D,OAAO,KAAK,CAAC;YACf,IAAI,OAAO,CAAC,KAAK,CAAC,MAAM,KAAK,CAAC,IAAI,OAAO,CAAC,SAAS,CAAC,MAAM,KAAK,CAAC;gBAC9D,OAAO,KAAK,CAAC;YAEf,IAAI,WAAW,GAAG,CAAC,EAAE,MAAM,GAAG,CAAC,EAAE,YAAY,GAAG,CAAC,EAAE,YAAY,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,aAAa,GAAG,CAAC,CAAC;YAEhG,YAAY;YACZ,IAAI,sBAAsB,EAAE;gBAC1B,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;gBACf,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAC7C,WAAW,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAClD,GAAG,GAAG,MAAM,GAAG,WAAW,CAAC;gBAC3B,aAAa,GAAG,GAAG,GAAG,GAAG,CAAC;gBAC1B,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,GAAG,aAAa,CAAC;gBACzC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,aAAa,CAAC;aACpC;iBAAM;gBACL,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;gBAC1B,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;aACf;YAED,aAAa;YACb,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YACnD,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YAC7C,aAAa,GAAG,GAAG,GAAG,CAAC,YAAY,GAAG,MAAM,CAAC,CAAC;YAC9C,aAAa,IAAI,aAAa,CAAC;YAC/B,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,MAAM,GAAG,aAAa,CAAC;YAC3C,IAAI,CAAC,CAAC,CAAC,GAAI,GAAG,GAAG,CAAC,MAAM,GAAG,YAAY,CAAC,GAAG,aAAa,CAAC;YACzD,KAAK,CAAC,CAAC,CAAC,GAAG,YAAY,GAAG,YAAY,GAAG,aAAa,CAAC;YAEvD,WAAW;YACX,IAAI,oBAAoB,EAAE;gBACxB,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBACnD,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBACnD,GAAG,GAAG,YAAY,GAAG,YAAY,CAAC;gBAClC,aAAa,GAAG,GAAG,GAAG,GAAG,CAAC;gBAC1B,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,YAAY,GAAG,aAAa,CAAC;gBACzC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,YAAY,GAAG,GAAG,CAAC,GAAG,aAAa,CAAC;gBAC/C,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;aAChB;iBAAM;gBACL,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;gBAC1B,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;aACf;YACD,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;;;;WAOG;QACK,MAAM,CAAC,wBAAwB,CAAC,KAAe,EAAE,IAAc,EAAE,KAAe,EAAE,OAAoC,EAAE,sBAA+B,EAAE,oBAA6B;YAC5L,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK;gBAC7B,OAAO,KAAK,CAAC;YACf,IAAI,KAAK,CAAC,MAAM,KAAK,IAAI,CAAC,MAAM,IAAI,KAAK,CAAC,MAAM,KAAK,KAAK,CAAC,MAAM,IAAI,KAAK,CAAC,MAAM,KAAK,OAAO,CAAC,KAAK,CAAC,MAAM;gBACxG,OAAO,KAAK,CAAC;YACf,IAAI,OAAO,CAAC,KAAK,CAAC,MAAM,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM;gBACnD,OAAO,KAAK,CAAC;YAEf,MAAM,YAAY,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;YAClD,MAAM,kBAAkB,GAAG,YAAY,GAAG,CAAC,CAAC;YAC5C,IAAI,WAAW,GAAG,CAAC,EAAE,MAAM,GAAG,CAAC,EAAE,YAAY,GAAG,CAAC,EAAE,YAAY,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,aAAa,GAAG,CAAC,CAAC;YAEhG,IAAI,OAAO,CAAC,MAAM,EAAE;gBAClB,YAAY;gBACZ,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAC7C,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,CAAC;gBACzF,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,CAAC,CAAC;gBACzF,WAAW,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAClD,IAAI,CAAC,qBAAqB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,WAAW,EAAE,MAAM,EAAE,YAAY,EAAE,YAAY,CAAC,CAAC;gBAEnG,aAAa;gBACb,YAAY,GAAG,YAAY,CAAC;gBAC5B,YAAY,GAAG,MAAM,CAAC;gBACtB,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAC7C,WAAW,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAClD,IAAI,CAAC,qBAAqB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,WAAW,EAAE,MAAM,EAAE,YAAY,EAAE,YAAY,CAAC,CAAC;gBAEnG,qDAAqD;gBACrD,WAAW,GAAG,YAAY,CAAC;gBAC3B,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,CAAC,CAAC;gBACnF,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,CAAC;gBAC7F,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,CAAC;gBACzF,IAAI,CAAC,qBAAqB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,kBAAkB,EAAE,WAAW,EAAE,MAAM,EAAE,YAAY,EAAE,YAAY,CAAC,CAAC;aACrH;iBAAM,EAAE,OAAO;gBACd,YAAY;gBACZ,IAAI,sBAAsB,EAAE;oBAC1B,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;oBACf,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;oBAC7C,WAAW,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;oBAClD,GAAG,GAAG,MAAM,GAAG,WAAW,CAAC;oBAC3B,aAAa,GAAG,GAAG,GAAG,GAAG,CAAC;oBAC1B,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,GAAG,aAAa,CAAC;oBACzC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,aAAa,CAAC;iBACpC;qBAAM;oBACL,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;oBAC1B,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;iBACf;gBAED,aAAa;gBACb,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAC7C,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBACnD,YAAY,GAAG,GAAG,CAAC;gBACnB,WAAW,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAClD,IAAI,CAAC,qBAAqB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,WAAW,EAAE,MAAM,EAAE,YAAY,EAAE,YAAY,CAAC,CAAC;gBAEnG,kBAAkB;gBAClB,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,CAAC,CAAC;gBACnF,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,CAAC;gBACzF,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,kBAAkB,GAAG,CAAC,CAAC,CAAC;gBAC7F,WAAW,GAAG,GAAG,CAAC;gBAClB,IAAI,CAAC,qBAAqB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,kBAAkB,EAAE,WAAW,EAAE,MAAM,EAAE,YAAY,EAAE,YAAY,CAAC,CAAC;gBAEpH,WAAW;gBACX,IAAI,oBAAoB,EAAE;oBACxB,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,YAAY,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC;oBAC7E,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC;oBACjF,GAAG,GAAG,YAAY,GAAG,YAAY,CAAC;oBAClC,aAAa,GAAG,GAAG,GAAG,GAAG,CAAC;oBAC1B,KAAK,CAAC,YAAY,CAAC,GAAG,CAAC,YAAY,GAAG,aAAa,CAAC;oBACpD,IAAI,CAAC,YAAY,CAAC,GAAG,CAAC,YAAY,GAAG,GAAG,CAAC,GAAG,aAAa,CAAC;oBAC1D,KAAK,CAAC,YAAY,CAAC,GAAG,GAAG,CAAC;iBAC3B;qBAAM;oBACL,KAAK,CAAC,YAAY,CAAC,GAAG,KAAK,CAAC,YAAY,CAAC,GAAG,GAAG,CAAC;oBAChD,IAAI,CAAC,YAAY,CAAC,GAAG,GAAG,CAAC;iBAC1B;aACF;YAED,cAAc;YACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,kBAAkB,EAAE,EAAE,CAAC,EAAE;gBAC3C,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBACjD,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC3D,YAAY,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBACvD,WAAW,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC1D,IAAI,CAAC,qBAAqB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,WAAW,EAAE,MAAM,EAAE,YAAY,EAAE,YAAY,CAAC,CAAC;aACpG;YACD,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;;WAKG;QACK,MAAM,CAAC,WAAW,CAAC,KAAe,EAAE,IAAc,EAAE,KAAe,EAAE,OAAoC;YAC/G,IAAI,sBAAsB,GAAG,KAAK,CAAC;YACnC,IAAI,oBAAoB,GAAG,KAAK,CAAC;YACjC,IAAI,OAAO,CAAC,iBAAiB,IAAI,CAAC,OAAO,CAAC,MAAM,EAAE;gBAChD,sBAAsB,GAAG,CAAC,SAAS,KAAK,OAAO,CAAC,YAAY,CAAC,CAAC;gBAC9D,oBAAoB,GAAG,CAAC,SAAS,KAAK,OAAO,CAAC,UAAU,CAAC,CAAC;aAC3D;YAED,IAAI,SAAS,GAAG,KAAK,CAAC;YACtB,IAAI,CAAC,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM;gBAChC,SAAS,GAAG,IAAI,CAAC,kBAAkB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC;iBACrD,IAAI,CAAC,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM;gBACrC,SAAS,GAAG,IAAI,CAAC,kBAAkB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,OAAO,EAAE,sBAAsB,EAAE,oBAAoB,CAAC,CAAC;iBAC5G,IAAI,CAAC,IAAI,OAAO,CAAC,SAAS,CAAC,MAAM;gBACpC,SAAS,GAAG,IAAI,CAAC,wBAAwB,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,OAAO,EAAE,sBAAsB,EAAE,oBAAoB,CAAC,CAAC;YAEvH,OAAO,SAAS,CAAC;QACnB,CAAC;QAED;;;;;WAKG;QACK,MAAM,CAAC,qBAAqB,CAAC,OAAkB,EAAE,OAAoC,EAAE,OAAgB;YAC7G,IAAI,OAAO,CAAC,MAAM,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC;gBACjD,OAAO,KAAK,CAAC;YACf,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK;gBAC7B,OAAO,KAAK,CAAC;YAEf,MAAM,KAAK,GAAG,GAAG,GAAC,GAAG,CAAC;YACtB,MAAM,YAAY,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;YAElD,IAAI,CAAC,KAAK,YAAY,EAAE,EAAE,gBAAgB;gBACxC,IAAI,OAAO;oBACT,OAAO,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;;oBAEtD,OAAO,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;gBACxD,OAAO,IAAI,CAAC;aACb;YAED,IAAI,CAAC,KAAK,YAAY,EAAE;gBACtB,MAAM,KAAK,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC5F,MAAM,IAAI,GAAG,GAAG,GAAG,KAAK,CAAC;gBACzB,MAAM,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,KAAK,GAAG,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;gBAC1F,IAAI,OAAO;oBACT,OAAO,CAAC,gBAAgB,CAAC,IAAI,EAAE,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,WAAW,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;;oBAElH,OAAO,CAAC,gBAAgB,CAAC,IAAI,EAAE,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,WAAW,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;gBAClH,OAAO,IAAI,CAAC;aACb;YAED,mBAAmB;YACnB,IAAI,OAAO,EAAE;gBACT,MAAM,KAAK,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC5F,MAAM,IAAI,GAAG,GAAG,GAAG,KAAK,CAAC;gBACzB,MAAM,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,KAAK,GAAG,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;gBAC1F,OAAO,CAAC,gBAAgB,CAAC,IAAI,EAAE,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,WAAW,CAAC,KAAK,EAAE,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;aACrH;iBAAM;gBACH,MAAM,KAAK,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,YAAY,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,YAAY,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC;gBAChJ,MAAM,IAAI,GAAG,GAAG,GAAG,KAAK,CAAC;gBACzB,MAAM,IAAI,GAAG,OAAO,CAAC,YAAY,CAAC,CAAC,WAAW,CAAC,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,EAAE,CAAC,KAAK,GAAG,KAAK,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;gBACnI,OAAO,CAAC,gBAAgB,CAAC,IAAI,EAAE,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,WAAW,CAAC,KAAK,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC;aAChK;YACD,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;;;WAMG;QACK,MAAM,CAAC,sBAAsB,CAAC,OAAkB,EAAE,OAAoC,EAAE,OAAgB;YAC9G,IAAI,OAAO,CAAC,MAAM,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC;gBACjD,OAAO,KAAK,CAAC;YAEf,MAAM,YAAY,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;YAClD,IAAI,CAAC,KAAK,YAAY;gBACpB,OAAO,IAAI,CAAC,qBAAqB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;YAE/D,IAAI,OAAO;gBACT,OAAO,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;;gBAExB,OAAO,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;YAC5D,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;;;WAMG;QACK,MAAM,CAAC,gCAAgC,CAAC,OAAkB,EAAE,OAAoC,EAAE,OAAgB;YACxH,IAAI,OAAO,CAAC,MAAM,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC;gBACjD,OAAO,KAAK,CAAC;YAEf,MAAM,OAAO,GAAG,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC,YAAY,CAAC,CAAC,CAAC,OAAO,CAAC,UAAU,CAAC;YACpE,IAAI,SAAS,KAAK,OAAO;gBACvB,OAAO,KAAK,CAAC;YAEf,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,+BAA+B;YAC7C,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,2FAA2F;YACzG,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,oCAAoC;YAElD,MAAM,YAAY,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;YAClD,IAAI,CAAC,KAAK,YAAY,EAAE,EAAE,yBAAyB;gBACjD,IAAI,OAAO,EAAE;oBACX,IAAI,GAAG,CAAC,CAAC;oBACT,IAAI,GAAG,CAAC,CAAC;oBACT,IAAI,GAAG,CAAC,CAAC;iBACV;qBAAM;oBACL,IAAI,GAAG,CAAC,CAAC;oBACT,IAAI,GAAG,CAAC,CAAC;oBACT,IAAI,GAAG,CAAC,CAAC;iBACV;aACF;iBAAM;gBACL,IAAI,OAAO,EAAE;oBACX,IAAI,GAAG,CAAC,CAAC;oBACT,IAAI,GAAG,CAAC,CAAC;oBACT,IAAI,GAAG,CAAC,CAAC;iBACV;qBAAM;oBACL,IAAI,GAAG,YAAY,GAAG,CAAC,CAAC;oBACxB,IAAI,GAAG,YAAY,GAAG,CAAC,CAAC;oBACxB,IAAI,GAAG,YAAY,CAAC;iBACrB;aACF;YAED,kCAAkC;YAClC,MAAM,WAAW,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC;YAC1D,OAAO,CAAC,IAAI,CAAC,CAAC,UAAU,CAAC,OAAO,EAAE,WAAW,GAAG,GAAG,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC;YACpE,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;;WAKG;QACK,MAAM,CAAC,iCAAiC,CAAC,OAAkB,EAAE,OAAoC,EAAE,OAAgB;YACzH,IAAI,OAAO,CAAC,MAAM,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC;gBACjD,OAAO,KAAK,CAAC;YAEf,MAAM,OAAO,GAAG,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC,YAAY,CAAC,CAAC,CAAC,OAAO,CAAC,UAAU,CAAC;YACpE,IAAI,SAAS,KAAK,OAAO;gBACvB,OAAO,KAAK,CAAC;YAEf,uCAAuC;YACvC,IAAI,CAAC,IAAI,CAAC,qBAAqB,CAAC,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC;gBACxD,OAAO,KAAK,CAAC;YAEf,MAAM,YAAY,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;YAClD,MAAM,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,+BAA+B;YAC5E,MAAM,IAAI,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,2FAA2F;YAExI,kFAAkF;YAClF,OAAO,CAAC,IAAI,CAAC,CAAC,UAAU,CAAC,OAAO,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC;YACxF,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;WAIG;QACK,MAAM,CAAC,+BAA+B,CAAC,OAAkB,EAAE,OAAoC;YACrG,MAAM,YAAY,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC;YAClD,IAAI,CAAC,OAAO,CAAC,kBAAkB;mBACxB,YAAY,IAAI,CAAC;mBACjB,CAAC,SAAS,KAAK,OAAO,CAAC,YAAY,IAAI,SAAS,KAAK,OAAO,CAAC,UAAU,CAAC;mBACxE,OAAO,CAAC,iBAAiB;mBACzB,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,EAAE;gBAC3D,OAAO,IAAI,CAAC;aACX;YACH,qFAAqF;YACrF,IAAI,SAAS,KAAK,OAAO,CAAC,YAAY,EAAE,EAAE,4DAA4D;gBACpG,MAAM,mBAAmB,GAAG,QAAQ,CAAC,cAAc,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;gBACxF,IAAI,SAAS,KAAK,mBAAmB,EAAE;oBACrC,MAAM,aAAa,GAAG,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC;oBACpF,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,UAAU,CAAC,mBAAmB,EAAE,aAAa,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC;iBACrG;aACF;iBAAM,IAAI,SAAS,KAAK,OAAO,CAAC,UAAU,EAAE,EAAG,4DAA4D;gBAC1G,MAAM,iBAAiB,GAAG,QAAQ,CAAC,cAAc,CAAC,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;gBACpH,IAAI,SAAS,KAAK,iBAAiB,EAAE;oBACnC,MAAM,eAAe,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;oBACxD,OAAO,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,iBAAiB,EAAE,eAAe,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;iBACvE;aACF;iBAAM,EAAG,8DAA8D;gBACtE,MAAM,aAAa,GAAG,QAAQ,CAAC,cAAc,CAAC,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;gBACjG,IAAI,SAAS,KAAK,aAAa,EAAE;oBAC/B,MAAM,eAAe,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;oBACxD,OAAO,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,aAAa,EAAE,eAAe,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;oBAClE,MAAM,aAAa,GAAG,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC;oBACpF,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,UAAU,CAAC,aAAa,EAAE,CAAC,aAAa,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC;iBAChG;aACF;YACD,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;WAIG;QACK,MAAM,CAAC,gBAAgB,CAAC,OAAkB,EAAE,OAAoC;YACtF,IAAI,OAAO,CAAC,MAAM,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM;gBAC7C,OAAO,KAAK,CAAC;YAEf,2CAA2C;YAC3C,MAAM,MAAM,GAAG,OAAO,CAAC,UAAU,EAAE,CAAC;YACpC,MAAM,MAAM,GAAG,OAAO,CAAC,UAAU,EAAE,CAAC;YACpC,OAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;YAC7B,OAAO,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;YAE9C,IAAI,SAAS,GAAG,KAAK,CAAC;YACtB,IAAI,SAAS,KAAK,OAAO,CAAC,YAAY,EAAE;gBACtC,IAAI,OAAO,CAAC,iBAAiB;oBAC3B,SAAS,GAAG,IAAI,CAAC,sBAAsB,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;;oBAEhE,SAAS,GAAG,IAAI,CAAC,qBAAqB,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;aAClE;iBAAM,EAAE,qBAAqB;gBAC5B,IAAI,OAAO,CAAC,kBAAkB;oBAC5B,SAAS,GAAG,IAAI,CAAC,gCAAgC,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;;oBAE1E,SAAS,GAAG,IAAI,CAAC,iCAAiC,CAAC,OAAO,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC;aAC9E;YAED,IAAI,SAAS,KAAK,OAAO,CAAC,UAAU,EAAE;gBACpC,IAAI,OAAO,CAAC,iBAAiB;oBAC3B,SAAS,GAAG,IAAI,CAAC,sBAAsB,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;;oBAEjE,SAAS,GAAG,IAAI,CAAC,qBAAqB,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;aACnE;iBAAM,EAAE,mBAAmB;gBAC1B,IAAI,OAAO,CAAC,kBAAkB;oBAC5B,SAAS,GAAG,IAAI,CAAC,gCAAgC,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;;oBAE3E,SAAS,GAAG,IAAI,CAAC,iCAAiC,CAAC,OAAO,EAAE,OAAO,EAAE,KAAK,CAAC,CAAC;aAC/E;YAED,IAAI,SAAS;gBACX,SAAS,GAAG,IAAI,CAAC,+BAA+B,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;YAErE,OAAO,SAAS,CAAC;QACnB,CAAC;QAED;;WAEG;QACK,MAAM,CAAC,oBAAoB,CAAC,MAAiB,EAAE,KAAe,EAAE,IAAc,EAAE,KAAe;YACrG,IAAI,KAAK,CAAC,MAAM,KAAK,IAAI,CAAC,MAAM,IAAI,KAAK,CAAC,MAAM,KAAK,KAAK,CAAC,MAAM,IAAI,KAAK,CAAC,MAAM,KAAK,MAAM,CAAC,MAAM;gBACjG,OAAO,SAAS,CAAC;YACnB,MAAM,KAAK,GAAc,EAAE,CAAC;YAC5B,MAAM,YAAY,GAAG,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC;YACvC,MAAM,OAAO,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAE,wBAAwB;YAC9D,IAAI,GAAG,GAAuB,GAAG,CAAC;YAElC,6BAA6B;YAC7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,YAAY,EAAE,EAAE,CAAC,EAAE;gBACrC,IAAI,SAAS,KAAK,CAAC,GAAG,GAAG,QAAQ,CAAC,yBAAyB,CAAC,CAAE,KAAK,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;oBACnF,OAAO,SAAS,CAAC;gBACnB,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC9B,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBAC9B,OAAO,CAAC,CAAC,CAAC,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;aAClD;YAED,8BAA8B;YAC9B,IAAI,SAAS,KAAK,CAAC,GAAG,GAAG,QAAQ,CAAC,yBAAyB,CAAC,GAAG,EAAE,IAAI,CAAC,YAAY,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC;gBACjH,OAAO,SAAS,CAAC;YACnB,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC;YAC/B,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC;YAC5C,KAAK,IAAI,CAAC,GAAG,YAAY,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,EAAE;gBAC1C,IAAI,SAAS,KAAK,CAAC,GAAG,GAAG,QAAQ,CAAC,yBAAyB,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;oBACxE,OAAO,SAAS,CAAC;gBACnB,OAAO,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAE,KAAK,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,YAAY,GAAG,CAAC,CAAC,EAAE,CAAE,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC;gBAChI,KAAK,CAAC,CAAC,CAAC,GAAG,CAAE,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,YAAY,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;aACnF;YAED,8BAA8B;YAC9B,IAAI,SAAS,KAAK,CAAC,GAAG,GAAG,QAAQ,CAAC,yBAAyB,CAAC,GAAG,EAAE,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC/E,OAAO,SAAS,CAAC;YACnB,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;YACjD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,YAAY,EAAE,EAAE,CAAC,EAAE;gBACrC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;aACzD;YACD,OAAO,KAAK,CAAC;QACf,CAAC;QAED;;WAEG;QACI,MAAM,CAAC,eAAe,CAAC,OAAoC;YAChE,OAAO,CAAC,KAAK,GAAG,CAAC,CAAC;YAElB,yGAAyG;YACzG,OAAO,CAAC,KAAK,GAAG,IAAI,CAAC,iCAAiC,CAAC,OAAO,CAAC,KAAK,EAAE,OAAO,CAAC,SAAS,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;YAEpG,oEAAoE;YACtE,IAAI,CAAC,IAAI,CAAC,wBAAwB,CAAC,OAAO,CAAC;gBACzC,OAAO,KAAK,CAAC;YAEf,uDAAuD;YACvD,IAAI,eAAe,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC,OAAO,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;YAC1G,IAAI,CAAC,KAAK,OAAO,CAAC,SAAS,CAAC,MAAM,IAAI,eAAe,EAAE;gBACrD,OAAO,CAAC,SAAS,CAAC,GAAG,EAAE,CAAC;gBACxB,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK;oBAC7B,OAAO,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC;gBACtB,eAAe,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC,OAAO,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;aACvG;YACD,IAAI,OAAO,CAAC,SAAS,CAAC,MAAM,IAAI,CAAC,EAAE;gBACjC,IAAI,eAAe;oBACjB,OAAO,KAAK,CAAC;gBACf,OAAO,CAAC,MAAM,GAAG,KAAK,CAAC;aACxB;YAED,kCAAkC;YAClC,IAAI,OAAO,CAAC,MAAM,EAAE;gBAClB,IAAI,CAAC,eAAe,EAAE;oBACpB,OAAO,CAAC,SAAS,CAAC,IAAI,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC;oBACrD,IAAI,SAAS,KAAK,OAAO,CAAC,KAAK,EAAE,EAAG,4BAA4B;wBAC9D,OAAO,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;qBACzJ;iBACF;gBACD,IAAI,OAAO,CAAC,SAAS,CAAC,MAAM,IAAI,CAAC;oBAC/B,OAAO,CAAC,MAAM,GAAG,KAAK,CAAC,CAAC,kDAAkD;aAC7E;YAED,IAAI,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC;gBAC9B,OAAO,KAAK,CAAC;YAEf,oCAAoC;YACpC,IAAI,SAAS,KAAK,OAAO,CAAC,YAAY,EAAE;gBACtC,IAAI,OAAO,CAAC,YAAY,CAAC,YAAY;oBACnC,OAAO,CAAC,YAAY,GAAG,SAAS,CAAC;;oBAEjC,OAAO,CAAC,YAAY,CAAC,gBAAgB,EAAE,CAAC;aAC3C;YACD,IAAI,SAAS,KAAK,OAAO,CAAC,UAAU,EAAE;gBACpC,IAAI,OAAO,CAAC,UAAU,CAAC,YAAY;oBACjC,OAAO,CAAC,UAAU,GAAG,SAAS,CAAC;;oBAE/B,OAAO,CAAC,UAAU,CAAC,gBAAgB,EAAE,CAAC;aACzC;YAED,OAAO,IAAI,CAAC;QACd,CAAC;QAED;;;;YAII;QACG,MAAM,CAAC,iCAAiC,CAAC,KAA2B,EAAE,YAAoB,EAAE,SAAmB;YACpH,IAAI,MAAM,GAAG,KAAK,aAAL,KAAK,uBAAL,KAAK,CAAE,KAAK,EAAE,CAAC;YAC5B,IAAI,SAAS,KAAK,MAAM,EAAE;gBACxB,MAAM,aAAa,GAAG,MAAM,CAAC,MAAM,GAAG,YAAY,CAAC;gBACnD,QAAQ,aAAa,EAAE;oBACrB,KAAK,CAAC,CAAC,CAAC,EAAI,yCAAyC;wBACnD,MAAM;qBACP;oBACD,KAAK,CAAC,CAAC,CAAK,0BAA0B;oBACtC,KAAK,CAAC,CAAC,CAAC,EAAI,0BAA0B;wBACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,aAAa,GAAG,CAAC,EAAE,EAAE,CAAC,EAAE;4BAC1C,MAAM,CAAC,GAAG,EAAE,CAAC;4BACb,MAAM,CAAC,KAAK,EAAE,CAAC;yBAChB;wBACD,MAAM;qBACP;oBACD,OAAO,CAAC,CAAC,EAAG,yCAAyC;wBACnD,MAAM,GAAG,SAAS,CAAC;wBACnB,MAAM;qBACP;iBACF;gBACD,IAAI,SAAS,IAAI,CAAC,IAAI,CAAC,cAAc,CAAC,MAAM,CAAC;oBAC3C,MAAM,GAAG,SAAS,CAAC;aACtB;YACD,OAAO,MAAM,CAAC;QAChB,CAAC;QAED;;;YAGI;QACG,MAAM,CAAC,iCAAiC,CAAC,MAA4B,EAAE,MAAgB,EAAE,MAAgB;YAC9G,MAAM,KAAK,GAAG,MAAM,aAAN,MAAM,uBAAN,MAAM,CAAE,KAAK,EAAE,CAAC;YAC9B,IAAI,SAAS,KAAK,KAAK,EAAE;gBACvB,MAAM,aAAa,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACrC,IAAI,MAAM,EAAE;oBACV,MAAM,KAAK,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC;oBAC/B,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,IAAI,aAAa,EAAE,KAAK,IAAI,CAAC,EAAE;wBACtD,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,CAAC,CAAE,oBAAoB;wBACxD,KAAK,CAAC,IAAI,CAAC,GAAG,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAK,0BAA0B;qBAC/D;iBACF;qBAAM;oBACL,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,aAAa,GAAG,CAAC,EAAE,EAAE,CAAC,EAAE;wBAC1C,KAAK,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;wBACnB,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;qBACjB;iBACF;aACF;YACD,OAAO,KAAK,CAAC;QACf,CAAC;QAED,sDAAsD;QAC/C,MAAM,CAAC,kBAAkB,CAAC,KAAgC;YAC/D,IAAI,SAAS,KAAK,KAAK,CAAC,KAAK,EAAE;gBAC7B,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,iCAAiC,CAAC,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,SAAS,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;gBACjG,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,iCAAiC,CAAC,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;aACvF;iBAAM;gBACL,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,gCAAgC,CAAC,YAAY,CAAC,iBAAiB,CAAC,KAAK,CAAC,SAAS,CAAC,EAAE,KAAK,CAAC,eAAe,EAAE,KAAK,CAAC,MAAM,CAAC,CAAC;gBAC1I,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,iCAAiC,CAAC,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;aACvF;QACH,CAAC;QAED;;WAEG;QACI,MAAM,CAAC,cAAc,CAAC,OAAoC;YAC/D,IAAI,CAAC,IAAI,CAAC,sBAAsB,CAAC,OAAO,CAAC,IAAI,CAAC,SAAS,KAAK,OAAO,CAAC,KAAK,CAAC;gBACxE,OAAO,SAAS,CAAC;YAEnB,MAAM,MAAM,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,CAAC;YACxC,MAAM,KAAK,GAAa,KAAK,CAAC,MAAM,CAAC,CAAC;YACtC,MAAM,IAAI,GAAa,KAAK,CAAC,MAAM,CAAC,CAAC;YACrC,MAAM,KAAK,GAAa,KAAK,CAAC,MAAM,CAAC,CAAC;YACtC,IAAI,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,OAAO,CAAC;gBAChD,OAAO,SAAS,CAAC;YAEnB,IAAI,KAAK,GAAyC,EAAE,CAAC;YACrD,IAAI,CAAC,OAAO,CAAC,MAAM,EAAE;gBACnB,MAAM,OAAO,GAAG,OAAO,CAAC,SAAS,CAAC,KAAK,EAAE,CAAC;gBAC1C,IAAI,CAAC,IAAI,CAAC,gBAAgB,CAAC,OAAO,EAAE,OAAO,CAAC;oBAC1C,OAAO,SAAS,CAAC;gBACnB,IAAI,OAAO,CAAC,MAAM,KAAK,MAAM,GAAG,CAAC;oBAC/B,OAAO,SAAS,CAAC,CAAC,+EAA+E;gBAEnG,iDAAiD;gBACjD,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;gBAC5C,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;gBACzC,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,WAAW,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,EAAE,IAAI,GAAG,MAAM,EAAE,EAAE,IAAI,EAAE;oBACvE,MAAM,CAAC,WAAW,EAAE,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC;oBACpC,MAAM,CAAC,WAAW,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC;oBACnC,MAAM,CAAC,WAAW,EAAE,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC;oBACpC,GAAG,CAAC,QAAQ,EAAE,CAAC,GAAG,OAAO,CAAC,IAAI,GAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACpC,GAAG,CAAC,QAAQ,EAAE,CAAC,GAAG,OAAO,CAAC,IAAI,GAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACpC,GAAG,CAAC,QAAQ,EAAE,CAAC,GAAG,OAAO,CAAC,IAAI,GAAC,CAAC,CAAC,CAAC,CAAC,CAAC;iBACrC;gBAED,MAAM,QAAQ,GAAG,YAAY,CAAC,4BAA4B,CAAC,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;gBACtF,IAAI,SAAS,KAAK,QAAQ;oBACxB,OAAO,SAAS,CAAC;gBAEnB,wCAAwC;gBACxC,KAAK,GAAG,IAAI,YAAY,CAAC,CAAC,GAAG,QAAQ,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;gBAClD,IAAI,MAAM,GAAG,CAAC,CAAC;gBACf,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACzC,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACzC,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACzC,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,QAAQ,CAAC,MAAM,GAAI;oBAC7C,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,QAAQ,CAAC,KAAK,EAAE,CAAC,CAAC;iBACrC;gBACD,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBACpE,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBACpE,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,OAAO,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;aACrE;iBAAM,EAAE,SAAS;gBAChB,IAAI,SAAS,KAAK,CAAC,KAAK,GAAG,IAAI,CAAC,oBAAoB,CAAC,OAAO,CAAC,SAAS,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC,EAAE;oBAC5F,IAAI,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE;wBACpB,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,GAAG,EAAG,CAAC,CAAC,CAAE,8CAA8C;wBAC5E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,KAAK,GAAG,CAAC,EAAE,EAAE,CAAC;4BACxC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,uCAAuC;qBACxE;iBACF;aACF;YACD,OAAO,KAAK,CAAC;QACf,CAAC;KACF;IAtwBY,0BAAU,aAswBtB,CAAA;AACH,CAAC,EA9wBgB,eAAe,KAAf,eAAe,QA8wB/B","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 Bspline\r\n */\r\n\r\n/* eslint-disable @typescript-eslint/naming-convention, no-empty, no-console*/\r\n\r\nimport { GrowableXYZArray } from \"../geometry3d/GrowableXYZArray\";\r\nimport { IndexedXYZCollection } from \"../geometry3d/IndexedXYZCollection\";\r\nimport { Geometry } from \"../Geometry\";\r\nimport { Point3d, Vector3d } from \"../geometry3d/Point3dVector3d\";\r\nimport { BandedSystem } from \"../numerics/BandedSystem\";\r\nimport { BSplineCurve3d } from \"./BSplineCurve\";\r\nimport { BSplineWrapMode, KnotVector } from \"./KnotVector\";\r\nimport { InterpolationCurve3dOptions, InterpolationCurve3dProps } from \"./InterpolationCurve3d\";\r\nimport { Point3dArray } from \"../geometry3d/PointHelpers\";\r\n\r\n/**\r\n * A class with static methods for creating B-spline curves.\r\n * @public\r\n */\r\nexport class BSplineCurveOps {\r\n  /**\r\n   * Greville knot algorithm: alternative to traditional c2 cubic algorithm allowing for any order, but no end conditions or periodicity.\r\n   * @param points pass-through points.\r\n   * @param order bspline order (1 more than degree)\r\n   */\r\n  public static createThroughPoints(points: IndexedXYZCollection | Point3d[], order: number): BSplineCurve3d | undefined {\r\n    const numPoints = points.length;\r\n    if (order > numPoints || order < 2)\r\n      return undefined;\r\n    const degree = order - 1;\r\n    const bw = 1 + 2 * degree;    // probably less than that . . just zeros at fringe.\r\n    const matrix = new Float64Array(bw * numPoints);\r\n    const basisFunctions = new Float64Array(order);\r\n    const rhs = new GrowableXYZArray();\r\n    const knots = KnotVector.createUniformClamped(numPoints, order - 1, 0.0, 1.0);\r\n    const xyz = Point3d.create();\r\n    for (let basePointIndex = 0; basePointIndex < numPoints; basePointIndex++) {\r\n      const u = knots.grevilleKnot(basePointIndex);\r\n      const spanIndex = knots.knotToLeftKnotIndex(u);\r\n      knots.evaluateBasisFunctions(spanIndex, u, basisFunctions);\r\n      // puzzlement .. how do the max points shift within the order spots?\r\n      let maxIndex = 0;\r\n      for (let i = 1; i < order; i++)\r\n        if (basisFunctions[i] > basisFunctions[maxIndex])\r\n          maxIndex = i;\r\n      const basisFunctionStartWithinRow = degree - maxIndex;\r\n      const rowStart = basePointIndex * bw;\r\n      for (let i = 0; i < order; i++) {\r\n        const realColumn = basePointIndex - degree + basisFunctionStartWithinRow + i;\r\n        if (rowStart + realColumn >= 0 && realColumn < numPoints)\r\n          matrix[rowStart + basisFunctionStartWithinRow + i] = basisFunctions[i];\r\n      }\r\n      if (points instanceof IndexedXYZCollection) {\r\n        rhs.push(points.getPoint3dAtUncheckedPointIndex(basePointIndex, xyz));\r\n      } else {\r\n        rhs.push(points[basePointIndex].clone());\r\n      }\r\n    }\r\n    const poles = BandedSystem.solveBandedSystemMultipleRHS(numPoints, bw, matrix, 3, rhs.float64Data());\r\n    if (poles) {\r\n      return BSplineCurve3d.create(poles, knots.knots, order);\r\n    }\r\n    return undefined;\r\n  }\r\n\r\n  /**\r\n   * Construct BSplineCurve3d that fit points using the C2 cubic algorithm.\r\n   * @param options curve definition, unmodified\r\n   */\r\n  public static createThroughPointsC2Cubic(options: InterpolationCurve3dOptions): BSplineCurve3d | undefined {\r\n    // Work on a copy rather than installing computed knots/tangents that could become stale when fit points change\r\n    // Knots/tangents that come in, however, are used without recomputation.\r\n    const validatedOptions = options.clone();\r\n    if (!this.C2CubicFit.validateOptions(validatedOptions))\r\n      return undefined;\r\n\r\n    const poles = this.C2CubicFit.constructPoles(validatedOptions);\r\n    if (undefined === poles)\r\n      return undefined;\r\n\r\n    const fullKnots = this.C2CubicFit.convertFitParamsToCubicKnotVector(validatedOptions.knots, validatedOptions.closed);\r\n    if (undefined === fullKnots)\r\n      return undefined;\r\n\r\n    const interpolant = BSplineCurve3d.create(poles, fullKnots, validatedOptions.order);\r\n\r\n    if (validatedOptions.closed)\r\n      interpolant?.setWrappable(BSplineWrapMode.OpenByAddingControlPoints);\r\n\r\n    return interpolant;\r\n  }\r\n}\r\n\r\n/**\r\n * Namespace for collecting curve fit API\r\n * @public\r\n */\r\nexport namespace BSplineCurveOps {\r\n  /**\r\n   * A helper class for creating C2 cubic fit curves.\r\n   * Knots herein are understood to be *interior* knots (including the start/end knot),\r\n   * so that there is one knot per fit point. In other words, the knots are fit parameters.\r\n   * @private\r\n   */\r\n  export class C2CubicFit {\r\n    /** Transform fit parameters to span [0,1]\r\n     * @param knots fit parameters, normalized in place\r\n     */\r\n    private static normalizeKnots(knots: number[] | undefined): boolean {\r\n      if (undefined === knots || knots.length < 2) {\r\n        knots = undefined;\r\n        return false;\r\n      }\r\n      const myKnots = KnotVector.create(knots, 1, false); // no extra knots at start/end\r\n      if (!myKnots.normalize()) {\r\n        knots = undefined;\r\n        return false;\r\n      }\r\n      for (let i = 0; i < knots.length; ++i)\r\n        knots[i] = myKnots.knots[i];\r\n      return true;\r\n    }\r\n\r\n    /** Compute chord-length fit parameters for C2 cubic fit algorithm */\r\n    private static constructChordLengthParameters(fitPoints: Point3d[]): number[] | undefined {\r\n      if (fitPoints.length < 2)\r\n        return undefined;\r\n      const params: number[] = [0.0];\r\n      for (let i = 1; i < fitPoints.length; ++i)\r\n        params[i] = params[i - 1] + fitPoints[i].distance(fitPoints[i - 1]);\r\n      if (!this.normalizeKnots(params))\r\n        return undefined;\r\n      return params;\r\n    }\r\n\r\n    /** Compute uniform fit parameters for C2 cubic fit algorithm */\r\n    private static constructUniformParameters(numParams: number): number[] | undefined {\r\n      if (numParams < 2)\r\n        return undefined;\r\n      const knots = KnotVector.createUniformClamped(numParams + 2, 3, 0.0, 1.0);\r\n      const params: number[] = [];\r\n      for (let i = knots.leftKnotIndex; i <= knots.rightKnotIndex; ++i)\r\n        params.push(knots.knots[i]);\r\n      return params;\r\n    }\r\n\r\n    /** Remove duplicate fit points, and their given knots in parallel */\r\n    private static removeDuplicateFitPoints(options: InterpolationCurve3dOptions): boolean {\r\n      if (undefined !== options.knots && options.knots.length !== options.fitPoints.length)\r\n        options.knots = undefined;\r\n\r\n      // get indices of duplicate points to be removed\r\n      const newPts = GrowableXYZArray.create(options.fitPoints);\r\n      const indices = newPts.findOrderedDuplicates();\r\n      newPts.clear();\r\n\r\n      // remove duplicate fit points\r\n      for (let iRead = 0, iIndex = 0; iRead < options.fitPoints.length; ++iRead) {\r\n        if (iRead === indices[iIndex])\r\n          ++iIndex; // skip the duplicate\r\n        else\r\n          newPts.push(options.fitPoints[iRead].clone());\r\n      }\r\n      options.fitPoints = newPts.getPoint3dArray();\r\n\r\n      // remove params corresponding to removed fit points\r\n      if (undefined !== options.knots) {\r\n        const newKnots: number[] = [];\r\n        for (let iRead = 0, iIndex = 0; iRead < options.knots.length; ++iRead) {\r\n          if (iRead === indices[iIndex])\r\n            ++iIndex; // skip\r\n          else\r\n            newKnots.push(options.knots[iRead]);\r\n        }\r\n        options.knots = newKnots.slice();\r\n      }\r\n      return true;\r\n    }\r\n\r\n    /** Construct fit parameters for the c2 cubic fit algorithm.\r\n     * @param fitPoints validated fit points (should not contain duplicates)\r\n     * @param isChordLength whether knots are computed using distances between successive fit points\r\n     * @param closed whether curve is periodically defined\r\n     * @return fit parameters, one per fit point\r\n    */\r\n     public static constructFitParametersFromPoints(fitPoints: Point3d[], isChordLength: number | undefined, closed: boolean | undefined): number[] | undefined {\r\n      let params: number[] | undefined;\r\n      if (isChordLength || !closed)\r\n        params = this.constructChordLengthParameters(fitPoints);\r\n      if (undefined === params)\r\n        params = this.constructUniformParameters(fitPoints.length);\r\n      return params;\r\n    }\r\n\r\n    /** Construct fit parameters for the c2 cubic fit algorithm, if they are missing.\r\n     * @param options validated as per validateOptions, possibly modified\r\n     * @return whether fit parameters are valid\r\n     */\r\n     public static constructFitParameters(options: InterpolationCurve3dOptions): boolean {\r\n      if (undefined === options.knots)\r\n        options.knots = this.constructFitParametersFromPoints(options.fitPoints, options.isChordLenKnots, options.closed);\r\n      return options.knots?.length === options.fitPoints.length;\r\n    }\r\n\r\n    /** Compute a row of the tridiagonal system matrix from Farin.\r\n     * @param alpha sub-diagonal, length = # fit points\r\n     * @param beta diagonal, length = # fit points\r\n     * @param gamma super-diagonal, length = # fit points\r\n     * @param index 0-based row index to set\r\n     */\r\n    private static computeAlphaBetaGamma(alpha: number[], beta: number[], gamma: number[], index: number, deltaIPlus1: number, deltaI: number, deltaIMinus1: number, deltaIMinus2: number) {\r\n      let denomReciprocal = 1.0 / (deltaIMinus2 + deltaIMinus1 + deltaI);\r\n      alpha[index] = deltaI * deltaI * denomReciprocal;\r\n      beta[index] = deltaI * (deltaIMinus2 + deltaIMinus1) * denomReciprocal;\r\n\r\n      denomReciprocal = 1.0 / (deltaIMinus1 + deltaI + deltaIPlus1);\r\n      beta[index] += deltaIMinus1 * (deltaI + deltaIPlus1) * denomReciprocal;\r\n      gamma[index] = deltaIMinus1 * deltaIMinus1 * denomReciprocal;\r\n\r\n      denomReciprocal = 1.0 / (deltaIMinus1 + deltaI);\r\n      alpha[index] *= denomReciprocal;\r\n      beta[index] *= denomReciprocal;\r\n      gamma[index] *= denomReciprocal;\r\n    }\r\n\r\n    /** Setup tridiagonal system for 2 fit points\r\n     * @param alpha sub-diagonal, length = 2\r\n     * @param beta diagonal, length = 2\r\n     * @param gamma super-diagonal, length = 2\r\n     */\r\n    private static setUpSystem2Points(alpha: number[], beta: number[], gamma: number[]): boolean {\r\n      if (alpha.length !== 2 || beta.length !== 2 || gamma.length !== 2)\r\n        return false;\r\n\r\n      // identity matrix\r\n      alpha[0] = alpha[1] = gamma[0] = gamma[1] = 0.0;\r\n      beta[0] = beta[1] = 1.0;\r\n      return true;\r\n    }\r\n\r\n    /** Setup tridiagonal system for 3 fit points\r\n     * @param alpha sub-diagonal, length = 3\r\n     * @param beta diagonal, length = 3\r\n     * @param gamma super-diagonal, length = 3\r\n     * @param options validated as per validateOptions, unmodified\r\n     * @param useNaturalStartTangent whether to bake the natural end condition into the first row\r\n     * @param useNaturalEndTangent whether to bake the natural end condition into the last row\r\n     */\r\n     private static setUpSystem3Points(alpha: number[], beta: number[], gamma: number[], options: InterpolationCurve3dOptions, useNaturalStartTangent: boolean, useNaturalEndTangent: boolean): boolean {\r\n      if (undefined === options.knots)\r\n        return false;\r\n      if (alpha.length !== 3 || beta.length !== 3 || gamma.length !== 3)\r\n        return false;\r\n      if (options.knots.length !== 3 || options.fitPoints.length !== 3)\r\n        return false;\r\n\r\n      let deltaIPlus1 = 0, deltaI = 0, deltaIMinus1 = 0, deltaIMinus2 = 0, sum = 0, sumReciprocal = 0;\r\n\r\n      // first row\r\n      if (useNaturalStartTangent) {\r\n        alpha[0] = 0.0;\r\n        deltaI = options.knots[1] - options.knots[0];\r\n        deltaIPlus1 = options.knots[2] - options.knots[1];\r\n        sum = deltaI + deltaIPlus1;\r\n        sumReciprocal = 1.0 / sum;\r\n        beta[0] = (deltaI + sum) * sumReciprocal;\r\n        gamma[0] = -deltaI * sumReciprocal;\r\n      } else {\r\n        alpha[0] = gamma[0] = 0.0;\r\n        beta[0] = 1.0;\r\n      }\r\n\r\n      // middle row\r\n      deltaIMinus1 = options.knots[1] - options.knots[0];\r\n      deltaI = options.knots[2] - options.knots[1];\r\n      sumReciprocal = 1.0 / (deltaIMinus1 + deltaI);\r\n      sumReciprocal *= sumReciprocal;\r\n      alpha[1] = deltaI * deltaI * sumReciprocal;\r\n      beta[1]  = 2.0 * (deltaI * deltaIMinus1) * sumReciprocal;\r\n      gamma[1] = deltaIMinus1 * deltaIMinus1 * sumReciprocal;\r\n\r\n      // last row\r\n      if (useNaturalEndTangent) {\r\n        deltaIMinus1 = options.knots[2] - options.knots[1];\r\n        deltaIMinus2 = options.knots[1] - options.knots[0];\r\n        sum = deltaIMinus2 + deltaIMinus1;\r\n        sumReciprocal = 1.0 / sum;\r\n        alpha[2] = -deltaIMinus1 * sumReciprocal;\r\n        beta[2] = (deltaIMinus1 + sum) * sumReciprocal;\r\n        gamma[2] = 0.0;\r\n      } else {\r\n        alpha[2] = gamma[2] = 0.0;\r\n        beta[2] = 1.0;\r\n      }\r\n      return true;\r\n    }\r\n\r\n    /** Setup tridiagonal system for 4 fit points or more\r\n     * @param alpha sub-diagonal, length = # fit points\r\n     * @param beta diagonal, length = # fit points\r\n     * @param gamma super-diagonal, length = # fit points\r\n     * @param options validated as per validateOptions, unmodified\r\n     * @param useNaturalStartTangent whether to bake the natural end condition into the first row\r\n     * @param useNaturalEndTangent whether to bake the natural end condition into the last row\r\n     */\r\n    private static setUpSystem4PointsOrMore(alpha: number[], beta: number[], gamma: number[], options: InterpolationCurve3dOptions, useNaturalStartTangent: boolean, useNaturalEndTangent: boolean): boolean {\r\n      if (undefined === options.knots)\r\n        return false;\r\n      if (alpha.length !== beta.length || alpha.length !== gamma.length || alpha.length !== options.knots.length)\r\n        return false;\r\n      if (options.knots.length !== options.fitPoints.length)\r\n        return false;\r\n\r\n      const numIntervals = options.fitPoints.length - 1;\r\n      const numIntervalsMinus1 = numIntervals - 1;\r\n      let deltaIPlus1 = 0, deltaI = 0, deltaIMinus1 = 0, deltaIMinus2 = 0, sum = 0, sumReciprocal = 0;\r\n\r\n      if (options.closed) {\r\n        // first row\r\n        deltaI = options.knots[1] - options.knots[0];\r\n        deltaIMinus2 = options.knots[numIntervalsMinus1] - options.knots[numIntervalsMinus1 - 1];\r\n        deltaIMinus1 = options.knots[numIntervalsMinus1 + 1] - options.knots[numIntervalsMinus1];\r\n        deltaIPlus1 = options.knots[2] - options.knots[1];\r\n        this.computeAlphaBetaGamma(alpha, beta, gamma, 0, deltaIPlus1, deltaI, deltaIMinus1, deltaIMinus2);\r\n\r\n        // second row\r\n        deltaIMinus2 = deltaIMinus1;\r\n        deltaIMinus1 = deltaI;\r\n        deltaI = options.knots[2] - options.knots[1];\r\n        deltaIPlus1 = options.knots[3] - options.knots[2];\r\n        this.computeAlphaBetaGamma(alpha, beta, gamma, 1, deltaIPlus1, deltaI, deltaIMinus1, deltaIMinus2);\r\n\r\n        // last row; there's one less equation than open case\r\n        deltaIPlus1 = deltaIMinus1;\r\n        deltaI = options.knots[numIntervalsMinus1 + 1] - options.knots[numIntervalsMinus1];\r\n        deltaIMinus2 = options.knots[numIntervalsMinus1 - 1] - options.knots[numIntervalsMinus1 - 2];\r\n        deltaIMinus1 = options.knots[numIntervalsMinus1] - options.knots[numIntervalsMinus1 - 1];\r\n        this.computeAlphaBetaGamma(alpha, beta, gamma, numIntervalsMinus1, deltaIPlus1, deltaI, deltaIMinus1, deltaIMinus2);\r\n      } else { // open\r\n        // first row\r\n        if (useNaturalStartTangent) {\r\n          alpha[0] = 0.0;\r\n          deltaI = options.knots[1] - options.knots[0];\r\n          deltaIPlus1 = options.knots[2] - options.knots[1];\r\n          sum = deltaI + deltaIPlus1;\r\n          sumReciprocal = 1.0 / sum;\r\n          beta[0] = (deltaI + sum) * sumReciprocal;\r\n          gamma[0] = -deltaI * sumReciprocal;\r\n        } else {\r\n          alpha[0] = gamma[0] = 0.0;\r\n          beta[0] = 1.0;\r\n        }\r\n\r\n        // second row\r\n        deltaI = options.knots[2] - options.knots[1];\r\n        deltaIMinus1 = options.knots[1] - options.knots[0];\r\n        deltaIMinus2 = 0.0;\r\n        deltaIPlus1 = options.knots[3] - options.knots[2];\r\n        this.computeAlphaBetaGamma(alpha, beta, gamma, 1, deltaIPlus1, deltaI, deltaIMinus1, deltaIMinus2);\r\n\r\n        // penultimate row\r\n        deltaI = options.knots[numIntervalsMinus1 + 1] - options.knots[numIntervalsMinus1];\r\n        deltaIMinus1 = options.knots[numIntervalsMinus1] - options.knots[numIntervalsMinus1 - 1];\r\n        deltaIMinus2 = options.knots[numIntervalsMinus1 - 1] - options.knots[numIntervalsMinus1 - 2];\r\n        deltaIPlus1 = 0.0;\r\n        this.computeAlphaBetaGamma(alpha, beta, gamma, numIntervalsMinus1, deltaIPlus1, deltaI, deltaIMinus1, deltaIMinus2);\r\n\r\n        // last row\r\n        if (useNaturalEndTangent) {\r\n          deltaIMinus1 = options.knots[numIntervals] - options.knots[numIntervals - 1];\r\n          deltaIMinus2 = options.knots[numIntervals - 1] - options.knots[numIntervals - 2];\r\n          sum = deltaIMinus2 + deltaIMinus1;\r\n          sumReciprocal = 1.0 / sum;\r\n          alpha[numIntervals] = -deltaIMinus1 * sumReciprocal;\r\n          beta[numIntervals] = (deltaIMinus1 + sum) * sumReciprocal;\r\n          gamma[numIntervals] = 0.0;\r\n        } else {\r\n          alpha[numIntervals] = gamma[numIntervals] = 0.0;\r\n          beta[numIntervals] = 1.0;\r\n        }\r\n      }\r\n\r\n      // middle rows\r\n      for (let i = 2; i < numIntervalsMinus1; ++i) {\r\n        deltaI = options.knots[i + 1] - options.knots[i];\r\n        deltaIMinus2 = options.knots[i - 1] - options.knots[i - 2];\r\n        deltaIMinus1 = options.knots[i] - options.knots[i - 1];\r\n        deltaIPlus1 = options.knots[i + 2] - options.knots[i + 1];\r\n        this.computeAlphaBetaGamma(alpha, beta, gamma, i, deltaIPlus1, deltaI, deltaIMinus1, deltaIMinus2);\r\n      }\r\n      return true;\r\n    }\r\n\r\n    /** Setup tridiagonal system\r\n     * @param alpha sub-diagonal, length = # fitPoints\r\n     * @param beta diagonal, length = # fitPoints\r\n     * @param gamma super-diagonal, length = # fitPoints\r\n     * @param options validated as per validateOptions, unmodified\r\n     */\r\n    private static setUpSystem(alpha: number[], beta: number[], gamma: number[], options: InterpolationCurve3dOptions): boolean {\r\n      let useNaturalStartTangent = false;\r\n      let useNaturalEndTangent = false;\r\n      if (options.isNaturalTangents && !options.closed) {\r\n        useNaturalStartTangent = (undefined === options.startTangent);\r\n        useNaturalEndTangent = (undefined === options.endTangent);\r\n      }\r\n\r\n      let succeeded = false;\r\n      if (2 === options.fitPoints.length)\r\n        succeeded = this.setUpSystem2Points(alpha, beta, gamma);\r\n      else if (3 === options.fitPoints.length)\r\n        succeeded = this.setUpSystem3Points(alpha, beta, gamma, options, useNaturalStartTangent, useNaturalEndTangent);\r\n      else if (4 <= options.fitPoints.length)\r\n        succeeded = this.setUpSystem4PointsOrMore(alpha, beta, gamma, options, useNaturalStartTangent, useNaturalEndTangent);\r\n\r\n      return succeeded;\r\n    }\r\n\r\n    /** Set the Bessel end condition for the linear system.\r\n     * @param dataPts array whose middle is the system rhs (augmented with first/last fitPoint at beginning/end);\r\n     *                2nd or penultimate point is set by this function.\r\n     * @param options validated as per validateOptions, unmodified\r\n     * @param atStart whether end condition is for start of curve (false: end of curve)\r\n     */\r\n    private static setBesselEndCondition(dataPts: Point3d[], options: InterpolationCurve3dOptions, atStart: boolean): boolean {\r\n      if (dataPts.length !== options.fitPoints.length + 2)\r\n        return false;\r\n      if (undefined === options.knots)\r\n        return false;\r\n\r\n      const scale = 1.0/3.0;\r\n      const numIntervals = options.fitPoints.length - 1;\r\n\r\n      if (1 === numIntervals) { // linear Bezier\r\n        if (atStart)\r\n          dataPts[0].interpolate(scale, dataPts[3], dataPts[1]);\r\n        else\r\n          dataPts[3].interpolate(scale, dataPts[0], dataPts[2]);\r\n        return true;\r\n      }\r\n\r\n      if (2 === numIntervals) {\r\n        const alpha = (options.knots[2] - options.knots[1]) / (options.knots[2] - options.knots[0]);\r\n        const beta = 1.0 - alpha;\r\n        const temp = dataPts[2].plus2Scaled(dataPts[0], -alpha * alpha, dataPts[4], -beta * beta);\r\n        if (atStart)\r\n          Point3d.createAdd2Scaled(temp, 1.0 / (2.0 * alpha), dataPts[0], alpha).interpolate(scale, dataPts[0], dataPts[1]);\r\n        else\r\n          Point3d.createAdd2Scaled(temp, 1.0 / (2.0 * beta), dataPts[4], beta).interpolate(scale, dataPts[4], dataPts[3]);\r\n        return true;\r\n      }\r\n\r\n      // numIntervals > 2\r\n      if (atStart) {\r\n          const alpha = (options.knots[2] - options.knots[1]) / (options.knots[2] - options.knots[0]);\r\n          const beta = 1.0 - alpha;\r\n          const temp = dataPts[2].plus2Scaled(dataPts[0], -alpha * alpha, dataPts[3], -beta * beta);\r\n          Point3d.createAdd2Scaled(temp, 1.0 / (2.0 * alpha), dataPts[0], alpha).interpolate(scale, dataPts[0], dataPts[1]);\r\n      } else {\r\n          const alpha = (options.knots[numIntervals] - options.knots[numIntervals - 1]) / (options.knots[numIntervals] - options.knots[numIntervals - 2]);\r\n          const beta = 1.0 - alpha;\r\n          const temp = dataPts[numIntervals].plus2Scaled(dataPts[numIntervals - 1], -alpha * alpha, dataPts[numIntervals + 2], -beta * beta);\r\n          Point3d.createAdd2Scaled(temp, 1.0 / (2.0 * beta), dataPts[numIntervals + 2], beta).interpolate(scale, dataPts[numIntervals + 2], dataPts[numIntervals + 1]);\r\n      }\r\n      return true;\r\n    }\r\n\r\n    /** Set the natural end condition for the linear system.\r\n     *  This is the end condition used by ADSK for fit-splines with a given zero tangent.\r\n     * @param dataPts array whose middle is the system rhs (augmented with first/last fitPoint at beginning/end);\r\n     *                2nd or penultimate point is set by this function.\r\n     * @param options validated as per validateOptions, unmodified\r\n     * @param atStart whether end condition is for start of curve (false: end of curve)\r\n     */\r\n    private static setNaturalEndCondition(dataPts: Point3d[], options: InterpolationCurve3dOptions, atStart: boolean): boolean {\r\n      if (dataPts.length !== options.fitPoints.length + 2)\r\n        return false;\r\n\r\n      const numIntervals = options.fitPoints.length - 1;\r\n      if (1 === numIntervals)\r\n        return this.setBesselEndCondition(dataPts, options, atStart);\r\n\r\n      if (atStart)\r\n        dataPts[1] = dataPts[0];\r\n      else\r\n        dataPts[dataPts.length - 2] = dataPts[dataPts.length - 1];\r\n      return true;\r\n    }\r\n\r\n    /** Set the end condition for the linear system to the given tangent, scaled by chord length.\r\n     *  This is the end condition used by ADSK for fit-splines with a given nonzero tangent.\r\n     * @param dataPts array whose middle is the system rhs (augmented with first/last fitPoint at beginning/end);\r\n     *                2nd or penultimate point is set by this function.\r\n     * @param options validated as per validateOptions, unmodified\r\n     * @param atStart whether end condition is for start of curve (false: end of curve)\r\n     */\r\n    private static setChordLengthScaledEndCondition(dataPts: Point3d[], options: InterpolationCurve3dOptions, atStart: boolean): boolean {\r\n      if (dataPts.length !== options.fitPoints.length + 2)\r\n        return false;\r\n\r\n      const tangent = atStart ? options.startTangent : options.endTangent;\r\n      if (undefined === tangent)\r\n        return false;\r\n\r\n      let iExt = 0; // index of first/last fitPoint\r\n      let iSet = 0; // index of 2nd/penultimate Bezier point to set (determines start/end tangent of the curve)\r\n      let iInt = 0; // index of 2nd/penultimate fitPoint\r\n\r\n      const numIntervals = options.fitPoints.length - 1;\r\n      if (1 === numIntervals) { // no interior fit points\r\n        if (atStart) {\r\n          iExt = 0;\r\n          iSet = 1;\r\n          iInt = 3;\r\n        } else {\r\n          iExt = 3;\r\n          iSet = 2;\r\n          iInt = 0;\r\n        }\r\n      } else {\r\n        if (atStart) {\r\n          iExt = 0;\r\n          iSet = 1;\r\n          iInt = 2;\r\n        } else {\r\n          iExt = numIntervals + 2;\r\n          iSet = numIntervals + 1;\r\n          iInt = numIntervals;\r\n        }\r\n      }\r\n\r\n      // NOTE: tangent points INTO curve\r\n      const chordLength = dataPts[iInt].distance(dataPts[iExt]);\r\n      dataPts[iExt].plusScaled(tangent, chordLength / 3.0, dataPts[iSet]);\r\n      return true;\r\n    }\r\n\r\n    /** Set the end condition for the linear system to the given tangent, scaled by bessel length.\r\n     * @param dataPts array whose middle is the system rhs (augmented with first/last fitPoint at beginning/end);\r\n     *                2nd or penultimate point is set by this function.\r\n     * @param options validated as per validateOptions, unmodified\r\n     * @param atStart whether end condition is for start of curve (false: end of curve)\r\n     */\r\n    private static setBesselLengthScaledEndCondition(dataPts: Point3d[], options: InterpolationCurve3dOptions, atStart: boolean): boolean {\r\n      if (dataPts.length !== options.fitPoints.length + 2)\r\n        return false;\r\n\r\n      const tangent = atStart ? options.startTangent : options.endTangent;\r\n      if (undefined === tangent)\r\n        return false;\r\n\r\n      // temporarily set bessel end condition\r\n      if (!this.setBesselEndCondition(dataPts, options, atStart))\r\n        return false;\r\n\r\n      const numIntervals = options.fitPoints.length - 1;\r\n      const iExt = atStart ? 0 : numIntervals + 2; // index of first/last fitPoint\r\n      const iSet = atStart ? 1 : numIntervals + 1; // index of 2nd/penultimate Bezier point to set (determines start/end tangent of the curve)\r\n\r\n      // reset end condition with our tangent, but scaled to the bessel tangent's length\r\n      dataPts[iExt].plusScaled(tangent, dataPts[iExt].distance(dataPts[iSet]), dataPts[iSet]);\r\n      return true;\r\n    }\r\n\r\n    /** Set the end condition for a physically closed (non-periodic) interpolant.\r\n     * @param dataPts array whose middle is the system rhs (augmented with first/last fitPoint at beginning/end);\r\n     *                2nd or penultimate point is set by this function.\r\n     * @param options validated as per validateOptions, unmodified\r\n     */\r\n    private static setPhysicallyClosedEndCondition(dataPts: Point3d[], options: InterpolationCurve3dOptions): boolean {\r\n      const numIntervals = options.fitPoints.length - 1;\r\n      if (!options.isColinearTangents\r\n          || numIntervals <= 2\r\n          || (undefined !== options.startTangent && undefined !== options.endTangent)\r\n          || options.isNaturalTangents\r\n          || !dataPts[0].isAlmostEqual(dataPts[numIntervals + 2])) {\r\n        return true;\r\n        }\r\n      // force parallel start/end tangents, using chord length scale for undefined tangents\r\n      if (undefined !== options.startTangent) { // start tangent is supplied; compute a parallel end tangent\r\n        const outwardStartTangent = Vector3d.createStartEnd(dataPts[1], dataPts[0]).normalize();\r\n        if (undefined !== outwardStartTangent) {\r\n          const endTangentMag = dataPts[numIntervals + 2].distance(dataPts[numIntervals + 1]);\r\n          dataPts[numIntervals + 2].plusScaled(outwardStartTangent, endTangentMag, dataPts[numIntervals + 1]);\r\n        }\r\n      } else if (undefined !== options.endTangent) {  // end tangent is supplied; compute a parallel start tangent\r\n        const outwardEndTangent = Vector3d.createStartEnd(dataPts[numIntervals + 1], dataPts[numIntervals + 2]).normalize();\r\n        if (undefined !== outwardEndTangent) {\r\n          const startTangentMag = dataPts[0].distance(dataPts[1]);\r\n          dataPts[0].plusScaled(outwardEndTangent, startTangentMag, dataPts[1]);\r\n        }\r\n      } else {  // neither tangent is supplied, compute both along same vector\r\n        const commonTangent = Vector3d.createStartEnd(dataPts[numIntervals + 1], dataPts[1]).normalize();\r\n        if (undefined !== commonTangent) {\r\n          const startTangentMag = dataPts[0].distance(dataPts[1]);\r\n          dataPts[0].plusScaled(commonTangent, startTangentMag, dataPts[1]);\r\n          const endTangentMag = dataPts[numIntervals + 2].distance(dataPts[numIntervals + 1]);\r\n          dataPts[numIntervals + 2].plusScaled(commonTangent, -endTangentMag, dataPts[numIntervals + 1]);\r\n        }\r\n      }\r\n      return true;\r\n    }\r\n\r\n    /** Set end conditions for the linear system to solve for the poles of the open interpolant, as per Farin 3e/4e.\r\n     * @param dataPts array whose interior is the system rhs and whose first/last entries are the first/last fitPoints;\r\n     *                points are inserted to become the 2nd and penultimate dataPts, the first/last rows of the system rhs.\r\n     * @param options validated as per validateOptions, unmodified\r\n     */\r\n    private static setEndConditions(dataPts: Point3d[], options: InterpolationCurve3dOptions): boolean {\r\n      if (dataPts.length !== options.fitPoints.length)\r\n        return false;\r\n\r\n      // insert dummy points to be computed below\r\n      const dummy0 = Point3d.createZero();\r\n      const dummy1 = Point3d.createZero();\r\n      dataPts.splice(1, 0, dummy0);\r\n      dataPts.splice(dataPts.length - 1, 0, dummy1);\r\n\r\n      let succeeded = false;\r\n      if (undefined === options.startTangent) {\r\n        if (options.isNaturalTangents)\r\n          succeeded = this.setNaturalEndCondition(dataPts, options, true);\r\n        else\r\n          succeeded = this.setBesselEndCondition(dataPts, options, true);\r\n      } else { // scale startTangent\r\n        if (options.isChordLenTangents)\r\n          succeeded = this.setChordLengthScaledEndCondition(dataPts, options, true);\r\n        else\r\n          succeeded = this.setBesselLengthScaledEndCondition(dataPts, options, true);\r\n      }\r\n\r\n      if (undefined === options.endTangent) {\r\n        if (options.isNaturalTangents)\r\n          succeeded = this.setNaturalEndCondition(dataPts, options, false);\r\n        else\r\n          succeeded = this.setBesselEndCondition(dataPts, options, false);\r\n      } else { // scale endTangent\r\n        if (options.isChordLenTangents)\r\n          succeeded = this.setChordLengthScaledEndCondition(dataPts, options, false);\r\n        else\r\n          succeeded = this.setBesselLengthScaledEndCondition(dataPts, options, false);\r\n      }\r\n\r\n      if (succeeded)\r\n        succeeded = this.setPhysicallyClosedEndCondition(dataPts, options);\r\n\r\n      return succeeded;\r\n    }\r\n\r\n    /** Solve the near tridiagonal system for a periodic C2 cubic interpolant.\r\n     * Alpha, beta, gamma are computed by setUpSystem, have same length as fitPts, and are overwritten.\r\n     */\r\n    private static solveNearTridiagonal(fitPts: Point3d[], alpha: number[], beta: number[], gamma: number[]): Point3d[] | undefined {\r\n      if (alpha.length !== beta.length || alpha.length !== gamma.length || alpha.length !== fitPts.length)\r\n        return undefined;\r\n      const poles: Point3d[] = [];\r\n      const numIntervals = fitPts.length - 1;\r\n      const leftPts = fitPts.slice(0, -1);  // last fitPt is ignored\r\n      let tmp: number | undefined = 0.0;\r\n\r\n      // first forward substitution\r\n      for (let i = 1; i < numIntervals; ++i) {\r\n        if (undefined === (tmp = Geometry.conditionalDivideFraction(- alpha[i], beta[i - 1])))\r\n          return undefined;\r\n        beta[i] += tmp * gamma[i - 1];\r\n        alpha[i] = tmp * alpha[i - 1];\r\n        leftPts[i].addScaledInPlace(leftPts[i - 1], tmp);\r\n      }\r\n\r\n      // first backward substitution\r\n      if (undefined === (tmp = Geometry.conditionalDivideFraction(1.0, beta[numIntervals - 1] + alpha[numIntervals - 1])))\r\n        return undefined;\r\n      gamma[numIntervals - 1] *= tmp;\r\n      leftPts[numIntervals - 1].scaleInPlace(tmp);\r\n      for (let i = numIntervals - 2; i >= 0; --i) {\r\n        if (undefined === (tmp = Geometry.conditionalDivideFraction(1.0, beta[i])))\r\n          return undefined;\r\n        Point3d.createScale(leftPts[i].plus2Scaled(leftPts[i + 1], - gamma[i], leftPts[numIntervals - 1], - alpha[i]), tmp, leftPts[i]);\r\n        gamma[i] = - (gamma[i] * gamma[i + 1] + alpha[i] * gamma[numIntervals - 1]) * tmp;\r\n      }\r\n\r\n      // second forward substitution\r\n      if (undefined === (tmp = Geometry.conditionalDivideFraction(1.0, 1.0 + gamma[0])))\r\n        return undefined;\r\n      poles.push(Point3d.createScale(leftPts[0], tmp));\r\n      for (let i = 1; i < numIntervals; ++i) {\r\n        poles.push(leftPts[i].plusScaled(poles[0], - gamma[i]));\r\n      }\r\n      return poles;\r\n    }\r\n\r\n    /** Adjust options by correcting invalid combinations\r\n     * @param options curve definition, possibly modified\r\n     */\r\n    public static validateOptions(options: InterpolationCurve3dOptions): boolean {\r\n      options.order = 4;\r\n\r\n      // remove relevant exterior knots so knots and fit points align *before* we start compressing fit points.\r\n      options.knots = this.convertCubicKnotVectorToFitParams(options.knots, options.fitPoints.length, true);\r\n\r\n        // compress out duplicate fit points (and their corresponding knots)\r\n      if (!this.removeDuplicateFitPoints(options))\r\n        return false;\r\n\r\n      // if only 2 unique points, then must create open curve\r\n      let hasClosurePoint = options.fitPoints[0].isAlmostEqual(options.fitPoints[options.fitPoints.length - 1]);\r\n      if (3 === options.fitPoints.length && hasClosurePoint) {\r\n        options.fitPoints.pop();\r\n        if (undefined !== options.knots)\r\n          options.knots.pop();\r\n        hasClosurePoint = options.fitPoints[0].isAlmostEqual(options.fitPoints[options.fitPoints.length - 1]);\r\n      }\r\n      if (options.fitPoints.length <= 2) {\r\n        if (hasClosurePoint)\r\n          return false;\r\n        options.closed = false;\r\n      }\r\n\r\n      // append closure point if missing\r\n      if (options.closed) {\r\n        if (!hasClosurePoint) {\r\n          options.fitPoints.push(options.fitPoints[0].clone());\r\n          if (undefined !== options.knots) {  // best guess: uniform knots\r\n            options.knots.push(options.knots[options.knots.length - 1] + (options.knots[options.knots.length - 1] - options.knots[0]) / (options.knots.length - 1));\r\n          }\r\n        }\r\n        if (options.fitPoints.length <= 4)\r\n          options.closed = false; // can't fit cubic closed curve to 3 unique points\r\n      }\r\n\r\n      if (options.fitPoints.length < 2)\r\n        return false;\r\n\r\n      // ASSUME: tangents point INTO curve\r\n      if (undefined !== options.startTangent) {\r\n        if (options.startTangent.isAlmostZero)\r\n          options.startTangent = undefined;\r\n        else\r\n          options.startTangent.normalizeInPlace();\r\n      }\r\n      if (undefined !== options.endTangent) {\r\n        if (options.endTangent.isAlmostZero)\r\n          options.endTangent = undefined;\r\n        else\r\n          options.endTangent.normalizeInPlace();\r\n      }\r\n\r\n      return true;\r\n    }\r\n\r\n    /** Converts a full cubic knot vector of expected length into fit parameters, by removing extraneous exterior knots.\r\n     * @param knots cubic knot vector, unmodified\r\n     * @param numFitPoints number of fit points\r\n     * @return fit parameters, or undefined if unexpected input\r\n     **/\r\n    public static convertCubicKnotVectorToFitParams(knots: number[] | undefined, numFitPoints: number, normalize?: boolean): number[] | undefined {\r\n      let params = knots?.slice();\r\n      if (undefined !== params) {\r\n        const numExtraKnots = params.length - numFitPoints;\r\n        switch (numExtraKnots) {\r\n          case 0: {   // ASSUME caller passed in interior knots\r\n            break;\r\n          }\r\n          case 4:     // modern full cubic knots\r\n          case 6: {   // legacy full cubic knots\r\n            for (let i = 0; i < numExtraKnots / 2; ++i) {\r\n              params.pop();\r\n              params.shift();\r\n            }\r\n            break;\r\n          }\r\n          default: {  // other knot configurations are unusable\r\n            params = undefined;\r\n            break;\r\n          }\r\n        }\r\n        if (normalize && !this.normalizeKnots(params))\r\n          params = undefined;\r\n      }\r\n      return params;\r\n    }\r\n\r\n    /** Return fit parameters augmented to a full cubic knot vector.\r\n     * @param params fit parameters, unmodified\r\n     * @param legacy whether to create a legacy (DGN) full knot vector, or modern vector with two less knots\r\n     **/\r\n    public static convertFitParamsToCubicKnotVector(params: number[] | undefined, closed?: boolean, legacy?: boolean): number[] | undefined {\r\n      const knots = params?.slice();\r\n      if (undefined !== knots) {\r\n        const numExtraKnots = legacy ? 6 : 4;\r\n        if (closed) {\r\n          const iTail = knots.length - 2;\r\n          for (let iHead = 2; iHead <= numExtraKnots; iHead += 2) {\r\n            knots.unshift(knots[iTail] - 1.0);  // index is constant\r\n            knots.push(1.0 + knots[iHead]);     // index increments by two\r\n          }\r\n        } else {\r\n          for (let i = 0; i < numExtraKnots / 2; ++i) {\r\n            knots.unshift(0.0);\r\n            knots.push(1.0);\r\n          }\r\n        }\r\n      }\r\n      return knots;\r\n    }\r\n\r\n    /** Ensure full legacy knot vector for JSON export **/\r\n    public static convertToJsonKnots(props: InterpolationCurve3dProps) {\r\n      if (undefined !== props.knots) {\r\n        props.knots = this.convertCubicKnotVectorToFitParams(props.knots, props.fitPoints.length, false);\r\n        props.knots = this.convertFitParamsToCubicKnotVector(props.knots, props.closed, true);\r\n      } else {\r\n        props.knots = this.constructFitParametersFromPoints(Point3dArray.clonePoint3dArray(props.fitPoints), props.isChordLenKnots, props.closed);\r\n        props.knots = this.convertFitParamsToCubicKnotVector(props.knots, props.closed, true);\r\n      }\r\n    }\r\n\r\n    /** Construct the control points for the c2 cubic fit algorithm\r\n     * @param options validated as per validateOptions, possibly modified\r\n     */\r\n    public static constructPoles(options: InterpolationCurve3dOptions): Point3d[] | Float64Array | undefined {\r\n      if (!this.constructFitParameters(options) || (undefined === options.knots))\r\n        return undefined;\r\n\r\n      const numRow = options.fitPoints.length;\r\n      const alpha: number[] = Array(numRow);\r\n      const beta: number[] = Array(numRow);\r\n      const gamma: number[] = Array(numRow);\r\n      if (!this.setUpSystem(alpha, beta, gamma, options))\r\n        return undefined;\r\n\r\n      let poles: Point3d[] | Float64Array | undefined = [];\r\n      if (!options.closed) {\r\n        const dataPts = options.fitPoints.slice();\r\n        if (!this.setEndConditions(dataPts, options))\r\n          return undefined;\r\n        if (dataPts.length !== numRow + 2)\r\n          return undefined; // sanity check: we added 2nd/penultimate points, rhs is middle numRow entries.\r\n\r\n        // construct tridiagonal banded system components\r\n        const matrix = new Float64Array(numRow * 3);\r\n        const rhs = new Float64Array(numRow * 3);\r\n        for (let iRow = 0, iMatrixRead = 0, iRhsRead = 0; iRow < numRow; ++iRow) {\r\n          matrix[iMatrixRead++] = alpha[iRow];\r\n          matrix[iMatrixRead++] = beta[iRow];\r\n          matrix[iMatrixRead++] = gamma[iRow];\r\n          rhs[iRhsRead++] = dataPts[iRow+1].x;\r\n          rhs[iRhsRead++] = dataPts[iRow+1].y;\r\n          rhs[iRhsRead++] = dataPts[iRow+1].z;\r\n        }\r\n\r\n        const solution = BandedSystem.solveBandedSystemMultipleRHS(numRow, 3, matrix, 3, rhs);\r\n        if (undefined === solution)\r\n          return undefined;\r\n\r\n        // pre/append first/last poles/fitPoints\r\n        poles = new Float64Array(3 + solution.length + 3);\r\n        let iWrite = 0;\r\n        poles[iWrite++] = options.fitPoints[0].x;\r\n        poles[iWrite++] = options.fitPoints[0].y;\r\n        poles[iWrite++] = options.fitPoints[0].z;\r\n        for (let iRead = 0; iRead < solution.length; ) {\r\n          poles[iWrite++] = solution[iRead++];\r\n        }\r\n        poles[iWrite++] = options.fitPoints[options.fitPoints.length - 1].x;\r\n        poles[iWrite++] = options.fitPoints[options.fitPoints.length - 1].y;\r\n        poles[iWrite++] = options.fitPoints[options.fitPoints.length - 1].z;\r\n      } else { // closed\r\n        if (undefined !== (poles = this.solveNearTridiagonal(options.fitPoints, alpha, beta, gamma))) {\r\n          if (poles.length > 2) {\r\n            poles.unshift(poles.pop()!);  // shift poles right to line up with the knots\r\n            for (let i = 0; i < options.order - 1; ++i)\r\n              poles.push(poles[i].clone()); // periodically extend (the modern way)\r\n          }\r\n        }\r\n      }\r\n      return poles;\r\n    }\r\n  }\r\n}\r\n"]}

package/lib/esm/bspline/Bezier1dNd.d.ts

@@ -64,42 +64,39 @@ export declare class Bezier1dNd {
      * @param poleIndexB second pole index
      */
     interpolatePoleInPlace(poleIndexA: number, fraction: number, poleIndexB: number): void;
-    private static _knotTolerance;
     /**
      * Compute new control points to "clamp" bspline unsaturated support to saturated form.
      * * At input time, the control points are associated with the input knots (unsaturated)
-     * * At output, they control points are modified by repeated knot insertion to be fully clamped.
+     * * At output, the control points are modified by repeated knot insertion to be fully clamped.
      * @param knots knot values for the current (unsaturated) pole set
      * @param spanIndex index of span whose (unsaturated) poles are in the bezier.
-     * @param optional function for `setInterval (knotA, knotB)` call to announce knot limits.
      */
     saturateInPlace(knots: KnotVector, spanIndex: number): boolean;
     /**
      * Saturate a univariate bspline coefficient array in place
-     * * On input, the array is the coefficients one span of a bspline, packed in an array of `(knots.order)` values.
+     * * On input, the array is the coefficients of one span of a bspline, packed in an array of `(knots.order)` values.
      * * These are modified in place, and on return are a bezier for the same knot interval.
      * @param coffs input as bspline coefficients, returned as bezier coefficients
      * @param knots knot vector
      * @param spanIndex index of span whose (unsaturated) poles are in the coefficients.
-     * @param optional function for `setInterval (knotA, knotB)` call to announce knot limits.
      */
     static saturate1dInPlace(coffs: Float64Array, knots: KnotVector, spanIndex: number): boolean;
     /**
-     * Apply deCasteljou interpolations to isolate a smaller bezier polygon, representing interval 0..fraction of the original
+     * Apply deCasteljau interpolations to isolate a smaller bezier polygon, representing interval 0..fraction of the original
      * @param fraction "end" fraction for split.
-     * @returns false if fraction is 0 -- no changes applied.
+     * @returns false if fraction is 0 or 1 -- no changes applied.
      */
     subdivideInPlaceKeepLeft(fraction: number): boolean;
     /**
-     * Apply deCasteljou interpolations to isolate a smaller bezier polygon, representing interval 0..fraction of the original
-     * @param fraction "end" fraction for split.
-     * @returns false if fraction is 0 -- no changes applied.
+     * Apply deCasteljau interpolations to isolate a smaller bezier polygon, representing interval fraction..1 of the original
+     * @param fraction "start" fraction for split.
+     * @returns false if fraction is 0 or 1 -- no changes applied.
      */
     subdivideInPlaceKeepRight(fraction: number): boolean;
     /**
      * Saturate a univariate bspline coefficient array in place
      * @param fraction0 fraction for first split.   This is the start of the output polygon
-     * @param fraction1 fraction for first split.   This is the start of the output polygon
+     * @param fraction1 fraction for second split.   This is the end of the output polygon
      * @return false if fractions are (almost) identical.
      */
     subdivideToIntervalInPlace(fraction0: number, fraction1: number): boolean;

package/lib/esm/bspline/Bezier1dNd.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"Bezier1dNd.d.ts","sourceRoot":"","sources":["../../../src/bspline/Bezier1dNd.ts"],"names":[],"mappings":"AAIA;;GAEG;AAGH,OAAO,EAAE,OAAO,EAAE,MAAM,+BAA+B,CAAC;AACxD,OAAO,EAAE,OAAO,EAAE,MAAM,+BAA+B,CAAC;AAExD,OAAO,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AACpD,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAEhD,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAE1C;;;;;GAKG;AACH,qBAAa,UAAU;IACrB,OAAO,CAAC,WAAW,CAAe;IAClC,OAAO,CAAC,MAAM,CAAS;IACvB,OAAO,CAAC,UAAU,CAAS;IAC3B,OAAO,CAAC,MAAM,CAAc;gBAET,SAAS,EAAE,MAAM,EAAE,OAAO,EAAE,YAAY;IAM3D,uCAAuC;IAChC,YAAY,CAAC,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IASxD,8BAA8B;IAC9B,IAAW,KAAK,WAA0B;IAC1C,uEAAuE;IACvE,IAAW,UAAU,iBAA+B;IACpD,2FAA2F;WAC7E,MAAM,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,OAAO,EAAE,GAAG,OAAO,EAAE,GAAG,UAAU,GAAG,SAAS;IAiCrF;;OAEG;IACI,QAAQ,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAG/D;;OAEG;IACI,kBAAkB,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAGzE,4DAA4D;IACrD,eAAe,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY,GAAG,SAAS;IAWlF,4DAA4D;IACrD,eAAe,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,EAAE,YAAY;IAOtD;;;OAGG;IACI,aAAa,CAAC,IAAI,EAAE,YAAY,EAAE,SAAS,EAAE,MAAM;IAK1D;;;;;OAKG;IACI,uBAAuB,CAAC,IAAI,EAAE,YAAY,EAAE,aAAa,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM;IAU3G,+FAA+F;IACxF,kBAAkB,IAAI,GAAG,EAAE;IAGlC,iDAAiD;IAC1C,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IAgBzC,8BAA8B;IACvB,cAAc;IAcrB;;;;;;OAMG;IACI,sBAAsB,CAAC,UAAU,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM;IAQtF,OAAO,CAAC,MAAM,CAAC,cAAc,CAAU;IACvC;;;;;;;OAOG;IACI,eAAe,CAAC,KAAK,EAAE,UAAU,EAAE,SAAS,EAAE,MAAM,GAAG,OAAO;IAsCrE;;;;;;;;OAQG;WACW,iBAAiB,CAAC,KAAK,EAAE,YAAY,EAAE,KAAK,EAAE,UAAU,EAAE,SAAS,EAAE,MAAM,GAAG,OAAO;IAqCnG;;;;OAIG;IACI,wBAAwB,CAAC,QAAQ,EAAE,MAAM,GAAG,OAAO;IAe1D;;;;OAIG;IACI,yBAAyB,CAAC,QAAQ,EAAE,MAAM,GAAG,OAAO;IAa3D;;;;;OAKG;IACI,0BAA0B,CAAC,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,GAAG,OAAO;IAahF,uDAAuD;IAChD,QAAQ,CAAC,EAAE,SAAS,CAAC;IAC5B,oDAAoD;IAC7C,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM;IAGvC,0CAA0C;IACnC,wBAAwB,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;CAC1D"}
+{"version":3,"file":"Bezier1dNd.d.ts","sourceRoot":"","sources":["../../../src/bspline/Bezier1dNd.ts"],"names":[],"mappings":"AAIA;;GAEG;AAGH,OAAO,EAAE,OAAO,EAAE,MAAM,+BAA+B,CAAC;AACxD,OAAO,EAAE,OAAO,EAAE,MAAM,+BAA+B,CAAC;AAExD,OAAO,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AACpD,OAAO,EAAE,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAEhD,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAE1C;;;;;GAKG;AACH,qBAAa,UAAU;IACrB,OAAO,CAAC,WAAW,CAAe;IAClC,OAAO,CAAC,MAAM,CAAS;IACvB,OAAO,CAAC,UAAU,CAAS;IAC3B,OAAO,CAAC,MAAM,CAAc;gBAET,SAAS,EAAE,MAAM,EAAE,OAAO,EAAE,YAAY;IAM3D,uCAAuC;IAChC,YAAY,CAAC,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IASxD,8BAA8B;IAC9B,IAAW,KAAK,WAA0B;IAC1C,uEAAuE;IACvE,IAAW,UAAU,iBAA+B;IACpD,2FAA2F;WAC7E,MAAM,CAAC,IAAI,EAAE,OAAO,EAAE,GAAG,OAAO,EAAE,GAAG,OAAO,EAAE,GAAG,UAAU,GAAG,SAAS;IAiCrF;;OAEG;IACI,QAAQ,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAG/D;;OAEG;IACI,kBAAkB,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY;IAGzE,4DAA4D;IACrD,eAAe,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,EAAE,YAAY,GAAG,YAAY,GAAG,SAAS;IAWlF,4DAA4D;IACrD,eAAe,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,EAAE,YAAY;IAOtD;;;OAGG;IACI,aAAa,CAAC,IAAI,EAAE,YAAY,EAAE,SAAS,EAAE,MAAM;IAK1D;;;;;OAKG;IACI,uBAAuB,CAAC,IAAI,EAAE,YAAY,EAAE,aAAa,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM;IAU3G,+FAA+F;IACxF,kBAAkB,IAAI,GAAG,EAAE;IAGlC,iDAAiD;IAC1C,aAAa,CAAC,KAAK,EAAE,GAAG,GAAG,OAAO;IAgBzC,8BAA8B;IACvB,cAAc;IAcrB;;;;;;OAMG;IACI,sBAAsB,CAAC,UAAU,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM;IAStF;;;;;;OAMG;IACI,eAAe,CAAC,KAAK,EAAE,UAAU,EAAE,SAAS,EAAE,MAAM,GAAG,OAAO;IAsCrE;;;;;;;OAOG;WACW,iBAAiB,CAAC,KAAK,EAAE,YAAY,EAAE,KAAK,EAAE,UAAU,EAAE,SAAS,EAAE,MAAM,GAAG,OAAO;IAuCnG;;;;OAIG;IACI,wBAAwB,CAAC,QAAQ,EAAE,MAAM,GAAG,OAAO;IAe1D;;;;OAIG;IACI,yBAAyB,CAAC,QAAQ,EAAE,MAAM,GAAG,OAAO;IAa3D;;;;;OAKG;IACI,0BAA0B,CAAC,SAAS,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,GAAG,OAAO;IAahF,uDAAuD;IAChD,QAAQ,CAAC,EAAE,SAAS,CAAC;IAC5B,oDAAoD;IAC7C,WAAW,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM;IAGvC,0CAA0C;IACnC,wBAAwB,CAAC,QAAQ,EAAE,MAAM,GAAG,MAAM;CAC1D"}

package/lib/esm/bspline/Bezier1dNd.js

@@ -12,6 +12,7 @@ import { Point3dArray } from "../geometry3d/PointHelpers";
 import { Segment1d } from "../geometry3d/Segment1d";
 import { Point4d } from "../geometry4d/Point4d";
 import { UnivariateBezier } from "../numerics/BezierPolynomials";
+import { KnotVector } from "./KnotVector";
 /**
  * Shared implementation details for derived bezier curve classes
  * * BezierCurve3d implements with blockSize 3.
@@ -168,7 +169,6 @@ export class Bezier1dNd {
             }
         }
     }
-    //
     /**
      * interpolate at `fraction` between poleA and poleB.
      * * Data is left "in place" in poleIndexA
@@ -187,10 +187,9 @@ export class Bezier1dNd {
     /**
      * Compute new control points to "clamp" bspline unsaturated support to saturated form.
      * * At input time, the control points are associated with the input knots (unsaturated)
-     * * At output, they control points are modified by repeated knot insertion to be fully clamped.
+     * * At output, the control points are modified by repeated knot insertion to be fully clamped.
      * @param knots knot values for the current (unsaturated) pole set
      * @param spanIndex index of span whose (unsaturated) poles are in the bezier.
-     * @param optional function for `setInterval (knotA, knotB)` call to announce knot limits.
      */
     saturateInPlace(knots, spanIndex) {
         const degree = knots.degree;
@@ -202,7 +201,7 @@ export class Bezier1dNd {
         const knotA = knotArray[kA];
         const knotB = knotArray[kB];
         this.setInterval(knotA, knotB);
-        if (knotB <= knotA + Bezier1dNd._knotTolerance)
+        if (knotB <= knotA + KnotVector.knotTolerance)
             return false;
         for (let numInsert = degree - 1; numInsert > 0; numInsert--) {
             //  left numInsert poles are pulled forward
@@ -232,12 +231,11 @@ export class Bezier1dNd {
     }
     /**
      * Saturate a univariate bspline coefficient array in place
-     * * On input, the array is the coefficients one span of a bspline, packed in an array of `(knots.order)` values.
+     * * On input, the array is the coefficients of one span of a bspline, packed in an array of `(knots.order)` values.
      * * These are modified in place, and on return are a bezier for the same knot interval.
      * @param coffs input as bspline coefficients, returned as bezier coefficients
      * @param knots knot vector
      * @param spanIndex index of span whose (unsaturated) poles are in the coefficients.
-     * @param optional function for `setInterval (knotA, knotB)` call to announce knot limits.
      */
     static saturate1dInPlace(coffs, knots, spanIndex) {
         const degree = knots.degree;
@@ -248,6 +246,8 @@ export class Bezier1dNd {
         const knotArray = knots.knots;
         const knotA = knotArray[kA];
         const knotB = knotArray[kB];
+        if (knotB <= knotA + KnotVector.knotTolerance)
+            return false;
         for (let numInsert = degree - 1; numInsert > 0; numInsert--) {
             //  left numInsert poles are pulled forward
             let k0 = kA - numInsert;
@@ -267,7 +267,7 @@ export class Bezier1dNd {
             if (knotArray[k2] > knotB) {
                 for (let i = 0; i < numInsert; i++, k2--) {
                     const knot2 = knotArray[k2]; // right side of moving window
-                    // left side of window ia always the (previously saturated) knotA
+                    // left side of window is always the (previously saturated) knotA
                     const fraction = (knotB - knot2) / (knotA - knot2);
                     k = degree - i;
                     coffs[k] += fraction * (coffs[k - 1] - coffs[k]);
@@ -277,9 +277,9 @@ export class Bezier1dNd {
         return true;
     }
     /**
-     * Apply deCasteljou interpolations to isolate a smaller bezier polygon, representing interval 0..fraction of the original
+     * Apply deCasteljau interpolations to isolate a smaller bezier polygon, representing interval 0..fraction of the original
      * @param fraction "end" fraction for split.
-     * @returns false if fraction is 0 -- no changes applied.
+     * @returns false if fraction is 0 or 1 -- no changes applied.
      */
     subdivideInPlaceKeepLeft(fraction) {
         if (Geometry.isAlmostEqualNumber(fraction, 1.0))
@@ -296,9 +296,9 @@ export class Bezier1dNd {
         return true;
     }
     /**
-     * Apply deCasteljou interpolations to isolate a smaller bezier polygon, representing interval 0..fraction of the original
-     * @param fraction "end" fraction for split.
-     * @returns false if fraction is 0 -- no changes applied.
+     * Apply deCasteljau interpolations to isolate a smaller bezier polygon, representing interval fraction..1 of the original
+     * @param fraction "start" fraction for split.
+     * @returns false if fraction is 0 or 1 -- no changes applied.
      */
     subdivideInPlaceKeepRight(fraction) {
         if (Geometry.isAlmostEqualNumber(fraction, 0.0))
@@ -315,7 +315,7 @@ export class Bezier1dNd {
     /**
      * Saturate a univariate bspline coefficient array in place
      * @param fraction0 fraction for first split.   This is the start of the output polygon
-     * @param fraction1 fraction for first split.   This is the start of the output polygon
+     * @param fraction1 fraction for second split.   This is the end of the output polygon
      * @return false if fractions are (almost) identical.
      */
     subdivideToIntervalInPlace(fraction0, fraction1) {
@@ -337,5 +337,4 @@ export class Bezier1dNd {
     /** map a fraction to the parent space. */
     fractionToParentFraction(fraction) { return this.interval ? this.interval.fractionToPoint(fraction) : fraction; }
 }
-Bezier1dNd._knotTolerance = 1.0e-8;
 //# sourceMappingURL=Bezier1dNd.js.map