@sdeverywhere/runtime 0.2.8 → 0.2.9

package/README.md

@@ -150,7 +150,7 @@ The `emcc` command line options should be similar to the following:
 
 ```
 $ emcc \
-build/<mymodel>.c build/macros.c build/model.c build/vensim.c \
+build/<mymodel>.c build/macros.c build/model.c build/vensim.c build/allocation.c \
 -Ibuild -o ./output/<mymodel>.js -Wall -Os \
 -s STRICT=1 -s MALLOC=emmalloc -s FILESYSTEM=0 -s MODULARIZE=1 \
 -s EXPORTED_FUNCTIONS="['_malloc','_free','_getInitialTime','_getFinalTime','_getSaveper','_setLookup','_runModelWithBuffers']" \

package/dist/index.cjs

@@ -1245,6 +1245,8 @@ function getJsModelFunctions() {
   let ctx;
   const cachedVectors = /* @__PURE__ */ new Map();
   const cachedSortVectors = /* @__PURE__ */ new Map();
+  const cachedInvertWorkMatrices = /* @__PURE__ */ new Map();
+  const cachedInvertResultMatrices = /* @__PURE__ */ new Map();
   return {
     setContext(context) {
       ctx = context;
@@ -1279,6 +1281,74 @@ function getJsModelFunctions() {
     INTEGER(x) {
       return Math.trunc(x);
     },
+    INVERT_MATRIX(matrix, n) {
+      let work = cachedInvertWorkMatrices.get(n);
+      if (work === void 0) {
+        work = Array(n);
+        for (let i = 0; i < n; i++) {
+          work[i] = Array(2 * n);
+        }
+        cachedInvertWorkMatrices.set(n, work);
+      }
+      let result = cachedInvertResultMatrices.get(n);
+      if (result === void 0) {
+        result = Array(n);
+        for (let i = 0; i < n; i++) {
+          result[i] = Array(n);
+        }
+        cachedInvertResultMatrices.set(n, result);
+      }
+      for (let i = 0; i < n; i++) {
+        for (let j = 0; j < n; j++) {
+          work[i][j] = matrix[i][j];
+          work[i][n + j] = i === j ? 1 : 0;
+        }
+      }
+      for (let k = 0; k < n; k++) {
+        let pivot = k;
+        let max = Math.abs(work[k][k]);
+        for (let i = k + 1; i < n; i++) {
+          const v = Math.abs(work[i][k]);
+          if (v > max) {
+            max = v;
+            pivot = i;
+          }
+        }
+        if (max === 0) {
+          for (let i = 0; i < n; i++) {
+            for (let j = 0; j < n; j++) {
+              result[i][j] = _NA_;
+            }
+          }
+          return result;
+        }
+        if (pivot !== k) {
+          const tmpRow = work[k];
+          work[k] = work[pivot];
+          work[pivot] = tmpRow;
+        }
+        const p = work[k][k];
+        for (let j = 0; j < 2 * n; j++) {
+          work[k][j] /= p;
+        }
+        for (let i = 0; i < n; i++) {
+          if (i !== k) {
+            const f = work[i][k];
+            if (f !== 0) {
+              for (let j = 0; j < 2 * n; j++) {
+                work[i][j] -= f * work[k][j];
+              }
+            }
+          }
+        }
+      }
+      for (let i = 0; i < n; i++) {
+        for (let j = 0; j < n; j++) {
+          result[i][j] = work[i][n + j];
+        }
+      }
+      return result;
+    },
     LN(x) {
       return Math.log(x);
     },