Modified and Checked the code with v 9.5.1.

diff --git a/Crystal_Growth_Phase_Field_Model/CMakeLists.txt b/Crystal_Growth_Phase_Field_Model/CMakeLists.txt
index 5bef9c9706c6562a3527889785a20ab17e88014b..f54ca8772627e119cfb021ebdd42e917d8b1f24e 100644 (file)
--- a/Crystal_Growth_Phase_Field_Model/CMakeLists.txt
+++ b/Crystal_Growth_Phase_Field_Model/CMakeLists.txt
@@ -17,7 +17,7 @@ SET(TARGET_SRC
         get_random_number.cpp
         )
 
-FIND_PACKAGE(deal.II 9.2.0 QUIET
+FIND_PACKAGE(deal.II 9.5.0 QUIET
         HINTS ${deal.II_DIR} ${DEAL_II_DIR} ../ ../../ $ENV{DEAL_II_DIR}
         )
 IF(NOT ${deal.II_FOUND})

diff --git a/Crystal_Growth_Phase_Field_Model/applying_bc.cpp b/Crystal_Growth_Phase_Field_Model/applying_bc.cpp
index b111f3ddad4c77dc6c0a4c228819c69a5a85a47d..44b805182ab8ee1d85b2801079acf48cae445aef 100644 (file)
--- a/Crystal_Growth_Phase_Field_Model/applying_bc.cpp
+++ b/Crystal_Growth_Phase_Field_Model/applying_bc.cpp
@@ -17,7 +17,7 @@ void PhaseFieldSolver::applying_bc(){
     // Prescribing p=1 at the left face (this will be maintained in the subsequent iterations when zero BC is applied in the Newton-Raphson iterations)
     VectorTools::interpolate_boundary_values (dof_handler,
                                               1,
-                                              ConstantFunction<2>(1., 2),
+                                              Functions::ConstantFunction<2>(1., 2),
                                               boundary_values,p_mask);
 
     // To apply the boundary values only to the solution vector without the Jacobian Matrix and RHS Vector

diff --git a/Crystal_Growth_Phase_Field_Model/run.cpp b/Crystal_Growth_Phase_Field_Model/run.cpp
index c775ec09e9c42da8a40d826f37bc0ecc35415697..d942a986f31edac090d6b92b5b0db177587a1417 100644 (file)
--- a/Crystal_Growth_Phase_Field_Model/run.cpp
+++ b/Crystal_Growth_Phase_Field_Model/run.cpp
@@ -46,7 +46,7 @@ void PhaseFieldSolver::run() {
             old_solution_np = old_solution;
             //Initialise the delta solution as zero
             VectorTools::interpolate(dof_handler,
-                                     ZeroFunction<2>(2),
+                                     Functions::ZeroFunction<2>(2),
                                      solution_update);
             solution_update.compress(VectorOperation::insert);
             //Assemble Jacobian and Residual