prepare some vectorized function signatures (to be implemented later)

Author: philippwindischhofer

include/Eisvogel/WeightingField.hh

@@ -60,6 +60,22 @@ public:
 
   template <typename KernelT = DefaultKernel>
   scalar_t E_phi(CoordVector pos);
+
+  // // Vectorized accessors for field components
+  // template <typename KernelT = DefaultKernel>
+  // std::vector<scalar_t> E_r(std::vector<CoordVector>& pos);
+
+  // template <typename KernelT = DefaultKernel>
+  // std::vector<scalar_t> E_z(std::vector<CoordVector>& pos);
+
+  // template <typename KernelT = DefaultKernel>
+  // std::vector<scalar_t> E_phi(std::vector<CoordVector>& pos);
+  
+  // template <typename KernelT = DefaultKernel>
+  // scalar_t Eip(CoordVector& pos, CoordVector& vec);
+
+  // template <typename KernelT = DefaultKernel>
+  // std::vector<scalar_t> Eip(std::vector<CoordVector>& pos, std::vector<CoordVector>& vec);
 
   DeltaVector GetSamplingIntervals() const;
   CoordVector GetStartCoords() const;

src/Integrator.hh

@@ -10,6 +10,9 @@
 template <class WeightingFieldT, typename KernelT = KeysCubicInterpolationKernel>
 scalar_t integrate(WeightingFieldT& wf, scalar_t t, const Current0D& curr, scalar_t os_factor = 1.0);
 
+// template <class WeightingFieldT, typename KernelT = KeysCubicInterpolationKernel>
+// std::vector<scalar_t> integrate(WeightingFieldT& wf, std::vector<scalar_t>& ts, const Current0D& curr, scalar_t os_factor = 1.0);
+
 template <class WeightingFieldT, typename KernelT = KeysCubicInterpolationKernel>
 scalar_t integrate(WeightingFieldT& wf, scalar_t t, const SparseCurrentDensity3D& current_distribution);
 

src/Integrator.hxx

@@ -66,6 +66,13 @@ scalar_t integrate(WeightingFieldT& wf, scalar_t t, const Current0D& curr, scala
   return signal;
 }
 
+// template <class WeightingFieldT, typename KernelT>
+// std::vector<scalar_t> integrate(WeightingFieldT& wf, std::vector<scalar_t>& ts, const Current0D& curr, scalar_t os_factor = 1.0) {
+
+  
+  
+// }
+
 template <class WeightingFieldT, typename KernelT>
 scalar_t integrate(WeightingFieldT& wf, scalar_t t, const SparseCurrentDensity3D& current_distribution) {
 

src/WeightingField.hxx

@@ -160,6 +160,14 @@ scalar_t WeightingField<FieldIndexerT, FieldStorageT>::E_phi(CoordVector pos) {
   return eval<KernelT>(pos, [&](auto& arg){return m_field_storage -> E_phi(arg);});
 };
 
+// template <class FieldIndexerT, class FieldSTorageT>
+// template <typename KernelT>
+// std::vector<scalar_t> WeightingField<FieldIndexerT, FieldStorageT>::E_r(std::vector<CoordVector>& pos) {
+//   std::vector<scalar_t> retval;
+  
+//   return retval;
+// }
+
 template <class FieldIndexerT, class FieldStorageT>
 template <typename KernelT, typename GetterT, typename ValueT>
 ValueT WeightingField<FieldIndexerT, FieldStorageT>::eval(CoordVector pos, GetterT getter) {