blocksparse/sparsematrix.h at master · onto/blocksparse · 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
45
46
47
48
49
50
51
52
53
54
55
56
#ifndef SPARSEMATRIX_H
#define SPARSEMATRIX_H

#include <vector>
#include <cstdlib>
#include <fstream>
#include <iostream>
#include <iomanip>
#include <ios>
#include <omp.h>

typedef std::vector<size_t> VectSizet;
typedef std::vector<size_t>::iterator VectSizetIterator;
typedef std::vector<double> VectDouble;
typedef std::vector<bool> VectBool;
typedef std::vector< VectSizet > VectVectSizet;

#define SPARSE_END 0

class SparseMatrix
{
public:
    SparseMatrix();
    SparseMatrix(size_t h, size_t w);
    SparseMatrix(const SparseMatrix &S);
    SparseMatrix(const char *file);

    VectDouble V;  //values
    VectSizet N;   //next in row
    VectSizet C;   //column
    VectSizet F;   //first in row

    size_t W; //width
    size_t H; //height

    void set(size_t row, size_t col, double value);
    double get(size_t row, size_t col) const;
    void rm(size_t row, size_t col);
    void add(size_t row, size_t col, double value);

    void swapCol(size_t c1, size_t c2);
    void swapRow(size_t r1, size_t r2);

    void permute(std::vector<size_t>& P);
    void permute(std::vector<size_t>& Pr, std::vector<size_t>& Pc);

    void print() const;
    void printStructure() const;

    void save2file(const char *file) const;
    void save2fileold(const char *file) const;

    void save2fileTransposed(const char *file) const;
};

#endif // SPARSEMATRIX_H