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

Go to the documentation of this file.

#ifndef OPENTISSUE_CORE_MATH_OPTIMIZATION_COMPUTE_INDEX_SETS_H
#define OPENTISSUE_CORE_MATH_OPTIMIZATION_COMPUTE_INDEX_SETS_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>
#include <OpenTissue/core/math/math_value_traits.h>
#include <cmath>
#include <cassert>

namespace OpenTissue
{
  namespace math
  {
    namespace optimization
    {
        template < typename T, typename bound_function_type>
        inline void compute_index_sets(
            ublas::vector<T> const & y
          , ublas::vector<T> const & x
          , bound_function_type const & l
          , bound_function_type const & u
          , ublas::vector<size_t> & bitmask
          , size_t  & cnt_active
          , size_t  & cnt_inactive
          )
        {
          using std::min;
          using std::max;

          typedef OpenTissue::math::ValueTraits<T>  value_traits;

          size_t n = x.size();
          bitmask.resize(n);

          cnt_active   = 0;
          cnt_inactive = 0;
          size_t cnt_lower = 0;
          size_t cnt_upper = 0;

          for (size_t i = 0; i < n; ++i)
          {
            if(y(i) > (x(i) - l(x,i)))
            {
              bitmask(i) = IN_LOWER;
              ++cnt_lower;
              ++cnt_inactive;
            }
            else if(y(i) < (x(i) - u(x,i)))
            {
              bitmask(i) = IN_UPPER;
              ++cnt_upper;
              ++cnt_inactive;
            }
            else
            {
              // y(i) <= (x(i) - lo(i)) &&  y(i) >= (x(i) - hi(i))
              bitmask(i) = IN_ACTIVE;
              ++cnt_active;
            }
          }
          assert( (cnt_active + cnt_inactive) == n        || !"compute_index_sets(): index sets were inconsistent");
          assert( (cnt_lower + cnt_upper) == cnt_inactive || !"compute_index_sets(): index sets were inconsistent");
        }

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

// OPENTISSUE_CORE_MATH_OPTIMIZATION_COMPUTE_INDEX_SETS_H
#endif