fix domeigestimator callback

Author: balos1

bindings/pysundials/sundials/pysundials_domeigestimator.cpp

@@ -55,10 +55,10 @@ void bind_sundomeigestimator(nb::module_& m)
 
       if (ATimes)
       {
-        return SUNDomEigEstimator_SetATimes(dee, nullptr,
+        return SUNDomEigEstimator_SetATimes(dee, fntable,
                                             sundomeigestimator_atimes_wrapper);
       }
-      else { return SUNDomEigEstimator_SetATimes(dee, nullptr, nullptr); }
+      else { return SUNDomEigEstimator_SetATimes(dee, fntable, nullptr); }
     },
     nb::arg("DEE"), nb::arg("ATimes").none());
 }

bindings/pysundials/sundials/pysundials_domeigestimator_usersupplied.hpp

@@ -46,11 +46,12 @@ inline SUNDomEigEstimatorFunctionTable* SUNDomEigEstimatorFunctionTable_Alloc()
   return fn_table;
 }
 
-SUNErrCode sundomeigestimator_atimes_wrapper(void* A_data, N_Vector v, N_Vector z)
+template<typename... Args>
+SUNErrCode sundomeigestimator_atimes_wrapper(Args... args)
 {
   return pysundials::user_supplied_fn_caller<
     std::remove_pointer_t<SUNATimesFn>, SUNDomEigEstimatorFunctionTable,
-    3>(&SUNDomEigEstimatorFunctionTable::atimes, A_data, v, z);
+    3>(&SUNDomEigEstimatorFunctionTable::atimes, std::forward<Args>(args)...);
 }
 
 #endif