Merge lp:~als88/fenics-plasticity/plasticity_update_ghost into lp:fenics-plasticity

[Antony Siahaan]

Adding a method (Strain::update_ghost_displacement in IPData.cpp) which updates the ghost value of displacement in order for GenericVector::get_local to catch the the correct values.

=== modified file 'src/IPData.cpp'
--- src/IPData.cpp	2013-03-04 13:06:28 +0000
+++ src/IPData.cpp	2013-04-04 21:31:23 +0000
@@ -166,7 +166,7 @@
     dolfin::error("Only 2D and 3D is supported");
 }
 //-----------------------------------------------------------------------------
-StrainData::StrainData(const dolfin::Function& U, const dolfin::FunctionSpace& Q) :
+StrainData::StrainData(dolfin::Function& U, const dolfin::FunctionSpace& Q) :
   _u_comp_elem(*U[0].function_space()->element()),
   _u_dofmap(*U.function_space()->dofmap()),
   _q_dofmap(*Q[0]->dofmap()),
@@ -194,6 +194,11 @@
 {
   delete [] _derivatives;
 }
+
+void StrainData::update_ghost_disp()
+{
+  _vector.update_ghost_values();  
+}
 //-----------------------------------------------------------------------------
 void StrainData::update(const dolfin::Cell& cell)
 {
 
=== modified file 'src/IPData.h'
--- src/IPData.h	2013-03-04 13:06:28 +0000
+++ src/IPData.h	2013-04-04 21:31:23 +0000
@@ -142,11 +142,14 @@
     public:
 
       /// Constructor
-      StrainData(const dolfin::Function& U, const dolfin::FunctionSpace& Q);
+      StrainData(dolfin::Function& U, const dolfin::FunctionSpace& Q);
 
       /// Destructor
       ~StrainData();
 
+
+      void update_ghost_disp();
+
       // Update ip coordinates and expansion coefficients for displacement
       // function for given cell.
       void update(const dolfin::Cell& cell);
@@ -167,7 +170,7 @@
       const dolfin::GenericDofMap& _q_dofmap;
 
       // Vector of expansion coefficients
-      const dolfin::GenericVector& _vector;
+      dolfin::GenericVector& _vector;
 
       // Number of integration points
       unsigned int _num_ip;
 
=== modified file 'src/PlasticityProblem.cpp'
--- src/PlasticityProblem.cpp	2013-02-26 11:20:02 +0000
+++ src/PlasticityProblem.cpp	2013-04-04 21:31:23 +0000
@@ -76,6 +76,8 @@
   // Elastic tangent
   const arma::mat& elastic_tangent = plastic_model.elastic_tangent;
 
+  strain_data.update_ghost_disp();
+
   // Zero local ip data
   cons_tangent.zeros();
   trial_stress.zeros();