@@ -4,6 +4,7 @@ members = [
   "nac3ast",
   "nac3parser",
   "nac3core",
+  "nac3core/src/codegen/externfns",

@@ -1874,6 +1876,37 @@ impl<'a> BuiltinBuilder<'a> {
         }
     }
 
+    /// Build the functions `try_invert_to` and `wilkinson_shift`

@@ -263,6 +265,8 @@ impl PrimDef {
             PrimDef::FunNpLdExp => fun("np_ldexp", None),
             PrimDef::FunNpHypot => fun("np_hypot", None),
             PrimDef::FunNpNextAfter => fun("np_nextafter", None),
+            PrimDef::FunTryInvertTo => fun("try_invert_to", None),

@@ -265,2 +267,4 @@
             PrimDef::FunNpNextAfter => fun("np_nextafter", None),
+            PrimDef::FunTryInvertTo => fun("try_invert_to", None),
+            PrimDef::FunWilkinsonShift => fun("wilkinson_shift", None),
             PrimDef::FunSome => fun("Some", None),

@@ -144,0 +145,4 @@
+        try:
+            y = np.linalg.inv(x)
+            x[:] = y
+        except np.linalg.LinAlgError:

@@ -144,0 +149,4 @@
+            return False
+        return True
+
+    def wilkinson_shift(x):

@@ -55,3 +58,3 @@
 
   clang -c -std=gnu11 -Wall -Wextra -O3 -o demo.o demo.c
-  clang -lm -o demo module.o demo.o
+  clang -lm -o demo module.o demo.o $externfns

@@ -0,0 +56,4 @@
+    let matrix2 = DMatrix::from_row_slice(dim2_0, dim2_1, data_slice2);
+    let mut result = DMatrix::<f64>::zeros(dim3_0, dim3_1);
+
+    matrix1.mul_to(&matrix2, &mut result);

@@ -1838,0 +1879,4 @@
+                            "",
+                        )
+                        .unwrap(),
+                    "0:ValueError",

@@ -1838,0 +1865,4 @@
+            // The following constraints must be satisfied:
+            // * Input must be 2D
+            // * number of rows should equal number of columns (square matrix)
+            if cfg!(debug_assertions) {

@@ -1838,0 +1883,4 @@
+
+        let (BasicTypeEnum::FloatType(_), BasicTypeEnum::FloatType(_)) = (n1_elem_ty, n2_elem_ty)
+        else {
+            unimplemented!("{FN_NAME} operates on float type NdArrays only");

@@ -1838,0 +2033,4 @@
+    }
+}
+
+/// Invokes the `np_linalg_svd` using `nalgebra` crate

@@ -0,0 +31,4 @@
+        raise_exn!($name, $fn_name, $message, 0, 0, 0)
+    }};
+}
+pub struct InputMatrix {

@@ -0,0 +8,4 @@
+use core::slice;
+use nalgebra::DMatrix;
+
+macro_rules! raise_exn {

@@ -0,0 +30,4 @@
+    }
+}
+
+pub unsafe fn raise(exception: *const Exception) -> ! {

@@ -45,2 +45,4 @@
+    linalg_externfns=../../target/debug/deps/liblinalg_externfns-?*.a
 elif [ -e ../../target/release/nac3standalone ]; then
     nac3standalone=../../target/release/nac3standalone
+    linalg_externfns=../../target/release/deps/liblinalg_externfns-?*.a

@@ -17,9 +17,7 @@ demo="$1"
 echo -n "Checking $demo... "
 ./interpret_demo.py "$demo" > interpreted.log
 ./run_demo.sh --out run.log "${nac3args[@]}" "$demo"
-./run_demo.sh --lli --out run_lli.log "${nac3args[@]}" "$demo"

@@ -63,0 +66,4 @@
+   
+   # Compile linalg crate to provide functions compatible with i386 architecture
+   if [ ! -f ../../target/i686-unknown-linux-gnu/release/liblinalg.a ]; then
+     cd linalg && nix-shell -p rustup --command "RUSTFLAGS=\"-C target-cpu=i386 -C target-feature=+sse2\" cargo build -q --release --target=i686-unknown-linux-gnu" && cd ..

@@ -1839,0 +1867,4 @@
+    out_ptr
+}
+
+/// Invokes the `np_linalg_matmul` linalg function

@@ -102,0 +114,4 @@
+    FunNpLinalgDet,
+    FunSpLinalgLu,
+    FunSpLinalgSchur,
+    FunSpLinalgHessenberg,

@@ -102,0 +104,4 @@
+
+    // Linalg functions
+    FunNpDot,
+    FunNpLinalgMatmul,

@@ -8,6 +8,7 @@ edition = "2021"
 parking_lot = "0.12"
 nac3parser = { path = "../nac3parser" }
 nac3core = { path = "../nac3core" }
+linalg = { path = "./demo/linalg" }

@@ -5,6 +5,7 @@ import importlib.util
 import importlib.machinery
 import math
 import numpy as np
+import scipy as sp

@@ -64,2 +64,2 @@
-   clang -m32 -c -std=gnu11 -Wall -Wextra -O3 -msse2 -o demo.o demo.c
-   clang -m32 -lm -Wl,--no-warn-search-mismatch -o demo module.o demo.o
+   
+  cd linalg && cargo build --release --target i686-unknown-linux-gnu -q && cd ..