utest-mesh-face-cell-connectivity.cpp

Go to the documentation of this file.

// Copyright (C) 2010-2013 von Karman Institute for Fluid Dynamics, Belgium
//
// This software is distributed under the terms of the
// GNU Lesser General Public License version 3 (LGPLv3).
// See doc/lgpl.txt and doc/gpl.txt for the license text.

#define BOOST_TEST_DYN_LINK
#define BOOST_TEST_MODULE "Tests cf3::mesh::FaceCellConnectivity"

#include <boost/test/unit_test.hpp>
#include <boost/assign/list_of.hpp>

#include "common/Log.hpp"
#include "common/OptionList.hpp"
#include "common/Core.hpp"
#include "common/FindComponents.hpp"

#include "Tools/Testing/TimedTestFixture.hpp"

#include "mesh/Mesh.hpp"
#include "mesh/Elements.hpp"
#include "mesh/Dictionary.hpp"
#include "mesh/Region.hpp"
#include "mesh/MeshReader.hpp"
#include "mesh/SimpleMeshGenerator.hpp"
#include "mesh/FaceCellConnectivity.hpp"
#include "mesh/ConnectivityData.hpp"

using namespace boost;
using namespace boost::assign;
using namespace cf3;
using namespace cf3::mesh;
using namespace cf3::common;
using namespace cf3::Tools;


struct FaceCellConnectivity_Fixture //: public Testing::TimedTestFixture
{
  FaceCellConnectivity_Fixture()
  {
     // uncomment if you want to use arguments to the test executable
     //int*    argc = &boost::unit_test::framework::master_test_suite().argc;
     //char*** argv = &boost::unit_test::framework::master_test_suite().argv;

  }

  ~FaceCellConnectivity_Fixture()
  {
  }



  static Handle< Mesh > m_mesh;
};

Handle< Mesh > FaceCellConnectivity_Fixture::m_mesh;

BOOST_FIXTURE_TEST_SUITE( FaceCellConnectivity_TestSuite, FaceCellConnectivity_Fixture )


BOOST_AUTO_TEST_CASE( Constructors )
{
  boost::shared_ptr<FaceCellConnectivity> c = allocate_component<FaceCellConnectivity>("faces_to_cells");
  BOOST_CHECK_EQUAL(c->name(),"faces_to_cells");
  BOOST_CHECK_EQUAL(FaceCellConnectivity::type_name(), "FaceCellConnectivity");
}


BOOST_AUTO_TEST_CASE( create_mesh )
{

 // create meshreader
 // boost::shared_ptr< MeshReader > meshreader = build_component_abstract_type<MeshReader>("cf3.mesh.neu.Reader","meshreader");
 // boost::filesystem::path fp_source ("quadtriag.neu");
 // m_mesh = meshreader->create_mesh_from(fp_source);

  m_mesh = Core::instance().root().create_component<Mesh>("mesh");
  Uint scale = 2;
  std::vector<Real> lengths  = list_of(4.)(2.);
  std::vector<Uint> nb_cells = list_of(scale*2u)(scale*2u);
  SimpleMeshGenerator& mesh_gen = *Core::instance().root().create_component<SimpleMeshGenerator>("mesh_gen");
  mesh_gen.options().set("mesh",m_mesh->uri());
  mesh_gen.options().set("lengths",lengths);
  mesh_gen.options().set("nb_cells",nb_cells);
  mesh_gen.execute();
  BOOST_CHECK(true);
}

BOOST_AUTO_TEST_CASE( face_elem_connectivity )
{

  // create and setup node to elements connectivity
  Handle<FaceCellConnectivity> c = m_mesh->create_component<FaceCellConnectivity>("face_cell_connectivity");
  c->setup( find_component<Region>(*m_mesh) );

  BOOST_CHECK_EQUAL(c->connectivity().size() , 40u);

  CFinfo << "nodes of face 0 : ";
  boost_foreach(const Uint node, c->face_nodes(0) )
    CFinfo << "  " << node;
  CFinfo << CFendl;
  // Output whole node to elements connectivity
  //CFinfo << c->connectivity() << CFendl;
}


BOOST_AUTO_TEST_SUITE_END()