Firing-Cost/ParallelDP.jl at master · pforero/Firing-Cost · 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
104
105
106
107
108
109
110
111
112
113
114
115
116
module ParallelDP

export LiLarge!, FeasLoopP, PossLoopP, Rows!

function LiLarge!(Li::BitArray{1},L1::Array{Int64,1},L2::Array{Int64,1})

    """Return a boolean vector saying if the change of
    portfolio is feasible (ie if tomorrow number of
    workers is larger than todays

    Parameters
    ----------
    Li : Ready Boolean Vector
    L1 : Todays Labour Portfolio
    L2 : Vector with all possible portfolios

    Output
    ----------
    Li : Filled Boolean Vector
    """

    for i=1:length(Li)

        Li[i]=L1[rem(i-1,length(L1))+1]>=L2[i];

    end

end

function FeasLoopP(L::Array{Any,1},L2::Array{Int64,1})

    """Find the indexes of the Feasible Portfolios
    for each current Portfolio

    Parameters
    ----------
    L : Labour Portfolio
    L2 : Vector with all possible portfolios

    Output
    ----------
    C : List with all Feasible Portfolio indexes for each Portfolio
    """

    C=Array(Array{Int64,1},length(L));

    Li=BitArray(length(L2));

    for i=1:length(L)

        LiLarge!(Li,vcat(L[i][1:end-2],sum(L[i][end-1:end])),L2);

        C[i]=find(prod(reshape(Li,length(L[1])-1,div(length(Li),length(L[1])-1)),1));

    end

    return C

end

function PossLoopP(L::Array{Any,1},L2::Array{Int64,1})

    """Find the indexes of the Possible Portfolio QnR
    for each current Portfolio

    Parameters
    ----------
    L : Labour Portfolio
    L2 : Vector with all possible portfolios

    Output
    ----------
    C : List with all Feasible Portfolio indexes for each Portfolio
    """

    C=Array(Array{Int64,1},length(L));

    Li=BitArray(length(L2));

    for i=1:length(L)

        LiLarge!(Li,L[i],L2);

        C[i]=find(prod(reshape(Li,length(L[1]),div(length(Li),length(L[1]))),1));

    end

    return C

end

function Rows!(R::Array{Array{Int64,1},1})

    """Give the Indexes for Todays Portfolio
    for each feasible Portfolio Change. It does
    it by filling the Feasible portfolio Matrix
    with its coloumn value

    Parameters
    ----------
    R : Empty Array of Indexes of Today Portfolio by the Number of Feasible Portfolios

    Output
    ----------
    R : Array of Indexes of Today Portfolio by the Number of Feasible Portfolios
    """

    for i=1:length(R)

        fill!(R[i],i);

    end

end

end