copper3d 1.15.14 → 1.15.16

package/dist/Utils/utils.js

@@ -138,39 +138,113 @@ export function switchEraserSize(size, urls) {
     }
     return url;
 }
-export function createFpsCap(loop, fps = 60) {
-    let targetFps = 0, fpsInterval = 0;
-    let lastTime = 0, lastOverTime = 0, prevOverTime = 0, deltaTime = 0;
-    function updateFps(value) {
-        targetFps = value;
-        fpsInterval = 1000 / targetFps;
-    }
-    updateFps(fps);
-    return {
-        // the targeted frame rate
-        get fps() {
-            return targetFps;
-        },
-        set fps(value) {
-            updateFps(value);
-        },
-        // the frame-capped loop function
-        loop: function (_scope, time) {
-            if (!!time) {
-                console.log(time, lastTime);
-                console.log(deltaTime);
-                deltaTime = time - lastTime;
-                console.log(deltaTime);
-                if (deltaTime < fpsInterval) {
-                    return;
-                }
-                prevOverTime = lastOverTime;
-                lastOverTime = deltaTime % fpsInterval;
-                lastTime = time - lastOverTime;
-                deltaTime -= prevOverTime;
+/**
+ * Cubic-Lagrange basis function
+ * @param x
+ * @returns
+ */
+export function L3(x) {
+    let L1 = 1 - x;
+    let L2 = x;
+    let sc = 9 / 2;
+    return [
+        0.5 * L1 * (3 * L1 - 1) * (3 * L1 - 2),
+        sc * L1 * L2 * (3 * L1 - 1),
+        sc * L1 * L2 * (3 * L2 - 1),
+        0.5 * L2 * (3 * L2 - 1) * (3 * L2 - 2),
+    ];
+}
+/**
+ * Cubic-Hermite basis function.
+ * @param x
+ * @returns
+ */
+export function H3(x) {
+    let x2 = x * x;
+    return [
+        1 - 3 * x2 + 2 * x * x2,
+        x * (x - 1) * (x - 1),
+        x2 * (3 - 2 * x),
+        x2 * (x - 1),
+    ];
+}
+/**
+ * To calculate the weights for each element of Xi using the cubic Lagrange basis functions.
+ * @param Xi
+ * @returns
+ */
+export function getWightsL3L3L3(Xi) {
+    let W0, W1, W2;
+    W0 = L3(Xi[0]);
+    W1 = L3(Xi[1]);
+    W2 = L3(Xi[2]);
+    let w = [];
+    for (let k2 in W2) {
+        for (let k1 in W1) {
+            for (let k0 in W0) {
+                w.push(W0[k0] * W1[k1] * W2[k2]);
             }
-            return loop(_scope, deltaTime);
-        },
-    };
+        }
+    }
+    return w;
+}
+export function getWightsH3H3H3(Xi) {
+    let mixIdx = [
+        [0, 0],
+        [1, 0],
+        [0, 1],
+        [1, 1],
+        [2, 0],
+        [3, 0],
+        [2, 1],
+        [3, 1],
+        [0, 2],
+        [1, 2],
+        [0, 3],
+        [1, 3],
+        [2, 2],
+        [3, 2],
+        [2, 3],
+        [3, 3],
+    ];
+    let weights = [];
+    let W0, W1;
+    for (let idx in Xi) {
+        W0 = H3(Xi[idx][0]);
+        W1 = H3(Xi[idx][1]);
+        let w = [];
+        for (let k in mixIdx) {
+            w.push(W0[mixIdx[k][0]] * W1[mixIdx[k][1]]);
+        }
+        weights.push(w);
+    }
+    return weights;
+}
+/**
+ * Euclidean distance n dimensions
+ * @param x0
+ * @param x1
+ * @returns
+ */
+export function calcDistance(x0, x1) {
+    let dx2 = 0;
+    for (let i = 0; i < x0.length; i++) {
+        let dx = x0[i] - x1[i];
+        dx2 += dx * dx;
+    }
+    return Math.sqrt(dx2);
+}
+export function perturbRandom(x, dx) {
+    let xp = [];
+    for (let i = 0; i < x.length; i++) {
+        xp.push(x[i] + 2 * dx * (Math.random() - 0.5));
+        if (xp[i] < 0) {
+            xp[i] = 0;
+        }
+        else if (xp[i] > 1) {
+            xp[i] = 1;
+        }
+    }
+    return xp;
 }
 //# sourceMappingURL=utils.js.map

package/dist/Utils/utils.js.map

@@ -1 +1 @@
-{"version":3,"file":"utils.js","sourceRoot":"","sources":["../../src/Utils/utils.ts"],"names":[],"mappings":"AAAA,MAAM,UAAU,mBAAmB,CACjC,SAA6C;IAE7C,MAAM,iBAAiB,GACrB,QAAQ,CAAC,iBAAiB,IAAK,QAAgB,CAAC,uBAAuB,CAAC;IAC1E,IAAI,CAAC,iBAAiB,EAAE;QACtB,IAAI,SAAS,CAAC,iBAAiB,EAAE;YAC/B,SAAS,CAAC,iBAAiB,EAAE,CAAC;SAC/B;aAAM,IAAK,SAAiB,CAAC,uBAAuB,EAAE;YACpD,SAAiB,CAAC,uBAAuB,EAAE,CAAC;SAC9C;KACF;SAAM;QACL,IAAI,QAAQ,CAAC,cAAc,EAAE;YAC3B,QAAQ,CAAC,cAAc,EAAE,CAAC;SAC3B;aAAM,IAAK,QAAgB,CAAC,oBAAoB,EAAE;YAChD,QAAgB,CAAC,oBAAoB,EAAE,CAAC;SAC1C;KACF;AACH,CAAC;AAED,MAAM,UAAU,KAAK;IACnB,OAAO,CACL;QACE,gBAAgB;QAChB,kBAAkB;QAClB,gBAAgB;QAChB,MAAM;QACN,QAAQ;QACR,MAAM;KACP,CAAC,QAAQ,CAAC,SAAS,CAAC,QAAQ,CAAC;QAC9B,CAAC,SAAS,CAAC,SAAS,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,YAAY,IAAI,QAAQ,CAAC,CAClE,CAAC;AACJ,CAAC;AAED,MAAM,UAAU,iBAAiB,CAC/B,MAAmB,EACnB,QAAiC;IAEjC,MAAM,CAAC,QAAQ,CAAC,CAAC,IAAI,EAAE,EAAE;QACvB,IAAI,CAAE,IAAmB,CAAC,MAAM;YAAE,OAAO;QACzC,IAAI,KAAK,CAAC,OAAO,CAAE,IAAmB,CAAC,QAAQ,CAAC,EAAE;YAChD,QAAQ,CAAE,IAAmB,CAAC,QAAQ,CAAC,CAAC;SACzC;aAAM;YACL,CAAE,IAAmB,CAAC,QAAQ,CAAC,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC;SACnD;IACH,CAAC,CAAC,CAAC;AACL,CAAC;AAED,MAAM,UAAU,OAAO;IACrB,IAAI,gBAAgB,GAAG,QAAQ,CAAC,aAAa,CAAC,KAAK,CAAC,CAAC;IACrD,MAAM,UAAU,GAAG,QAAQ,CAAC,aAAa,CAAC,KAAK,CAAC,CAAC;IACjD,IAAI,QAAQ,GAAG,QAAQ,CAAC,aAAa,CAAC,KAAK,CAAC,CAAC;IAC7C,gBAAgB,CAAC,KAAK,CAAC,QAAQ,GAAG,UAAU,CAAC;IAC7C,sCAAsC;IACtC,uCAAuC;IACvC,gBAAgB,CAAC,KAAK,CAAC,OAAO,GAAG,MAAM,CAAC;IACxC,QAAQ,CAAC,KAAK,CAAC,KAAK,GAAG,MAAM,CAAC;IAC9B,gBAAgB,CAAC,KAAK,CAAC,MAAM,GAAG,MAAM,CAAC;IACvC,gBAAgB,CAAC,SAAS,GAAG,sBAAsB,CAAC;IACpD,UAAU,CAAC,SAAS,GAAG,6BAA6B,CAAC;IACrD,QAAQ,CAAC,SAAS,GAAG,2BAA2B,CAAC;IAEjD,UAAU,CAAC,SAAS,GAAG,oNAAoN,CAAC;IAE5O,gBAAgB,CAAC,WAAW,CAAC,UAAU,CAAC,CAAC;IACzC,gBAAgB,CAAC,WAAW,CAAC,QAAQ,CAAC,CAAC;IAEvC,OAAO,EAAE,gBAAgB,EAAE,QAAQ,EAAE,CAAC;AACxC,CAAC;AAED,MAAM,UAAU,gBAAgB,CAAC,IAAY,EAAE,IAAe;IAC5D,IAAI,GAAG,GAAG,EAAE,CAAC;IACb,IAAI,CAAC,CAAC,IAAI,IAAI,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE;QAC7B,IAAI,IAAI,IAAI,CAAC,EAAE;YACb,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,kBAAkB,CAAC;SACxC;aAAM,IAAI,CAAC,GAAG,IAAI,IAAI,IAAI,IAAI,CAAC,EAAE;YAChC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,kBAAkB,CAAC;SACxC;aAAM,IAAI,CAAC,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YACjC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,EAAE,CAAC,oBAAoB,CAAC;SAC3C;aAAM;YACL,GAAG,GAAG,OAAO,IAAI,CAAC,EAAE,CAAC,oBAAoB,CAAC;SAC3C;KACF;SAAM;QACL,IAAI,IAAI,IAAI,CAAC,EAAE;YACb,GAAG,GAAG,2HAA2H,CAAC;SACnI;aAAM,IAAI,CAAC,GAAG,IAAI,IAAI,IAAI,IAAI,CAAC,EAAE;YAChC,GAAG,GAAG,2HAA2H,CAAC;SACnI;aAAM,IAAI,CAAC,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YACjC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM;YACL,GAAG,GAAG,8HAA8H,CAAC;SACtI;KACF;IAED,OAAO,GAAG,CAAC;AACb,CAAC;AAED,MAAM,UAAU,YAAY,CAAC,IAAc,EAAE,GAAG,GAAG,EAAE;IACnD,IAAI,SAAS,GAAG,CAAC,EACf,WAAW,GAAG,CAAC,CAAC;IAClB,IAAI,QAAQ,GAAG,CAAC,EACd,YAAY,GAAG,CAAC,EAChB,YAAY,GAAG,CAAC,EAChB,SAAS,GAAG,CAAC,CAAC;IAEhB,SAAS,SAAS,CAAC,KAAa;QAC9B,SAAS,GAAG,KAAK,CAAC;QAClB,WAAW,GAAG,IAAI,GAAG,SAAS,CAAC;IACjC,CAAC;IAED,SAAS,CAAC,GAAG,CAAC,CAAC;IAEf,OAAO;QACL,0BAA0B;QAC1B,IAAI,GAAG;YACL,OAAO,SAAS,CAAC;QACnB,CAAC;QACD,IAAI,GAAG,CAAC,KAAK;YACX,SAAS,CAAC,KAAK,CAAC,CAAC;QACnB,CAAC;QAED,iCAAiC;QACjC,IAAI,EAAE,UAAU,MAAW,EAAE,IAAY;YACvC,IAAI,CAAC,CAAC,IAAI,EAAE;gBACV,OAAO,CAAC,GAAG,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;gBAC5B,OAAO,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;gBAEvB,SAAS,GAAG,IAAI,GAAG,QAAQ,CAAC;gBAC5B,OAAO,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;gBAEvB,IAAI,SAAS,GAAG,WAAW,EAAE;oBAC3B,OAAO;iBACR;gBAED,YAAY,GAAG,YAAY,CAAC;gBAC5B,YAAY,GAAG,SAAS,GAAG,WAAW,CAAC;gBACvC,QAAQ,GAAG,IAAI,GAAG,YAAY,CAAC;gBAE/B,SAAS,IAAI,YAAY,CAAC;aAC3B;YAED,OAAO,IAAI,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC;QACjC,CAAC;KACF,CAAC;AACJ,CAAC"}
+{"version":3,"file":"utils.js","sourceRoot":"","sources":["../../src/Utils/utils.ts"],"names":[],"mappings":"AAAA,MAAM,UAAU,mBAAmB,CACjC,SAA6C;IAE7C,MAAM,iBAAiB,GACrB,QAAQ,CAAC,iBAAiB,IAAK,QAAgB,CAAC,uBAAuB,CAAC;IAC1E,IAAI,CAAC,iBAAiB,EAAE;QACtB,IAAI,SAAS,CAAC,iBAAiB,EAAE;YAC/B,SAAS,CAAC,iBAAiB,EAAE,CAAC;SAC/B;aAAM,IAAK,SAAiB,CAAC,uBAAuB,EAAE;YACpD,SAAiB,CAAC,uBAAuB,EAAE,CAAC;SAC9C;KACF;SAAM;QACL,IAAI,QAAQ,CAAC,cAAc,EAAE;YAC3B,QAAQ,CAAC,cAAc,EAAE,CAAC;SAC3B;aAAM,IAAK,QAAgB,CAAC,oBAAoB,EAAE;YAChD,QAAgB,CAAC,oBAAoB,EAAE,CAAC;SAC1C;KACF;AACH,CAAC;AAED,MAAM,UAAU,KAAK;IACnB,OAAO,CACL;QACE,gBAAgB;QAChB,kBAAkB;QAClB,gBAAgB;QAChB,MAAM;QACN,QAAQ;QACR,MAAM;KACP,CAAC,QAAQ,CAAC,SAAS,CAAC,QAAQ,CAAC;QAC9B,CAAC,SAAS,CAAC,SAAS,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,YAAY,IAAI,QAAQ,CAAC,CAClE,CAAC;AACJ,CAAC;AAED,MAAM,UAAU,iBAAiB,CAC/B,MAAmB,EACnB,QAAiC;IAEjC,MAAM,CAAC,QAAQ,CAAC,CAAC,IAAI,EAAE,EAAE;QACvB,IAAI,CAAE,IAAmB,CAAC,MAAM;YAAE,OAAO;QACzC,IAAI,KAAK,CAAC,OAAO,CAAE,IAAmB,CAAC,QAAQ,CAAC,EAAE;YAChD,QAAQ,CAAE,IAAmB,CAAC,QAAQ,CAAC,CAAC;SACzC;aAAM;YACL,CAAE,IAAmB,CAAC,QAAQ,CAAC,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC;SACnD;IACH,CAAC,CAAC,CAAC;AACL,CAAC;AAED,MAAM,UAAU,OAAO;IACrB,IAAI,gBAAgB,GAAG,QAAQ,CAAC,aAAa,CAAC,KAAK,CAAC,CAAC;IACrD,MAAM,UAAU,GAAG,QAAQ,CAAC,aAAa,CAAC,KAAK,CAAC,CAAC;IACjD,IAAI,QAAQ,GAAG,QAAQ,CAAC,aAAa,CAAC,KAAK,CAAC,CAAC;IAC7C,gBAAgB,CAAC,KAAK,CAAC,QAAQ,GAAG,UAAU,CAAC;IAC7C,sCAAsC;IACtC,uCAAuC;IACvC,gBAAgB,CAAC,KAAK,CAAC,OAAO,GAAG,MAAM,CAAC;IACxC,QAAQ,CAAC,KAAK,CAAC,KAAK,GAAG,MAAM,CAAC;IAC9B,gBAAgB,CAAC,KAAK,CAAC,MAAM,GAAG,MAAM,CAAC;IACvC,gBAAgB,CAAC,SAAS,GAAG,sBAAsB,CAAC;IACpD,UAAU,CAAC,SAAS,GAAG,6BAA6B,CAAC;IACrD,QAAQ,CAAC,SAAS,GAAG,2BAA2B,CAAC;IAEjD,UAAU,CAAC,SAAS,GAAG,oNAAoN,CAAC;IAE5O,gBAAgB,CAAC,WAAW,CAAC,UAAU,CAAC,CAAC;IACzC,gBAAgB,CAAC,WAAW,CAAC,QAAQ,CAAC,CAAC;IAEvC,OAAO,EAAE,gBAAgB,EAAE,QAAQ,EAAE,CAAC;AACxC,CAAC;AAED,MAAM,UAAU,gBAAgB,CAAC,IAAY,EAAE,IAAe;IAC5D,IAAI,GAAG,GAAG,EAAE,CAAC;IACb,IAAI,CAAC,CAAC,IAAI,IAAI,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE;QAC7B,IAAI,IAAI,IAAI,CAAC,EAAE;YACb,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,kBAAkB,CAAC;SACxC;aAAM,IAAI,CAAC,GAAG,IAAI,IAAI,IAAI,IAAI,CAAC,EAAE;YAChC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,kBAAkB,CAAC;SACxC;aAAM,IAAI,CAAC,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YACjC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC,CAAC,oBAAoB,CAAC;SAC1C;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,OAAO,IAAI,CAAC,EAAE,CAAC,oBAAoB,CAAC;SAC3C;aAAM;YACL,GAAG,GAAG,OAAO,IAAI,CAAC,EAAE,CAAC,oBAAoB,CAAC;SAC3C;KACF;SAAM;QACL,IAAI,IAAI,IAAI,CAAC,EAAE;YACb,GAAG,GAAG,2HAA2H,CAAC;SACnI;aAAM,IAAI,CAAC,GAAG,IAAI,IAAI,IAAI,IAAI,CAAC,EAAE;YAChC,GAAG,GAAG,2HAA2H,CAAC;SACnI;aAAM,IAAI,CAAC,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YACjC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM,IAAI,EAAE,GAAG,IAAI,IAAI,IAAI,IAAI,EAAE,EAAE;YAClC,GAAG,GAAG,8HAA8H,CAAC;SACtI;aAAM;YACL,GAAG,GAAG,8HAA8H,CAAC;SACtI;KACF;IAED,OAAO,GAAG,CAAC;AACb,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,EAAE,CAAC,CAAS;IAC1B,IAAI,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,IAAI,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,OAAO;QACL,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QACtC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QAC3B,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QAC3B,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;KACvC,CAAC;AACJ,CAAC;AAED;;;;GAIG;AAEH,MAAM,UAAU,EAAE,CAAC,CAAS;IAC1B,IAAI,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACf,OAAO;QACL,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE;QACvB,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QACrB,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAChB,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;KACb,CAAC;AACJ,CAAC;AAED;;;;GAIG;AAEH,MAAM,UAAU,eAAe,CAAC,EAAY;IAC1C,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;IACf,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACf,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACf,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACf,IAAI,CAAC,GAAG,EAAE,CAAC;IACX,KAAK,IAAI,EAAE,IAAI,EAAE,EAAE;QACjB,KAAK,IAAI,EAAE,IAAI,EAAE,EAAE;YACjB,KAAK,IAAI,EAAE,IAAI,EAAE,EAAE;gBACjB,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;aAClC;SACF;KACF;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,UAAU,eAAe,CAAC,EAAc;IAC5C,IAAI,MAAM,GAAG;QACX,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;QACN,CAAC,CAAC,EAAE,CAAC,CAAC;KACP,CAAC;IAEF,IAAI,OAAO,GAAG,EAAE,CAAC;IACjB,IAAI,EAAE,EAAE,EAAE,CAAC;IACX,KAAK,IAAI,GAAG,IAAI,EAAE,EAAE;QAClB,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpB,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpB,IAAI,CAAC,GAAG,EAAE,CAAC;QACX,KAAK,IAAI,CAAC,IAAI,MAAM,EAAE;YACpB,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;SAC7C;QACD,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;KACjB;IACD,OAAO,OAAO,CAAC;AACjB,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,YAAY,CAAC,EAAY,EAAE,EAAY;IACrD,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QAClC,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QACvB,GAAG,IAAI,EAAE,GAAG,EAAE,CAAC;KAChB;IACD,OAAO,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;AACxB,CAAC;AAED,MAAM,UAAU,aAAa,CAAC,CAAW,EAAE,EAAU;IACnD,IAAI,EAAE,GAAG,EAAE,CAAC;IACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACjC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,CAAC,MAAM,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC;QAC/C,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE;YACb,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;SACX;aAAM,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE;YACpB,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;SACX;KACF;IACD,OAAO,EAAE,CAAC;AACZ,CAAC"}

package/dist/bundle.esm.js

@@ -38478,6 +38478,69 @@ function switchEraserSize(size, urls) {
         }
     }
     return url;
+}
+/**
+ * Cubic-Lagrange basis function
+ * @param x
+ * @returns
+ */
+function L3(x) {
+    let L1 = 1 - x;
+    let L2 = x;
+    let sc = 9 / 2;
+    return [
+        0.5 * L1 * (3 * L1 - 1) * (3 * L1 - 2),
+        sc * L1 * L2 * (3 * L1 - 1),
+        sc * L1 * L2 * (3 * L2 - 1),
+        0.5 * L2 * (3 * L2 - 1) * (3 * L2 - 2),
+    ];
+}
+/**
+ * To calculate the weights for each element of Xi using the cubic Lagrange basis functions.
+ * @param Xi
+ * @returns
+ */
+function getWightsL3L3L3(Xi) {
+    let W0, W1, W2;
+    W0 = L3(Xi[0]);
+    W1 = L3(Xi[1]);
+    W2 = L3(Xi[2]);
+    let w = [];
+    for (let k2 in W2) {
+        for (let k1 in W1) {
+            for (let k0 in W0) {
+                w.push(W0[k0] * W1[k1] * W2[k2]);
+            }
+        }
+    }
+    return w;
+}
+/**
+ * Euclidean distance n dimensions
+ * @param x0
+ * @param x1
+ * @returns
+ */
+function calcDistance(x0, x1) {
+    let dx2 = 0;
+    for (let i = 0; i < x0.length; i++) {
+        let dx = x0[i] - x1[i];
+        dx2 += dx * dx;
+    }
+    return Math.sqrt(dx2);
+}
+function perturbRandom(x, dx) {
+    let xp = [];
+    for (let i = 0; i < x.length; i++) {
+        xp.push(x[i] + 2 * dx * (Math.random() - 0.5));
+        if (xp[i] < 0) {
+            xp[i] = 0;
+        }
+        else if (xp[i] > 1) {
+            xp[i] = 1;
+        }
+    }
+    return xp;
 }
 
 /* eslint no-bitwise: 0 */
@@ -56552,10 +56615,6 @@ class copperRenderer extends baseRenderer {
         this.interval = 1 / this.fps;
         this.preRenderCallbackFunctions = [];
         this.animate = (time) => {
-            // if (this.start) {
-            //   this.fpsCap = createFpsCap(this.render, 10);
-            //   this.start = false;
-            // }
             var _a, _b, _c;
             switch ((_a = this.options) === null || _a === void 0 ? void 0 : _a.fpsMode) {
                 case "1":
@@ -68864,6 +68923,16 @@ class nrrd_tools {
                 else if (this.gui_states.sphere) {
                     let mouseX = e.offsetX;
                     let mouseY = e.offsetY;
+                    //  record mouseX,Y, and enable crosshair function
+                    this.nrrd_states.sphereOrigin[this.axis] = [
+                        mouseX,
+                        mouseY,
+                        this.nrrd_states.currentIndex,
+                    ];
+                    this.setUpSphereOrigins(mouseX, mouseY);
+                    this.nrrd_states.cursorPageX = mouseX;
+                    this.nrrd_states.cursorPageY = mouseY;
+                    this.enableCrosshair();
                     // draw circle setup width/height for sphere canvas
                     this.drawSphere(mouseX, mouseY, this.nrrd_states.sphereRadius);
                     this.drawingCanvas.addEventListener("wheel", this.drawingPrameters.handleSphereWheel, true);
@@ -68976,18 +69045,8 @@ class nrrd_tools {
                 else if (this.gui_states.sphere &&
                     !this.nrrd_states.enableCursorChoose) {
                     // let { ctx, canvas } = this.setCurrentLayer();
-                    let mouseX = e.offsetX;
-                    let mouseY = e.offsetY;
-                    this.nrrd_states.sphereOrigin[this.axis] = [
-                        mouseX,
-                        mouseY,
-                        this.nrrd_states.currentIndex,
-                    ];
-                    /************ */
-                    this.setUpSphereOrigins(mouseX, mouseY);
-                    this.nrrd_states.cursorPageX = mouseX;
-                    this.nrrd_states.cursorPageY = mouseY;
-                    this.enableCrosshair();
+                    // let mouseX = e.offsetX;
+                    // let mouseY = e.offsetY;
                     // plan B
                     // findout all index in the sphere radius range in Axial view
                     if (this.nrrd_states.spherePlanB) {
@@ -69227,9 +69286,10 @@ class nrrd_tools {
             }
             // limited the radius max and min
             this.nrrd_states.sphereRadius = Math.max(1, Math.min(this.nrrd_states.sphereRadius, 50));
+            console.log(this.nrrd_states.sphereOrigin[this.axis][0], this.nrrd_states.sphereOrigin[this.axis][1]);
             // get mouse position
-            const mouseX = e.offsetX;
-            const mouseY = e.offsetY;
+            const mouseX = this.nrrd_states.sphereOrigin[this.axis][0];
+            const mouseY = this.nrrd_states.sphereOrigin[this.axis][1];
             this.drawSphere(mouseX, mouseY, this.nrrd_states.sphereRadius);
         };
         return sphereEvent;
@@ -70100,8 +70160,214 @@ class nrrd_tools {
     }
 }
 
+class Node {
+    constructor(id, p) {
+        this.elements = [];
+        this.id = id;
+        this.p = p;
+    }
+}
+class Element {
+    constructor(id, basis, nodes) {
+        this.id = id;
+        this.basis = basis;
+        this.nodes = nodes;
+    }
+}
+class MeshNodeTool {
+    constructor() {
+        this.nodes = {};
+        this.elements = {};
+    }
+    addNode(id, p) {
+        this.nodes[id] = new Node(id, p);
+    }
+    addElement(id, basis, nodeIds) {
+        let nodes = [];
+        for (let nid in nodeIds) {
+            let enid = nodeIds[nid];
+            nodes.push(this.nodes[enid]);
+        }
+        let element = new Element(id, basis, nodes);
+        this.elements[id] = element;
+        for (let idx in element.nodes) {
+            element.nodes[idx].elements.push(element);
+        }
+    }
+    loadMesh(json) {
+        for (let nid in json.nodes) {
+            this.addNode(nid, json.nodes[nid]);
+        }
+        for (let eid in json.elements) {
+            let elem = json.elements[eid];
+            this.addElement(eid, elem.basis, elem.nodes);
+        }
+    }
+    evaluate(elementId, xi) {
+        let w = getWightsL3L3L3(xi);
+        return evaluateElement(this.elements[elementId], w);
+    }
+    find(point, startingNodeId) {
+        let pt = [point.x, point.y, point.z];
+        let nodeMaterialPoints = this.getNodeMaterialPoints(startingNodeId);
+        for (let i in nodeMaterialPoints) {
+            let mpt = nodeMaterialPoints[i];
+            let x = this.evaluate(mpt.elementId, mpt.xi);
+            let r = calcDistance(x, pt);
+            mpt.error = r;
+        }
+        let t0 = new Date().getTime();
+        let n = 0;
+        for (let step = 0; step < 10000; step++) {
+            for (let i in nodeMaterialPoints) {
+                let mpt = nodeMaterialPoints[i];
+                let xi = perturbRandom(mpt.xi, 0.1);
+                let x = this.evaluate(mpt.elementId, xi);
+                let r = calcDistance(x, pt);
+                if (r < mpt.error) {
+                    mpt.error = r;
+                    mpt.xi = xi;
+                    if (r < 1.0) {
+                        console.log("Time: (" +
+                            n +
+                            ") " +
+                            (new Date().getTime() - t0) +
+                            "ms Error: " +
+                            r);
+                        return mpt;
+                    }
+                }
+            }
+            n = step;
+        }
+        let bestIndex = 0;
+        let bestError = nodeMaterialPoints[0].error;
+        for (let i in nodeMaterialPoints) {
+            let mpt = nodeMaterialPoints[i];
+            if (mpt.error < bestError) {
+                bestIndex = Number(i);
+                bestError = mpt.error;
+            }
+        }
+        let mpt = nodeMaterialPoints[bestIndex];
+        console.log("Time: (" +
+            n +
+            ") " +
+            (new Date().getTime() - t0) +
+            "ms Error: " +
+            mpt.error);
+        return mpt;
+    }
+    search(point, startingNodeId, tol) {
+        let pt = [point.x, point.y, point.z];
+        let nodeMaterialPoints = this.getNodeMaterialPoints(startingNodeId);
+        for (let i in nodeMaterialPoints) {
+            let mpt = nodeMaterialPoints[i];
+            let x = this.evaluate(mpt.elementId, mpt.xi);
+            mpt.error = calcDistance(x, pt);
+        }
+        let t0 = new Date().getTime();
+        let xiStepSize = [0.1, 0.01, 0.001];
+        for (let i in nodeMaterialPoints) {
+            let mpt = nodeMaterialPoints[i];
+            let numDims = mpt.xi.length;
+            let nIter = 0;
+            let go = true;
+            while (go || nIter < 10) {
+                let minErr = 1 * mpt.error;
+                let minXi;
+                for (let xiIdx in xiStepSize) {
+                    let dXi = xiStepSize[xiIdx];
+                    for (let dim = 0; dim < numDims; dim++) {
+                        let xi = mpt.xi.slice();
+                        xi[dim] += dXi;
+                        if (xi[dim] >= 0 && xi[dim] <= 1) {
+                            let x = this.evaluate(mpt.elementId, mpt.xi);
+                            let err = calcDistance(x, pt);
+                            if (err < minErr) {
+                                minErr = err;
+                                minXi = xi;
+                                if (err < tol) {
+                                    mpt.error = err;
+                                    mpt.xi = xi;
+                                    console.log("Time 1: " +
+                                        (new Date().getTime() - t0) +
+                                        "ms Error: " +
+                                        mpt.error);
+                                    return mpt;
+                                }
+                            }
+                        }
+                        xi = mpt.xi.slice();
+                        xi[dim] -= dXi;
+                        if (xi[dim] >= 0 && xi[dim] <= 1) {
+                            let x = this.evaluate(mpt.elementId, mpt.xi);
+                            let err = calcDistance(x, pt);
+                            if (err < minErr) {
+                                minErr = err;
+                                minXi = xi;
+                                if (err < tol) {
+                                    mpt.error = err;
+                                    mpt.xi = xi;
+                                    console.log("Time 2: " +
+                                        (new Date().getTime() - t0) +
+                                        "ms Error: " +
+                                        mpt.error);
+                                    return mpt;
+                                }
+                            }
+                        }
+                    }
+                }
+                if (minXi === null || minXi === undefined) {
+                    go = false;
+                }
+                else {
+                    mpt.xi = minXi;
+                    mpt.error = minErr;
+                }
+                nIter++;
+            }
+        }
+        let bestIndex = 0;
+        let bestError = nodeMaterialPoints[0].error;
+        for (let i in nodeMaterialPoints) {
+            let mpt = nodeMaterialPoints[i];
+            if (mpt.error < bestError) {
+                bestIndex = Number(i);
+                bestError = mpt.error;
+            }
+        }
+        let mpt = nodeMaterialPoints[bestIndex];
+        console.log("Time 3: " + (new Date().getTime() - t0) + "ms Error: " + mpt.error);
+        return mpt;
+    }
+    getNodeMaterialPoints(nodeId) {
+        let materialPoints = [];
+        let node = this.nodes[nodeId];
+        for (let eid in node.elements) {
+            let element = node.elements[eid];
+            let nodeIdx = element.nodes.indexOf(node);
+            let xi0 = (nodeIdx % 4) / 3;
+            let xi1 = Math.floor((nodeIdx % 16) / 4) / 3;
+            let xi2 = Math.floor(nodeIdx / 16) / 3;
+            materialPoints.push({ elementId: element.id, xi: [xi0, xi1, xi2] });
+        }
+        return materialPoints;
+    }
+}
+function evaluateElement(element, weights) {
+    let x = [0, 0, 0];
+    for (let f = 0; f < 3; f++) {
+        for (let i = 0; i < weights.length; i++) {
+            x[f] += weights[i] * element.nodes[i].p[f];
+        }
+    }
+    return x;
+}
+
 // import * as kiwrious from "copper3d_plugin_heart_k";
-const REVISION = "v1.15.14";
-console.log("%cCopper3D Visualisation %cBeta:v1.15.14", "padding: 3px;color:white; background:#023047", "padding: 3px;color:white; background:#f50a25");
+const REVISION = "v1.15.16";
+console.log("%cCopper3D Visualisation %cBeta:v1.15.16", "padding: 3px;color:white; background:#023047", "padding: 3px;color:white; background:#f50a25");
 
-export { CameraViewPoint, Copper3dTrackballControls, REVISION, addBoxHelper, addLabelToScene, configKiwriousHeart, convert3DPostoScreenPos, convertScreenPosto3DPos, copperMScene, copperMSceneRenderer, copperRenderer, copperRendererOnDemond, copperScene, copperSceneOnDemond, createTexture2D_NRRD, fullScreenListenner, kiwrious, loading, nrrd_tools, setHDRFilePath };
+export { CameraViewPoint, Copper3dTrackballControls, MeshNodeTool, REVISION, addBoxHelper, addLabelToScene, configKiwriousHeart, convert3DPostoScreenPos, convertScreenPosto3DPos, copperMScene, copperMSceneRenderer, copperRenderer, copperRendererOnDemond, copperScene, copperSceneOnDemond, createTexture2D_NRRD, fullScreenListenner, kiwrious, loading, nrrd_tools, setHDRFilePath };