From 15abe1962a1d3433ec319f0767893ea2815302b9 Mon Sep 17 00:00:00 2001
From: Jean-Paul Pelteret <jppelteret@gmail.com>
Date: Sat, 23 May 2020 18:19:53 +0200
Subject: [PATCH] Update to deal.II 9.1:
 MultipointFluxMixedFiniteElementMethods

---
 MultipointFluxMixedFiniteElementMethods/CMakeLists.txt | 2 +-
 MultipointFluxMixedFiniteElementMethods/mfmfe.cc       | 5 +++--
 2 files changed, 4 insertions(+), 3 deletions(-)

diff --git a/MultipointFluxMixedFiniteElementMethods/CMakeLists.txt b/MultipointFluxMixedFiniteElementMethods/CMakeLists.txt
index 36be816..055911f 100644
--- a/MultipointFluxMixedFiniteElementMethods/CMakeLists.txt
+++ b/MultipointFluxMixedFiniteElementMethods/CMakeLists.txt
@@ -24,7 +24,7 @@ SET(TARGET_SRC
 
 CMAKE_MINIMUM_REQUIRED(VERSION 2.8.12)
 
-FIND_PACKAGE(deal.II 9.0.0 QUIET
+FIND_PACKAGE(deal.II 9.1 QUIET
   HINTS ${deal.II_DIR} ${DEAL_II_DIR} ../ ../../ $ENV{DEAL_II_DIR}
   )
 IF(NOT ${deal.II_FOUND})
diff --git a/MultipointFluxMixedFiniteElementMethods/mfmfe.cc b/MultipointFluxMixedFiniteElementMethods/mfmfe.cc
index 9fdca37..a551288 100644
--- a/MultipointFluxMixedFiniteElementMethods/mfmfe.cc
+++ b/MultipointFluxMixedFiniteElementMethods/mfmfe.cc
@@ -17,6 +17,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/base/logstream.h>
 #include <deal.II/base/timer.h>
+#include <deal.II/base/utilities.h>
 #include <deal.II/base/work_stream.h>
 
 #include <deal.II/lac/full_matrix.h>
@@ -417,7 +418,7 @@ namespace MFMFE
     const FEValuesExtractors::Vector velocity (0);
     const FEValuesExtractors::Scalar pressure (dim);
 
-    const unsigned int n_vel = dim*pow(degree+1,dim);
+    const unsigned int n_vel = dim*Utilities::pow(degree+1,dim);
     std::unordered_map<unsigned int, std::unordered_map<unsigned int, double>> div_map;
 
     // One, we need to be able to assemble the communication between velocity and
@@ -818,7 +819,7 @@ namespace MFMFE
 
     pres_solution.reinit(n_p);
 
-    SolverControl solver_control (2.0*n_p, 1e-10);
+    SolverControl solver_control (static_cast<int>(2.0*n_p), 1e-10);
     SolverCG<> solver (solver_control);
 
     PreconditionIdentity identity;
-- 
2.39.5