[Feature request] Non-const access to PreBasis
I would like to have some possibility to either
- access the PreBasis of a basis in a non-const way, or
- have a more general method for updates than
update(gridView)
.
For the former there is const PreBasis& preBasis() const
in the DefaultGlobalBasis
but I can't update anything through that due to the const-ness.
For the latter something like
template<typename... Args> void update(Args&&... args) { preBasis_.update(std::forward<Args>(args)...); preBasis_.initializeIndices(); }
in the DefaultGlobalBasis
comes to mind. That way, if I implement a custom basis through a custom PreBasis, I could forward any update information to my PreBasis. I'm not so sure about the preBasis_.initializeIndices()
part there, though. On the one hand it's consistent with the existing update(gv)
behaviour, on the other hand maybe one does not always want to trigger the index initialization when using the update method.