Add tensorProduct for DiscreteFunctionP0

The tensor product is performed cell by cell

Also give access to the tensor product in the language. In this special case vectors have to be of the same dimension since non-square matrices are not allowed there

Also add related tests

src/language/modules/MathFunctionRegisterForVh.cpp

@@ -212,6 +212,76 @@ MathFunctionRegisterForVh::MathFunctionRegisterForVh(SchemeModule& scheme_module
 
                                       ));
 
+  scheme_module._addBuiltinFunction("tensorProduct", std::function(
+
+                                                       [](std::shared_ptr<const DiscreteFunctionVariant> a,
+                                                          std::shared_ptr<const DiscreteFunctionVariant> b)
+                                                         -> std::shared_ptr<const DiscreteFunctionVariant> {
+                                                         return tensorProduct(a, b);
+                                                       }
+
+                                                       ));
+
+  scheme_module._addBuiltinFunction("tensorProduct",
+                                    std::function(
+
+                                      [](std::shared_ptr<const DiscreteFunctionVariant> a,
+                                         const TinyVector<1> b) -> std::shared_ptr<const DiscreteFunctionVariant> {
+                                        return tensorProduct(a, b);
+                                      }
+
+                                      ));
+
+  scheme_module._addBuiltinFunction("tensorProduct",
+                                    std::function(
+
+                                      [](std::shared_ptr<const DiscreteFunctionVariant> a,
+                                         const TinyVector<2> b) -> std::shared_ptr<const DiscreteFunctionVariant> {
+                                        return tensorProduct(a, b);
+                                      }
+
+                                      ));
+
+  scheme_module._addBuiltinFunction("tensorProduct",
+                                    std::function(
+
+                                      [](std::shared_ptr<const DiscreteFunctionVariant> a,
+                                         const TinyVector<3>& b) -> std::shared_ptr<const DiscreteFunctionVariant> {
+                                        return tensorProduct(a, b);
+                                      }
+
+                                      ));
+
+  scheme_module._addBuiltinFunction("tensorProduct",
+                                    std::function(
+
+                                      [](const TinyVector<1> a, std::shared_ptr<const DiscreteFunctionVariant> b)
+                                        -> std::shared_ptr<const DiscreteFunctionVariant> {
+                                        return tensorProduct(a, b);
+                                      }
+
+                                      ));
+
+  scheme_module._addBuiltinFunction("tensorProduct",
+                                    std::function(
+
+                                      [](const TinyVector<2> a, std::shared_ptr<const DiscreteFunctionVariant> b)
+                                        -> std::shared_ptr<const DiscreteFunctionVariant> {
+                                        return tensorProduct(a, b);
+                                      }
+
+                                      ));
+
+  scheme_module._addBuiltinFunction("tensorProduct",
+                                    std::function(
+
+                                      [](const TinyVector<3>& a, std::shared_ptr<const DiscreteFunctionVariant> b)
+                                        -> std::shared_ptr<const DiscreteFunctionVariant> {
+                                        return tensorProduct(a, b);
+                                      }
+
+                                      ));
+
   scheme_module._addBuiltinFunction("det", std::function(
 
                                              [](std::shared_ptr<const DiscreteFunctionVariant> A)