blocksparse/matrixcontainers.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
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
94
95
96
97
98
99
100
101
102
103
#ifndef LUCONTAINER_H
#define LUCONTAINER_H

#include <vector>

#include "matrix.h"
#include "sparsematrix.h"
#include "decompositor.h"

struct LUPS
{
    SparseMatrix U;
    std::vector<size_t> LF;
    std::vector<size_t> P;
    std::vector<size_t> Pt;

    void printL() const
    {
        for (size_t i = 1; i <= U.H; ++i)
        {
            size_t j = 1, jj;
            for (size_t q = LF[i]; q != SPARSE_END; q = U.N[q])
            {
                jj = U.C[q];
                while (j < jj) { std::cout << "0" << '\t'; ++j; }
                std::cout << U.V[q] << '\t'; ++j;
            }
            while (j <= U.W) { std::cout << "0" << '\t'; ++j; }
            std::cout << ";\n";
        }
    }

    void printU() const
    {
        U.print();
    }

    void save2fileL(const char *file)
    {
        std::ofstream out(file);

        out << U.H << " " << U.W << '\n';

        for (size_t i = 1; i <= U.H; ++i)
        {
            for (size_t q = LF[i]; q != SPARSE_END; q = U.N[q])
            {
                out << U.C[q] << " " << U.V[q] << '\t';
            }
            out << "0\n";
        }
        out.close();
    }

    void save2fileU(const char *file)
    {
        U.save2file(file);
    }
};

struct LUPM
{
    Matrix M;
    std::vector<size_t> P;
    std::vector<size_t> Pt;

//    void printL() const
//    {
//        for (size_t i = 1; i <= M.H; ++i)
//        {
//            for (size_t j = 1; j < i; ++j)
//                std::cout << M(i,j) << '\t';
//            for (size_t j = i; j <= M.H; ++j)
//                std::cout << "*\t";
//            std::cout << '\n';
//        }
//    }

//    void printU() const
//    {
//        for (size_t i = 1; i <= M.H; ++i)
//        {
//            for (size_t j = 1; j < i; ++j)
//                std::cout << "*\t";
//            for (size_t j = i; j <= M.H; ++j)
//                std::cout << M(i,j) << '\t';
//            std::cout << '\n';
//        }
//    }
};

struct FactorizedBBDSparseMatrix
{
    std::vector< LUPS > Alu;
    std::vector< Matrix > Bh;
    std::vector< SparseMatrix > C;
    LUPM H;
    BBDStruct BBDS;

    size_t N, Nb;
};

#endif // LUCONTAINER_H