Reduce wave::operator[] cost

common/wave/element.cpp

@@ -17,8 +17,12 @@ void WaveMultiField<T_Physic>::initializeField(const WaveField<T_Physic>& field)
 
 template<Physic T_Physic>
 const WaveField<T_Physic>& WaveMultiField<T_Physic>::operator[](const unsigned int index) const
+{
+    if(index!=_index)
     {
         _field.setAllVector(_dofValues[index], gmshfem::dofs::RawOrder::Hash);
+        _index=index;
+    }
     return _field;
 }
 

common/wave/element.h

@@ -19,8 +19,9 @@ private:
     mutable WaveField<T_Physic> _field;
     bool _isInitialized;
     std::vector< gmshfem::algebra::Vector< std::complex< double > > > _dofValues;
+    mutable unsigned int _index;
 public:
-    WaveMultiField(unsigned int n, const std::string& name = "") : _field(), _isInitialized(false), _dofValues(n) {_field.name(name);};
+    WaveMultiField(unsigned int n, const std::string& name = "") : _field(), _isInitialized(false), _dofValues(n), _index(n+1) {_field.name(name);};
 
     void initializeField(const WaveField<T_Physic>& field);
     bool isInitialized() const {return _isInitialized;};