/home/hauberg/Dokumenter/Capture/humim-tracker-0.1/src/OpenTissue/OpenTissue/dynamics/fem/fem_init.h

Go to the documentation of this file.

#ifndef OPENTISSUE_DYNAMICS_FEM_FEM_INIT_H
#define OPENTISSUE_DYNAMICS_FEM_FEM_INIT_H
//
// OpenTissue Template Library
// - A generic toolbox for physics-based modeling and simulation.
// Copyright (C) 2008 Department of Computer Science, University of Copenhagen.
//
// OTTL is licensed under zlib: http://opensource.org/licenses/zlib-license.php
//
#include <OpenTissue/configuration.h>

#include <OpenTissue/dynamics/fem/fem_uniform_young.h>
#include <OpenTissue/dynamics/fem/fem_uniform_poisson.h>
#include <OpenTissue/dynamics/fem/fem_uniform_density.h>
#include <OpenTissue/dynamics/fem/fem_compute_mass.h>
#include <OpenTissue/dynamics/fem/fem_initialize_stiffness_elements.h>
#include <OpenTissue/dynamics/fem/fem_clear_stiffness_assembly.h>
#include <OpenTissue/dynamics/fem/fem_initialize_plastic.h>

namespace OpenTissue
{
  namespace fem
  {
    template < typename fem_mesh,typename real_type >
    inline void init(
      fem_mesh & mesh,
      real_type const & young,
      real_type const & poisson,
      real_type const & density,
      real_type const & c_yield,
      real_type const & c_creep,
      real_type const & c_max
      )
    {
      //--- Assign material parameters
      detail::uniform_young( mesh.tetrahedron_begin(), mesh.tetrahedron_end(), young );
      detail::uniform_poisson( mesh.tetrahedron_begin(), mesh.tetrahedron_end(), poisson );
      detail::uniform_density(mesh.tetrahedron_begin(),mesh.tetrahedron_end(),density);
      //--- Compute stiffness and mass matrices
      detail::initialize_stiffness_elements(mesh.tetrahedron_begin(),mesh.tetrahedron_end());
      detail::clear_stiffness_assembly(mesh.node_begin(),mesh.node_end());
      detail::compute_mass(mesh);
      detail::initialize_plastic(mesh.tetrahedron_begin(),mesh.tetrahedron_end(),c_yield,c_creep,c_max);
    }

  } // namespace fem
} // namespace OpenTissue

//OPENTISSUE_DYNAMICS_FEM_FEM_INIT_H
#endif