added divide. now for the basic arithmetic operations dervis can be c…

…alculated
Konrad1991 · Jun 3, 2024 · 77db7cd · 77db7cd
1 parent 97dfe77
commit 77db7cd
Show file tree

Hide file tree

Showing 2 changed files with 21 additions and 0 deletions.
diff --git a/include/etr_bits/Vector/DerivativeCalc.hpp b/include/etr_bits/Vector/DerivativeCalc.hpp
@@ -60,6 +60,19 @@ static constexpr auto walkTD() {
                               TimesDerivTrait>();
 }
 
+template <typename TRaw>
+  requires IsDivide<TRaw>
+static constexpr auto walkTD() {
+  using TDDivision = std::remove_const_t<std::remove_reference_t<TRaw>>;
+  constexpr auto LT = produceVariableType<typename TDDivision::typeTraitL>();
+  constexpr auto RT = produceVariableType<typename TDDivision::typeTraitR>();
+  constexpr auto LDeriv = walkTD<typename TDDivision::typeTraitL>();
+  constexpr auto RDeriv = walkTD<typename TDDivision::typeTraitR>();
+  return produceQuarternyType<decltype(LT), decltype(RT), decltype(LDeriv),
+                              decltype(RDeriv), QuarternaryTrait,
+                              DivideDerivTrait>();
+}
+
 template <typename TD>
   requires IsAddition<TD>
 static constexpr auto walkTD() {

diff --git a/tests/Derivatives_Tests.cpp b/tests/Derivatives_Tests.cpp
@@ -66,5 +66,13 @@ int main() {
     print(b, av);
     print(get_derivs(a));
     print(get_derivs(b));
+
+    std::cout << "\n"
+              << "a = a / b" << std::endl;
+    a = a / b;
+    print(a, av);
+    print(b, av);
+    print(get_derivs(a));
+    print(get_derivs(b));
   }
 }