diff --git a/templates/semi_analytic_solutions.cpp.in b/templates/semi_analytic_solutions.cpp.in
index bb76219f1..d2e336e54 100644
--- a/templates/semi_analytic_solutions.cpp.in
+++ b/templates/semi_analytic_solutions.cpp.in
@@ -99,7 +99,7 @@ void @ModelName@_semi_analytic_solutions::set_to_boundary_values(
 }
 
 void @ModelName@_semi_analytic_solutions::evaluate_solutions(
-   @ModelName@_soft_parameters& model)
+   @ModelName@_soft_parameters& model) const
 {
 @evaluateSemiAnalyticSolns@
 }
diff --git a/templates/semi_analytic_solutions.hpp.in b/templates/semi_analytic_solutions.hpp.in
index 33b71083d..2290eb2bd 100644
--- a/templates/semi_analytic_solutions.hpp.in
+++ b/templates/semi_analytic_solutions.hpp.in
@@ -64,7 +64,7 @@ public:
 @boundaryValueSetters@
 @coefficientGetters@
    void calculate_coefficients(const @ModelName@_soft_parameters&);
-   void evaluate_solutions(@ModelName@_soft_parameters&);
+   void evaluate_solutions(@ModelName@_soft_parameters&) const;
 
 private:
    struct Boundary_values {