← Back to branch summary

~njansson/dolfin/hpc

« back to all changes in this revision

Viewing changes to dolfin/fem/Assembler.cpp

  • Committer: Anders Logg
  • Date: 2008-05-28 08:13:15 UTC
  • mto: (2668.1.38 trunk)
  • mto: This revision was merged to the branch mainline in revision 2670.
  • Revision ID: logg@simula.no-20080528081315-fgsvgvhtxt8h54rv
Remove vertex_map and cell_map from BoundaryMesh interface.
They are now instead available as mesh data named "vertex map" and "cell map".

Show diffs side-by-side

added added

removed removed

Lines of Context:
dolfin/fem/Assembler.cpp
5
5
// Modified by Ola Skavhaug, 2007
6
6
//
7
7
// First added:  2007-01-17
8
 
// Last changed: 2008-04-10
 
8
// Last changed: 2008-05-28
9
9
 
10
10
#include <dolfin/log/dolfin_log.h>
11
11
#include <dolfin/common/Array.h>
15
15
#include <dolfin/mesh/Mesh.h>
16
16
#include <dolfin/mesh/Cell.h>
17
17
#include <dolfin/mesh/Facet.h>
 
18
#include <dolfin/mesh/MeshData.h>
18
19
#include <dolfin/mesh/BoundaryMesh.h>
19
20
#include <dolfin/mesh/MeshFunction.h>
20
21
#include <dolfin/mesh/SubDomain.h>
215
216
  ufc::exterior_facet_integral* integral = ufc.exterior_facet_integrals[0];
216
217
 
217
218
  // Create boundary mesh
218
 
  MeshFunction<uint> vertex_map;
219
 
  MeshFunction<uint> cell_map;
220
 
  BoundaryMesh boundary(mesh, vertex_map, cell_map);
221
 
  
 
219
  BoundaryMesh boundary(mesh);
 
220
  MeshFunction<uint>* cell_map = boundary.data().meshFunction("cell map");
 
221
  dolfin_assert(cell_map);
 
222
 
222
223
  // Assemble over exterior facets (the cells of the boundary)
223
224
  Progress p(progressMessage(A.rank(), "exterior facets"), boundary.numCells());
224
225
  for (CellIterator boundary_cell(boundary); !boundary_cell.end(); ++boundary_cell)
225
226
  {
226
227
    // Get mesh facet corresponding to boundary cell
227
 
    Facet mesh_facet(mesh, cell_map(*boundary_cell));
 
228
    Facet mesh_facet(mesh, (*cell_map)(*boundary_cell));
228
229
 
229
230
    // Get integral for sub domain (if any)
230
231
    if (domains && domains->size() > 0)