utilities.hpp

// This file is part of the imaging2 class library.
//
// University of Innsbruck, Infmath Imaging, 2009.
// http://infmath.uibk.ac.at
//
// All rights reserved.


#ifndef FEM_UTILITIES_H
#define FEM_UTILITIES_H

#include <fem/Grid.hpp>
#include <fem/fem_2d_triangle_types.hpp>
#include <fem/fem_1d_types.hpp>

namespace imaging
{
  void uniform_grid(float_t lower_bound, float_t upper_bound, std::size_t n_elements, Grid<fem_1d_types> & grid, ublas::compressed_matrix<float_t> & stiffness_matrix_prototype, std::size_t system_size = 1);
  
  void circle_grid(float_t radius, std::size_t n_rings, Grid<fem_2d_triangle_types> & grid, ublas::compressed_matrix<float_t> & stiffness_matrix_prototype, std::size_t system_size = 1);
  
  void ellipse_grid(float_t a, float_t b, std::size_t n_rings, Grid<fem_2d_triangle_types> & grid, ublas::compressed_matrix<float_t> & stiffness_matrix_prototype, std::size_t system_size = 1);
  
  void triangulate_shape(const BoundaryDiscretizer<2> & shape_discretizer, float_t max_triangle_area, Grid<fem_2d_triangle_types> & grid, ublas::compressed_matrix<float_t> & stiffness_matrix_prototype, std::size_t system_size = 1);
}

#endif

---