/home/hauberg/Dokumenter/Capture/humim-tracker-0.1/src/OpenTissue/OpenTissue/dynamics/mbd/math/mbd_math_compute_f.h

Go to the documentation of this file.

#ifndef OPENTISSUE_DYNAMICS_MBD_UTIL_MATH_COMPUTE_F_H
#define OPENTISSUE_DYNAMICS_MBD_UTIL_MATH_COMPUTE_F_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/core/math/math_is_number.h>

#include <OpenTissue/dynamics/mbd/math/mbd_math_update_f.h>

#include <cassert>

namespace OpenTissue
{
  namespace mbd
  {
    namespace math
    {

      template<typename matrix_type, typename vector_type>
      void compute_f(matrix_type const & WJT,vector_type const & x,vector_type & f)
      {
        typedef typename matrix_type::size_type  size_type;

        assert( WJT.size1() >0          || !"compute_f(): WJT was empty?");
        assert( WJT.size2() >0          || !"compute_f(): WJT was empty?");
        assert( WJT.size1() == x.size() || !"compute_f(): Incorrect dimensions");

        f.resize( WJT.size2(), false );
        f.clear();

        for (size_type i = 0;i < x.size() ; ++i)
          update_f(WJT,i,x(i),f);
      }


    } // end  namespace math
  } // end  namespace mbd
} // end namespace OpenTissue
// OPENTISSUE_DYNAMICS_MBD_UTIL_MATH_COMPUTE_F_H
#endif