matrix/mioarray.h at master · MAYHEM-Lab/matrix · GitHub

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
#ifndef ARRAY_H
#define ARRAY_H

#include "mio.h"

struct array_stc_2d
{
        int xdim;
        int ydim;
	MIO *mio;
        double *data;
};

typedef struct array_stc_2d Array2D;
Array2D *MakeArray2D(int xdim, int ydim);
Array2D *MakeArray2DFromMIO(MIO *a_mio);
void FreeArray2D(Array2D *a);
void PrintArray2D(Array2D *a);
Array2D *CopyArray2D(Array2D *a);
Array2D *TransposeArray2D(Array2D *a);
Array2D *MultiplyArray2D(Array2D *a, Array2D *b);
Array2D *InvertArray2D(Array2D *a);
Array2D *AddArray2D(Array2D *a, Array2D *b);
Array2D *SubtractArray2D(Array2D *a, Array2D *b);
Array2D *NormalizeRowsArray2D(Array2D *a);
Array2D *NormalizeColsArray2D(Array2D *a);
double NormArray1D(Array2D *a);

Array2D *EigenVectorArray2D(Array2D *a);

/*
 * choose 1D to be column vectors
 */

#define Array1D Array2D
#define MakeArray1D(dim) MakeArray2D((dim),1)
#define MakeArray1DFromMIO(dim) MakeArray2DFromMIO((dim),1)
#define FreeArray1D(a) FreeArray2D(a)
#define PrintArray1D(a) PrintArray2D(a)
#define CopyArray1D(a) CopyArray2D(a)
#define TransposeArray1D(a) TransposeArray2D(a)

Array1D *EigenValueArray2D(Array2D *a);
#endif