ginac/utils.h

   1 /** @file utils.h
   2  *
   3  *  Interface to several small and furry utilities. */
   4
   5 /*
   6  *  GiNaC Copyright (C) 1999 Johannes Gutenberg University Mainz, Germany
   7  *
   8  *  This program is free software; you can redistribute it and/or modify
   9  *  it under the terms of the GNU General Public License as published by
  10  *  the Free Software Foundation; either version 2 of the License, or
  11  *  (at your option) any later version.
  12  *
  13  *  This program is distributed in the hope that it will be useful,
  14  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
  15  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  16  *  GNU General Public License for more details.
  17  *
  18  *  You should have received a copy of the GNU General Public License
  19  *  along with this program; if not, write to the Free Software
  20  *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  21  */
  22
  23 #ifndef __GINAC_UTILS_H__
  24 #define __GINAC_UTILS_H__
  25
  26 #include <strstream>
  27 #include <string>
  28 #include "config.h"
  29
  30 namespace GiNaC {
  31
  32 template<class T>
  33 string ToString(T const & t)
  34 {
  35     char buf[256];
  36     ostrstream(buf,sizeof(buf)) << t << ends;
  37     return buf;
  38 }
  39
  40 unsigned log2(unsigned n);
  41
  42 int compare_pointers(void const * a, void const * b);
  43
  44 // modified from stl_algo.h: always do com(*first1,*first2) instead of comp(*first2,*first1)
  45 template <class InputIterator1, class InputIterator2, class OutputIterator,
  46           class Compare>
  47 OutputIterator mymerge(InputIterator1 first1, InputIterator1 last1,
  48                        InputIterator2 first2, InputIterator2 last2,
  49                        OutputIterator result, Compare comp) {
  50   while (first1 != last1 && first2 != last2) {
  51     if (comp(*first1, *first2)) {
  52       *result = *first1;
  53       ++first1;
  54     }
  55     else {
  56       *result = *first2;
  57       ++first2;
  58     }
  59     ++result;
  60   }
  61   return copy(first2, last2, copy(first1, last1, result));
  62 }
  63
  64 // like merge(), but three lists with *last2<*first3
  65 template <class InputIterator1, class InputIterator2, class InputIterator3,
  66           class OutputIterator, class Compare>
  67 OutputIterator mymerge3(InputIterator1 first1, InputIterator1 last1,
  68                         InputIterator2 first2, InputIterator2 last2,
  69                         InputIterator3 first3, InputIterator3 last3,
  70                         OutputIterator result, Compare comp) {
  71   while (first1 != last1 && first2 != last2) {
  72     if (comp(*first1, *first2)) {
  73       *result = *first1;
  74       ++first1;
  75     }
  76     else {
  77       *result = *first2;
  78       ++first2;
  79     }
  80     ++result;
  81   }
  82
  83   if (first1==last1) {
  84     // list1 empty, copy rest of list2, then list3
  85     return copy(first3, last3, copy(first2, last2, result));
  86   } else {
  87     // list2 empty, merge rest of list1 with list3
  88     return mymerge(first1,last1,first3,last3,result,comp);
  89   }
  90 }
  91
  92 } // namespace GiNaC
  93
  94 #endif // ndef __GINAC_UTILS_H__