implemented exponential function operator

Author: LukasMolzberger

CMakeLists.txt

@@ -20,6 +20,7 @@ set(SOURCES
         src/fields/subtraction.cpp
         src/fields/multiplication.cpp
         src/fields/division.cpp
+        src/fields/exponential_function.cpp
         src/fields/field_link_definition.cpp
         src/fields/obj.cpp
         src/fields/type.cpp

include/fields/exponential_function.h

@@ -0,0 +1,19 @@
+#ifndef EXPONENTIAL_FUNCTION_H
+#define EXPONENTIAL_FUNCTION_H
+
+#include "fields/abstract_function_definition.h"
+#include "fields/type.h"
+#include "fields/obj.h"
+
+class ExponentialFunction : public AbstractFunctionDefinition {
+public:
+    // Constructor
+    ExponentialFunction(Type* ref, const std::string& name);
+
+    void initializeField(Field* field) override;
+    
+    // Overridden method from AbstractFunctionDefinition
+    double computeUpdate(Obj* obj, FieldLinkDefinition* fl, double u) override;
+};
+
+#endif // EXPONENTIAL_FUNCTION_H 

install_manifest.txt

@@ -5,6 +5,7 @@
 /Users/lukasmolzberger/CLionProjects/aika-cpp/.venv/lib/python3.12/site-packages/aika/include/fields/null_terminated_array.h
 /Users/lukasmolzberger/CLionProjects/aika-cpp/.venv/lib/python3.12/site-packages/aika/include/fields/relation.h
 /Users/lukasmolzberger/CLionProjects/aika-cpp/.venv/lib/python3.12/site-packages/aika/include/fields/field_definition.h
+/Users/lukasmolzberger/CLionProjects/aika-cpp/.venv/lib/python3.12/site-packages/aika/include/fields/exponential_function.h
 /Users/lukasmolzberger/CLionProjects/aika-cpp/.venv/lib/python3.12/site-packages/aika/include/fields/flattened_type_relation.h
 /Users/lukasmolzberger/CLionProjects/aika-cpp/.venv/lib/python3.12/site-packages/aika/include/fields/test_object.h
 /Users/lukasmolzberger/CLionProjects/aika-cpp/.venv/lib/python3.12/site-packages/aika/include/fields/test_type.h

python-tests/exponential-test.py

@@ -0,0 +1,48 @@
+import unittest
+import sys
+import os
+
+# Add the project root to Python's module search path
+sys.path.append(os.path.abspath(os.path.join(os.path.dirname(__file__), "..")))
+
+import aika
+
+class ExponentialFunctionTestCase(unittest.TestCase):
+
+    def testExponentialFunction(self):
+        print("Module 'aika' was loaded from:", aika.__file__)
+
+        TEST_RELATION_FROM = aika.RelationOne(1, "TEST_FROM")
+        TEST_RELATION_TO = aika.RelationOne(2, "TEST_TO")
+        TEST_RELATION_TO.setReversed(TEST_RELATION_FROM)
+        TEST_RELATION_FROM.setReversed(TEST_RELATION_TO)
+
+        registry = aika.TypeRegistry()
+
+        typeA = aika.TestType(registry, "A")
+        typeB = aika.TestType(registry, "B")
+
+        a = typeA.inputField("a")
+
+        exp_func = typeB.exp("exp_func")
+
+        exp_func.input(TEST_RELATION_FROM, a, 0)
+
+        registry.flattenTypeHierarchy()
+
+        oa = typeA.instantiate()
+        ob = typeB.instantiate()
+
+        aika.TestObj.linkObjects(oa, ob)
+        ob.initFields()
+
+        oa.setFieldValue(a, 5.0)  
+
+        expected_value = 148.413159102576603
+        actual_value = ob.getFieldValue(exp_func)
+        print("ob.getFieldValue(exp_func):", actual_value)
+
+        self.assertAlmostEqual(expected_value, actual_value, places=5)
+
+if __name__ == '__main__':
+    unittest.main() 

src/fields/division.cpp

@@ -24,15 +24,15 @@ void Division::initializeField(Field* field) {
 // based on which argument is being updated.
 double Division::computeUpdate(Obj* obj, FieldLinkDefinition* fl, double u) {
     if (fl->getArgument() == 0) {
-        double divisor = getInputs()[1]->getInputValue(obj);
+        double divisor = inputs[1]->getInputValue(obj);
         if (divisor == 0.0) {
             // Handle division by zero if necessary
             return 0.0;
         }
         return u / divisor;
     } else {
-        double dividend = getInputs()[0]->getInputValue(obj);
-        double divisor = getInputs()[1]->getUpdatedInputValue(obj);
+        double dividend = inputs[0]->getInputValue(obj);
+        double divisor = inputs[1]->getUpdatedInputValue(obj);
         if (divisor == 0.0) {
             // Handle division by zero if necessary
             return 0.0;

src/fields/exponential_function.cpp

@@ -0,0 +1,19 @@
+#include "fields/exponential_function.h"
+#include <cmath>
+
+// Constructor
+ExponentialFunction::ExponentialFunction(Type* ref, const std::string& name)
+    : AbstractFunctionDefinition(ref, name, 1, 0.0) {}
+
+void ExponentialFunction::initializeField(Field* field) {
+    Obj* toObj = field->getObject();
+    double valueArg0 = inputs[0]->getInputValue(toObj);
+
+    field->setValue(std::exp(valueArg0));
+}
+
+// Overridden computeUpdate method
+double ExponentialFunction::computeUpdate(Obj* obj, FieldLinkDefinition* fl, double u) {
+    double valueArg0 = inputs[0]->getUpdatedInputValue(obj);
+    return std::exp(valueArg0) - obj->getFieldValue(this);
+} 

src/fields/multiplication.cpp

@@ -6,8 +6,6 @@ Multiplication::Multiplication(Type* ref, const std::string& name)
 
 void Multiplication::initializeField(Field* field) {
     Obj* toObj = field->getObject();
-    FieldDefinition* fieldDef = field->getFieldDefinition();
-    std::vector<FieldLinkDefinition*> inputs = fieldDef->getInputs();
     double valueArg0 = inputs[0]->getInputValue(toObj);
     double valueArg1 = inputs[1]->getInputValue(toObj);
 

src/fields/python_bindings.cpp

@@ -14,6 +14,7 @@
 #include "fields/addition.h"
 #include "fields/multiplication.h"
 #include "fields/division.h"
+#include "fields/exponential_function.h"
 
 
 // ----------------
@@ -61,6 +62,8 @@ PYBIND11_MODULE(aika, m)
       // Bind Division (inherits from AbstractFunctionDefinition)
       py::class_<Division, AbstractFunctionDefinition>(m, "Division");
 
+      py::class_<ExponentialFunction, AbstractFunctionDefinition>(m, "ExponentialFunction");
+
       py::class_<InputField, FieldDefinition>(m, "InputField")
             .def(py::init<Type*, const std::string &>())
             .def("__str__", [](const InputField &f) {
@@ -101,6 +104,12 @@ PYBIND11_MODULE(aika, m)
                         const_cast<Type*>(&ref),
                         name
                   );
+            }, py::return_value_policy::reference_internal)
+            .def("exp", [](const Type &ref, const std::string &name) {
+                  return new ExponentialFunction(
+                        const_cast<Type*>(&ref),
+                        name
+                  );
             }, py::return_value_policy::reference_internal);
 
       py::class_<Obj>(m, "Obj")