program.cpp 4.19 KB
Newer Older
Paul's avatar
Paul committed
1
#include <rtg/program.hpp>
Paul's avatar
Paul committed
2
#include <rtg/stringutils.hpp>
Paul's avatar
Paul committed
3
#include <rtg/instruction.hpp>
Paul's avatar
Paul committed
4
#include <iostream>
Paul's avatar
Paul committed
5
6
7
8
#include <algorithm>

namespace rtg {

Paul's avatar
Paul committed
9
10
11
12
13
14
struct program_impl
{
    // A list is used to keep references to an instruction stable
    std::list<instruction> instructions;
};

Paul's avatar
Paul committed
15
program::program() : impl(std::make_unique<program_impl>()) {}
Paul's avatar
Paul committed
16

Paul's avatar
Paul committed
17
program::program(program&&) noexcept = default;
Paul's avatar
Paul committed
18
19
program& program::operator=(program&&) noexcept = default;
program::~program() noexcept                    = default;
Paul's avatar
Paul committed
20

Paul's avatar
Paul committed
21
instruction_ref program::add_instruction(operation op, std::vector<instruction_ref> args)
Paul's avatar
Paul committed
22
23
24
{
    return insert_instruction(impl->instructions.end(), std::move(op), std::move(args));
}
Paul's avatar
Paul committed
25
26
instruction_ref
program::insert_instruction(instruction_ref ins, operation op, std::vector<instruction_ref> args)
Paul's avatar
Paul committed
27
{
Paul's avatar
Paul committed
28
29
30
    assert(std::all_of(
               args.begin(), args.end(), [&](instruction_ref x) { return has_instruction(x); }) &&
           "Argument is not an exisiting instruction");
Paul's avatar
Paul committed
31
32
    std::vector<shape> shapes(args.size());
    std::transform(
Paul's avatar
Paul committed
33
        args.begin(), args.end(), shapes.begin(), [](instruction_ref i) { return i->result; });
Paul's avatar
Paul committed
34
    shape r     = op.compute_shape(shapes);
Paul's avatar
Paul committed
35
36
    auto result = impl->instructions.insert(ins, {op, r, args});
    assert(result->arguments == args);
Paul's avatar
Paul committed
37
38
    for(auto&& arg : args)
        arg->output.push_back(result);
Paul's avatar
Paul committed
39
    return result;
Paul's avatar
Paul committed
40
41
}

Paul's avatar
Paul committed
42
instruction_ref program::add_literal(literal l)
Paul's avatar
Paul committed
43
44
{
    impl->instructions.emplace_back(std::move(l));
Paul's avatar
Paul committed
45
    return std::prev(impl->instructions.end());
Paul's avatar
Paul committed
46
47
}

Paul's avatar
Paul committed
48
instruction_ref program::add_parameter(std::string name, shape s)
Paul's avatar
Paul committed
49
50
{
    impl->instructions.push_back({builtin::param{std::move(name)}, s, {}});
Paul's avatar
Paul committed
51
    return std::prev(impl->instructions.end());
Paul's avatar
Paul committed
52
53
}

Paul's avatar
Paul committed
54
bool program::has_instruction(instruction_ref ins) const
Paul's avatar
Paul committed
55
{
Paul's avatar
Paul committed
56
57
58
59
    return std::find_if(
               impl->instructions.begin(), impl->instructions.end(), [&](const instruction& x) {
                   return std::addressof(*ins) == std::addressof(x);
               }) != impl->instructions.end();
Paul's avatar
Paul committed
60
61
}

62
63
64
65
66
67
68
69
70
instruction_ref program::begin()
{
    return impl->instructions.begin();
}
instruction_ref program::end()
{
    return impl->instructions.end();
}

Paul's avatar
Paul committed
71
literal program::eval(std::unordered_map<std::string, argument> params) const
Paul's avatar
Paul committed
72
73
74
{
    std::unordered_map<const instruction*, argument> results;
    argument result;
Paul's avatar
Paul committed
75
    for(auto& ins : impl->instructions)
Paul's avatar
Paul committed
76
    {
Paul's avatar
Paul committed
77
        if(ins.op.name() == "@literal")
Paul's avatar
Paul committed
78
79
80
        {
            result = ins.lit.get_argument();
        }
Paul's avatar
Paul committed
81
        else if(starts_with(ins.op.name(), "@param"))
Paul's avatar
Paul committed
82
        {
Paul's avatar
Paul committed
83
            result = params.at(ins.op.name().substr(7));
Paul's avatar
Paul committed
84
        }
Paul's avatar
Paul committed
85
86
87
        else
        {
            std::vector<argument> values(ins.arguments.size());
Paul's avatar
Paul committed
88
89
90
            std::transform(ins.arguments.begin(),
                           ins.arguments.end(),
                           values.begin(),
Paul's avatar
Paul committed
91
                           [&](instruction_ref i) { return results.at(std::addressof(*i)); });
Paul's avatar
Paul committed
92
            result = ins.op.compute(values);
Paul's avatar
Paul committed
93
94
95
        }
        results.emplace(std::addressof(ins), result);
    }
Paul's avatar
Paul committed
96
    return literal{result.get_shape(), result.data()};
Paul's avatar
Paul committed
97
98
}

Paul's avatar
Paul committed
99
std::ostream& operator<<(std::ostream& os, const program& p)
Paul's avatar
Paul committed
100
101
102
103
{
    std::unordered_map<const instruction*, std::string> names;
    int count = 0;

Paul's avatar
Paul committed
104
    for(auto& ins : p.impl->instructions)
Paul's avatar
Paul committed
105
106
107
108
109
110
111
    {
        std::string var_name = "@" + std::to_string(count);
        if(starts_with(ins.op.name(), "@param"))
        {
            var_name = ins.op.name().substr(7);
        }

Paul's avatar
Paul committed
112
        os << var_name << " = ";
Paul's avatar
Paul committed
113

Paul's avatar
Paul committed
114
        os << ins.op.name();
Paul's avatar
Paul committed
115
116
117

        if(ins.op.name() == "@literal")
        {
Paul's avatar
Paul committed
118
            if(ins.lit.get_shape().elements() > 10)
Paul's avatar
Paul committed
119
                os << "{ ... }";
Paul's avatar
Paul committed
120
            else
Paul's avatar
Paul committed
121
                os << "{" << ins.lit << "}";
Paul's avatar
Paul committed
122
123
124
125
126
        }

        if(!ins.arguments.empty())
        {
            char delim = '(';
Paul's avatar
Paul committed
127
            for(auto&& arg : ins.arguments)
Paul's avatar
Paul committed
128
            {
Paul's avatar
Paul committed
129
130
                assert(p.has_instruction(arg) && "Instruction not found");
                os << delim << names.at(std::addressof(*arg));
Paul's avatar
Paul committed
131
132
                delim = ',';
            }
Paul's avatar
Paul committed
133
            os << ")";
Paul's avatar
Paul committed
134
135
        }

Paul's avatar
Paul committed
136
        os << " -> " << ins.result;
Paul's avatar
Paul committed
137

Paul's avatar
Paul committed
138
        os << std::endl;
Paul's avatar
Paul committed
139
140

        names.emplace(std::addressof(ins), var_name);
Paul's avatar
Paul committed
141
        count++;
Paul's avatar
Paul committed
142
    }
Paul's avatar
Paul committed
143
    return os;
Paul's avatar
Paul committed
144
145
}

Paul's avatar
Paul committed
146
} // namespace rtg