From f18fed1f711c69d44930e5a900fb3f1d3e916882 Mon Sep 17 00:00:00 2001
From: Denis Davydov <davydden@gmail.com>
Date: Sat, 27 Feb 2016 20:07:22 +0100
Subject: [PATCH] fix apply_boundary_values() for PETSc

---
 source/numerics/matrix_tools.cc | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/source/numerics/matrix_tools.cc b/source/numerics/matrix_tools.cc
index 07a54505c5..28f784333d 100644
--- a/source/numerics/matrix_tools.cc
+++ b/source/numerics/matrix_tools.cc
@@ -2399,7 +2399,7 @@ namespace MatrixTools
             // figure out which rows of the matrix we
             // have to eliminate on this processor
             std::vector<types::global_dof_index> constrained_rows;
-            for (std::map<types::global_dof_index,double>::const_iterator
+            for (std::map<types::global_dof_index,PetscScalar>::const_iterator
                  dof  = boundary_values.begin();
                  dof != boundary_values.end();
                  ++dof)
@@ -2422,7 +2422,7 @@ namespace MatrixTools
 
             std::vector<types::global_dof_index> indices;
             std::vector<PetscScalar>  solution_values;
-            for (std::map<types::global_dof_index,double>::const_iterator
+            for (std::map<types::global_dof_index,PetscScalar>::const_iterator
                  dof  = boundary_values.begin();
                  dof != boundary_values.end();
                  ++dof)
@@ -2522,7 +2522,7 @@ namespace MatrixTools
               block++;
             }
           const types::global_dof_index index = dof->first - offset;
-          block_boundary_values[block].insert(std::pair<types::global_dof_index, double> (index,dof->second));
+          block_boundary_values[block].insert(std::pair<types::global_dof_index, PetscScalar> (index,dof->second));
         }
     }
 
-- 
2.39.5