Instead of loading the volume load as a fixed (x,y,z) vector, load it as a...

Instead of loading the volume load as a fixed (x,y,z) vector, load it as a space-dependent Python function

src/cosserat-continuum.cc

@@ -403,15 +403,14 @@ int main (int argc, char *argv[]) try
       return nV;
     };
 
-    FieldVector<double,3> volumeLoadValues {0,0,0};
-    if (parameterSet.hasKey("volumeLoad"))
-        volumeLoadValues = parameterSet.get<FieldVector<double,3> >("volumeLoad");
-
-    auto volumeLoad = [&]( FieldVector<double,dimworld>) {
-      auto vL = volumeLoadValues;
-      vL *= (-homotopyParameter);
-      return vL;
-    };
+    Python::Reference volumeLoadClass = Python::import(parameterSet.get<std::string>("volumeLoadPythonFunction", "zero-volume-load"));
+    Python::Callable volumeLoadCallable = volumeLoadClass.get("VolumeLoad");
+
+    // Call a constructor
+    Python::Reference volumeLoadPythonObject = volumeLoadCallable(homotopyParameter);
+
+    // Extract object member functions as Dune functions
+    auto volumeLoad = Python::make_function<FieldVector<double,3> > (volumeLoadPythonObject.get("volumeLoad"));
 
     if (mpiHelper.rank() == 0) {
         std::cout << "Material parameters:" << std::endl;