Oups.
[pysvrt.git] / misc.h
1 /*
2  *  svrt is the ``Synthetic Visual Reasoning Test'', an image
3  *  generator for evaluating classification performance of machine
4  *  learning systems, humans and primates.
5  *
6  *  Copyright (c) 2009 Idiap Research Institute, http://www.idiap.ch/
7  *  Written by Francois Fleuret <francois.fleuret@idiap.ch>
8  *
9  *  This file is part of svrt.
10  *
11  *  svrt is free software: you can redistribute it and/or modify it
12  *  under the terms of the GNU General Public License version 3 as
13  *  published by the Free Software Foundation.
14  *
15  *  svrt is distributed in the hope that it will be useful, but
16  *  WITHOUT ANY WARRANTY; without even the implied warranty of
17  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
18  *  General Public License for more details.
19  *
20  *  You should have received a copy of the GNU General Public License
21  *  along with svrt.  If not, see <http://www.gnu.org/licenses/>.
22  *
23  */
24
25 #ifndef MISC_H
26 #define MISC_H
27
28 #include <iostream>
29 #include <cmath>
30 #include <fstream>
31 #include <cfloat>
32 #include <stdlib.h>
33 #include <string.h>
34
35 using namespace std;
36
37 typedef double scalar_t;
38 // typedef float scalar_t;
39
40 const int buffer_size = 1024;
41
42 using namespace std;
43
44 #ifdef DEBUG
45 #define ASSERT(x) if(!(x)) { \
46   std::cerr << "ASSERT FAILED IN " << __FILE__ << ":" << __LINE__ << endl; \
47   abort(); \
48 }
49 #else
50 #define ASSERT(x)
51 #endif
52
53 template<class T>
54 T **allocate_array(int a, int b) {
55   T *tmp = new T[a * b];
56   T **array = new T *[a];
57   for(int k = 0; k < a; k++) {
58     array[k] = tmp;
59     tmp += b;
60   }
61   return array;
62 }
63
64 template<class T>
65 void deallocate_array(T **array) {
66   delete[] array[0];
67   delete[] array;
68 }
69
70 template <class T>
71 void write_var(ostream *os, const T *x) { os->write((char *) x, sizeof(T)); }
72
73 template <class T>
74 void read_var(istream *is, T *x) { is->read((char *) x, sizeof(T)); }
75
76 template <class T>
77 inline T sq(T x) {
78   return x * x;
79 }
80
81 inline scalar_t log2(scalar_t x) {
82   return log(x)/log(2.0);
83 }
84
85 template <class T>
86 void grow(int *nb_max, int nb, T** current, int factor) {
87   ASSERT(*nb_max > 0);
88   if(nb == *nb_max) {
89     T *tmp = new T[*nb_max * factor];
90     memcpy(tmp, *current, *nb_max * sizeof(T));
91     delete[] *current;
92     *current = tmp;
93     *nb_max *= factor;
94   }
95 }
96
97 #endif