main.cpp

Dolfin provide a high-performance linear algebra library

// Copyright (C) 2006-2007 Anders Logg.
// Licensed under the GNU LGPL Version 2.1.
//
// Modified by Johan Hoffman 2006.
//
// First added:  2006-10-26
// Last changed: 2007-06-31

#include <dolfin.h>

using namespace dolfin;

int main()
{
  // Create mesh of unit square
  UnitSquare mesh(5, 5);
  plot(mesh);

  // Uniform refinement
  mesh.refine();
  plot(mesh);

  // Refine mesh close to x = (0.5, 0.5)
  Point p(0.5, 0.5);
  for (unsigned int i = 0; i < 5; i++)
  {
    // Mark cells for refinement
    MeshFunction<bool> cell_markers(mesh, mesh.topology().dim());
    cell_markers = false;
    for (CellIterator c(mesh); !c.end(); ++c)
    {
      if (c->midpoint().distance(p) < 0.1)
        cell_markers(*c) = true;
    }
    
    // Refine mesh
    mesh.refine(cell_markers);
    
    // Smooth mesh
    mesh.smooth();

    plot(mesh);
  }

  return 0;
}