fix test.constructionModelesLassoMLE; TODO: det(rho[,,r]) ? numerical errors ? ....
[valse.git] / src / sources / utils.h
CommitLineData
8e92c49c
BA
1#ifndef valse_utils_H
2#define valse_utils_H
3
9ff729fb
BA
4//#include <stdint.h>
5
6/********
7 * Types
8 *******/
9
c3bc4705 10typedef double Real;
9ff729fb
BA
11//typedef uint32_t UInt;
12//typedef int32_t Int;
13
8e92c49c
BA
14/*******************
15 * tune parallelism
16 *******************/
17
18// Number of OpenMP threads
19#define OMP_NUM_THREADS 8
20
21// CHUNK_SIZE = number of lambda values to be treated sequentially by a single core
22#define CHUNK_SIZE 1
23
24/*******************************
25 * Matrix and arrays indexation
26 *******************************/
27
28// Matrix Index ; TODO? ncol unused
29#define mi(i,j,nrow,ncol)\
30 j*nrow + i
31
32// Array Index ; TODO? d3 unused
33#define ai(i,j,k,d1,d2,d3)\
34 k*d1*d2 + j*d1 + i
35
36// Array4 Index ; TODO? ...
37#define ai4(i,j,k,m,d1,d2,d3,d4)\
38 m*d1*d2*d3 + k*d1*d2 + j*d1 + i
39
8e92c49c
BA
40/*************************
41 * Array copy & "zeroing"
42 ************************/
43
44// Fill an array with zeros
45#define zeroArray(array, size)\
46{\
47 for (int u=0; u<size; u++)\
48 array[u] = 0;\
49}
50
51// Copy an 1D array
52#define copyArray(array, copy, size)\
53{\
54 for (int u=0; u<size; u++)\
55 copy[u] = array[u];\
56}
57
58#endif