pMinMax_impl.h 9.77 KB
Newer Older
Pierre Aubert's avatar
Pierre Aubert committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17

/***************************************
	Auteur : Pierre Aubert
	Mail : aubertp7@gmail.com
	Licence : CeCILL-C
****************************************/

#ifndef __PMIN_MAX_IMPL_H__
#define __PMIN_MAX_IMPL_H__

#include "pMinMax.h"

///fonction qui renvoie la plus petite valeurs entre a et b
/**	@param a : valeur
 * 	@param b : valeur
 * 	@return petite valeurs entre a et b
*/
Pierre Aubert's avatar
Pierre Aubert committed
18
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
19 20 21 22 23 24 25 26 27
T pmin(const T & a, const T & b){
	return (a < b)? a:b;
}

///fonction qui renvoie la plus grande valeurs entre a et b
/**	@param a : valeur
 * 	@param b : valeur
 * 	@return grande valeurs entre a et b
*/
Pierre Aubert's avatar
Pierre Aubert committed
28
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
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
T pmax(const T & a, const T & b){
	return (a > b)? a:b;
}

///renvoie le minimum entre x, y, z et t
/**	@param x : nombre
	@param y : nombre
	@param z : nombre
	@param t : nombre
	@return minimum entre x, y, z et t
*/
template <class T>
T pmin4(const T & x, const T & y, const T & z, const T & t){
	if(x < y){
		if(x < z){
			if(x < t) return x;
			else return t;
		}else{
			if(z < t) return z;
			else return t;
		}
	}else{
		if(y < z){
			if(y < t) return y;
			else return t;
		}else{
			if(z < t) return z;
			else return t;
		}
	}
}

///renvoie le maximum entre x, y, z et t
/**	@param x : nombre
	@param y : nombre
	@param z : nombre
	@param t : nombre
	@return maximum entre x, y, z et t
*/
template <class T>
T pmax4(const T & x, const T & y, const T & z, const T & t){
	if(x > y){
		if(x > z){
			if(x > t) return x;
			else return t;
		}else{
			if(z > t) return z;
			else return t;
		}
	}else{
		if(y > z){
			if(y > t) return y;
			else return t;
		}else{
			if(z > t) return z;
			else return t;
		}
	}
}

///fonction qui renvoie le minimum d'un tableau
/**	@param tab : tableau dont on veut le minimum
 * 	@param size : taille du tableau
 * 	@return minimum du tableau
*/
Pierre Aubert's avatar
Pierre Aubert committed
94
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110
T pminTab(const T* tab, size_t size){
	if(tab == NULL || size == 0lu) return 0.0;
	T res = tab[0];
	if(size == 1lu) return res;
	size_t i;
	for(i = 1lu; i < size; ++i){
		if(tab[i] < res) res = tab[i];
	}
	return res;
}

///fonction qui renvoie la position du minimum d'un tableau
/**	@param tab : tableau dont on veut le minimum
 * 	@param size : taille du tableau
 * 	@return position du minimum du tableau
*/
Pierre Aubert's avatar
Pierre Aubert committed
111
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131
size_t pminTabPos(const T* tab, size_t size){
	if(tab == NULL || size == 0lu) return 0lu;
	T res = tab[0];
	if(size == 1lu) return 0lu;
	size_t i, pos(0lu);
	for(i = 1lu; i < size; ++i){
		if(tab[i] < res){
			res = tab[i];
			pos = i;
		}
	}
	return pos;
}

///fonction qui renvoie la position du minimum d'un tableau
/**	@param[out] minVal : minimum value
 * 	@param tab : tableau dont on veut le minimum
 * 	@param size : taille du tableau
 * 	@return position du minimum du tableau
*/
Pierre Aubert's avatar
Pierre Aubert committed
132
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151
size_t pminTabPos(T & minVal, const T* tab, size_t size){
	if(tab == NULL || size == 0lu) return 0lu;
	minVal = tab[0];
	if(size == 1lu) return 0lu;
	size_t i, pos(0lu);
	for(i = 1lu; i < size; ++i){
		if(tab[i] < minVal){
			minVal = tab[i];
			pos = i;
		}
	}
	return pos;
}

///fonction qui renvoie le maximum d'un tableau
/**	@param tab : tableau dont on veut le maximum
 * 	@param size : taille du tableau
 * 	@return maximum du tableau
*/
Pierre Aubert's avatar
Pierre Aubert committed
152
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168
T pmaxTab(const T* tab, size_t size){
	if(tab == NULL || size == 0lu) return 0.0;
	T res = tab[0];
	if(size == 1lu) return res;
	size_t i;
	for(i = 1lu; i < size; ++i){
		if(tab[i] > res) res = tab[i];
	}
	return res;
}

///fonction qui renvoie la position du maximum d'un tableau
/**	@param tab : tableau dont on veut le maximum
 * 	@param size : taille du tableau
 * 	@return position du maximum du tableau
*/
Pierre Aubert's avatar
Pierre Aubert committed
169
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189
size_t pmaxTabPos(const T* tab, size_t size){
	if(tab == NULL || size == 0lu) return 0lu;
	T res = tab[0];
	if(size == 1lu) return 0lu;
	size_t i, pos(0lu);
	for(i = 1lu; i < size; ++i){
		if(tab[i] > res){
			res = tab[i];
			pos = i;
		}
	}
	return pos;
}

///fonction qui renvoie la position du maximum d'un tableau
/**	@param[out] maxVal : minimum value
 * 	@param tab : tableau dont on veut le minimum
 * 	@param size : taille du tableau
 * 	@return position du maximum du tableau
*/
Pierre Aubert's avatar
Pierre Aubert committed
190
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210
size_t pmaxTabPos(T & maxVal, const T* tab, size_t size){
	if(tab == NULL || size == 0lu) return 0lu;
	maxVal = tab[0];
	if(size == 1lu) return 0lu;
	size_t i, pos(0lu);
	for(i = 1lu; i < size; ++i){
		if(tab[i] > maxVal){
			maxVal = tab[i];
			pos = i;
		}
	}
	return pos;
}

///Gets the minimum and maximum values in a table
/**	@param[out] xmin : minimum value in the table
 * 	@param[out] xmax : maximum value in the table
 * 	@param tab : table we want min and max
 * 	@param size : number of elements in the table
*/
Pierre Aubert's avatar
Pierre Aubert committed
211
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
212 213 214 215 216 217 218 219 220 221
void pminmaxTab(T & xmin, T & xmax, const T* tab, size_t size){
	if(tab == NULL || size == 0lu) return;
	xmin = tab[0];
	xmax = xmin;
	if(size == 1lu) return;
	size_t i;
	T tmp;
	for(i = 1; i < size; ++i){
		tmp = tab[i];
		if(tmp > xmax) xmax = tmp;
Pierre Aubert's avatar
Pierre Aubert committed
222
		else{xmin = tmp;}
Pierre Aubert's avatar
Pierre Aubert committed
223 224 225 226 227 228 229 230 231
	}
}

///Compute the miminum of two tables
/**	@param[out] vecMin : vector of the minimum element of tab1 and tab2
 * 	@param tab1 : table
 * 	@param tab2 : table
 * 	@param size : number of element in each table
*/
Pierre Aubert's avatar
Pierre Aubert committed
232
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
233 234 235 236 237 238 239 240 241 242 243 244
void pminVec(T * vecMin, const T* tab1, const T * tab2, size_t size){
	for(size_t i(0lu); i < size; ++i){
		vecMin[i] = pmin(tab1[i], tab2[i]);
	}
}

///Compute the miminum of two tables
/**	@param[out] vecMax : vector of the maximum element of tab1 and tab2
 * 	@param tab1 : table
 * 	@param tab2 : table
 * 	@param size : number of element in each table
*/
Pierre Aubert's avatar
Pierre Aubert committed
245
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
246 247 248 249 250 251 252 253 254 255 256 257 258
void pmaxVec(T * vecMax, const T* tab1, const T * tab2, size_t size){
	for(size_t i(0lu); i < size; ++i){
		vecMax[i] = pmax(tab1[i], tab2[i]);
	}
}

///Gets the minimum and maximum values in a matrix
/**	@param[out] tabMin : minimum value in the matrix for each column
 * 	@param[out] tabMax : maximum value in the matrix for each column
 * 	@param nbCol : number of columns in the matrix (and number of elements in the tabXmin and tabXmax table)
 * 	@param matrix : table we want min and max
 * 	@param nbRow : number of rows in the matrix
*/
Pierre Aubert's avatar
Pierre Aubert committed
259
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282
void pminmaxMatrix(T * tabMin, T * tabMax, size_t nbCol, const T* matrix, size_t nbRow){
	if(tabMax == NULL || tabMin == NULL || nbCol == 0lu || matrix == NULL || nbRow == 0lu) return;
	memcpy(tabMin, matrix, sizeof(T)*nbCol);
	memcpy(tabMax, tabMin, sizeof(T)*nbCol);
	if(nbRow == 1lu) return;
	T tmp;
	for(size_t i(0lu); i < nbRow; ++i){
		for(size_t j(0lu); j < nbCol; ++j){
			tmp = matrix[i*nbCol + j];
			if(tmp > tabMax[j]) tabMax[j] = tmp;
			else if(tmp < tabMin[j]) tabMin[j] = tmp;
		}
	}
}

///Get the minimum and maximum value in a Matrix pitch
/**	@param [out] xmin : minimum value in the Matrix pitch
 * 	@param[out] xmax : maximum value in the Matrix pitch
 * 	@param mat : input matrix
 * 	@param nbRow : number of rows in the matrix pitch
 * 	@param nbCol : number of columns in the matrix pitch
 * 	@param pitch : number of extra columns used to aligned the first element of each line of the matrix (so it is a matrix pitch)
*/
Pierre Aubert's avatar
Pierre Aubert committed
283
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305
void pminmaxMatrixPitch(T & xmin, T & xmax, const T* mat, size_t nbRow, size_t nbCol, size_t pitch){
	if(mat == NULL || nbRow == 0lu || nbCol == 0lu) return;
	xmin = mat[0];
	xmax = xmin;
	T tmp;
	for(size_t i(0lu); i < nbRow; ++i){
		for(size_t j(0lu); j < nbCol; ++j){
			tmp = mat[i*(nbCol + pitch) + j];
			if(tmp > xmax) xmax = tmp;
			else if(tmp < xmin) xmin = tmp;
		}
	}
}

///Fonction qui calcule le minimum et le maximum des paires de valeurs dans un tableaux à une dimension
/**	@param xmin : valeur minimale de x
 * 	@param xmax : valeur maximale de x
 * 	@param ymin : valeur minimale de y
 * 	@param ymax : valeur maximale de y
 * 	@param tab : tableau de valeurs de XY
 * 	@param sizeNbPair : nombre de paires (moitié de la taille du tableau de XY totale)
*/
Pierre Aubert's avatar
Pierre Aubert committed
306
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327
void pminmaxTabXY(T & xmin, T & xmax, T & ymin, T & ymax, const T * tab, size_t sizeNbPair){
	if(tab == NULL || sizeNbPair == 0lu) return;
	xmin = tab[0lu];
	xmax = tab[0lu];
	ymin = tab[1lu];
	ymax = tab[1lu];
	if(sizeNbPair > 1lu){
		for(size_t i(0lu); i < sizeNbPair; ++i){
			if(xmin > tab[2lu*i]) xmin = tab[2lu*i];
			if(xmax < tab[2lu*i]) xmax = tab[2lu*i];
			if(ymin > tab[2lu*i + 1lu]) ymin = tab[2lu*i + 1lu];
			if(ymax < tab[2lu*i + 1lu]) ymax = tab[2lu*i + 1lu];
		}
	}
}

///Converts a table of couple (x,y) in a table of (x-y) values
/**	@param tabX : table of (x-y) values
 * 	@param tabXY : table of couple (x,y)
 * 	@param sizeTabX : size of the tabX table and half-size of the tabXY one
*/
Pierre Aubert's avatar
Pierre Aubert committed
328
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
329 330 331 332 333 334 335 336 337 338 339 340
void tabXYtoTabXmY(T* tabX, const T* tabXY, size_t sizeTabX){
	if(tabX == NULL || tabXY == NULL || sizeTabX == 0lu) return;
	for(size_t i(0lu); i < sizeTabX; ++i){
		tabX[i] = tabXY[2lu*i] - tabXY[2lu*i + 1lu];
	}
}

///Converts a table of couple (x,y) in a table of (x-y) values
/**	@param tabX : table of (x-y) values
 * 	@param tabXY : table of couple (x,y)
 * 	@param sizeTabX : size of the tabX table and half-size of the tabXY one
*/
Pierre Aubert's avatar
Pierre Aubert committed
341
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
342 343 344 345 346 347 348 349 350 351 352 353 354 355
void tabXYtoTabXmYonX(T* tabX, const T* tabXY, size_t sizeTabX){
	if(tabX == NULL || tabXY == NULL || sizeTabX == 0lu) return;
	for(size_t i(0lu); i < sizeTabX; ++i){
		if(fabs(tabXY[2lu*i]) > 0.00001f) tabX[i] = (tabXY[2lu*i] - tabXY[2lu*i + 1lu])/tabXY[2lu*i];
		else tabX[i] = 0.0f;
	}
}

///Get the number of diffrents value in a table
/**	@param tab : table of value
 * 	@param size : number of elements in the table
 * 	@param threshold : considers the minimal distance between to values
 * 	@return number of diffrents value in a table
*/
Pierre Aubert's avatar
Pierre Aubert committed
356
template<typename T>
Pierre Aubert's avatar
Pierre Aubert committed
357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372
size_t getNumberOfDiffrentValues(const T * tab, size_t size, float threshold){
	if(tab == NULL || size == 0lu) return 0lu;
	size_t nb(1lu);
	float max(tab[0]);
	for(size_t i(1lu); i < size; ++i){
		if(max + threshold < tab[i]){
			++nb;
			max = tab[i];
		}
	}
	return nb;
}

#endif