dataset.h

Go to the documentation of this file.

/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
/*                                                                           */
/*                  This file is part of the class library                   */
/*       SoPlex --- the Sequential object-oriented simPlex.                  */
/*                                                                           */
/*  Copyright (c) 1996-2024 Zuse Institute Berlin (ZIB)                      */
/*                                                                           */
/*  Licensed under the Apache License, Version 2.0 (the "License");          */
/*  you may not use this file except in compliance with the License.         */
/*  You may obtain a copy of the License at                                  */
/*                                                                           */
/*      http://www.apache.org/licenses/LICENSE-2.0                           */
/*                                                                           */
/*  Unless required by applicable law or agreed to in writing, software      */
/*  distributed under the License is distributed on an "AS IS" BASIS,        */
/*  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. */
/*  See the License for the specific language governing permissions and      */
/*  limitations under the License.                                           */
/*                                                                           */
/*  You should have received a copy of the Apache-2.0 license                */
/*  along with SoPlex; see the file LICENSE. If not email to soplex@zib.de.  */
/*                                                                           */
/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
 
/**@file  dataset.h
 * @brief Set of data objects.
 */
#ifndef _DATASET_H_
#define _DATASET_H_
 
 
#include <string.h>
#include <assert.h>
 
#include "soplex/array.h"
#include "soplex/dataarray.h"
#include "soplex/datakey.h"
#include "soplex/spxalloc.h"
#include "soplex/exceptions.h"
 
namespace soplex
{
/**@class DataSet
   @brief   Set of data objects.
   @ingroup Elementary
 
   Class DataSet manages of sets of data objects  of a
   template type DATA. For constructing a DataSet the maximum number
   of entries must be given. The current maximum number may be inquired
   with method max().
 
   Adding more then max() elements to a DataSet will core dump. However,
   method reMax() allows to reset max() without loss of elements currently
   in the DataSet. The current number of elements in a DataSet is returned
   by method num().
 
   Adding elements to a DataSet is done via methods add() or create(),
   while remove() removes elements from a DataSet. When adding an element
   to a DataSet the new element is assigned a DataKey. DataKeys serve to
   access DATA elements in a set via a version of the subscript
   operator[](DataKey).
 
   For convenience all elements in a DataSet are implicitely numbered
   from 0 through num()-1 and can be accessed with these numbers
   using a 2nd subscript operator[](int). The reason for providing
   DataKeys to access elements of a DataSet is that the Key of an
   element remains unchanged as long as the element is a member of the
   DataSet, while the numbers will change in an undefined way, if
   other elements are added to or removed from the DataSet.
 
   The elements in a DataSet and their DataKeys are stored in two arrays:
   - theitem keeps the elements data along with their number stored in item.
   - thekey  keeps the DataKey::idx's of the elements in a DataSet.
 
   Both arrays have size themax.
 
   In #thekey only elements 0 thru thenum-1 contain DataKey::idx%'s of
   valid elements, i.e., elements currently in the DataSet.
   The current number of elements in the DataSet is counted in thenum.
 
   In #theitem only elements 0 thru thesize-1 are used, but only some of
   them actually contain real data elements of the DataSet. They are
   recognized by having info >= 0, which gives the number of that
   element. Otherwise info < 0 indicates an unused element. Unused
   elements are linked in a single linked list: starting with element
   <tt>-firstfree-1</tt>, the next free element is given by
   <tt>-info-1.</tt> The last free element in the list is marked by
   <tt>info == -themax-1.</tt> Finally all elements in theitem with
   <tt>index >= thesize</tt> are unused as well.
 
   @warning malloc/realloc and memcpy are used to handle the members
      of the set. If you use DataSet with something that is not
      a \ref DataObjects "Data Object" you will be in severe trouble.
*/
template<class DATA>
class DataSet
{
   static_assert(std::is_trivially_copyable<DATA>::value,
                 "Only trivially copyable types are allowed with DataSet, since it does memcopy");
protected:
 
   //-----------------------------------
   /**@name Types */
   ///@{
   ///
   struct Item
   {
      DATA data;       ///< data element
      int  info;       ///< element number. info \f$\in\f$ [0,thesize-1]
      ///< iff element is used
   }* theitem;         ///< array of elements in the DataSet
   ///@}
 
   //-----------------------------------
   /**@name Data */
   ///@{
   DataKey* thekey;    ///< DataKey::idx's of elements
   int themax;         ///< length of arrays theitem and thekey
   int thesize;        ///< highest used element in theitem
   int thenum;         ///< number of elements in DataSet
   int firstfree;      ///< first unused element in theitem
   ///@}
 
public:
 
   //-----------------------------------
   /**@name Extension
    *  Whenever a new element is added to a DataSet, the latter assigns it a
    *  DataKey. For this all methods that extend a DataSet by one ore more
    *  elements are provided with two signatures, one of them having a
    *  parameter for returning the assigned DataKey(s).
    */
   ///@{
   /// adds an element.
   void add(DataKey& newkey, const DATA& item)
   {
      DATA* data = create(newkey);
 
      assert(data != nullptr);
 
      *data = item;
   }
   /// adds element \p item.
   void add(const DATA& item)
   {
      DATA* data = create();
 
      assert(data != nullptr);
 
      *data = item;
   }
 
   /// add several items.
   void add(DataKey newkey[], const DATA* item, int n)
   {
      assert(n         >= 0);
      assert(num() + n <= max());
 
      for(int i = 0; i < n; ++i)
         add(newkey[i], item[i]);
   }
 
   /// adds \p n elements from \p items.
   void add(const DATA* items, int n)
   {
      assert(n         >= 0);
      assert(num() + n <= max());
 
      for(int i = 0; i < n; ++i)
         add(items[i]);
   }
 
   /// adds several new items.
   void add(DataKey newkey[], const DataSet < DATA >& set)
   {
      assert(num() + set.num() <= max());
 
      for(int i = 0; i < set.num(); ++i)
         add(newkey[i], set[i]);
   }
 
   /// adds all elements of \p set.
   void add(const DataSet < DATA >& set)
   {
      assert(num() + set.num() <= max());
 
      for(int i = 0; i < set.num(); ++i)
         add(set[i]);
   }
 
   /// creates new data element in DataSet.
   /**@return Pointer to the newly created element.
    */
   DATA* create(DataKey& newkey)
   {
      assert(num() < max());
 
      if(firstfree != -themax - 1)
      {
         newkey.idx = -firstfree - 1;
         firstfree = theitem[newkey.idx].info;
      }
      else
         newkey.idx = thesize++;
 
      thekey[thenum] = newkey;
      theitem[newkey.idx].info = thenum;
      ++thenum;
 
      return &(theitem[newkey.idx].data);
   }
   /// creates new (uninitialized) data element in DataSet.
   /**@return Pointer to the newly created element.
    */
   DATA* create()
   {
      DataKey tmp;
      return create(tmp);
   }
   ///@}
 
   //-----------------------------------
   /**@name Shrinkage
    * When elements are removed from a DataSet, the remaining ones are
    * renumbered from 0 through the new size()-1. How this renumbering is
    * performed will not be revealed, since it might be target of future
    * changes. However, some methods provide a parameter
    * <tt>int* perm</tt>, which
    * returns the new order after the removal: If <tt>perm[i] < 0</tt>,
    * the element numbered i prior to the removal operation has been removed
    * from the set. Otherwise, <tt>perm[i] = j >= 0</tt> means that the
    * element with number i prior to the removal operation has been
    * renumbered to j. Removing a single element from a DataSet yields a
    * simple renumbering of the elements: The last element in the set
    * (i.e., element num()-1) is moved to the index of the removed element.
    */
   ///@{
   /// removes the \p removenum 'th element.
   void remove(int removenum)
   {
      if(has(removenum))
      {
         int idx = thekey[removenum].idx;
 
         theitem[idx].info = firstfree;
         firstfree = -idx - 1;
 
         while(-firstfree == thesize)
         {
            firstfree = theitem[ -firstfree - 1].info;
            --thesize;
         }
 
         --thenum;
 
         if(removenum != thenum)
         {
            thekey[removenum] = thekey[thenum];
            theitem[thekey[removenum].idx].info = removenum;
         }
      }
   }
 
   /// removes element with key \p removekey.
   void remove(const DataKey& removekey)
   {
      remove(number(removekey));
   }
 
   /// remove multiple elements.
   /** This method removes all elements for the DataSet with an
    *  index i such that \p perm[i] < 0. Upon completion, \p perm contains
    *  the new numbering of elements.
    */
   void remove(int perm[])
   {
      int k, j, first = -1;
 
      // setup permutation and remove items
      for(k = j = 0; k < num(); ++k)
      {
         if(perm[k] >= 0)       // j has not been removed ...
            perm[k] = j++;
         else
         {
            int idx = thekey[k].idx;
            theitem[idx].info = firstfree;
            firstfree = -idx - 1;
 
            if(first < 0)
               first = k;
         }
      }
 
      if(first >= 0)         // move remaining items
      {
         for(k = first, j = num(); k < j; ++k)
         {
            if(perm[k] >= 0)
            {
               thekey[perm[k]] = thekey[k];
               theitem[thekey[k].idx].info = perm[k];
               thekey[k].idx = -1;
            }
            else
               --thenum;
         }
      }
   }
 
   /// remove \p n elements given by \p keys and \p perm.
   void remove(const DataKey* keys, int n, int* perm)
   {
      assert(perm != nullptr);
 
      for(int i = num() - 1; i >= 0; --i)
         perm[i] = i;
 
      while(--n >= 0)
         perm[number(keys[n])] = -1;
 
      remove(perm);
   }
   /// remove \p n elements given by \p keys.
   void remove(const DataKey* keys, int n)
   {
      DataArray<int> perm(num());
      remove(keys, n, perm.get_ptr());
   }
   /// remove \p n elements given by \p nums and \p perm.
   void remove(const int* nums, int n, int* perm)
   {
      assert(perm != nullptr);
 
      for(int i = num() - 1; i >= 0; --i)
         perm[i] = i;
 
      while(--n >= 0)
         perm[nums[n]] = -1;
 
      remove(perm);
   }
   /// remove \p n elements with numbers \p nums.
   void remove(const int* nums, int n)
   {
      DataArray<int> perm(num());
      remove(nums, n, perm.get_ptr());
   }
 
   /// remove all elements.
   void clear()
   {
      thesize = 0;
      thenum = 0;
      firstfree = -themax - 1;
   }
   ///@}
 
   //-----------------------------------
   /**@name Access
    * When accessing elements from a DataSet with one of the index
    * operators, it must be ensured that the index is valid for the
    * DataSet. If this is not known afore, it is the programmers
    * responsability to ensure this using the inquiry methods below.
    */
   ///@{
   ///
   DATA& operator[](int n)
   {
      assert(n >= 0 && n < thenum);
      return theitem[thekey[n].idx].data;
   }
   /// returns element number \p n.
   const DATA& operator[](int n) const
   {
      assert(n >= 0 && n < thenum);
      return theitem[thekey[n].idx].data;
   }
 
   ///
   DATA& operator[](const DataKey& k)
   {
      assert(k.idx < thesize);
      return theitem[k.idx].data;
   }
   /// returns element with DataKey \p k.
   const DATA& operator[](const DataKey& k) const
   {
      assert(k.idx < thesize);
      return theitem[k.idx].data;
   }
   ///@}
 
   //-----------------------------------
   /**@name Inquiry */
   ///@{
   /// returns maximum number of elements that would fit into DataSet.
   int max() const
   {
      return themax;
   }
 
   /// returns number of elements currently in DataSet.
   int num() const
   {
      return thenum;
   }
 
   /// returns the maximum DataKey::idx currently in DataSet.
   int size() const
   {
      return thesize;
   }
 
   /// returns DataKey of \p n 'th element in DataSet.
   DataKey key(int n) const
   {
      assert(n >= 0 && n < num());
      return thekey[n];
   }
 
   /// returns DataKey of element \p item in DataSet.
   DataKey key(const DATA* item) const
   {
      assert(number(item) >= 0);
      return thekey[number(item)];
   }
 
   /// returns the number of the element with DataKey \p k in DataSet or -1,
   /// if it doesn't exist.
   int number(const DataKey& k) const
   {
      if(k.idx < 0 || k.idx >= size())
         throw SPxException("Invalid index");
 
      return theitem[k.idx].info;
   }
 
   /**@todo Please check whether this is correctly implemented! */
   /// returns the number of element \p item in DataSet,
   /// throws exception if it doesn't exist.
   int number(const DATA* item) const
   {
      ptrdiff_t idx = reinterpret_cast<const struct Item*>(item) - theitem;
 
      if(idx < 0 || idx >= size())
         throw SPxException("Invalid index");
 
      return theitem[idx].info;
   }
 
   /// Is \p k a valid DataKey of an element in DataSet?
   bool has(const DataKey& k) const
   {
      return theitem[k.idx].info >= 0;
   }
 
   /// Is \p n a valid number of an element in DataSet?
   bool has(int n) const
   {
      return (n >= 0 && n < num());
   }
 
   /// Does \p item belong to DataSet?
   bool has(const DATA* item) const
   {
      int n;
 
      try
      {
         n = number(item);
      }
      catch(...)
      {
         return false;
      }
 
      return n >= 0;
   }
   ///@}
 
   //-----------------------------------
   /**@name Miscellaneous */
   ///@{
   /// resets max() to \p newmax.
   /** This method will not succeed if \p newmax < size(), in which case
    *  \p newmax == size() will be taken. As generally this method involves
    *  copying the #DataSet%s elements in memory, reMax() returns the
    *  number of bytes the addresses of elements in the DataSet have been
    *  moved. Note, that this is identical for all elements in the
    *  DataSet.
    */
   ptrdiff_t reMax(int newmax = 0)
   {
      struct Item* old_theitem = theitem;
      newmax = (newmax < size()) ? size() : newmax;
 
      int* lastfree = &firstfree;
 
      while(*lastfree != -themax - 1)
         lastfree = &(theitem[ -1 - *lastfree].info);
 
      *lastfree = -newmax - 1;
 
      themax = newmax;
      spx_realloc(theitem, themax);
      spx_realloc(thekey,  themax);
 
      return reinterpret_cast<char*>(theitem)
             - reinterpret_cast<char*>(old_theitem);
   }
 
   /// consistency check.
   bool isConsistent() const
   {
#ifdef ENABLE_CONSISTENCY_CHECKS
 
      if(theitem == 0 || thekey == 0)
         return SPX_MSG_INCONSISTENT("DataSet");
 
      if(thesize > themax || thenum > themax || thenum > thesize)
         return SPX_MSG_INCONSISTENT("DataSet");
 
      if(thesize == thenum && firstfree != -themax - 1)
         return SPX_MSG_INCONSISTENT("DataSet");
 
      if(thesize != thenum && firstfree == -themax - 1)
         return SPX_MSG_INCONSISTENT("DataSet");
 
      for(int i = 0; i < thenum; ++i)
         if(theitem[thekey[i].idx].info != i)
            return SPX_MSG_INCONSISTENT("DataSet");
 
#endif
 
      return true;
   }
   ///@}
 
   //-----------------------------------
   /**@name Constructors / Destructors */
   ///@{
   /// default constructor.
   explicit
   DataSet(int pmax = 8)
      : theitem(nullptr)
      , thekey(nullptr)
      , themax(pmax < 1 ? 8 : pmax)
      , thesize(0)
      , thenum(0)
 
   {
      firstfree = -themax - 1;
 
      spx_alloc(theitem, themax);
 
      try
      {
         spx_alloc(thekey, themax);
      }
      catch(const SPxMemoryException& x)
      {
         spx_free(theitem);
         throw x;
      }
 
      assert(isConsistent());
   }
 
   /// copy constructor.
   DataSet(const DataSet& old)
      : theitem(0)
      , thekey(nullptr)
      , themax(old.themax)
      , thesize(old.thesize)
      , thenum(old.thenum)
   {
      firstfree = (old.firstfree == -old.themax - 1)
                  ? -themax - 1
                  : old.firstfree;
 
      spx_alloc(theitem, themax);
 
      try
      {
         spx_alloc(thekey, themax);
      }
      catch(const SPxMemoryException& x)
      {
         spx_free(theitem);
         throw x;
      }
 
      memcpy(theitem, old.theitem, themax * sizeof(*theitem));
      memcpy(thekey,  old.thekey,  themax * sizeof(*thekey));
 
      assert(isConsistent());
   }
 
   /// assignment operator.
   /** The assignment operator involves #reMax()%ing the lvalue DataSet
    *  to the size needed for copying all elements of the rvalue. After the
    *  assignment all DataKeys from the lvalue are valid for the rvalue as
    *  well. They refer to a copy of the corresponding data elements.
    */
   DataSet < DATA >& operator=(const DataSet < DATA >& rhs)
   {
      if(this != &rhs)
      {
         int i;
 
         if(rhs.size() > max())
            reMax(rhs.size());
 
         clear();
 
         for(i = 0; i < rhs.size(); ++i)
            memcpy(&theitem[i], &rhs.theitem[i], sizeof(*theitem));
 
         for(i = 0; i < rhs.num(); ++i)
            thekey[i] = rhs.thekey[i];
 
         if(rhs.firstfree == -rhs.themax - 1)
            firstfree = -themax - 1;
         else
         {
            firstfree = rhs.firstfree;
            i = rhs.firstfree;
 
            while(rhs.theitem[ -i - 1].info != -rhs.themax - 1)
               i = rhs.theitem[ -i - 1].info;
 
            theitem[ -i - 1].info = -themax - 1;
         }
 
         thenum = rhs.thenum;
         thesize = rhs.thesize;
 
         assert(isConsistent());
      }
 
      return *this;
   }
 
   /// destructor.
   ~DataSet()
   {
      if(theitem)
         spx_free(theitem);
 
      if(thekey)
         spx_free(thekey);
   }
   ///@}
};
 
} // namespace soplex
#endif // _DATASET_H_