/home/hauberg/Dokumenter/Capture/humim-tracker-0.1/src/OpenTissue/OpenTissue/core/math/optimization/optimization_agglomerate_vector.h

Go to the documentation of this file.

#ifndef OPENTISSUE_CORE_MATH_OPTIMIZATION_AGGLOMERATE_VECTOR_H
#define OPENTISSUE_CORE_MATH_OPTIMIZATION_AGGLOMERATE_VECTOR_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/big/big_types.h>
#include <stdexcept>

namespace OpenTissue
{
  namespace math
  {
    namespace optimization
    {
      template<typename T>
      inline void agglomerate_vector(
        ublas::vector<T> const & x_a
        , ublas::vector<T> const & x_b
        , ublas::vector<size_t> const & new2old
        , ublas::vector<T> & x
        )
      {
        size_t n   = x.size();
        size_t n_a = x_a.size();
        size_t n_b = x_b.size();

        if ( (n_a + n_b) != n )
          throw std::invalid_argument( "Incompatible dimensions" );

        for(size_t i = 0; i< n_a;++i)
          x( new2old( i ) ) = x_a( i );
        for(size_t i = 0; i< n_b;++i)
          x( new2old( i + n_a ) ) = x_b( i );
      }

    } // namespace optimization
  } // namespace math
} // namespace OpenTissue

// OPENTISSUE_CORE_MATH_OPTIMIZATION_AGGLOMERATE_VECTOR_H
#endif