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

Go to the documentation of this file.

#ifndef OPENTISSUE_CORE_MATH_OPTIMIZATION_PARTITION_VECTOR_H
#define OPENTISSUE_CORE_MATH_OPTIMIZATION_PARTITION_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 <OpenTissue/core/math/optimization/optimization_constants.h>

namespace OpenTissue
{
  namespace math
  {
    namespace optimization
    {
      template<typename T >
      inline void  partition_vector( 
        ublas::vector<T> const & x
        , ublas::vector<size_t> const & bitmask
        , ublas::vector<size_t> const & old2new
        , size_t const & cnt_active
        , size_t const & cnt_inactive
        , ublas::vector<T>  & x_a
        , ublas::vector<T>  & x_b
        )
      {
        x_a.resize( cnt_active );
        x_b.resize( cnt_inactive );

        size_t n = x.size();

        for(size_t i_old = 0; i_old < n; ++i_old)
        {
          size_t i_new = old2new( i_old );
          if( bitmask(i_old) == IN_ACTIVE )
            x_a( i_new ) = x( i_old);
          else
            x_b( i_new-cnt_active ) = x( i_old);
        }
      }

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

// OPENTISSUE_CORE_MATH_OPTIMIZATION_PARTITION_VECTOR_H
#endif