Changed: Let IntegrandBase::setMode() allocate primsol for RECOVERY mode

Author: kmokstad

src/ASM/IntegrandBase.C

@@ -25,6 +25,18 @@
 #include <cstdio>
 
 
+/*!
+  The default implementation of this method allocates for one solution vector
+  in \ref primsol when the \a mode is set to SIM::RECOVERY or SIM::NORMS.
+  This is sufficient for most linear problems with a single right-hand side.
+*/
+
+void IntegrandBase::setMode (SIM::SolutionMode mode)
+{
+  primsol.resize((m_mode = mode) >= SIM::RECOVERY ? 1 : 0);
+}
+
+
 /*!
   Override this method if the integrand needs some patch-specific data
   to be initialized before performing the numerical integration.

src/ASM/IntegrandBase.h

@@ -57,7 +57,7 @@ class IntegrandBase : public Integrand
   // ===============================
 
   //! \brief Defines the solution mode before the element assembly is started.
-  virtual void setMode(SIM::SolutionMode mode) { m_mode = mode; }
+  virtual void setMode(SIM::SolutionMode mode);
   //! \brief Returns current solution mode.
   SIM::SolutionMode getMode() const override { return m_mode; }
   //! \brief Initializes an integration parameter for the integrand.