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

[Antony Siahaan]

MPI is now working in fenics-plasticity.

Change in IPData.cpp and IPData.h to enable parallel execution :
- using local_size() method instead of size() method for vector object.
- using get_local() method instead of gather() method for vector object.

=== modified file 'src/IPData.cpp'
--- src/IPData.cpp	2013-02-26 11:20:02 +0000
+++ src/IPData.cpp	2013-03-14 10:55:25 +0000
@@ -46,7 +46,7 @@
 void HistoryData::get_old_values(dolfin::GenericVector& vector)
 {
   // Check that size are equal becore setting values.
-  assert(vector.size() == _old_vals.size());
+  assert(vector.local_size() == _old_vals.size());
   vector.set(&_old_vals[0], _old_vals.size(), &_dofs[0]);
 }
 //-----------------------------------------------------------------------------
@@ -204,7 +204,8 @@
   _element_dofs = _u_dofmap.cell_dofs(cell.index());
 
   // Get displacements at nodes on element
-  _vector.gather(_element_vals, _element_dofs);
+  //_vector.gather(_element_vals, _element_dofs);
+  _vector.get_local(&(_element_vals[0]), _element_dofs.size(), &(_element_dofs[0]));
 
   // Tabulate coordinates of quadrature points on current cell
   _q_dofmap.tabulate_coordinates(_coordinates, _ufc_cell);
 
=== modified file 'src/IPData.h'
--- src/IPData.h	2012-12-30 12:20:58 +0000
+++ src/IPData.h	2013-03-14 10:55:25 +0000
@@ -89,13 +89,15 @@
     /// Update local armadillo data structure from associated function vector
     void update_arma(arma::mat& matrix) const
     {
-      _vector.gather(_ip_values, _ip_dofs);
+     // _vector.gather(_ip_values, _ip_dofs);
+      _vector.get_local(&(_ip_values[0]), _ip_dofs.size(), &(_ip_dofs[0]));
       local_to_arma(matrix);
     }
 
     void update_arma(arma::vec& vector) const
     {
-      _vector.gather(_ip_values, _ip_dofs);
+      //_vector.gather(_ip_values, _ip_dofs);
+      _vector.get_local(&(_ip_values[0]), _ip_dofs.size(), &(_ip_dofs[0]));
       local_to_arma(vector);
     }