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

Go to the documentation of this file.

#ifndef OPENTISSUE_CORE_MATH_OPTIMIZATION_COMPUTE_INDEX_REORDERING_H
#define OPENTISSUE_CORE_MATH_OPTIMIZATION_COMPUTE_INDEX_REORDERING_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 <cassert>

namespace OpenTissue
{
  namespace math
  {
    namespace optimization
    {

      inline void  compute_index_reordering( 
        ublas::vector<size_t> const & bitmask
        , ublas::vector<size_t>    & old2new
        , ublas::vector<size_t>    & new2old
        )
      {
        size_t n = bitmask.size();
        old2new.resize(n);
        new2old.resize(n);

        size_t r = 0;

        for (size_t i = 0; i < n; ++ i)
        {
          if(bitmask(i) == IN_ACTIVE)
          {
            old2new( i ) = r;
            new2old( r ) = i;
            ++r;
          }
        }

        for (size_t i = 0; i < n; ++ i)
        {
          if(bitmask(i) == IN_LOWER)
          {
            old2new( i ) = r;
            new2old( r ) = i;
            ++r;
          }
        }

        for (size_t i = 0; i < n; ++ i)
        {
          if(bitmask(i) == IN_UPPER)
          {
            old2new( i ) = r;
            new2old( r ) = i;
            ++r;
          }
        }

        assert( r==n || !"compute_index_reordering(): something went wrong");
      }


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

// OPENTISSUE_CORE_MATH_OPTIMIZATION_COMPUTE_INDEX_REORDERING_H
#endif