elements.cpp

#include <iostream>
#include "core/element.hpp"
#include "library/lib_shape.hpp"
 
void volume()
{
    typedef NiHu::volume_element<NiHu::quad_1_shape_set, double> Element;
 
    Element::coords_t coords;
    coords <<
        0.0, 1.0, 1.0, 0.0,
        0.0, 0.0, 1.0, 1.0;
 
    Element elem(coords);
 
    Element::xi_t xi;
    xi << 0.0, 0.0;
 
    std::cout << "location:\n" << elem.get_x(xi) << std::endl;
 
    std::cout << "gradient:\n" << elem.get_dx(xi) << std::endl;
    std::cout << "jacobian:\n" << elem.get_dx(xi).determinant() << std::endl;
 
    std::cout << "2nd derivative:\n" << elem.get_ddx(xi) << std::endl;
}
 
 
void surface()
{
    typedef NiHu::surface_element<NiHu::quad_1_shape_set, double> Element;
 
    Element::coords_t coords;
    coords <<
        0.0, 1.0, 1.0, 0.0,
        0.0, 0.0, 1.0, 1.0,
        0.0, 0.0, 0.0, 0.0;
 
    Element elem(coords);
 
    Element::xi_t xi;
    xi << 0.0, 0.0;
 
    std::cout << "location:\n" << elem.get_x(xi) << std::endl;
 
    std::cout << "gradient:\n" << elem.get_dx(xi) << std::endl;
    std::cout << "normal:\n" << elem.get_normal(xi) << std::endl;
    std::cout << "jacobian:\n" << elem.get_normal(xi).norm() << std::endl;
 
    std::cout << "2nd derivative:\n" << elem.get_ddx(xi) << std::endl;
}
 
int main()
{
    volume();
    surface();
    return 0;
}