b78d46a32f04c91cd5826f32f36bc02551ffb28c
[ginac.git] / 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 template<class T>
31 string ToString(T const & t)
32 {
33     char buf[256];
34     ostrstream(buf,sizeof(buf)) << t << ends;
35     return buf;
36 }
37
38 unsigned log2(unsigned n);
39
40 int compare_pointers(void const * a, void const * b);
41
42 #define DYNCONSTCAST(FINALTYPE,BASICTYPE,EXPRESSION) \
43     dynamic_cast<FINALTYPE>(const_cast<BASICTYPE>(EXPRESSION))
44
45 // modified from stl_algo.h: always do com(*first1,*first2) instead of comp(*first2,*first1)
46 template <class InputIterator1, class InputIterator2, class OutputIterator,
47           class Compare>
48 OutputIterator mymerge(InputIterator1 first1, InputIterator1 last1,
49                        InputIterator2 first2, InputIterator2 last2,
50                        OutputIterator result, Compare comp) {
51   while (first1 != last1 && first2 != last2) {
52     if (comp(*first1, *first2)) {
53       *result = *first1;
54       ++first1;
55     }
56     else {
57       *result = *first2;
58       ++first2;
59     }
60     ++result;
61   }
62   return copy(first2, last2, copy(first1, last1, result));
63 }
64
65 // like merge(), but three lists with *last2<*first3
66 template <class InputIterator1, class InputIterator2, class InputIterator3,
67           class OutputIterator, class Compare>
68 OutputIterator mymerge3(InputIterator1 first1, InputIterator1 last1,
69                         InputIterator2 first2, InputIterator2 last2,
70                         InputIterator3 first3, InputIterator3 last3,
71                         OutputIterator result, Compare comp) {
72   while (first1 != last1 && first2 != last2) {
73     if (comp(*first1, *first2)) {
74       *result = *first1;
75       ++first1;
76     }
77     else {
78       *result = *first2;
79       ++first2;
80     }
81     ++result;
82   }
83
84   if (first1==last1) {
85     // list1 empty, copy rest of list2, then list3
86     return copy(first3, last3, copy(first2, last2, result));
87   } else {
88     // list2 empty, merge rest of list1 with list3
89     return mymerge(first1,last1,first3,last3,result,comp);
90   }
91 }
92
93 #endif // ndef __GINAC_UTILS_H__