utest-mesh-lagrangep2b-triag2d.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 "Test module for ETYPE"

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

#include "common/Log.hpp"

#include "mesh/LagrangeP2B/Triag2D.hpp"

#include "Tools/Testing/Difference.hpp"

using namespace boost::assign;
using namespace cf3;
using namespace cf3::mesh;
using namespace cf3::mesh::LagrangeP2B;


typedef Triag2D ETYPE;

struct LagrangeP2BTriag2DFixture
{
  typedef ETYPE::NodesT NodesT;
  LagrangeP2BTriag2DFixture() :
    mapped_coords(1./6., 1./6.),
    nodes((NodesT() <<
           1.0, 0.0,
           0.0, 2.0,
           -1.0, 0.0,
           0.5, 1.0,
           -0.5, 1.0,
           0.0, 0.0,
           0.0, 2./3.).finished())
  {}

  ~LagrangeP2BTriag2DFixture()  {}

  const ETYPE::MappedCoordsT mapped_coords;
  const NodesT nodes;

};


BOOST_FIXTURE_TEST_SUITE( LagrangeP2BTriag2DSuite, LagrangeP2BTriag2DFixture )


BOOST_AUTO_TEST_CASE( ShapeFunction )
{
  ETYPE::SF::ValueT reference_result;
  reference_result << 5./18., -1./18., -1./18., 2./9., -1./9., 2./9., 0.5;
  ETYPE::SF::ValueT result;
  ETYPE::SF::compute_value(mapped_coords, result);

  std::cout<< reference_result << std::endl;
  std::cout<< reference_result.sum() << std::endl;

  std::cout<< result << std::endl;
  std::cout<< result.sum() << std::endl;


  BOOST_CHECK_CLOSE ( result.sum(), 1.0, 1e-8    );

}

BOOST_AUTO_TEST_CASE( Gradient )
{
//  const cf3::RealVector reference_result = list_of(5./18.)(-1./18.)(-1./18.)(2./9.)(-1./9.)(2./9.)(0.5);
  ETYPE::SF::GradientT result;
  ETYPE::SF::compute_gradient(mapped_coords, result);

//  std::cout<< reference_result << std::endl;
//  std::cout<< reference_result.sum() << std::endl;

  std::cout<< result << std::endl;

  std::cout<< result(XX,0) + result(XX,1) + result(XX,2) + result(XX,3) + result(XX,4) + result(XX,5) + result(XX,6)  << std::endl;
  std::cout<< result(YY,0) + result(YY,1) + result(YY,2) + result(YY,3) + result(YY,4) + result(YY,5) + result(YY,6)  << std::endl;


  const ETYPE::MappedCoordsT m1(1./6.,  5./6.);
  ETYPE::SF::compute_gradient(m1, result);

  std::cout<< result << std::endl;

  std::cout<< result(XX,0) + result(XX,1) + result(XX,2) + result(XX,3) + result(XX,4) + result(XX,5) + result(XX,6)  << std::endl;
  std::cout<< result(YY,0) + result(YY,1) + result(YY,2) + result(YY,3) + result(YY,4) + result(YY,5) + result(YY,6)  << std::endl;

  BOOST_CHECK_LE ( std::abs(result.sum()), 1e-14 );
}


BOOST_AUTO_TEST_SUITE_END()