unitinterval.cpp

来自「利用C」· C++ 代码 · 共 47 行

// Copyright (C) 2007 Kristian B. Oelgaard.
// Licensed under the GNU LGPL Version 2.1.
//
// First added:  2007-11-23
// Last changed: 2007-11-23

#include "MeshEditor.h"
#include "UnitInterval.h"

using namespace dolfin;

//-----------------------------------------------------------------------------
UnitInterval::UnitInterval(uint nx) : Mesh()
{
  if ( nx < 1 )
    error("Size of unit interval must be at least 1.");

  rename("mesh", "Mesh of the unit interval (0,1)");

  // Open mesh for editing
  MeshEditor editor;
  editor.open(*this, CellType::interval, 1, 1);

  // Create vertices and cells:
  editor.initVertices((nx+1));
  editor.initCells(nx);

  // Create main vertices:
  for (uint ix = 0; ix <= nx; ix++)
  {
    const real x = static_cast<real>(ix) / static_cast<real>(nx);
    editor.addVertex(ix, x);
  }

  // Create intervals
  for (uint ix = 0; ix < nx; ix++) {
  	const uint v0 = ix;
  	const uint v1 = v0 + 1;
  	editor.addCell(ix, v0, v1);
  }

  // Close mesh editor
  editor.close();

}
//-----------------------------------------------------------------------------