1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
|
/***************************************************************************
defs.hh - Commonly used definitions
-------------------
begin : Sun Jun 11 2000
copyright : (C) 2000 by Erik Schnetter
email : schnetter@astro.psu.edu
$Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/defs.hh,v 1.7 2002/03/11 13:17:13 schnetter Exp $
***************************************************************************/
/***************************************************************************
* *
* This program is free software; you can redistribute it and/or modify *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation; either version 2 of the License, or *
* (at your option) any later version. *
* *
***************************************************************************/
#ifndef DEFS_HH
#define DEFS_HH
#ifdef HAVE_CONFIG_H
#include <config.h>
#endif
#include <assert.h>
#include <algorithm>
#include <iostream>
#include <list>
#include <set>
#include <vector>
using namespace std;
// Stringification
#define STR(s) #s
// Fortranification
#define FORTRAN_NAME(x) x##_
// A general type
enum centering { vertex_centered, cell_centered };
// Useful helper
template<class T>
inline T square (const T& x) { return x*x; }
// Another useful helper
template<class T>
inline T ipow (const T& x, const int y) {
if (y<0) {
return T(1)/ipow(x,-y);
} else if (y==0) {
return T(1);
} else if (y%2) {
return x * ipow(x*x,y/2);
} else {
return ipow(x*x,y/2);
}
}
// Skip whitespace
void skipws (istream& is);
// Container input
template<class T> istream& input (istream& is, vector<T>& v);
template<class T>
inline istream& operator>> (istream& is, vector<T>& v) {
return input(is,v);
}
// Container output
template<class T> ostream& output (ostream& os, const list<T>& l);
template<class T> ostream& output (ostream& os, const set<T>& s);
template<class T> ostream& output (ostream& os, const vector<T>& v);
template<class T>
inline ostream& operator<< (ostream& os, const list<T>& l) {
return output(os,l);
}
template<class T>
inline ostream& operator<< (ostream& os, const set<T>& s) {
return output(os,s);
}
template<class T>
inline ostream& operator<< (ostream& os, const vector<T>& v) {
return output(os,v);
}
#if defined(TMPL_IMPLICIT)
# include "defs.cc"
#endif
#endif // DEFS_HH
|
