/home/hauberg/Dokumenter/Capture/humim-tracker-0.1/src/OpenTissue/OpenTissue/dynamics/sph/solvers/sph_buoyancy.h

Go to the documentation of this file.

#ifndef OPENTISSUE_DYNAMICS_SPH_SOLVERS_SPH_BUOYANCY_H
#define OPENTISSUE_DYNAMICS_SPH_SOLVERS_SPH_BUOYANCY_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/sph/sph_solver.h>
#include <OpenTissue/dynamics/sph/sph_particle.h>

namespace OpenTissue
{
  namespace sph
  {

    template< typename Types >
    class Buoyancy : Solver< Types, typename Types::vector >
    {
    public:
      typedef Solver<Types, typename Types::vector>  base_type;
      typedef typename base_type::value  value;
      typedef typename Types::real_type  real_type;
      typedef typename Types::vector  vector;
      typedef typename Types::particle  particle;
      typedef typename Types::particle_cptr_container::const_iterator  particle_cptr_container_citerator;

    public:
      Buoyancy(const real_type& buoyancy, const real_type& rest_density, const vector& gravity) : base_type()
        , m_b(buoyancy)
        , m_r0(rest_density)
        , m_g(gravity)
      {
      }

      ~Buoyancy()
      {
      }

      Buoyancy& operator=(const Buoyancy&)
      {
        return *this;
      }

    public:
      virtual value apply(const particle& par, particle_cptr_container_citerator, particle_cptr_container_citerator) const
      {
        return value(m_g*(m_b*(par.density()-m_r0)));
      }

      value apply(const particle& par, const particle&) const
      {
        return value(m_g*(m_b*(par.density()-m_r0)));
      }

    private:
      const real_type  m_b;    
      const real_type  m_r0;   
      const vector  m_g;  

    }; // End class Buoyancy

  } // namespace sph
} // namespace OpenTissue

// OPENTISSUE_DYNAMICS_SPH_SOLVERS_SPH_BUOYANCY_H
#endif