RubiscoSimulation/cornermove.hpp at master · binarybin/RubiscoSimulation · 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
//
//  cornermove.hpp
//  polymersim
//
//  Created by Bin Xu on 11/18/15.
//  Copyright © 2015 Bin Xu. All rights reserved.
//

#ifndef cornermove_hpp
#define cornermove_hpp

#include <tuple>
#include <vector>

using std::tuple;
using std::vector;
using std::make_tuple;

template <class S, class P> class CornerMove
{
private:
    S& space;
public:
    CornerMove(S& thespace) : space(thespace) {}

    void UpdatePolymer(Polymer<P>& poly, int pointid, P& newpoint)
    {
        poly.locs[pointid] = newpoint;
    }

    void UpdateReverseCheckingSpace(P& oldpoint, P& newpoint, Polymer<P>& poly)
    {
        space.RSpacePointMove(oldpoint, newpoint);
    }
    vector<tuple<int, int, P>> GetPossibleMoves(Polymer<P>& poly);

};

template <> vector<tuple<int, int, Pos2d2l>> CornerMove<Space2D2L, Pos2d2l>::GetPossibleMoves(Polymer<Pos2d2l>& poly)
{
    vector<tuple<int, int, Pos2d2l>> possible_moves;
    for (int pos = 0; pos < poly.locs.size() - 2; pos++)
    {
        int x1 = poly.locs[pos].x; int y1 = poly.locs[pos].y;
        int x2 = poly.locs[pos+1].x; int y2 = poly.locs[pos+1].y;
        int x3 = poly.locs[pos+2].x; int y3 = poly.locs[pos+2].y;
        bool siml = poly.locs[pos].siml;
        int layer = siml? 0:1;

        if (x1 != x3 && y1 != y3)
        {
            if (x1 == x2)
            {
                assert(y1 != y2 && y2 == y3 && x2 != x3);
                if (space.space[layer][x3][y1] == 0)
                    possible_moves.push_back(make_tuple(pos+1, pos+1, Pos2d2l(x3, y1, siml)));
            }
            else
            {
                assert(y1 == y2 && x2 == x3 && y2 != y3);
                if (space.space[layer][x1][y3] == 0)
                    possible_moves.push_back(make_tuple(pos+1, pos+1, Pos2d2l(x1, y3, siml)));
            }
        }
    }
    return possible_moves;
}


#endif /* cornermove_hpp */