Changeset 642e65 in git for kernel/F5cData.cc

Timestamp:

Jul 7, 2009, 8:54:56 PM (14 years ago)

Author:

Christian Eder

Branches:

(u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')

Children:

ab3db62fd566bf974e5dcf6da2405b18b464b0c4

Parents:

e9f5b23f3e4572bcb33e40a8f5debc2ca9bd11a6

Message:

implementation of classes Label and CPair in F5cData.*


git-svn-id: file:///usr/local/Singular/svn/trunk@11957 2c84dea3-7e68-4137-9b89-c4e89433aadc

File:

• kernel/F5cData.cc (modified) (2 diffs)

kernel/F5cData.cc

@@ +1 @@
+/*!
+ * \file F5cData.cc
+ * \author Christian Eder
+ * \brief This file includes the inlined methods of the classes CPair and Rule.
+ * \details This file include the inlined methods of the class CPair,
+ * representing critical pairs in the F5C Algorithm, and the class Rule,
+ * representing the labels in the F5C Algorithm. In this file the handling of
+ * these data types is defined.
+ */
+
 #include "mod2.h"
 
@@ -18 +28 @@
 #include "F5cData.h"
 #include "F5cLists.h"
+
 /*
-=====================
-everything is inlined 
-=====================
-*/
+ * CONSTRUCTOR / DESTRUCTOR of class Label
+ */
+
+Label::Label(int* expVec) 
+{
+  m_pExpVec       = expVec;
+  m_nShortExpVec  = computeShortExpVec(m_pExpVec);
+}
+
+Label::~Label() 
+{
+  delete &m_nShortExpVec;
+  delete [] m_pExpVec;
+}
+
+/*
+ * GETTER of class Label
+ */
+inline int* Label::getExpVec() 
+{
+  return m_pExpVec;
+}
+
+inline unsigned long Label::getShortExpVec() 
+{
+  return m_nShortExpVec;
+}
+
+/*
+ * SETTER of class Label
+ */
+
+
+/*
+ * MISC of class Label
+ */
+inline unsigned long Label::computeShortExpVec(int* expVec) 
+{
+  //if (p == NULL) return 0;
+  unsigned long ev = 0; // short exponent vector
+  unsigned int n = BIT_SIZEOF_LONG / currRing->N; // number of bits per exp
+  unsigned int m1; // highest bit which is filled with (n+1)
+  unsigned int i = 0, j=1;
+
+  if (n == 0)
+  {
+    if (currRing->N <2*BIT_SIZEOF_LONG)
+    {
+      n=1;
+      m1=0;
+    }
+    else
+    {
+      for (; j<=(unsigned long) currRing->N; j++)
+      {
+        if (expVec[j] > 0) i++;
+        if (i == BIT_SIZEOF_LONG) break;
+      }
+      if (i>0)
+      ev = ~((unsigned long)0) >> ((unsigned long) (BIT_SIZEOF_LONG - i));
+      return ev;
+    }
+  }
+  else
+  {
+    m1 = (n+1)*(BIT_SIZEOF_LONG - n*currRing->N);
+  }
+
+  n++;
+  while (i<m1)
+  {
+    ev |= getBitFields(expVec[j], i, n);
+    i += n;
+    j++;
+  }
+
+  n--;
+  while (i<BIT_SIZEOF_LONG)
+  {
+    ev |= getBitFields(expVec[j], i, n);
+    i += n;
+    j++;
+  }
+  return ev;
+}
+
+inline unsigned long Label::getBitFields(int e, unsigned int s, unsigned int n)
+{
+#define Sy_bit_L(x)     (((unsigned long)1L)<<(x))
+  unsigned int i = 0;
+  unsigned long  ev = 0L;
+  assume(n > 0 && s < BIT_SIZEOF_LONG);
+  do
+  {
+    assume(s+i < BIT_SIZEOF_LONG);
+    if (e > (int) i) ev |= Sy_bit_L(s+i);
+    else break;
+    i++;
+  }
+  while (i < n);
+  return ev;
+}
+
+
+
 #endif
-// HAVE_F5c
+// HAVE_F5C