Changeset dfed1c for src/base/helper.hpp

Timestamp:

Nov 22, 2011, 9:22:10 PM (14 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

facba0

Parents:

66f24d

Message:

Major vmg update.

git-svn-id: ​https://svn.version.fz-juelich.de/scafacos/trunk@1136 5161e1c8-67bf-11de-9fd5-51895aff932f

File:

• src/base/helper.hpp (modified) (6 diffs)

src/base/helper.hpp

@@ -3 +3 @@
  * @author Julian Iseringhausen <isering@ins.uni-bonn.de>
  * @date   Tue Apr  5 21:03:47 2011
- * 
+ *
  * @brief  Provides various helper functions.
- * 
- * 
+ *
  */
 
@@ -15 +14 @@
 #include <cstdio>
 #include <limits>
+#include <sstream>
+#include <string>
 
 namespace VMG
 {
 
+class Grid;
+class Vector;
+
 namespace Helper
 {
+  char* GetCharArray(const std::string& str);
+  std::string ReplaceWhitespaces(const char* buffer, const char* replace);
+
+  template <class T>
+  std::string ToString(const T& val)
+  {
+    std::stringstream str;
+    str << std::scientific << val;
+    return str.str();
+  }
+
+  template <class T>
+  T ToVal(const char* val_str)
+  {
+    T val;
+    std::stringstream str;
+    str << val_str;
+    str >> val;
+
+    return val;
+  }
+
+  template <class T>
+  T ToValWithDefault(const char* val_str, const T& def)
+  {
+    T val;
+    std::stringstream str(val_str);
+    str >> val;
+
+    if (str.fail() || str.bad() || !str.eof()) {
+#ifdef DEBUG_VERBOSE
+      std::printf("VMG::Helper::ToValWithDefault: Using default value.\n");
+#endif
+      val = def;
+    }
+
+    return val;
+  }
 
   /**
    * Checks a number for validity, i.e. it is neither nan nor inf.
-   * 
    */
   inline bool CheckNumber(const vmg_float& number)
@@ -50 +91 @@
   }
 
-  inline int intpow(unsigned int base, unsigned int power)
+  inline int intpow(int base, unsigned int power)
   {
-    assert (power >= 0);
-    unsigned int result = 1;
+    int result = 1;
     while (power != 0) {
       if (power & 1)
@@ -62 +102 @@
     return result;
   }
+
+  inline vmg_float pow(vmg_float base, unsigned int power)
+  {
+    vmg_float result = 1.0;
+    while (power != 0) {
+      if (power & 1)
+        result *= base;
+      base *= base;
+      power >>= 1;
+    }
+    return result;
+  }
+
+  inline unsigned int fact(unsigned int number)
+  {
+    unsigned int result = 1;
+    for (unsigned int i=2; i<=number; ++i)
+      result *= i;
+    return result;
+  }
+
+  inline vmg_float pow_3(vmg_float val)
+  {
+    return val*val*val;
+  }
+
 
   /**
@@ -72 +138 @@
     bool rval = (val1 == val2);
 
+#ifdef DEBUG_OUTPUT
     if (!rval)
       printf("WARNING: Values are not equal (%s)\n", name);
+#endif /* DEBUG_OUTPUT */
 
     assert(rval);
@@ -80 +148 @@
   }
 
-  inline void MM(int nn, vmg_float u, vmg_float *tbl)
-  {
-    int i, n, m;
-
-    assert(u >= 0.0 && u <= 1.0);
-
-    /* n = 1: */
-    /* i = 0: */
-    tbl[0] = 1.0;
-
-    for (n = 2; n <= nn; n++) {
-      m = n - 1;
-
-      /* i = n - 1: */
-      tbl[n-1] =  (1 - u) * tbl[n-2] / m;
-
-      for (i = n - 2; i >= 1; i--)
-        tbl[i] = ((u + i) * tbl[i] + (n - (u + i)) * tbl[i-1]) / m;
-
-      /* i = 0: */
-      tbl[0] *= u / m;
-    }
-  }
-
-  inline void MdM(int nn, vmg_float u,
-                  vmg_float * restrict tbl, vmg_float * restrict dtbl)
-  {
-    int i, n, m;
-
-    assert(nn >= 2);
-    assert(u >= 0.0 && u <= 1.0);
-
-    /* n = 1: */
-    /* i = 0: */
-    tbl[0] = 1.0;
-
-    for (n = 2; n <= nn; n++) {
-      m = n - 1;
-
-      /* differentiation */
-      if (n == nn) {
-        dtbl[n-1] = - tbl[n-2];
-
-        for (i = n - 2; i >= 1; i--)
-          dtbl[i] = tbl[i] - tbl[i-1];
-
-        dtbl[0] = tbl[0];
-      }
-
-      /* i = n - 1: */
-      tbl[n-1] =  (1 - u) * tbl[n-2] / m;
-
-      for (i = n - 2; i >= 1; i--)
-        tbl[i] = ((u + i) * tbl[i] + (n - (u + i)) * tbl[i-1]) / m;
-
-      /* i = 0: */
-      tbl[0] *= u / m;
-
-    }
-  }
+  vmg_float InterpolateTrilinear(const Vector& point, const Grid& grid);
 
 }