@itwin/core-geometry 3.0.0-dev.76 → 3.0.0-dev.80

package/lib/esm/geometry4d/Point4d.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Point4d.js","sourceRoot":"","sources":["../../../src/geometry4d/Point4d.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAE/F;;GAEG;AACH,OAAO,EAAmB,QAAQ,EAAE,MAAM,aAAa,CAAC;AACxD,OAAO,EAAE,4BAA4B,EAAE,MAAM,4CAA4C,CAAC;AAC1F,OAAO,EAAE,yBAAyB,EAAE,MAAM,yCAAyC,CAAC;AACpF,OAAO,EAAE,OAAO,EAAE,MAAM,+BAA+B,CAAC;AACxD,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,+BAA+B,CAAC;AAClE,OAAO,EAAE,KAAK,EAAE,MAAM,qBAAqB,CAAC;AAQ5C;;;;;;;;;GASG;AACH,SAAS,mBAAmB,CAAC,GAAW,EAAE,EAAU,EAAE,GAAW,EAC/D,CAAS,EAAE,EAAU,EAAE,IAAY;IACnC,OAAO,IAAI,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;AAChD,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,OAAO,OAAO;IA+BlB,kCAAkC;IAClC,YAAsB,IAAY,CAAC,EAAE,IAAY,CAAC,EAAE,IAAY,CAAC,EAAE,IAAY,CAAC;QAC9E,IAAI,CAAC,IAAI,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAChC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjB,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjB,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjB,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACnB,CAAC;IAnCD,kCAAkC;IAC3B,GAAG,CAAC,IAAY,CAAC,EAAE,IAAY,CAAC,EAAE,IAAY,CAAC,EAAE,IAAY,CAAC;QACnE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjB,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjB,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjB,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACjB,OAAO,IAAI,CAAC;IACd,CAAC;IACD;;OAEG;IACI,YAAY,CAAC,KAAa,EAAE,KAAa;QAC9C,IAAI,KAAK,IAAI,CAAC,IAAI,KAAK,GAAG,CAAC,EAAE;YAC3B,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,KAAK,CAAC;SAC1B;IACH,CAAC;IACD,uBAAuB;IACvB,IAAW,CAAC,KAAK,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACvC,IAAW,CAAC,CAAC,GAAW,IAAI,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IACjD,uBAAuB;IACvB,IAAW,CAAC,KAAK,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACvC,IAAW,CAAC,CAAC,GAAW,IAAI,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IACjD,uBAAuB;IACvB,IAAW,CAAC,KAAK,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACvC,IAAW,CAAC,CAAC,GAAW,IAAI,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IACjD,qCAAqC;IACrC,IAAW,CAAC,KAAK,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACvC,IAAW,CAAC,CAAC,GAAW,IAAI,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IASjD,8CAA8C;IACvC,MAAM,CAAC,MAAM,CAAC,IAAY,CAAC,EAAE,IAAY,CAAC,EAAE,IAAY,CAAC,EAAE,IAAY,CAAC,EAAE,MAAgB;QAC/F,OAAO,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACnE,CAAC;IACD,qCAAqC;IAC9B,OAAO,CAAC,KAAc;QAC3B,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC7B,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC7B,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC7B,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC7B,OAAO,IAAI,CAAC;IACd,CAAC;IACD,uBAAuB;IAChB,KAAK,CAAC,MAAgB;QAC3B,OAAO,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7G,CAAC;IACD,0DAA0D;IACnD,WAAW,CAAC,IAAmB;QACpC,IAAI,QAAQ,CAAC,aAAa,CAAC,IAAI,EAAE,CAAC,CAAC;YACjC,IAAI,CAAC,GAAG,CAAC,IAAK,CAAC,CAAC,CAAC,EAAE,IAAK,CAAC,CAAC,CAAC,EAAE,IAAK,CAAC,CAAC,CAAC,EAAE,IAAK,CAAC,CAAC,CAAC,CAAC,CAAC;;YAEjD,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACzB,CAAC;IAED,wEAAwE;IACjE,MAAM,CAAC,QAAQ,CAAC,IAAmB;QACxC,MAAM,MAAM,GAAG,IAAI,OAAO,EAAE,CAAC;QAC7B,MAAM,CAAC,WAAW,CAAC,IAAI,CAAC,CAAC;QACzB,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,6EAA6E;IACtE,aAAa,CAAC,KAAc;QACjC,OAAO,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC;eAC5C,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC;eAC1C,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC;eAC1C,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IAClD,CAAC;IACD;;;;;;OAMG;IACI,iBAAiB,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS;QACjE,OAAO,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC;eACtC,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC;eACpC,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC;eACpC,QAAQ,CAAC,gBAAgB,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5C,CAAC;IAED;;;OAGG;IACI,MAAM;QACX,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IAClE,CAAC;IACD;;OAEG;IACI,YAAY,CAAC,KAAc;QAChC,OAAO,QAAQ,CAAC,cAAc,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACzJ,CAAC;IACD;;OAEG;IACI,mBAAmB,CAAC,KAAc;QACvC,OAAO,QAAQ,CAAC,qBAAqB,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IAChK,CAAC;IACD;OACG;IACI,cAAc,CAAC,KAAc;QAClC,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC;QAClB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC;QACnB,IAAI,QAAQ,CAAC,qBAAqB,CAAC,EAAE,CAAC,IAAI,QAAQ,CAAC,qBAAqB,CAAC,EAAE,CAAC;YAC1E,OAAO,SAAS,CAAC;QACnB,OAAO,QAAQ,CAAC,YAAY,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,EAAE,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IAC/G,CAAC;IACD;;OAEG;IACI,OAAO,CAAC,KAAc;QAC3B,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAClL,CAAC;IACD,oEAAoE;IAC7D,MAAM;QACX,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAClH,CAAC;IACD,+DAA+D;IACxD,aAAa;QAClB,OAAO,QAAQ,CAAC,cAAc,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACzF,CAAC;IACD,wHAAwH;IACjH,mBAAmB;QACxB,OAAO,QAAQ,CAAC,oBAAoB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACjF,CAAC;IAED,wEAAwE;IACjE,KAAK,CAAC,KAAc,EAAE,MAAgB;QAC3C,OAAO,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;IACxJ,CAAC;IACD,sDAAsD;IAC/C,kBAAkB,CAAC,KAAc,EAAE,MAAiB;QACzD,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACxB,MAAM,EAAE,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACzB,OAAO,QAAQ,CAAC,MAAM,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;IACzJ,CAAC;IACD,kFAAkF;IAC3E,yBAAyB,CAAC,KAAc,EAAE,MAAiB;QAChE,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACxB,OAAO,QAAQ,CAAC,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;IACxH,CAAC;IACD,uEAAuE;IAChE,IAAI,CAAC,KAAc,EAAE,MAAgB;QAC1C,OAAO,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;IACxJ,CAAC;IACD,8CAA8C;IAC9C,IAAW,YAAY;QACrB,OAAO,QAAQ,CAAC,qBAAqB,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC,CAAC;IACvD,CAAC;IACD,mDAAmD;IAC5C,MAAM,CAAC,UAAU,KAAc,OAAO,IAAI,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACvE;;;;OAIG;IACI,MAAM,CAAC,sBAAsB,CAAC,MAAe,EAAE,MAAe,EAAE,MAAgB;QACrF,OAAO,OAAO,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;IACtK,CAAC;IACD;;;;OAIG;IACI,MAAM,CAAC,oBAAoB,CAAC,IAAkB,EAAE,SAAiB,CAAC,EAAE,MAAgB;QACzF,OAAO,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,EAAE,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;IACpG,CAAC;IACD,0FAA0F;IACnF,MAAM,CAAC,wBAAwB,CAAC,GAAW,EAAE,CAAS;QAC3D,OAAO,IAAI,OAAO,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC7C,CAAC;IACH;;;;;;OAMG;IAEM,MAAM,CAAC,eAAe,CAAC,KAA0C;QACtE,IAAI,KAAK,YAAY,OAAO;YAC1B,OAAO,IAAI,OAAO,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC7C,IAAI,KAAK,YAAY,OAAO;YAC1B,OAAO,IAAI,OAAO,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACnD,IAAI,KAAK,YAAY,OAAO;YAC1B,OAAO,KAAK,CAAC,KAAK,EAAE,CAAC;QACvB,yCAAyC;QACzC,IAAI,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE;YACxB,MAAM,EAAE,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;YAC7C,MAAM,EAAE,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;YAC7C,MAAM,EAAE,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;YAC7C,MAAM,EAAE,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;YAC7C,OAAO,IAAI,OAAO,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;SACpC;QACD,MAAM,CAAC,GAAG,KAAK,CAAC,cAAc,CAAE,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QACrD,MAAM,CAAC,GAAG,KAAK,CAAC,cAAc,CAAE,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QACrD,MAAM,CAAC,GAAG,KAAK,CAAC,cAAc,CAAE,GAAG,CAAC,CAAC,CAAC,CAAE,KAAa,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QAC9D,MAAM,CAAC,GAAG,KAAK,CAAC,cAAc,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,KAAa,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QAC7D,OAAO,IAAI,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAEjC,CAAC;IACD,uCAAuC;IAChC,UAAU,CAAC,MAAe,EAAE,WAAmB,EAAE,MAAgB;QACtE,OAAO,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,WAAW,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,WAAW,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,WAAW,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,WAAW,EAAE,MAAM,CAAC,CAAC;IACpN,CAAC;IACD;OACG;IACI,WAAW,CAAC,QAAgB,EAAE,MAAe,EAAE,MAAgB;QACpE,MAAM,CAAC,GAAG,GAAG,GAAG,QAAQ,CAAC;QACzB,OAAO,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,EAAE,MAAM,CAAC,CAAC;IACxN,CAAC;IACD,6DAA6D;IACtD,WAAW,CAAC,OAAgB,EAAE,OAAe,EAAE,OAAgB,EAAE,OAAe,EAAE,MAAgB;QACvG,OAAO,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,MAAM,CAAC,CAAC;IACxT,CAAC;IACD,iFAAiF;IAC1E,WAAW,CAAC,OAAgB,EAAE,OAAe,EAAE,OAAgB,EAAE,OAAe,EAAE,OAAgB,EAAE,OAAe,EAAE,MAAgB;QAC1I,OAAO,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,MAAM,CAAC,CAAC;IACxa,CAAC;IACD,6DAA6D;IACtD,MAAM,CAAC,gBAAgB,CAAC,OAAgB,EAAE,OAAe,EAAE,OAAgB,EAAE,OAAe,EAAE,MAAgB;QACnH,OAAO,OAAO,CAAC,MAAM,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,MAAM,CAAC,CAAC;IAC5P,CAAC;IACD,iFAAiF;IAC1E,MAAM,CAAC,gBAAgB,CAAC,OAAgB,EAAE,OAAe,EAAE,OAAgB,EAAE,OAAe,EAAE,OAAgB,EAAE,OAAe,EAAE,MAAgB;QACtJ,OAAO,OAAO,CAAC,MAAM,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,EAAE,MAAM,CAAC,CAAC;IAC5W,CAAC;IACD,kFAAkF;IAC3E,mBAAmB,CAAC,OAAgB,EAAE,OAAgB;QAC3D,OAAO,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACxE,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACnE,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACnE,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACxE,CAAC;IACD,+DAA+D;IACxD,UAAU,CAAC,KAAc;QAC9B,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACnI,CAAC;IACD,wDAAwD;IACjD,cAAc,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS;QAC9D,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACnF,CAAC;IACD,8FAA8F;IACvF,QAAQ,CAAC,KAAc;QAC5B,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACjG,CAAC;IACD,4EAA4E;IACrE,WAAW,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;QAChD,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC/E,CAAC;IACD,oGAAoG;IAC7F,gBAAgB,CAAC,KAAc;QACpC,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IAC3G,CAAC;IACD,mGAAmG;IAC5F,QAAQ,CAAC,MAAgB;QAC9B,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IACrF,CAAC;IACD,4EAA4E;IACrE,WAAW,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;QAChD,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAChE,CAAC;IACD,oBAAoB;IACb,MAAM,CAAC,KAAK,KAAc,OAAO,IAAI,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAClE,oBAAoB;IACb,MAAM,CAAC,KAAK,KAAc,OAAO,IAAI,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAClE,oBAAoB;IACb,MAAM,CAAC,KAAK,KAAc,OAAO,IAAI,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAClE,oBAAoB;IACb,MAAM,CAAC,KAAK,KAAc,OAAO,IAAI,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAClE,0EAA0E;IACnE,gBAAgB,CAAC,WAAmB,EAAE,MAAgB;QAC3D,IAAI,WAAW,KAAK,GAAG,EAAE;YACvB,OAAO,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,WAAW,EAAE,MAAM,CAAC,CAAC;SAC9C;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,2CAA2C;IACpC,KAAK,CAAC,KAAa,EAAE,MAAgB;QAC1C,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,OAAO,EAAE,CAAC;QACzC,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QACtC,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QACtC,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QACtC,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QACtC,OAAO,MAAM,CAAC;IAChB,CAAC;IACD,wCAAwC;IACjC,MAAM,CAAC,MAAgB;QAC5B,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,OAAO,EAAE,CAAC;QACzC,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC/B,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC/B,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC/B,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC/B,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;;;OAIG;IACI,eAAe,CAAC,MAAgB;QACrC,MAAM,GAAG,GAAG,QAAQ,CAAC,0BAA0B,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9D,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,OAAO,EAAE,CAAC;QACzC,OAAO,IAAI,CAAC,gBAAgB,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;IAC5C,CAAC;IACD;;;;OAIG;IACI,SAAS,CAAC,MAAgB;QAC/B,MAAM,GAAG,GAAG,QAAQ,CAAC,0BAA0B,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9D,IAAI,GAAG,KAAK,GAAG;YACb,OAAO,SAAS,CAAC;QACnB,MAAM,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,8DAA8D;QACnF,OAAO,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC,CAAC;IACtF,CAAC;IAED;;;OAGG;IACI,iBAAiB;QACtB,MAAM,GAAG,GAAG,QAAQ,CAAC,0BAA0B,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9D,IAAI,GAAG,KAAK,GAAG;YACb,OAAO,QAAQ,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QACjD,MAAM,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,8DAA8D;QACnF,OAAO,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC5D,CAAC;IACD;;;;;;;;OAQG;IACI,MAAM,CAAC,2BAA2B,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,MAAgB;QACpG,MAAM,GAAG,GAAG,QAAQ,CAAC,0BAA0B,CAAC,CAAC,CAAC,CAAC;QACnD,MAAM,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,8DAA8D;QACvG,OAAO,OAAO,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC,CAAC;IACrD,CAAC;IACD;;;;;;;;;;;;OAYG;IACI,MAAM,CAAC,mCAAmC,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,MAAc;QAC1J,MAAM,GAAG,GAAG,QAAQ,CAAC,0BAA0B,CAAC,CAAC,CAAC,CAAC;QACnD,qBAAqB;QACrB,qFAAqF;QACrF,MAAM,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,8DAA8D;QACvG,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,OAAO,KAAK,CAAC,YAAY,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,MAAM,CAAC,CAAC;IACjI,CAAC;IACD;;;;;;;;;;;;OAYG;IACI,MAAM,CAAC,uDAAuD,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,GAAW,EAAE,GAAW,EAAE,GAAW,EAAE,GAAW,EAAE,MAAkC;QACtP,MAAM,GAAG,GAAG,QAAQ,CAAC,0BAA0B,CAAC,CAAC,CAAC,CAAC;QACnD,qBAAqB;QACrB,qFAAqF;QACrF,MAAM,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,8DAA8D;QACvG,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,GAAG,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;QACnC,MAAM,GAAG,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;QACnC,MAAM,GAAG,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;QACnC,OAAO,yBAAyB,CAAC,yBAAyB,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,mBAAmB,CAAC,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,mBAAmB,CAAC,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,mBAAmB,CAAC,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;IAC7O,CAAC;IACD;;;OAGG;IACI,mBAAmB,CAAC,MAAgB;QACzC,MAAM,GAAG,GAAG,QAAQ,CAAC,0BAA0B,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9D,IAAI,GAAG,KAAK,GAAG;YACb,OAAO,OAAO,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;QACzC,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,OAAO,EAAE,CAAC;QACzC,MAAM,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC;QACpB,OAAO,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC,CAAC;IACtF,CAAC;IACD;;;;OAIG;IACI,aAAa,CAAC,MAAgB;QACnC,MAAM,GAAG,GAAG,QAAQ,CAAC,0BAA0B,CAAC,IAAI,CAAC,aAAa,EAAE,CAAC,CAAC;QACtE,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,OAAO,EAAE,CAAC;QACzC,OAAO,IAAI,CAAC,gBAAgB,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;IAC5C,CAAC;IAED;;OAEG;IACI,MAAM,CAAC,qBAAqB,CAAC,MAAe,EAAE,MAAe,EAAE,MAAe,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS;QACpH,OAAO,QAAQ,CAAC,aAAa,CAC3B,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAC9C,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAC9C,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACpD,CAAC;IACD;;;;;OAKG;IACI,MAAM,CAAC,yBAAyB,CAAC,MAAe,EAAE,MAAe,EAAE,MAAe;QACvF,OAAO,OAAO,CAAC,MAAM,CACnB,OAAO,CAAC,qBAAqB,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAC9D,CAAC,OAAO,CAAC,qBAAqB,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAC/D,OAAO,CAAC,qBAAqB,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAC9D,CAAC,OAAO,CAAC,qBAAqB,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACrE,CAAC;IACD,sFAAsF;IAC/E,8BAA8B,CAAC,MAAqC;QACzE,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,mBAAmB,EAAE,CAAC,CAAC;QAChD,MAAM,SAAS,GAAG,QAAQ,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QAC1D,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;QACjB,MAAM,IAAI,GAAG,QAAQ,CAAC,yBAAyB,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACxD,IAAI,IAAI,KAAK,SAAS,EAAE;YACtB,MAAM,WAAW,GAAG,IAAI,GAAG,IAAI,CAAC;YAChC,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,WAAW,CAAC;YAC3B,SAAS,CAAC,YAAY,CAAC,WAAW,CAAC,CAAC;YACpC,OAAO,4BAA4B,CAAC,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,SAAS,EAAE,MAAM,CAAC,CAAC;SACnH;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,6DAA6D;IACtD,mBAAmB;QACxB,MAAM,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QAEnG,IAAI,SAAS,GAAG,GAAG,EAAE;YACnB,MAAM,CAAC,GAAG,GAAG,GAAG,SAAS,CAAC;YAC1B,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC;YACZ,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC;YACZ,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC;YACZ,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC;SACb;QACD,OAAO,SAAS,CAAC;IACnB,CAAC;IACD,6EAA6E;IACtE,MAAM,CAAC,sBAAsB,CAAC,WAAoB,EAAE,iBAAyB,EAAE,WAAoB,EAAE,MAAgB;QAC1H,IAAI,CAAC,MAAM;YACT,MAAM,GAAG,IAAI,OAAO,EAAE,CAAC;QACzB,MAAM,aAAa,GAAG,MAAM,CAAC;QAE7B,8CAA8C;QAC9C,IAAI,GAAG,KAAK,iBAAiB,EAAE;YAC7B,MAAM,GAAG,WAAW,CAAC;YACrB,OAAO,MAAM,CAAC;SACf;QACD,IAAI,GAAG,KAAK,iBAAiB,EAAE;YAC7B,MAAM,GAAG,WAAW,CAAC;YACrB,OAAO,MAAM,CAAC;SACf;QACD,IAAI,GAAG,KAAK,iBAAiB,EAAE;YAC7B,WAAW,CAAC,IAAI,CAAC,WAAW,EAAE,MAAM,CAAC,CAAC;YACtC,MAAM,CAAC,mBAAmB,EAAE,CAAC;YAC7B,OAAO,MAAM,CAAC;SACf;QAED,MAAM,EAAE,GAAG,WAAW,CAAC,KAAK,EAAE,CAAC;QAC/B,MAAM,EAAE,GAAG,WAAW,CAAC,KAAK,EAAE,CAAC;QAC/B,IAAI,GAAG,GAAG,WAAW,CAAC,UAAU,CAAC,WAAW,CAAC,CAAC;QAE9C,qDAAqD;QACrD,IAAI,GAAG,GAAG,GAAG,EAAE;YACb,EAAE,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC;YACd,GAAG,GAAG,CAAC,GAAG,CAAC;SACZ;QAED,wDAAwD;QACxD,IAAI,GAAG,GAAG,aAAa,EAAE;YACvB,EAAE,CAAC,WAAW,CAAC,iBAAiB,EAAE,EAAE,EAAE,MAAM,CAAC,CAAC;YAC9C,MAAM,CAAC,mBAAmB,EAAE,CAAC;YAC7B,OAAO,MAAM,CAAC;SACf;QAED,eAAe;QACf,IAAI,GAAG,GAAG,CAAC,GAAG;YACZ,GAAG,GAAG,CAAC,GAAG,CAAC;aACR,IAAI,GAAG,GAAG,GAAG;YAChB,GAAG,GAAG,GAAG,CAAC;QAEZ,oCAAoC;QACpC,MAAM,EAAE,GAAG,IAAI,OAAO,EAAE,CAAC;QACzB,EAAE,CAAC,UAAU,CAAC,EAAE,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,CAAE,0DAA0D;QACxF,EAAE,CAAC,mBAAmB,EAAE,CAAC;QAEzB,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAC7B,MAAM,oBAAoB,GAAG,KAAK,GAAG,iBAAiB,CAAC;QACvD,MAAM,GAAG,OAAO,CAAC,gBAAgB,CAAC,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,oBAAoB,CAAC,EAAE,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,oBAAoB,CAAC,CAAC,CAAC;QAC1G,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;OAEG;IACI,oBAAoB,CAAC,KAAc;QACxC,MAAM,IAAI,GAAG,IAAI,CAAC,aAAa,EAAE,CAAC;QAClC,MAAM,IAAI,GAAG,KAAK,CAAC,aAAa,EAAE,CAAC;QACnC,MAAM,GAAG,GAAG,IAAI,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC,CAAE,+BAA+B;QACpE,MAAM,GAAG,GAAG,QAAQ,CAAC,yBAAyB,CAAC,GAAG,EAAE,IAAI,GAAG,IAAI,CAAC,CAAC;QACjE,IAAI,GAAG,KAAK,SAAS;YACnB,OAAO,SAAS,CAAC;QACnB,OAAO,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IACxB,CAAC;CACF","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n\r\n/** @packageDocumentation\r\n * @module Numerics\r\n */\r\nimport { BeJSONFunctions, Geometry } from \"../Geometry\";\r\nimport { Plane3dByOriginAndUnitNormal } from \"../geometry3d/Plane3dByOriginAndUnitNormal\";\r\nimport { Plane3dByOriginAndVectors } from \"../geometry3d/Plane3dByOriginAndVectors\";\r\nimport { Point2d } from \"../geometry3d/Point2dVector2d\";\r\nimport { Point3d, Vector3d } from \"../geometry3d/Point3dVector3d\";\r\nimport { Ray3d } from \"../geometry3d/Ray3d\";\r\nimport { XAndY, XYAndZ } from \"../geometry3d/XYZProps\";\r\n\r\n/**\r\n * 4d point packed in an array of 4 numbers.\r\n * @public\r\n */\r\nexport type Point4dProps = number[];\r\n/**\r\n *\r\n * @param ddg numerator second derivative\r\n * @param dh denominator derivative\r\n * @param ddh denominator second derivative\r\n * @param f primary function (g/h)\r\n * @param df derivative of (g/h)\r\n * @param divH = (1/h)\r\n * @internal\r\n */\r\nfunction quotientDerivative2(ddg: number, dh: number, ddh: number,\r\n  f: number, df: number, divH: number): number {\r\n  return divH * (ddg - 2.0 * df * dh - f * ddh);\r\n}\r\n\r\n/** 4 Dimensional point (x,y,z,w) used in perspective calculations.\r\n * * the coordinates are stored in a Float64Array of length 4.\r\n * * properties `x`, `y`, `z`, `w` access array members.\r\n * *\r\n * * The coordinates are physically stored as a single Float64Array with 4 entries. (w last)\r\n * *\r\n * @public\r\n */\r\nexport class Point4d implements BeJSONFunctions {\r\n  /** x,y,z,w are packed into a Float64Array */\r\n  public xyzw: Float64Array;\r\n  /** Set x,y,z,w of this point.  */\r\n  public set(x: number = 0, y: number = 0, z: number = 0, w: number = 0): Point4d {\r\n    this.xyzw[0] = x;\r\n    this.xyzw[1] = y;\r\n    this.xyzw[2] = z;\r\n    this.xyzw[3] = w;\r\n    return this;\r\n  }\r\n  /** Set a component by index.\r\n   * * No change if index is out of range.\r\n   */\r\n  public setComponent(index: number, value: number) {\r\n    if (index >= 0 && index < 4) {\r\n      this.xyzw[index] = value;\r\n    }\r\n  }\r\n  /** The x component. */\r\n  public get x() { return this.xyzw[0]; }\r\n  public set x(val: number) { this.xyzw[0] = val; }\r\n  /** The y component. */\r\n  public get y() { return this.xyzw[1]; }\r\n  public set y(val: number) { this.xyzw[1] = val; }\r\n  /** The z component. */\r\n  public get z() { return this.xyzw[2]; }\r\n  public set z(val: number) { this.xyzw[2] = val; }\r\n  /** The w component of this point. */\r\n  public get w() { return this.xyzw[3]; }\r\n  public set w(val: number) { this.xyzw[3] = val; }\r\n  /** Construct from coordinates. */\r\n  protected constructor(x: number = 0, y: number = 0, z: number = 0, w: number = 0) {\r\n    this.xyzw = new Float64Array(4);\r\n    this.xyzw[0] = x;\r\n    this.xyzw[1] = y;\r\n    this.xyzw[2] = z;\r\n    this.xyzw[3] = w;\r\n  }\r\n  /** Return a Point4d with specified x,y,z,w */\r\n  public static create(x: number = 0, y: number = 0, z: number = 0, w: number = 0, result?: Point4d): Point4d {\r\n    return result ? result.set(x, y, z, w) : new Point4d(x, y, z, w);\r\n  }\r\n  /** Copy coordinates from `other`. */\r\n  public setFrom(other: Point4d): Point4d {\r\n    this.xyzw[0] = other.xyzw[0];\r\n    this.xyzw[1] = other.xyzw[1];\r\n    this.xyzw[2] = other.xyzw[2];\r\n    this.xyzw[3] = other.xyzw[3];\r\n    return this;\r\n  }\r\n  /** Clone this point */\r\n  public clone(result?: Point4d): Point4d {\r\n    return result ? result.setFrom(this) : new Point4d(this.xyzw[0], this.xyzw[1], this.xyzw[2], this.xyzw[3]);\r\n  }\r\n  /** Set this point's xyzw from a json array `[x,y,z,w]` */\r\n  public setFromJSON(json?: Point4dProps) {\r\n    if (Geometry.isNumberArray(json, 4))\r\n      this.set(json![0], json![1], json![2], json![3]);\r\n    else\r\n      this.set(0, 0, 0, 0);\r\n  }\r\n\r\n  /** Create a new point with coordinates from a json array `[x,y,z,w]` */\r\n  public static fromJSON(json?: Point4dProps): Point4d {\r\n    const result = new Point4d();\r\n    result.setFromJSON(json);\r\n    return result;\r\n  }\r\n  /** Near-equality test, using `Geometry.isSameCoordinate` on all 4 x,y,z,w */\r\n  public isAlmostEqual(other: Point4d): boolean {\r\n    return Geometry.isSameCoordinate(this.x, other.x)\r\n      && Geometry.isSameCoordinate(this.y, other.y)\r\n      && Geometry.isSameCoordinate(this.z, other.z)\r\n      && Geometry.isSameCoordinate(this.w, other.w);\r\n  }\r\n  /**\r\n   * Test for same coordinate by direct x,y,z,w args\r\n   * @param x x to test\r\n   * @param y y to test\r\n   * @param z z to test\r\n   * @param w w to test\r\n   */\r\n  public isAlmostEqualXYZW(x: number, y: number, z: number, w: number): boolean {\r\n    return Geometry.isSameCoordinate(this.x, x)\r\n      && Geometry.isSameCoordinate(this.y, y)\r\n      && Geometry.isSameCoordinate(this.z, z)\r\n      && Geometry.isSameCoordinate(this.w, w);\r\n  }\r\n\r\n  /**\r\n   * Convert an Angle to a JSON object.\r\n   * @return {*} [x,y,z,w]\r\n   */\r\n  public toJSON(): Point4dProps {\r\n    return [this.xyzw[0], this.xyzw[1], this.xyzw[2], this.xyzw[3]];\r\n  }\r\n  /** Return the 4d distance from this point to other, with all 4 components squared into the hypotenuse.\r\n   * * x,y,z,w all participate without normalization.\r\n   */\r\n  public distanceXYZW(other: Point4d): number {\r\n    return Geometry.hypotenuseXYZW(other.xyzw[0] - this.xyzw[0], other.xyzw[1] - this.xyzw[1], other.xyzw[2] - this.xyzw[2], other.xyzw[3] - this.xyzw[3]);\r\n  }\r\n  /** Return the squared 4d distance from this point to other, with all 4 components squared into the hypotenuse.\r\n   * * x,y,z,w all participate without normalization.\r\n   */\r\n  public distanceSquaredXYZW(other: Point4d): number {\r\n    return Geometry.hypotenuseSquaredXYZW(other.xyzw[0] - this.xyzw[0], other.xyzw[1] - this.xyzw[1], other.xyzw[2] - this.xyzw[2], other.xyzw[3] - this.xyzw[3]);\r\n  }\r\n  /** Return the distance between the instance and other after normalizing by weights\r\n   */\r\n  public realDistanceXY(other: Point4d): number | undefined {\r\n    const wA = this.w;\r\n    const wB = other.w;\r\n    if (Geometry.isSmallMetricDistance(wA) || Geometry.isSmallMetricDistance(wB))\r\n      return undefined;\r\n    return Geometry.hypotenuseXY(other.xyzw[0] / wB - this.xyzw[0] / wA, other.xyzw[1] / wB - this.xyzw[1] / wA);\r\n  }\r\n  /** Return the largest absolute distance between corresponding components\r\n   * * x,y,z,w all participate without normalization.\r\n   */\r\n  public maxDiff(other: Point4d): number {\r\n    return Math.max(Math.abs(other.xyzw[0] - this.xyzw[0]), Math.abs(other.xyzw[1] - this.xyzw[1]), Math.abs(other.xyzw[2] - this.xyzw[2]), Math.abs(other.xyzw[3] - this.xyzw[3]));\r\n  }\r\n  /** Return the largest absolute entry of all 4 components x,y,z,w */\r\n  public maxAbs(): number {\r\n    return Math.max(Math.abs(this.xyzw[0]), Math.abs(this.xyzw[1]), Math.abs(this.xyzw[2]), Math.abs(this.xyzw[3]));\r\n  }\r\n  /** Returns the magnitude including all 4 components x,y,z,w */\r\n  public magnitudeXYZW(): number {\r\n    return Geometry.hypotenuseXYZW(this.xyzw[0], this.xyzw[1], this.xyzw[2], this.xyzw[3]);\r\n  }\r\n  /** Returns the magnitude of the leading xyz components.  w is ignored.  (i.e. the leading xyz are NOT divided by w.) */\r\n  public magnitudeSquaredXYZ(): number {\r\n    return Geometry.hypotenuseSquaredXYZ(this.xyzw[0], this.xyzw[1], this.xyzw[2]);\r\n  }\r\n\r\n  /** Return the difference (this-other) using all 4 components x,y,z,w */\r\n  public minus(other: Point4d, result?: Point4d): Point4d {\r\n    return Point4d.create(this.xyzw[0] - other.xyzw[0], this.xyzw[1] - other.xyzw[1], this.xyzw[2] - other.xyzw[2], this.xyzw[3] - other.xyzw[3], result);\r\n  }\r\n  /** Return `((other.w * this) -  (this.w * other))` */\r\n  public crossWeightedMinus(other: Point4d, result?: Vector3d): Vector3d {\r\n    const wa = this.xyzw[3];\r\n    const wb = other.xyzw[3];\r\n    return Vector3d.create(wb * this.xyzw[0] - wa * other.xyzw[0], wb * this.xyzw[1] - wa * other.xyzw[1], wb * this.xyzw[2] - wa * other.xyzw[2], result);\r\n  }\r\n  /** Return `((other.w * this) -  (this.w * other))`, with other.w known to be 1 */\r\n  public crossWeightedMinusPoint3d(other: Point3d, result?: Vector3d): Vector3d {\r\n    const wa = this.xyzw[3];\r\n    return Vector3d.create(this.xyzw[0] - wa * other.x, this.xyzw[1] - wa * other.y, this.xyzw[2] - wa * other.z, result);\r\n  }\r\n  /** Return the sum of this and other, using all 4 components x,y,z,w */\r\n  public plus(other: Point4d, result?: Point4d): Point4d {\r\n    return Point4d.create(this.xyzw[0] + other.xyzw[0], this.xyzw[1] + other.xyzw[1], this.xyzw[2] + other.xyzw[2], this.xyzw[3] + other.xyzw[3], result);\r\n  }\r\n  /** Test if all components are nearly zero. */\r\n  public get isAlmostZero(): boolean {\r\n    return Geometry.isSmallMetricDistance(this.maxAbs());\r\n  }\r\n  /** Create a point with zero in all coordinates. */\r\n  public static createZero(): Point4d { return new Point4d(0, 0, 0, 0); }\r\n  /**\r\n   * Create plane coefficients for the plane containing pointA, pointB, and 0010.\r\n   * @param pointA first point\r\n   * @param pointB second point\r\n   */\r\n  public static createPlanePointPointZ(pointA: Point4d, pointB: Point4d, result?: Point4d) {\r\n    return Point4d.create(pointA.y * pointB.w - pointA.w * pointB.y, pointA.w * pointB.x - pointA.x * pointB.w, 0.0, pointA.x * pointB.y - pointA.y * pointB.x, result);\r\n  }\r\n  /**\r\n   * extract 4 consecutive numbers from a Float64Array into a Point4d.\r\n   * @param data buffer of numbers\r\n   * @param xIndex first index for x,y,z,w sequence\r\n   */\r\n  public static createFromPackedXYZW(data: Float64Array, xIndex: number = 0, result?: Point4d): Point4d {\r\n    return Point4d.create(data[xIndex], data[xIndex + 1], data[xIndex + 2], data[xIndex + 3], result);\r\n  }\r\n  /** Create a `Point4d` with x,y,z from an `XYAndZ` input, and w from a separate number. */\r\n  public static createFromPointAndWeight(xyz: XYAndZ, w: number): Point4d {\r\n    return new Point4d(xyz.x, xyz.y, xyz.z, w);\r\n  }\r\n/** Create a `Point4d` from\r\n * * Point2d, Point3d, or Point4d\r\n * * other structure with members x,y and optional z,w\r\n * * array of numbers\r\n * * default z is 0.0\r\n * * default 2 is 1.0  (array[3] can replace)\r\n */\r\n\r\n  public static createFromPoint(point: XAndY | XYAndZ | Point4d | number[]): Point4d {\r\n    if (point instanceof Point2d)\r\n      return new Point4d(point.x, point.y, 0, 1);\r\n    if (point instanceof Point3d)\r\n      return new Point4d(point.x, point.y, point.z, 1);\r\n    if (point instanceof Point4d)\r\n      return point.clone();\r\n    // hm ... some flavor of x,y,z subset ...\r\n    if (Array.isArray(point)) {\r\n      const x1 = point.length > 0 ? point[0] : 0.0;\r\n      const y1 = point.length > 1 ? point[1] : 0.0;\r\n      const z1 = point.length > 2 ? point[2] : 0.0;\r\n      const w1 = point.length > 3 ? point[3] : 1.0;\r\n      return new Point4d(x1, y1, z1, w1);\r\n    }\r\n    const x = point.hasOwnProperty (\"x\") ? point.x : 0.0;\r\n    const y = point.hasOwnProperty (\"y\") ? point.y : 0.0;\r\n    const z = point.hasOwnProperty (\"z\") ? (point as any).z : 0.0;\r\n    const w = point.hasOwnProperty(\"w\") ? (point as any).w : 0.0;\r\n    return new Point4d(x, y, z, w);\r\n\r\n  }\r\n  /** Return `point + vector * scalar` */\r\n  public plusScaled(vector: Point4d, scaleFactor: number, result?: Point4d): Point4d {\r\n    return Point4d.create(this.xyzw[0] + vector.xyzw[0] * scaleFactor, this.xyzw[1] + vector.xyzw[1] * scaleFactor, this.xyzw[2] + vector.xyzw[2] * scaleFactor, this.xyzw[3] + vector.xyzw[3] * scaleFactor, result);\r\n  }\r\n  /** Return interpolation between instance and pointB at fraction\r\n   */\r\n  public interpolate(fraction: number, pointB: Point4d, result?: Point4d): Point4d {\r\n    const v = 1.0 - fraction;\r\n    return Point4d.create(this.xyzw[0] * v + pointB.xyzw[0] * fraction, this.xyzw[1] * v + pointB.xyzw[1] * fraction, this.xyzw[2] * v + pointB.xyzw[2] * fraction, this.xyzw[3] * v + pointB.xyzw[3] * fraction, result);\r\n  }\r\n  /** Return `point + vectorA * scalarA + vectorB * scalarB` */\r\n  public plus2Scaled(vectorA: Point4d, scalarA: number, vectorB: Point4d, scalarB: number, result?: Point4d): Point4d {\r\n    return Point4d.create(this.xyzw[0] + vectorA.xyzw[0] * scalarA + vectorB.xyzw[0] * scalarB, this.xyzw[1] + vectorA.xyzw[1] * scalarA + vectorB.xyzw[1] * scalarB, this.xyzw[2] + vectorA.xyzw[2] * scalarA + vectorB.xyzw[2] * scalarB, this.xyzw[3] + vectorA.xyzw[3] * scalarA + vectorB.xyzw[3] * scalarB, result);\r\n  }\r\n  /** Return `point + vectorA * scalarA + vectorB * scalarB + vectorC * scalarC` */\r\n  public plus3Scaled(vectorA: Point4d, scalarA: number, vectorB: Point4d, scalarB: number, vectorC: Point4d, scalarC: number, result?: Point4d): Point4d {\r\n    return Point4d.create(this.xyzw[0] + vectorA.xyzw[0] * scalarA + vectorB.xyzw[0] * scalarB + vectorC.xyzw[0] * scalarC, this.xyzw[1] + vectorA.xyzw[1] * scalarA + vectorB.xyzw[1] * scalarB + vectorC.xyzw[1] * scalarC, this.xyzw[2] + vectorA.xyzw[2] * scalarA + vectorB.xyzw[2] * scalarB + vectorC.xyzw[2] * scalarC, this.xyzw[3] + vectorA.xyzw[3] * scalarA + vectorB.xyzw[3] * scalarB + vectorC.xyzw[3] * scalarC, result);\r\n  }\r\n  /** Return `point + vectorA * scalarA + vectorB * scalarB` */\r\n  public static createAdd2Scaled(vectorA: Point4d, scalarA: number, vectorB: Point4d, scalarB: number, result?: Point4d): Point4d {\r\n    return Point4d.create(vectorA.xyzw[0] * scalarA + vectorB.xyzw[0] * scalarB, vectorA.xyzw[1] * scalarA + vectorB.xyzw[1] * scalarB, vectorA.xyzw[2] * scalarA + vectorB.xyzw[2] * scalarB, vectorA.xyzw[3] * scalarA + vectorB.xyzw[3] * scalarB, result);\r\n  }\r\n  /** Return `point + vectorA \\ scalarA + vectorB * scalarB + vectorC * scalarC` */\r\n  public static createAdd3Scaled(vectorA: Point4d, scalarA: number, vectorB: Point4d, scalarB: number, vectorC: Point4d, scalarC: number, result?: Point4d): Point4d {\r\n    return Point4d.create(vectorA.xyzw[0] * scalarA + vectorB.xyzw[0] * scalarB + vectorC.xyzw[0] * scalarC, vectorA.xyzw[1] * scalarA + vectorB.xyzw[1] * scalarB + vectorC.xyzw[1] * scalarC, vectorA.xyzw[2] * scalarA + vectorB.xyzw[2] * scalarB + vectorC.xyzw[2] * scalarC, vectorA.xyzw[3] * scalarA + vectorB.xyzw[3] * scalarB + vectorC.xyzw[3] * scalarC, result);\r\n  }\r\n  /** Return dot product of (4d) vectors from the instance to targetA and targetB */\r\n  public dotVectorsToTargets(targetA: Point4d, targetB: Point4d): number {\r\n    return (targetA.xyzw[0] - this.xyzw[0]) * (targetB.xyzw[0] - this.xyzw[0]) +\r\n      (targetA.xyzw[1] - this.xyzw[1]) * (targetB.xyzw[1] - this.xyzw[1]) +\r\n      (targetA.xyzw[2] - this.xyzw[2]) * (targetB.xyzw[2] - this.xyzw[2]) +\r\n      (targetA.xyzw[3] - this.xyzw[3]) * (targetB.xyzw[3] - this.xyzw[3]);\r\n  }\r\n  /** return (4d) dot product of the instance and other point. */\r\n  public dotProduct(other: Point4d): number {\r\n    return this.xyzw[0] * other.xyzw[0] + this.xyzw[1] * other.xyzw[1] + this.xyzw[2] * other.xyzw[2] + this.xyzw[3] * other.xyzw[3];\r\n  }\r\n  /** return (4d) dot product of the instance with xyzw */\r\n  public dotProductXYZW(x: number, y: number, z: number, w: number): number {\r\n    return this.xyzw[0] * x + this.xyzw[1] * y + this.xyzw[2] * z + this.xyzw[3] * w;\r\n  }\r\n  /** dotProduct with (point.x, point.y, point.z, 1) Used in PlaneAltitudeEvaluator interface */\r\n  public altitude(point: Point3d): number {\r\n    return this.xyzw[0] * point.x + this.xyzw[1] * point.y + this.xyzw[2] * point.z + this.xyzw[3];\r\n  }\r\n  /** dotProduct with (x, y, z, 1) Used in PlaneAltitudeEvaluator interface */\r\n  public altitudeXYZ(x: number, y: number, z: number): number {\r\n    return this.xyzw[0] * x + this.xyzw[1] * y + this.xyzw[2] * z + this.xyzw[3];\r\n  }\r\n  /** dotProduct with (point.x, point.y, point.z, point.w) Used in PlaneAltitudeEvaluator interface */\r\n  public weightedAltitude(point: Point4d): number {\r\n    return this.xyzw[0] * point.x + this.xyzw[1] * point.y + this.xyzw[2] * point.z + this.xyzw[3] * point.w;\r\n  }\r\n  /** dotProduct with (vector.x, vector.y, vector.z, 0).  Used in PlaneAltitudeEvaluator interface */\r\n  public velocity(vector: Vector3d): number {\r\n    return this.xyzw[0] * vector.x + this.xyzw[1] * vector.y + this.xyzw[2] * vector.z;\r\n  }\r\n  /** dotProduct with (x,y,z, 0).  Used in PlaneAltitudeEvaluator interface */\r\n  public velocityXYZ(x: number, y: number, z: number): number {\r\n    return this.xyzw[0] * x + this.xyzw[1] * y + this.xyzw[2] * z;\r\n  }\r\n  /** unit X vector */\r\n  public static unitX(): Point4d { return new Point4d(1, 0, 0, 0); }\r\n  /** unit Y vector */\r\n  public static unitY(): Point4d { return new Point4d(0, 1, 0, 0); }\r\n  /** unit Z vector */\r\n  public static unitZ(): Point4d { return new Point4d(0, 0, 1, 0); }\r\n  /** unit W vector */\r\n  public static unitW(): Point4d { return new Point4d(0, 0, 0, 1); }\r\n  /** Divide by denominator, but return undefined if denominator is zero. */\r\n  public safeDivideOrNull(denominator: number, result?: Point4d): Point4d | undefined {\r\n    if (denominator !== 0.0) {\r\n      return this.scale(1.0 / denominator, result);\r\n    }\r\n    return undefined;\r\n  }\r\n  /** scale all components (including w!!) */\r\n  public scale(scale: number, result?: Point4d): Point4d {\r\n    result = result ? result : new Point4d();\r\n    result.xyzw[0] = this.xyzw[0] * scale;\r\n    result.xyzw[1] = this.xyzw[1] * scale;\r\n    result.xyzw[2] = this.xyzw[2] * scale;\r\n    result.xyzw[3] = this.xyzw[3] * scale;\r\n    return result;\r\n  }\r\n  /** Negate components (including w!!) */\r\n  public negate(result?: Point4d): Point4d {\r\n    result = result ? result : new Point4d();\r\n    result.xyzw[0] = -this.xyzw[0];\r\n    result.xyzw[1] = -this.xyzw[1];\r\n    result.xyzw[2] = -this.xyzw[2];\r\n    result.xyzw[3] = -this.xyzw[3];\r\n    return result;\r\n  }\r\n  /**\r\n   * If `this.w` is nonzero, return a 4d point `(x/w,y/w,z/w, 1)`\r\n   * If `this.w` is zero, return undefined.\r\n   * @param result optional result\r\n   */\r\n  public normalizeWeight(result?: Point4d): Point4d | undefined {\r\n    const mag = Geometry.correctSmallMetricDistance(this.xyzw[3]);\r\n    result = result ? result : new Point4d();\r\n    return this.safeDivideOrNull(mag, result);\r\n  }\r\n  /**\r\n   * If `this.w` is nonzero, return a 3d point `(x/w,y/w,z/w)`\r\n   * If `this.w` is zero, return undefined.\r\n   * @param result optional result\r\n   */\r\n  public realPoint(result?: Point3d): Point3d | undefined {\r\n    const mag = Geometry.correctSmallMetricDistance(this.xyzw[3]);\r\n    if (mag === 0.0)\r\n      return undefined;\r\n    const a = 1.0 / mag; // in zero case everything multiplies right back to true zero.\r\n    return Point3d.create(this.xyzw[0] * a, this.xyzw[1] * a, this.xyzw[2] * a, result);\r\n  }\r\n\r\n  /** Convert the homogeneous point to a (strongly typed) point or vector.\r\n   * * If `this.w` is nonzero, return a Point3d `(x/w,y/w,z/w)`\r\n   * * If `this.w` is zero, return a Vector3d `(x,y,z)`\r\n   */\r\n  public realPointOrVector(): Point3d | Vector3d {\r\n    const mag = Geometry.correctSmallMetricDistance(this.xyzw[3]);\r\n    if (mag === 0.0)\r\n      return Vector3d.create(this.x, this.y, this.z);\r\n    const a = 1.0 / mag; // in zero case everything multiplies right back to true zero.\r\n    return Point3d.create(this.x * a, this.y * a, this.z * a);\r\n  }\r\n  /**\r\n   * * If w is nonzero, return Point3d with x/w,y/w,z/w.\r\n   * * If w is zero, return 000\r\n   * @param x x coordinate\r\n   * @param y y coordinate\r\n   * @param z z coordinate\r\n   * @param w w coordinate\r\n   * @param result optional result\r\n   */\r\n  public static createRealPoint3dDefault000(x: number, y: number, z: number, w: number, result?: Point3d): Point3d {\r\n    const mag = Geometry.correctSmallMetricDistance(w);\r\n    const a = mag === 0 ? 0.0 : (1.0 / mag); // in zero case everything multiplies right back to true zero.\r\n    return Point3d.create(x * a, y * a, z * a, result);\r\n  }\r\n  /**\r\n   * * If w is nonzero, return Vector3d which is the derivative of the projected xyz with given w and 4d derivatives.\r\n   * * If w is zero, return 000\r\n   * @param x x coordinate\r\n   * @param y y coordinate\r\n   * @param z z coordinate\r\n   * @param w w coordinate\r\n   * @param dx x coordinate of derivative\r\n   * @param dy y coordinate of derivative\r\n   * @param dz z coordinate of derivative\r\n   * @param dw w coordinate of derivative\r\n   * @param result optional result\r\n   */\r\n  public static createRealDerivativeRay3dDefault000(x: number, y: number, z: number, w: number, dx: number, dy: number, dz: number, dw: number, result?: Ray3d): Ray3d {\r\n    const mag = Geometry.correctSmallMetricDistance(w);\r\n    // real point is X/w.\r\n    // real derivative is (X' * w - X *w) / ww, and weight is always 0 by cross products.\r\n    const a = mag === 0 ? 0.0 : (1.0 / mag); // in zero case everything multiplies right back to true zero.\r\n    const aa = a * a;\r\n    return Ray3d.createXYZUVW(x * a, y * a, z * a, (dx * w - dw * x) * aa, (dy * w - dw * y) * aa, (dz * w - dw * z) * aa, result);\r\n  }\r\n  /**\r\n   * * If w is nonzero, return Vector3d which is the derivative of the projected xyz with given w and 4d derivatives.\r\n   * * If w is zero, return 000\r\n   * @param x x coordinate\r\n   * @param y y coordinate\r\n   * @param z z coordinate\r\n   * @param w w coordinate\r\n   * @param dx x coordinate of derivative\r\n   * @param dy y coordinate of derivative\r\n   * @param dz z coordinate of derivative\r\n   * @param dw w coordinate of derivative\r\n   * @param result optional result\r\n   */\r\n  public static createRealDerivativePlane3dByOriginAndVectorsDefault000(x: number, y: number, z: number, w: number, dx: number, dy: number, dz: number, dw: number, ddx: number, ddy: number, ddz: number, ddw: number, result?: Plane3dByOriginAndVectors): Plane3dByOriginAndVectors {\r\n    const mag = Geometry.correctSmallMetricDistance(w);\r\n    // real point is X/w.\r\n    // real derivative is (X' * w - X *w) / ww, and weight is always 0 by cross products.\r\n    const a = mag === 0 ? 0.0 : (1.0 / mag); // in zero case everything multiplies right back to true zero.\r\n    const aa = a * a;\r\n    const fx = x * a;\r\n    const fy = y * a;\r\n    const fz = z * a;\r\n    const dfx = (dx * w - dw * x) * aa;\r\n    const dfy = (dy * w - dw * y) * aa;\r\n    const dfz = (dz * w - dw * z) * aa;\r\n    return Plane3dByOriginAndVectors.createOriginAndVectorsXYZ(fx, fy, fz, dfx, dfy, dfz, quotientDerivative2(ddx, dw, ddw, fx, dfx, a), quotientDerivative2(ddy, dw, ddw, fy, dfy, a), quotientDerivative2(ddz, dw, ddw, fz, dfz, a), result);\r\n  }\r\n  /**\r\n   * * If this.w is nonzero, return Point3d with x/w,y/w,z/w.\r\n   * * If this.w is zero, return 000\r\n   */\r\n  public realPointDefault000(result?: Point3d): Point3d {\r\n    const mag = Geometry.correctSmallMetricDistance(this.xyzw[3]);\r\n    if (mag === 0.0)\r\n      return Point3d.create(0, 0, 0, result);\r\n    result = result ? result : new Point3d();\r\n    const a = 1.0 / mag;\r\n    return Point3d.create(this.xyzw[0] * a, this.xyzw[1] * a, this.xyzw[2] * a, result);\r\n  }\r\n  /** divide all components (x,y,z,w) by the 4d magnitude.\r\n   *\r\n   * * This is appropriate for normalizing a quaternion\r\n   * * Use normalizeWeight to divide by the w component.\r\n   */\r\n  public normalizeXYZW(result?: Point4d): Point4d | undefined {\r\n    const mag = Geometry.correctSmallMetricDistance(this.magnitudeXYZW());\r\n    result = result ? result : new Point4d();\r\n    return this.safeDivideOrNull(mag, result);\r\n  }\r\n\r\n  /**\r\n   * Return the determinant of the 3x3 matrix using components i,j,k of the 3 inputs.\r\n   */\r\n  public static determinantIndexed3X3(pointA: Point4d, pointB: Point4d, pointC: Point4d, i: number, j: number, k: number) {\r\n    return Geometry.tripleProduct(\r\n      pointA.xyzw[i], pointA.xyzw[j], pointA.xyzw[k],\r\n      pointB.xyzw[i], pointB.xyzw[j], pointB.xyzw[k],\r\n      pointC.xyzw[i], pointC.xyzw[j], pointC.xyzw[k]);\r\n  }\r\n  /**\r\n   * Return a Point4d perpendicular to all 3 inputs. (A higher level cross product concept)\r\n   * @param pointA first point\r\n   * @param pointB second point\r\n   * @param pointC third point\r\n   */\r\n  public static perpendicularPoint4dPlane(pointA: Point4d, pointB: Point4d, pointC: Point4d): Point4d {\r\n    return Point4d.create(\r\n      Point4d.determinantIndexed3X3(pointA, pointB, pointC, 1, 2, 3),\r\n      -Point4d.determinantIndexed3X3(pointA, pointB, pointC, 2, 3, 0),\r\n      Point4d.determinantIndexed3X3(pointA, pointB, pointC, 3, 0, 1),\r\n      -Point4d.determinantIndexed3X3(pointA, pointB, pointC, 0, 1, 2));\r\n  }\r\n  /** Treating this Point4d as plane coefficients, convert to origin and normal form. */\r\n  public toPlane3dByOriginAndUnitNormal(result?: Plane3dByOriginAndUnitNormal): Plane3dByOriginAndUnitNormal | undefined {\r\n    const a = Math.sqrt(this.magnitudeSquaredXYZ());\r\n    const direction = Vector3d.create(this.x, this.y, this.z);\r\n    const w = this.w;\r\n    const divA = Geometry.conditionalDivideFraction(1.0, a);\r\n    if (divA !== undefined) {\r\n      const divASquared = divA * divA;\r\n      const b = -w * divASquared;\r\n      direction.scaleInPlace(divASquared);\r\n      return Plane3dByOriginAndUnitNormal.create(Point3d.create(this.x * b, this.y * b, this.z * b), direction, result);\r\n    }\r\n    return undefined;\r\n  }\r\n  /** Normalize so sum of squares of all 4 coordinates is 1. */\r\n  public normalizeQuaternion() {\r\n    const magnitude = Math.sqrt(this.x * this.x + this.y * this.y + this.z * this.z + this.w * this.w);\r\n\r\n    if (magnitude > 0.0) {\r\n      const f = 1.0 / magnitude;\r\n      this.x *= f;\r\n      this.y *= f;\r\n      this.z *= f;\r\n      this.w *= f;\r\n    }\r\n    return magnitude;\r\n  }\r\n  /** Return a (normalized) quaternion interpolated between two quaternions. */\r\n  public static interpolateQuaternions(quaternion0: Point4d, fractionParameter: number, quaternion1: Point4d, result?: Point4d): Point4d {\r\n    if (!result)\r\n      result = new Point4d();\r\n    const maxSafeCosine = 0.9995;\r\n\r\n    // return exact quaternions for special values\r\n    if (0.0 === fractionParameter) {\r\n      result = quaternion0;\r\n      return result;\r\n    }\r\n    if (1.0 === fractionParameter) {\r\n      result = quaternion1;\r\n      return result;\r\n    }\r\n    if (0.5 === fractionParameter) {\r\n      quaternion0.plus(quaternion1, result);\r\n      result.normalizeQuaternion();\r\n      return result;\r\n    }\r\n\r\n    const q0 = quaternion0.clone();\r\n    const q1 = quaternion1.clone();\r\n    let dot = quaternion0.dotProduct(quaternion1);\r\n\r\n    // prevent interpolation through the longer great arc\r\n    if (dot < 0.0) {\r\n      q1.negate(q1);\r\n      dot = -dot;\r\n    }\r\n\r\n    // if nearly parallel, use interpolate and renormalize .\r\n    if (dot > maxSafeCosine) {\r\n      q0.interpolate(fractionParameter, q1, result);\r\n      result.normalizeQuaternion();\r\n      return result;\r\n    }\r\n\r\n    // safety check\r\n    if (dot < -1.0)\r\n      dot = -1.0;\r\n    else if (dot > 1.0)\r\n      dot = 1.0;\r\n\r\n    // create orthonormal basis {q0, q2}\r\n    const q2 = new Point4d();\r\n    q1.plusScaled(q0, -dot, q2);  //  bsiDPoint4d_addScaledDPoint4d(& q2, & q1, & q0, -dot);\r\n    q2.normalizeQuaternion();\r\n\r\n    const angle = Math.acos(dot);\r\n    const angleOfInterpolation = angle * fractionParameter;\r\n    result = Point4d.createAdd2Scaled(q0, Math.cos(angleOfInterpolation), q2, Math.sin(angleOfInterpolation));\r\n    return result;\r\n  }\r\n  /** Measure the \"angle\" between two points, using all 4 components in the dot product that\r\n   * gives the cosine of the angle.\r\n   */\r\n  public radiansToPoint4dXYZW(other: Point4d): number | undefined {\r\n    const magA = this.magnitudeXYZW();\r\n    const magB = other.magnitudeXYZW();\r\n    const dot = this.dotProduct(other);  // == cos (theta) * magA * magB\r\n    const cos = Geometry.conditionalDivideFraction(dot, magA * magB);\r\n    if (cos === undefined)\r\n      return undefined;\r\n    return Math.acos(cos);\r\n  }\r\n}\r\n"]}

package/lib/esm/numerics/BandedSystem.d.ts

@@ -0,0 +1,39 @@
+/**
+ * `BandedSystem` is a class with static methods for solving banded linear systems, such as in computing
+ * Bspline poles for pass-through points
+ * @internal
+ */
+export declare class BandedSystem {
+    /** apply LU decomposition to a banded system */
+    static decomposeLU(numRow: number, bw: number, /* band width */ data: Float64Array): boolean;
+    /**
+     *
+     * @param sum evolving sum.  sum.length
+     * @param source data being added
+     * @param sourceRow row in source.  Plain offset is sourceRow * sum.length
+     * @param scale scale factor to apply.
+     */
+    private static arrayAddScaledBlock;
+    private static blockAssignBlockMinusArray;
+    private static blockSumOfScaledBlockScaledArray;
+    /**
+     * Solve a linear system A*X=B where
+     * * A is nominally an `numRow*numRow` matrix, but is stored in banded row-major form
+     * * The band storage is `bw` numbers per row, with the middle value being the diagonal of that row.
+     *    * Hence rows near top and bottom have band values `outside` the matrix.
+     * * The right hand side is an `numRow*numRHS` matrix in row-major order.
+     * @param numRow number of rows (and columns) of the nominal full matrix.
+     * @param bw total bandwidth (diagonal + equal number of values to left and right)
+     * @param matrix the banded matrix, as packed row-major
+     * @param numRHS the number of right hand sides.
+     * @param rhs the right hand sides
+     */
+    static solveBandedSystemMultipleRHS(numRow: number, bw: number, /* band width */ matrix: Float64Array, numRHS: number, // number of components in each RHS row.
+    rhs: Float64Array): Float64Array | undefined;
+    /**
+     * Multiply a banded numRow*numRow matrix times a full numRow*numRHS, return as new matrix
+     */
+    static multiplyBandedTimesFull(numRow: number, bw: number, /* band width */ bandedMatrix: Float64Array, numRHS: number, // number of components in each RHS row.
+    rhs: Float64Array): Float64Array;
+}
+//# sourceMappingURL=BandedSystem.d.ts.map

package/lib/esm/numerics/BandedSystem.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"BandedSystem.d.ts","sourceRoot":"","sources":["../../../src/numerics/BandedSystem.ts"],"names":[],"mappings":"AAKA;;;;GAIG;AACH,qBAAa,YAAY;IACvB,gDAAgD;WAClC,WAAW,CACvB,MAAM,EAAE,MAAM,EACd,EAAE,EAAE,MAAM,EAAI,gBAAgB,CAC9B,IAAI,EAAE,YAAY,GACjB,OAAO;IAwCV;;;;;;OAMG;IACH,OAAO,CAAC,MAAM,CAAC,mBAAmB;IAQlC,OAAO,CAAC,MAAM,CAAC,0BAA0B;IASzC,OAAO,CAAC,MAAM,CAAC,gCAAgC;IAS/C;;;;;;;;;;;OAWG;WACW,4BAA4B,CACxC,MAAM,EAAE,MAAM,EACd,EAAE,EAAE,MAAM,EAAI,gBAAgB,CAC9B,MAAM,EAAE,YAAY,EACpB,MAAM,EAAE,MAAM,EAAE,wCAAwC;IACxD,GAAG,EAAE,YAAY,GAChB,YAAY,GAAG,SAAS;IAyC3B;;OAEG;WACW,uBAAuB,CACnC,MAAM,EAAE,MAAM,EACd,EAAE,EAAE,MAAM,EAAI,gBAAgB,CAC9B,YAAY,EAAE,YAAY,EAC1B,MAAM,EAAE,MAAM,EAAE,wCAAwC;IACxD,GAAG,EAAE,YAAY,GAChB,YAAY;CAsBhB"}

package/lib/esm/numerics/BandedSystem.js

@@ -0,0 +1,156 @@
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+import { Geometry } from "../Geometry";
+/**
+ * `BandedSystem` is a class with static methods for solving banded linear systems, such as in computing
+ * Bspline poles for pass-through points
+ * @internal
+ */
+export class BandedSystem {
+    /** apply LU decomposition to a banded system */
+    static decomposeLU(numRow, bw, /* band width */ data) {
+        const n = numRow - 1;
+        const sbw = Math.floor(bw / 2); // ASSUMES bw is odd?
+        let sum;
+        // Phase 1:
+        //   [A b C]
+        //   [d q f]
+        //   [G h I]
+        // q is a diagonal (pivot
+        // d, f are row vectors
+        // A,C,I are blocks
+        // b,h are column vectors.
+        // Phase 1:  [q,f] -= d * [b,C]
+        // Phase 2: h = (h- G*b)/q
+        // This is standard gaussian elimination, but in row-think rather than the usual column-think
+        for (let i = 0; i <= n; i++) {
+            const jh = Math.min(n, i + sbw);
+            for (let j = i; j <= jh; j++) {
+                const kl = Math.max(0, j - sbw);
+                sum = 0.0;
+                for (let k = kl; k < i; k++)
+                    sum += data[i * bw + k - i + sbw] * data[k * bw + j - k + sbw];
+                data[i * bw + j - i + sbw] -= sum;
+            }
+            for (let j = i + 1; j <= jh; j++) {
+                const kl = Math.max(0, j - sbw);
+                sum = 0.0;
+                for (let k = kl; k < i; k++)
+                    sum += data[j * bw + k - j + sbw] * data[k * bw + i - k + sbw];
+                if (Math.abs(data[i * bw + sbw]) < 1e-9) // TODO -- tolerance !!!
+                    return false;
+                data[j * bw + i - j + sbw] = (data[j * bw + i - j + sbw] - sum) / data[i * bw + sbw];
+            }
+        }
+        return true;
+    }
+    /**
+     *
+     * @param sum evolving sum.  sum.length
+     * @param source data being added
+     * @param sourceRow row in source.  Plain offset is sourceRow * sum.length
+     * @param scale scale factor to apply.
+     */
+    static arrayAddScaledBlock(sum, source, sourceRow, scale) {
+        const n = sum.length;
+        let k = n * sourceRow;
+        for (let i = 0; i < n; i++, k++) {
+            sum[i] += source[k] * scale;
+        }
+    }
+    //   dest[destRow][*] = sourceA[sourceRow][*] - sourceB[*]
+    static blockAssignBlockMinusArray(dest, destRow, sourceA, sourceARow, sourceB) {
+        const n = sourceB.length;
+        let destIndex = destRow * n;
+        let sourceIndex = sourceARow * n;
+        for (let i = 0; i < n; i++, sourceIndex++, destIndex++) {
+            dest[destIndex] = sourceA[sourceIndex] - sourceB[i];
+        }
+    }
+    //   dest[destRow][*] = sourceA[sourceBRow][*] * scaleA - sourceB[*] * scaleB
+    static blockSumOfScaledBlockScaledArray(dest, destRow, sourceA, sourceBRow, scaleA, sourceB, scaleB) {
+        const n = sourceB.length;
+        let destIndex = destRow * n;
+        let sourceBIndex = sourceBRow * n;
+        for (let i = 0; i < n; i++, sourceBIndex++, destIndex++) {
+            dest[destIndex] = sourceA[sourceBIndex] * scaleA + sourceB[i] * scaleB;
+        }
+    }
+    /**
+     * Solve a linear system A*X=B where
+     * * A is nominally an `numRow*numRow` matrix, but is stored in banded row-major form
+     * * The band storage is `bw` numbers per row, with the middle value being the diagonal of that row.
+     *    * Hence rows near top and bottom have band values `outside` the matrix.
+     * * The right hand side is an `numRow*numRHS` matrix in row-major order.
+     * @param numRow number of rows (and columns) of the nominal full matrix.
+     * @param bw total bandwidth (diagonal + equal number of values to left and right)
+     * @param matrix the banded matrix, as packed row-major
+     * @param numRHS the number of right hand sides.
+     * @param rhs the right hand sides
+     */
+    static solveBandedSystemMultipleRHS(numRow, bw, /* band width */ matrix, numRHS, // number of components in each RHS row.
+    rhs) {
+        if (!this.decomposeLU(numRow, bw, matrix))
+            return undefined;
+        const n = numRow - 1;
+        const sbw = Math.floor(bw / 2);
+        const rhsRowS = new Float64Array(numRHS);
+        /* Compute solution vector */
+        // Z is solution vector . . .
+        const reducedRHS = new Float64Array(numRHS * numRow);
+        const result = new Float64Array(numRHS * numRow);
+        for (let i = 0; i <= n; i++) {
+            rhsRowS.fill(0);
+            const jl = Math.max(0, i - sbw);
+            for (let j = jl; j < i; j++) {
+                this.arrayAddScaledBlock(rhsRowS, reducedRHS, j, matrix[i * bw + j - i + sbw]);
+                // S.SumOf(S, Z[j], data[i * bw + j - i + sbw]);
+            }
+            this.blockAssignBlockMinusArray(reducedRHS, i, rhs, i, rhsRowS);
+        }
+        for (let i = n; i >= 0; i--) {
+            const fact = Geometry.conditionalDivideCoordinate(1.0, matrix[i * bw + sbw]);
+            if (fact === undefined)
+                return undefined;
+            rhsRowS.fill(0);
+            const jh = Math.min(n, i + sbw);
+            for (let j = i + 1; j <= jh; j++) {
+                // S.SumOf(S, Q[j], data[i * bw + j - i + sbw]);
+                this.arrayAddScaledBlock(rhsRowS, result, j, matrix[i * bw + j - i + sbw]);
+            }
+            this.blockSumOfScaledBlockScaledArray(result, i, reducedRHS, i, fact, rhsRowS, -fact);
+            // Q[i].SumOf(O, Z[i], fact, S, -fact);
+        }
+        return result;
+    }
+    /**
+     * Multiply a banded numRow*numRow matrix times a full numRow*numRHS, return as new matrix
+     */
+    static multiplyBandedTimesFull(numRow, bw, /* band width */ bandedMatrix, numRHS, // number of components in each RHS row.
+    rhs) {
+        const result = new Float64Array(rhs.length);
+        const halfBandWidth = Math.floor(bw / 2);
+        let sum, k0, k1, kRef;
+        for (let pivot = 0; pivot < numRow; pivot++) {
+            // "k" vars are nominal column indices in the matrix, and nominal row indices in the rhs.
+            k0 = pivot - halfBandWidth;
+            if (k0 < 0)
+                k0 = 0;
+            k1 = pivot + halfBandWidth + 1;
+            if (k1 > numRow)
+                k1 = numRow;
+            kRef = halfBandWidth + pivot * (bw - 1);
+            for (let m = 0; m < numRHS; m++) {
+                sum = 0;
+                for (let k = k0; k < k1; k++) {
+                    sum += bandedMatrix[kRef + k] * rhs[k * numRHS + m];
+                }
+                result[pivot * numRHS + m] = sum;
+            }
+        }
+        return result;
+    }
+}
+//# sourceMappingURL=BandedSystem.js.map

package/lib/esm/numerics/BandedSystem.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"BandedSystem.js","sourceRoot":"","sources":["../../../src/numerics/BandedSystem.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F,OAAO,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AACvC;;;;GAIG;AACH,MAAM,OAAO,YAAY;IACvB,gDAAgD;IACzC,MAAM,CAAC,WAAW,CACvB,MAAc,EACd,EAAU,EAAI,gBAAgB,CAC9B,IAAkB;QAElB,MAAM,CAAC,GAAG,MAAM,GAAG,CAAC,CAAC;QACrB,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,qBAAqB;QACrD,IAAI,GAAG,CAAC;QACR,WAAW;QACX,YAAY;QACZ,YAAY;QACZ,YAAY;QACZ,yBAAyB;QACzB,uBAAuB;QACvB,mBAAmB;QACnB,0BAA0B;QAC1B,+BAA+B;QAC/B,0BAA0B;QAC1B,6FAA6F;QAC7F,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;YAC3B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC;YAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE;gBAC5B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC;gBAChC,GAAG,GAAG,GAAG,CAAC;gBACV,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBACzB,GAAG,IAAI,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;gBACjE,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC;aACnC;YAED,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE;gBAChC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC;gBAChC,GAAG,GAAG,GAAG,CAAC;gBACV,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBACzB,GAAG,IAAI,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;gBAEjE,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,EAAI,wBAAwB;oBACjE,OAAO,KAAK,CAAC;gBAEf,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC;aACtF;SACF;QACD,OAAO,IAAI,CAAC;IACd,CAAC;IAED;;;;;;OAMG;IACK,MAAM,CAAC,mBAAmB,CAAC,GAAiB,EAAE,MAAoB,EAAE,SAAiB,EAAE,KAAa;QAC1G,MAAM,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC;QACrB,IAAI,CAAC,GAAG,CAAC,GAAG,SAAS,CAAC;QACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE;YAC/B,GAAG,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;SAC7B;IACH,CAAC;IACD,0DAA0D;IAClD,MAAM,CAAC,0BAA0B,CAAC,IAAkB,EAAE,OAAe,EAAE,OAAqB,EAAE,UAAkB,EAAE,OAAqB;QAC7I,MAAM,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC;QACzB,IAAI,SAAS,GAAG,OAAO,GAAG,CAAC,CAAC;QAC5B,IAAI,WAAW,GAAG,UAAU,GAAG,CAAC,CAAC;QACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,WAAW,EAAE,EAAE,SAAS,EAAE,EAAE;YACtD,IAAI,CAAC,SAAS,CAAC,GAAG,OAAO,CAAC,WAAW,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;SACrD;IACH,CAAC;IACD,6EAA6E;IACrE,MAAM,CAAC,gCAAgC,CAAC,IAAkB,EAAE,OAAe,EAAE,OAAqB,EAAE,UAAkB,EAAE,MAAc,EAAE,OAAqB,EAAE,MAAc;QACnL,MAAM,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC;QACzB,IAAI,SAAS,GAAG,OAAO,GAAG,CAAC,CAAC;QAC5B,IAAI,YAAY,GAAG,UAAU,GAAG,CAAC,CAAC;QAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,YAAY,EAAE,EAAE,SAAS,EAAE,EAAE;YACvD,IAAI,CAAC,SAAS,CAAC,GAAG,OAAO,CAAC,YAAY,CAAC,GAAG,MAAM,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;SACxE;IACH,CAAC;IAED;;;;;;;;;;;OAWG;IACI,MAAM,CAAC,4BAA4B,CACxC,MAAc,EACd,EAAU,EAAI,gBAAgB,CAC9B,MAAoB,EACpB,MAAc,EAAE,wCAAwC;IACxD,GAAiB;QAEjB,IAAI,CAAC,IAAI,CAAC,WAAW,CAAC,MAAM,EAAE,EAAE,EAAE,MAAM,CAAC;YACvC,OAAO,SAAS,CAAC;QAEnB,MAAM,CAAC,GAAG,MAAM,GAAG,CAAC,CAAC;QACrB,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAC/B,MAAM,OAAO,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,CAAC;QAEzC,6BAA6B;QAC7B,6BAA6B;QAC7B,MAAM,UAAU,GAAG,IAAI,YAAY,CAAC,MAAM,GAAG,MAAM,CAAC,CAAC;QACrD,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,MAAM,GAAG,MAAM,CAAC,CAAC;QAEjD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;YAC3B,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAEhB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC;YAChC,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;gBAC3B,IAAI,CAAC,mBAAmB,CAAC,OAAO,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;gBAC/E,gDAAgD;aACjD;YACD,IAAI,CAAC,0BAA0B,CAAC,UAAU,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;SACjE;QACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;YAC3B,MAAM,IAAI,GAAG,QAAQ,CAAC,2BAA2B,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC;YAC7E,IAAI,IAAI,KAAK,SAAS;gBACpB,OAAO,SAAS,CAAC;YAEnB,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAEhB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC;YAChC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE;gBAChC,gDAAgD;gBAChD,IAAI,CAAC,mBAAmB,CAAC,OAAO,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;aAC5E;YACD,IAAI,CAAC,gCAAgC,CAAC,MAAM,EAAE,CAAC,EAAE,UAAU,EAAE,CAAC,EAAE,IAAI,EAAE,OAAO,EAAE,CAAC,IAAI,CAAC,CAAC;YACtF,uCAAuC;SACxC;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;IACD;;OAEG;IACI,MAAM,CAAC,uBAAuB,CACnC,MAAc,EACd,EAAU,EAAI,gBAAgB,CAC9B,YAA0B,EAC1B,MAAc,EAAE,wCAAwC;IACxD,GAAiB;QAEjB,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;QAC5C,MAAM,aAAa,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QACzC,IAAI,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC;QAEtB,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,MAAM,EAAE,KAAK,EAAE,EAAE;YAC3C,yFAAyF;YACzF,EAAE,GAAG,KAAK,GAAG,aAAa,CAAC;YAC3B,IAAI,EAAE,GAAG,CAAC;gBAAE,EAAE,GAAG,CAAC,CAAC;YACnB,EAAE,GAAG,KAAK,GAAG,aAAa,GAAG,CAAC,CAAC;YAC/B,IAAI,EAAE,GAAG,MAAM;gBAAE,EAAE,GAAG,MAAM,CAAC;YAC7B,IAAI,GAAG,aAAa,GAAG,KAAK,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;YACxC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE;gBAC/B,GAAG,GAAG,CAAC,CAAC;gBACR,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE;oBAC5B,GAAG,IAAI,YAAY,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,MAAM,GAAG,CAAC,CAAC,CAAC;iBACrD;gBACD,MAAM,CAAC,KAAK,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;aAClC;SACF;QACD,OAAO,MAAM,CAAC;IAChB,CAAC;CACF","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\nimport { Geometry } from \"../Geometry\";\r\n/**\r\n * `BandedSystem` is a class with static methods for solving banded linear systems, such as in computing\r\n * Bspline poles for pass-through points\r\n * @internal\r\n */\r\nexport class BandedSystem {\r\n  /** apply LU decomposition to a banded system */\r\n  public static decomposeLU(\r\n    numRow: number,\r\n    bw: number,   /* band width */\r\n    data: Float64Array,\r\n  ): boolean {\r\n    const n = numRow - 1;\r\n    const sbw = Math.floor(bw / 2); // ASSUMES bw is odd?\r\n    let sum;\r\n    // Phase 1:\r\n    //   [A b C]\r\n    //   [d q f]\r\n    //   [G h I]\r\n    // q is a diagonal (pivot\r\n    // d, f are row vectors\r\n    // A,C,I are blocks\r\n    // b,h are column vectors.\r\n    // Phase 1:  [q,f] -= d * [b,C]\r\n    // Phase 2: h = (h- G*b)/q\r\n    // This is standard gaussian elimination, but in row-think rather than the usual column-think\r\n    for (let i = 0; i <= n; i++) {\r\n      const jh = Math.min(n, i + sbw);\r\n      for (let j = i; j <= jh; j++) {\r\n        const kl = Math.max(0, j - sbw);\r\n        sum = 0.0;\r\n        for (let k = kl; k < i; k++)\r\n          sum += data[i * bw + k - i + sbw] * data[k * bw + j - k + sbw];\r\n        data[i * bw + j - i + sbw] -= sum;\r\n      }\r\n\r\n      for (let j = i + 1; j <= jh; j++) {\r\n        const kl = Math.max(0, j - sbw);\r\n        sum = 0.0;\r\n        for (let k = kl; k < i; k++)\r\n          sum += data[j * bw + k - j + sbw] * data[k * bw + i - k + sbw];\r\n\r\n        if (Math.abs(data[i * bw + sbw]) < 1e-9)   // TODO -- tolerance !!!\r\n          return false;\r\n\r\n        data[j * bw + i - j + sbw] = (data[j * bw + i - j + sbw] - sum) / data[i * bw + sbw];\r\n      }\r\n    }\r\n    return true;\r\n  }\r\n\r\n  /**\r\n   *\r\n   * @param sum evolving sum.  sum.length\r\n   * @param source data being added\r\n   * @param sourceRow row in source.  Plain offset is sourceRow * sum.length\r\n   * @param scale scale factor to apply.\r\n   */\r\n  private static arrayAddScaledBlock(sum: Float64Array, source: Float64Array, sourceRow: number, scale: number) {\r\n    const n = sum.length;\r\n    let k = n * sourceRow;\r\n    for (let i = 0; i < n; i++, k++) {\r\n      sum[i] += source[k] * scale;\r\n    }\r\n  }\r\n  //   dest[destRow][*] = sourceA[sourceRow][*] - sourceB[*]\r\n  private static blockAssignBlockMinusArray(dest: Float64Array, destRow: number, sourceA: Float64Array, sourceARow: number, sourceB: Float64Array) {\r\n    const n = sourceB.length;\r\n    let destIndex = destRow * n;\r\n    let sourceIndex = sourceARow * n;\r\n    for (let i = 0; i < n; i++, sourceIndex++, destIndex++) {\r\n      dest[destIndex] = sourceA[sourceIndex] - sourceB[i];\r\n    }\r\n  }\r\n  //   dest[destRow][*] = sourceA[sourceBRow][*] * scaleA - sourceB[*] * scaleB\r\n  private static blockSumOfScaledBlockScaledArray(dest: Float64Array, destRow: number, sourceA: Float64Array, sourceBRow: number, scaleA: number, sourceB: Float64Array, scaleB: number) {\r\n    const n = sourceB.length;\r\n    let destIndex = destRow * n;\r\n    let sourceBIndex = sourceBRow * n;\r\n    for (let i = 0; i < n; i++, sourceBIndex++, destIndex++) {\r\n      dest[destIndex] = sourceA[sourceBIndex] * scaleA + sourceB[i] * scaleB;\r\n    }\r\n  }\r\n\r\n  /**\r\n   * Solve a linear system A*X=B where\r\n   * * A is nominally an `numRow*numRow` matrix, but is stored in banded row-major form\r\n   * * The band storage is `bw` numbers per row, with the middle value being the diagonal of that row.\r\n   *    * Hence rows near top and bottom have band values `outside` the matrix.\r\n   * * The right hand side is an `numRow*numRHS` matrix in row-major order.\r\n   * @param numRow number of rows (and columns) of the nominal full matrix.\r\n   * @param bw total bandwidth (diagonal + equal number of values to left and right)\r\n   * @param matrix the banded matrix, as packed row-major\r\n   * @param numRHS the number of right hand sides.\r\n   * @param rhs the right hand sides\r\n   */\r\n  public static solveBandedSystemMultipleRHS(\r\n    numRow: number,\r\n    bw: number,   /* band width */\r\n    matrix: Float64Array,\r\n    numRHS: number, // number of components in each RHS row.\r\n    rhs: Float64Array, // RHS data, packed, overwritten by solution\r\n  ): Float64Array | undefined {\r\n    if (!this.decomposeLU(numRow, bw, matrix))\r\n      return undefined;\r\n\r\n    const n = numRow - 1;\r\n    const sbw = Math.floor(bw / 2);\r\n    const rhsRowS = new Float64Array(numRHS);\r\n\r\n    /* Compute solution vector */\r\n    // Z is solution vector . . .\r\n    const reducedRHS = new Float64Array(numRHS * numRow);\r\n    const result = new Float64Array(numRHS * numRow);\r\n\r\n    for (let i = 0; i <= n; i++) {\r\n      rhsRowS.fill(0);\r\n\r\n      const jl = Math.max(0, i - sbw);\r\n      for (let j = jl; j < i; j++) {\r\n        this.arrayAddScaledBlock(rhsRowS, reducedRHS, j, matrix[i * bw + j - i + sbw]);\r\n        // S.SumOf(S, Z[j], data[i * bw + j - i + sbw]);\r\n      }\r\n      this.blockAssignBlockMinusArray(reducedRHS, i, rhs, i, rhsRowS);\r\n    }\r\n    for (let i = n; i >= 0; i--) {\r\n      const fact = Geometry.conditionalDivideCoordinate(1.0, matrix[i * bw + sbw]);\r\n      if (fact === undefined)\r\n        return undefined;\r\n\r\n      rhsRowS.fill(0);\r\n\r\n      const jh = Math.min(n, i + sbw);\r\n      for (let j = i + 1; j <= jh; j++) {\r\n        // S.SumOf(S, Q[j], data[i * bw + j - i + sbw]);\r\n        this.arrayAddScaledBlock(rhsRowS, result, j, matrix[i * bw + j - i + sbw]);\r\n      }\r\n      this.blockSumOfScaledBlockScaledArray(result, i, reducedRHS, i, fact, rhsRowS, -fact);\r\n      // Q[i].SumOf(O, Z[i], fact, S, -fact);\r\n    }\r\n\r\n    return result;\r\n  }\r\n  /**\r\n   * Multiply a banded numRow*numRow matrix times a full numRow*numRHS, return as new matrix\r\n   */\r\n  public static multiplyBandedTimesFull(\r\n    numRow: number,\r\n    bw: number,   /* band width */\r\n    bandedMatrix: Float64Array,\r\n    numRHS: number, // number of components in each RHS row.\r\n    rhs: Float64Array, // RHS data, packed, overwritten by solution\r\n  ): Float64Array {\r\n    const result = new Float64Array(rhs.length);\r\n    const halfBandWidth = Math.floor(bw / 2);\r\n    let sum, k0, k1, kRef;\r\n\r\n    for (let pivot = 0; pivot < numRow; pivot++) {\r\n      // \"k\" vars are nominal column indices in the matrix, and nominal row indices in the rhs.\r\n      k0 = pivot - halfBandWidth;\r\n      if (k0 < 0) k0 = 0;\r\n      k1 = pivot + halfBandWidth + 1;\r\n      if (k1 > numRow) k1 = numRow;\r\n      kRef = halfBandWidth + pivot * (bw - 1);\r\n      for (let m = 0; m < numRHS; m++) {\r\n        sum = 0;\r\n        for (let k = k0; k < k1; k++) {\r\n          sum += bandedMatrix[kRef + k] * rhs[k * numRHS + m];\r\n        }\r\n        result[pivot * numRHS + m] = sum;\r\n      }\r\n    }\r\n    return result;\r\n  }\r\n}\r\n"]}