Package: GrAL Cartesian2D       

---

Gral/Grids/Cartesian2D/test-cartesian-geometry2d.h

Go to the documentation of this file.

#ifndef GRAL_GB_CARTESIAN2D_TEST_CARTESIAN_GEOMETRY_2D_H
#define GRAL_GB_CARTESIAN2D_TEST_CARTESIAN_GEOMETRY_2D_H


/* ------------------------------------------------------------

    Copyright (C) 1997 - 2009 Guntram Berti

    This file is part of the Grid Algorithms Library (GrAL),
    available at http://gral.berti-cmm.de

    GrAL is distributed under the MIT license,
    see the file LICENSE or http://gral.berti-cmm.de/license

--------------------------------------------------------------- */


#include "Gral/Base/common-grid-basics.h"

namespace GrAL {

template<class GEOM>
void test_cartesian2d_geometry(GEOM & Geo, ::std::ostream& out)
{
  typedef grid_types<typename GEOM::grid_type> gt;

  for(typename gt::VertexIterator v(Geo.TheGrid()); !v.IsDone(); ++v)
    out << "v=" << (*v).index() << " @ " << Geo.coord(*v) << "\n";
  for(typename gt::EdgeIterator e(Geo.TheGrid()); ! e.IsDone(); ++e)
    out << "e=[" << (*e).V1().index() << " -> " << (*e).V2().index() << "] "
        << " length(e)=" << Geo.length(*e) 
        << " center(e)=" << Geo.center(*e)
        << " normal(e)=" << Geo.normal(*e)
        << "\n";
  for(typename gt::CellIterator c(Geo.TheGrid()); !c.IsDone(); ++c) {
    out << "c=" << (*c).index() 
        << " volume(c)=" << Geo.volume(*c) 
        << " center(c)=" << Geo.center(*c) 
        << "\n";
    for(typename gt::FacetOnCellIterator fc(*c); !fc.IsDone(); ++fc)
      out << "fc=[" << fc.V1().index() << " -> " << fc.V2().index() << "] "
          << " outer_area_normal(fc)=" << Geo.outer_area_normal(fc) 
          << " outer_normal(fc)="      << Geo.outer_normal(fc) 
          << "\n";
  }

}

} // namespace GrAL 

#endif

---

©  Guntram Berti 1997-2009. See the GrAL Homepage for up-to-date information.

---