operators.hpp 17.3 KB
Newer Older
Paul's avatar
Paul committed
1
2
#ifndef MIGRAPH_GUARD_OPERATORS_HPP
#define MIGRAPH_GUARD_OPERATORS_HPP
Paul's avatar
Paul committed
3

4
#include <array>
Paul's avatar
Paul committed
5
#include <migraph/operation.hpp>
Paul's avatar
Paul committed
6
#include <migraph/check_shapes.hpp>
Paul's avatar
Paul committed
7
8
#include <migraph/stringutils.hpp>
#include <migraph/streamutils.hpp>
Paul's avatar
Paul committed
9
#include <cmath>
Paul's avatar
Paul committed
10
#include <utility>
Paul's avatar
Paul committed
11

Paul's avatar
Paul committed
12
namespace migraph {
Paul's avatar
Paul committed
13

Paul's avatar
Paul committed
14
15
struct not_computable
{
Paul's avatar
Paul committed
16
    argument compute(context&, const shape&, const std::vector<argument>&) const
Paul's avatar
Paul committed
17
18
19
    {
        MIGRAPH_THROW("not computable");
    }
Paul's avatar
Paul committed
20
21
};

22
23
struct batch_norm_inference
{
24
25
    float epsilon  = 1.0e-6f;
    float momentum = 0.9f;
26
27
28

    std::string name() const { return "batch_norm_inference"; }

29
30
31
32
33
34
35
36
    enum bn_infer_mode_t
    {
        per_activation,
        spatial,
    };

    bn_infer_mode_t bn_mode = spatial;

37
38
    bool is_test = false;

39
40
41
42
43
44
    shape compute_shape(std::vector<shape> inputs) const
    {
        check_shapes{inputs, *this}.has(5);
        return inputs.front();
    }

Paul's avatar
Paul committed
45
    argument compute(context&, const shape&, const std::vector<argument>&) const
46
47
48
49
50
    {
        MIGRAPH_THROW("not computable");
    }
};

Paul's avatar
Paul committed
51
struct convolution
Paul's avatar
Paul committed
52
{
Paul's avatar
Paul committed
53
54
55
    std::array<std::size_t, 2> padding  = {{0, 0}};
    std::array<std::size_t, 2> stride   = {{1, 1}};
    std::array<std::size_t, 2> dilation = {{1, 1}};
Paul's avatar
Paul committed
56
57
58
59
60
61
62
    enum padding_mode_t
    {
        default_, // NOLINT
        same,
        valid
    };
    padding_mode_t padding_mode = default_;
Paul's avatar
Paul committed
63
    std::string name() const { return "convolution"; }
Paul's avatar
Paul committed
64
65
    shape compute_shape(std::vector<shape> inputs) const
    {
Paul's avatar
Paul committed
66
        check_shapes{inputs, *this}.has(2).same_type().same_ndims().only_dims(4);
Paul's avatar
Paul committed
67

Paul's avatar
Paul committed
68
        const shape& input   = inputs.at(0);
Paul's avatar
Paul committed
69
        const shape& weights = inputs.at(1);
Paul's avatar
Paul committed
70
        auto t               = input.type();
Paul's avatar
Paul committed
71
72
        if(padding_mode == default_)
        {
Paul's avatar
Paul committed
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
            return {t,
                    {
                        input.lens()[0],
                        weights.lens()[0],
                        std::size_t(std::max<std::ptrdiff_t>(
                            1,
                            (input.lens()[2] - (1 + dilation[0] * (weights.lens()[2] - 1)) +
                             2 * padding[0]) /
                                    stride[0] +
                                1)),
                        std::size_t(std::max<std::ptrdiff_t>(
                            1,
                            (input.lens()[3] - (1 + dilation[1] * (weights.lens()[3] - 1)) +
                             2 * padding[1]) /
                                    stride[1] +
                                1)),
                    }};
Paul's avatar
Paul committed
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
        }
        else if(padding_mode == same)
        {
            return {t,
                    {input.lens()[0],
                     weights.lens()[0],
                     static_cast<std::size_t>(
                         std::ceil(static_cast<double>(input.lens()[2]) / stride[0])),
                     static_cast<std::size_t>(
                         std::ceil(static_cast<double>(input.lens()[3]) / stride[1]))}};
        }
        else if(padding_mode == valid)
        {
            return {
                t,
                {input.lens()[0],
                 weights.lens()[0],
                 static_cast<std::size_t>(std::ceil(
                     static_cast<double>(input.lens()[2] - weights.lens()[2] + 1) / stride[0])),
                 static_cast<std::size_t>(std::ceil(
                     static_cast<double>(input.lens()[3] - weights.lens()[3] + 1) / stride[1]))}};
        }
        else
        {
Paul's avatar
Paul committed
114
            MIGRAPH_THROW("Invalid padding mode");
Paul's avatar
Paul committed
115
        }
Paul's avatar
Paul committed
116
    }
Paul's avatar
Paul committed
117

Paul's avatar
Paul committed
118
    argument compute(context&, const shape&, const std::vector<argument>&) const
Paul's avatar
Paul committed
119
120
121
    {
        MIGRAPH_THROW("not computable");
    }
Paul's avatar
Paul committed
122

Paul's avatar
Paul committed
123
    friend std::ostream& operator<<(std::ostream& os, const convolution& op)
Paul's avatar
Paul committed
124
    {
Paul's avatar
Paul committed
125
126
127
128
129
        os << op.name() << "[";
        os << "padding={" << stream_range(op.padding) << "}, ";
        os << "stride={" << stream_range(op.stride) << "}, ";
        os << "dilation={" << stream_range(op.dilation) << "}";
        os << "]";
Paul's avatar
Paul committed
130
131
        return os;
    }
Paul's avatar
Paul committed
132
133
};

Scott Thornton's avatar
Scott Thornton committed
134
135
struct im2col
{
Scott Thornton's avatar
Scott Thornton committed
136
137
138
139
140
141
142
143
144
145
146
147
148
149
    std::array<std::size_t, 2> padding  = {{0, 0}};
    std::array<std::size_t, 2> stride   = {{1, 1}};
    std::array<std::size_t, 2> dilation = {{1, 1}};
    enum padding_mode_t
    {
        default_, // NOLINT
        same,
        valid
    };

    std::string name() const { return "im2col"; }

    shape compute_shape(std::vector<shape> inputs) const
    {
Scott Thornton's avatar
Scott Thornton committed
150
151
152
        auto input          = inputs[0];
        auto weights        = inputs[1];
        auto batch_size     = input.lens()[0];
Scott Thornton's avatar
Scott Thornton committed
153
        auto input_channels = weights.lens()[1];
Scott Thornton's avatar
Scott Thornton committed
154
155
        auto kernel_height  = weights.lens()[2];
        auto kernel_width   = weights.lens()[3];
Scott Thornton's avatar
Scott Thornton committed
156
        check_shapes{inputs, *this}.has(2);
Scott Thornton's avatar
Scott Thornton committed
157
158
        if(batch_size != 1)
            MIGRAPH_THROW("im2col only support batch_size 1");
Scott Thornton's avatar
Scott Thornton committed
159
        auto output_height = std::size_t(std::max<std::ptrdiff_t>(
Scott Thornton's avatar
Scott Thornton committed
160
161
162
            1,
            (input.lens()[2] - (1 + dilation[0] * (kernel_height - 1)) + 2 * padding[0]) /
                    stride[0] +
Scott Thornton's avatar
Scott Thornton committed
163
                1));
Scott Thornton's avatar
Scott Thornton committed
164
165
166
167
        auto output_width  = std::size_t(std::max<std::ptrdiff_t>(
            1,
            (input.lens()[3] - (1 + dilation[1] * (kernel_width - 1)) + 2 * padding[1]) /
                    stride[1] +
Scott Thornton's avatar
Scott Thornton committed
168
                1));
Scott Thornton's avatar
Scott Thornton committed
169
170
        auto channels_col  = kernel_height * kernel_width * input_channels;
        return {input.type(), {output_height * output_width, channels_col}};
Scott Thornton's avatar
Scott Thornton committed
171
172
173
174
175
176
177
178
    }

    argument compute(context&, const shape&, const std::vector<argument>&) const
    {
        MIGRAPH_THROW("not computable");
    }
};

Paul's avatar
Paul committed
179
struct pooling
Paul's avatar
Paul committed
180
{
Paul's avatar
Paul committed
181
    std::string mode                   = "average";
Paul's avatar
Paul committed
182
183
184
    std::array<std::size_t, 2> padding = {{0, 0}};
    std::array<std::size_t, 2> stride  = {{1, 1}};
    std::array<std::size_t, 2> lengths = {{1, 1}};
Paul's avatar
Paul committed
185
    std::string name() const { return "pooling"; }
Scott Thornton's avatar
Scott Thornton committed
186

Paul's avatar
Paul committed
187
188
    shape compute_shape(std::vector<shape> inputs) const
    {
Paul's avatar
Paul committed
189
        check_shapes{inputs, *this}.has(1).only_dims(4);
Paul's avatar
Paul committed
190

Paul's avatar
Paul committed
191
        const shape& input = inputs.at(0);
Paul's avatar
Paul committed
192
        auto t             = input.type();
Paul's avatar
Paul committed
193

Paul's avatar
Paul committed
194
195
        assert(lengths[0] <= (input.lens()[2] + 2 * padding[0]));
        assert(lengths[1] <= (input.lens()[3] + 2 * padding[1]));
Paul's avatar
Paul committed
196

Scott Thornton's avatar
Scott Thornton committed
197
198
199
200
201
202
        return {t,
                {
                    input.lens()[0],
                    input.lens()[1],
                    std::size_t(std::max<std::ptrdiff_t>(
                        1,
Paul's avatar
Paul committed
203
                        std::ptrdiff_t(std::floor((input.lens()[2] + 2 * padding[0] - lengths[0]) /
Paul's avatar
Paul committed
204
                                                  static_cast<float>(stride[0]))) +
Scott Thornton's avatar
Scott Thornton committed
205
206
207
                            1)),
                    std::size_t(std::max<std::ptrdiff_t>(
                        1,
Paul's avatar
Paul committed
208
                        std::ptrdiff_t(std::floor((input.lens()[3] + 2 * padding[1] - lengths[1]) /
Paul's avatar
Paul committed
209
                                                  static_cast<float>(stride[1]))) +
Scott Thornton's avatar
Scott Thornton committed
210
211
                            1)),
                }};
Paul's avatar
Paul committed
212
    }
Paul's avatar
Paul committed
213

Paul's avatar
Paul committed
214
    argument compute(context&, const shape&, const std::vector<argument>&) const
Paul's avatar
Paul committed
215
216
217
    {
        MIGRAPH_THROW("not computable");
    }
Paul's avatar
Paul committed
218

Paul's avatar
Paul committed
219
    friend std::ostream& operator<<(std::ostream& os, const pooling& op)
Paul's avatar
Paul committed
220
    {
Paul's avatar
Paul committed
221
222
223
224
225
        os << op.name() << "[";
        os << "padding={" << stream_range(op.padding) << "}, ";
        os << "stride={" << stream_range(op.stride) << "}, ";
        os << "lengths={" << stream_range(op.lengths) << "}";
        os << "]";
Paul's avatar
Paul committed
226
227
        return os;
    }
Paul's avatar
Paul committed
228
229
};

Paul's avatar
Paul committed
230
struct activation
Paul's avatar
Paul committed
231
232
{
    std::string mode;
Paul's avatar
Paul committed
233
    std::string name() const { return "activation"; }
Paul's avatar
Paul committed
234
235
    shape compute_shape(std::vector<shape> inputs) const
    {
Paul's avatar
Paul committed
236
        check_shapes{inputs, *this}.has(1);
Paul's avatar
Paul committed
237
238
        return inputs.front();
    }
Paul's avatar
Paul committed
239

Paul's avatar
Paul committed
240
    argument compute(context&, const shape&, const std::vector<argument>&) const
Paul's avatar
Paul committed
241
242
243
    {
        MIGRAPH_THROW("not computable");
    }
Paul's avatar
Paul committed
244
    friend std::ostream& operator<<(std::ostream& os, const activation& op)
Paul's avatar
Paul committed
245
    {
Paul's avatar
Paul committed
246
        os << op.name() << ":" << op.mode;
Paul's avatar
Paul committed
247
248
        return os;
    }
Paul's avatar
Paul committed
249
250
};

251
252
253
254
255
256
struct transpose
{
    std::vector<int64_t> dims;
    std::string name() const { return "transpose"; }
    shape compute_shape(std::vector<shape> inputs) const
    {
Paul's avatar
Paul committed
257
        check_shapes{inputs, *this}.has(1);
258
        auto input         = inputs.at(0);
259
        auto input_lens    = input.lens();
260
261
        auto input_strides = input.strides();
        auto t             = input.type();
Paul's avatar
Paul committed
262
263
        if(dims.size() != input_lens.size())
        {
Paul's avatar
Paul committed
264
            MIGRAPH_THROW("Permutation has wrong number of axes");
265
266
267
        }
        std::vector<int64_t> axes(dims.size());
        std::iota(axes.begin(), axes.end(), 0);
Paul's avatar
Paul committed
268
269
        if(!std::is_permutation(axes.begin(), axes.end(), dims.begin()))
        {
Paul's avatar
Paul committed
270
            MIGRAPH_THROW("Invalid permutation");
271
        }
272
273
        std::vector<size_t> output_lens(input_lens.size());
        std::vector<size_t> output_strides(input_lens.size());
Paul's avatar
Paul committed
274
275
276
        for(int i = 0; i < output_lens.size(); i++)
        {
            output_lens[i]    = input_lens[dims[i]];
277
278
            output_strides[i] = input_strides[dims[i]];
        }
279
        return {t, output_lens, output_strides};
280
    }
Paul's avatar
Paul committed
281
    argument compute(context&, shape output_shape, std::vector<argument> args) const
Paul's avatar
Paul committed
282
    {
Paul's avatar
Paul committed
283
        return {std::move(output_shape), std::move(args.front().data)};
Paul's avatar
Paul committed
284
    }
Paul's avatar
Paul committed
285
286
287
288
289
290
291
    friend std::ostream& operator<<(std::ostream& os, const transpose& op)
    {
        os << op.name() << "[";
        os << "dims={" << stream_range(op.dims) << "}";
        os << "]";
        return os;
    }
292
293
};

Paul's avatar
Paul committed
294
struct contiguous
295
296
297
298
{
    std::string name() const { return "contiguous"; }
    shape compute_shape(std::vector<shape> inputs) const
    {
Paul's avatar
Paul committed
299
        check_shapes{inputs, *this}.has(1);
Paul's avatar
Paul committed
300
301
302
303
        auto lens = inputs.at(0).lens();
        auto t    = inputs.at(0).type();
        if(lens.size() < 2)
        {
Paul's avatar
Paul committed
304
            MIGRAPH_THROW("Number of dimensions should exceed 1");
305
306
307
        }
        return {t, lens};
    }
Paul's avatar
Paul committed
308
    argument compute(context&, const shape&, const std::vector<argument>&) const
Paul's avatar
Paul committed
309
310
311
    {
        MIGRAPH_THROW("not computable");
    }
312
313
};

Paul's avatar
Paul committed
314
315
316
struct reshape
{
    std::vector<int64_t> dims;
Paul's avatar
Paul committed
317
    std::string name() const { return "reshape"; }
Paul's avatar
Paul committed
318
319
    shape compute_shape(std::vector<shape> inputs) const
    {
Paul's avatar
Paul committed
320
        check_shapes{inputs, *this}.has(1);
Paul's avatar
Paul committed
321
322
        auto&& idims = inputs.front().lens();
        std::vector<std::size_t> rdims(dims.begin(), dims.end());
323
324
325
        auto n_neg_dims = std::count(dims.begin(), dims.end(), -1);
        if(n_neg_dims > 1)
            MIGRAPH_THROW("Dimensions for reshape can only have one -1 dim");
Paul's avatar
Paul committed
326
        for(std::size_t i = 0; i < dims.size(); i++)
Paul's avatar
Paul committed
327
328
329
330
        {
            if(dims[i] == 0)
                rdims[i] = idims[i];
        }
331
332
333
334
335
336
337
338
339
340
341
        if(n_neg_dims > 0)
        {
            size_t missing_dim =
                -inputs.front().elements() /
                std::accumulate(rdims.begin(), rdims.end(), 1, std::multiplies<int64_t>());
            for(std::size_t i = 0; i < rdims.size(); i++)
            {
                if(dims[i] == -1)
                    rdims[i] = missing_dim;
            }
        }
Paul's avatar
Paul committed
342
343
344
        if(dims.back() == -1)
        {
            rdims.pop_back();
Paul's avatar
Paul committed
345
            std::copy(idims.begin() + rdims.size(), idims.end(), std::back_inserter(rdims));
Paul's avatar
Paul committed
346
        }
Scott Thornton's avatar
Scott Thornton committed
347
        shape s{inputs.front().type(), rdims};
Paul's avatar
Paul committed
348
        if(s.elements() != inputs.front().elements())
Paul's avatar
Paul committed
349
            MIGRAPH_THROW("Wrong number of elements for reshape");
Scott Thornton's avatar
Scott Thornton committed
350
        return s;
Paul's avatar
Paul committed
351
352
    }

Paul's avatar
Paul committed
353
    argument compute(context&, shape output_shape, std::vector<argument> args) const
Paul's avatar
Paul committed
354
    {
Paul's avatar
Paul committed
355
        return {std::move(output_shape), std::move(args.front().data)};
Paul's avatar
Paul committed
356
    }
Paul's avatar
Paul committed
357

Paul's avatar
Paul committed
358
    friend std::ostream& operator<<(std::ostream& os, const reshape& op)
Paul's avatar
Paul committed
359
    {
Paul's avatar
Paul committed
360
        os << op.name() << "[";
Paul's avatar
Paul committed
361
        os << "dims={" << stream_range(op.dims) << "}";
Paul's avatar
Paul committed
362
        os << "]";
Paul's avatar
Paul committed
363
364
        return os;
    }
Paul's avatar
Paul committed
365
366
};

367
368
struct gemm
{
Paul's avatar
Paul committed
369
    float alpha = 1.0;
Paul's avatar
Paul committed
370
    float beta  = 0.0;
371
    std::string name() const { return "gemm"; }
372
373
    shape compute_shape(std::vector<shape> inputs) const
    {
Paul's avatar
Paul committed
374
        check_shapes{inputs, *this}.has(2).same_type();
375
376
        const shape& a = inputs.at(0);
        const shape& b = inputs.at(1);
Scott Thornton's avatar
Scott Thornton committed
377
        auto t         = a.type();
378

379
        if(a.lens()[1] != b.lens()[0])
Paul's avatar
Paul committed
380
381
            MIGRAPH_THROW("Inner dimensions do not match: {" + to_string_range(a.lens()) + "} x {" +
                          to_string_range(b.lens()) + "}");
Scott Thornton's avatar
Scott Thornton committed
382
        return {t, {a.lens()[0], b.lens()[1]}};
383
    }
384

Paul's avatar
Paul committed
385
    argument compute(context&, const shape&, const std::vector<argument>&) const
Paul's avatar
Paul committed
386
387
388
    {
        MIGRAPH_THROW("not computable");
    }
389
390

    friend std::ostream& operator<<(std::ostream& os, const gemm& op)
391
392
    {
        os << op.name() << "[";
393
        os << "]";
Scott Thornton's avatar
Scott Thornton committed
394
        return os;
395
396
397
    }
};

398
struct unary
Scott Thornton's avatar
Scott Thornton committed
399
{
400
401
    shape compute_shape(std::vector<shape> inputs) const
    {
402
403
        check_shapes{inputs}.has(1);
        return inputs.at(0);
404
    }
Paul's avatar
Paul committed
405
    argument compute(context&, const shape&, const std::vector<argument>&) const
Paul's avatar
Paul committed
406
407
408
    {
        MIGRAPH_THROW("not computable");
    }
Scott Thornton's avatar
Scott Thornton committed
409
410
};

411
412
struct identity : unary
{
413
    std::string name() const { return "identity"; }
414
415
416
};

struct abs : unary
Scott Thornton's avatar
Scott Thornton committed
417
{
418
    std::string name() const { return "abs"; }
Scott Thornton's avatar
Scott Thornton committed
419
420
};

421
struct exp : unary
Scott Thornton's avatar
Scott Thornton committed
422
{
423
    std::string name() const { return "exp"; }
Scott Thornton's avatar
Scott Thornton committed
424
425
};

426
struct sin : unary
Scott Thornton's avatar
Scott Thornton committed
427
{
428
    std::string name() const { return "sin"; }
Scott Thornton's avatar
Scott Thornton committed
429
430
};

431
struct cos : unary
Scott Thornton's avatar
Scott Thornton committed
432
{
433
    std::string name() const { return "cos"; }
Scott Thornton's avatar
Scott Thornton committed
434
435
};

436
struct tan : unary
Scott Thornton's avatar
Scott Thornton committed
437
{
438
    std::string name() const { return "tan"; }
Scott Thornton's avatar
Scott Thornton committed
439
440
};

441
struct asin : unary
Scott Thornton's avatar
Scott Thornton committed
442
{
443
    std::string name() const { return "asin"; }
Scott Thornton's avatar
Scott Thornton committed
444
445
};

446
struct acos : unary
Scott Thornton's avatar
Scott Thornton committed
447
{
448
    std::string name() const { return "acos"; }
Scott Thornton's avatar
Scott Thornton committed
449
450
};

451
struct atan : unary
Scott Thornton's avatar
Scott Thornton committed
452
{
453
    std::string name() const { return "atan"; }
Scott Thornton's avatar
Scott Thornton committed
454
455
};

456
struct softmax : unary
Scott Thornton's avatar
Scott Thornton committed
457
{
458
    std::string name() const { return "softmax"; }
Scott Thornton's avatar
Scott Thornton committed
459
460
};

461
struct tanh : unary
Scott Thornton's avatar
Scott Thornton committed
462
{
463
    std::string name() const { return "tanh"; }
Scott Thornton's avatar
Scott Thornton committed
464
465
};

466
struct sigmoid : unary
Scott Thornton's avatar
Scott Thornton committed
467
{
468
    std::string name() const { return "sigmoid"; }
Scott Thornton's avatar
Scott Thornton committed
469
470
};

471
struct neg : unary
Scott Thornton's avatar
Scott Thornton committed
472
{
473
    std::string name() const { return "neg"; }
Scott Thornton's avatar
Scott Thornton committed
474
475
};

476
struct flatten
Scott Thornton's avatar
Scott Thornton committed
477
{
Paul's avatar
Paul committed
478
    uint64_t axis = 0;
Scott Thornton's avatar
Scott Thornton committed
479
    std::string name() const { return "flatten"; }
Paul's avatar
Paul committed
480
481
482
    shape compute_shape(std::vector<shape> inputs) const
    {
        check_shapes{inputs}.has(1);
Paul's avatar
Paul committed
483
484
        auto&& lens = inputs.front().lens();

Paul's avatar
Paul committed
485
        if(axis > lens.size())
Paul's avatar
Paul committed
486
        {
Paul's avatar
Paul committed
487
            MIGRAPH_THROW("axis for flatten must be less than tensor rank");
Paul's avatar
Paul committed
488
        }
Paul's avatar
Paul committed
489
490
491
492
        auto x =
            std::accumulate(lens.begin(), lens.begin() + axis, std::size_t{1}, std::multiplies<>{});
        auto y =
            std::accumulate(lens.begin() + axis, lens.end(), std::size_t{1}, std::multiplies<>{});
493
        return {inputs.at(0).type(), {x, y}};
Paul's avatar
Paul committed
494
495
496
    }
    argument compute(context&, shape output_shape, std::vector<argument> args) const
    {
Paul's avatar
Paul committed
497
        return {std::move(output_shape), std::move(args.front().data)};
Paul's avatar
Paul committed
498
    }
Paul's avatar
Paul committed
499
500
501
502
503
504
505
    friend std::ostream& operator<<(std::ostream& os, const flatten& op)
    {
        os << op.name() << "[";
        os << "axis=" << op.axis;
        os << "]";
        return os;
    }
Scott Thornton's avatar
Scott Thornton committed
506
};
507
508
509
510
511
512
struct broadcast
{
    uint64_t axis = 0;
    std::string name() const { return "broadcast"; }
    shape compute_shape(std::vector<shape> inputs) const
    {
Paul's avatar
Paul committed
513
514
515
516
        auto t      = inputs.at(0).type();
        auto result = inputs.at(0);
        auto input  = inputs.at(1);

Paul's avatar
Paul committed
517
        std::vector<size_t> bcast_strides(result.lens().size(), 0);
518

Paul's avatar
Paul committed
519
520
        if(std::all_of(
               result.lens().cbegin(), result.lens().cend(), [&](auto x) { return x == 1; }))
521
        {
Scott Thornton's avatar
Scott Thornton committed
522
            if(axis != 0)
Paul's avatar
Paul committed
523
                MIGRAPH_THROW("when broadcasting tensor of size 1, axis should be 0");
Paul's avatar
Paul committed
524
            return {t, result.lens(), std::move(bcast_strides)};
525
526
527
        }
        else
        {
Paul's avatar
Paul committed
528
529
            assert(result.lens().size() - axis >= input.lens().size());
            if(!std::equal(input.lens().begin(), input.lens().end(), result.lens().begin() + axis))
Paul's avatar
Paul committed
530
                MIGRAPH_THROW("when broadcasting success sizes must match");
Paul's avatar
Paul committed
531
            std::copy(input.strides().begin(), input.strides().end(), bcast_strides.begin() + axis);
Paul's avatar
Paul committed
532
            return {t, result.lens(), std::move(bcast_strides)};
533
534
        }
    }
Paul's avatar
Paul committed
535
    argument compute(context&, shape output_shape, std::vector<argument> args) const
Scott Thornton's avatar
Scott Thornton committed
536
    {
Paul's avatar
Paul committed
537
        return {std::move(output_shape), std::move(args.at(1).data)};
Scott Thornton's avatar
Scott Thornton committed
538
    }
Paul's avatar
Paul committed
539
540
541
542
543
544
545
    friend std::ostream& operator<<(std::ostream& os, const broadcast& op)
    {
        os << op.name() << "[";
        os << "axis=" << op.axis;
        os << "]";
        return os;
    }
546
547
};

548
struct binary
Scott Thornton's avatar
Scott Thornton committed
549
{
550
    uint64_t broadcast = 0;
551
552
    shape compute_shape(std::vector<shape> inputs) const
    {
553
554
        check_shapes{inputs}.has(2).same_type().same_dims();
        return inputs.at(0);
555
    }
Paul's avatar
Paul committed
556
    argument compute(context&, const shape&, const std::vector<argument>&) const
Paul's avatar
Paul committed
557
558
559
    {
        MIGRAPH_THROW("not computable");
    }
Scott Thornton's avatar
Scott Thornton committed
560
561
};

562
563
564
565
566
567
struct add : binary
{
    std::string name() const { return "add"; }
};

struct sub : binary
Scott Thornton's avatar
Scott Thornton committed
568
569
570
571
{
    std::string name() const { return "sub"; }
};

572
struct mul : binary
Scott Thornton's avatar
Scott Thornton committed
573
574
575
576
{
    std::string name() const { return "mul"; }
};

577
struct div : binary
Scott Thornton's avatar
Scott Thornton committed
578
579
580
581
{
    std::string name() const { return "div"; }
};

Paul's avatar
Paul committed
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
struct load
{
    shape s;
    std::size_t offset = 0;
    std::string name() const { return "load"; }
    shape compute_shape(const std::vector<shape>& inputs) const
    {
        check_shapes{inputs}.has(1);
        return s;
    }
    argument compute(context&, const shape&, const std::vector<argument>& args) const
    {
        return {s, args[0].data() + offset};
    }
};

Paul's avatar
Paul committed
598
struct outline
Scott Thornton's avatar
Scott Thornton committed
599
{
Paul's avatar
Paul committed
600
601
    shape s;
    std::string name() const { return "outline"; }
Paul's avatar
Paul committed
602
    shape compute_shape(const std::vector<shape>& inputs) const
Paul's avatar
Paul committed
603
    {
Paul's avatar
Paul committed
604
        check_shapes{inputs, *this}.has(0);
Paul's avatar
Paul committed
605
606
        return s;
    }
Paul's avatar
Paul committed
607
608
609
610
    argument compute(context&, const shape&, const std::vector<argument>&) const
    {
        return {s, nullptr};
    }
Scott Thornton's avatar
Scott Thornton committed
611
612
};

Paul's avatar
Paul committed
613
} // namespace migraph
Paul's avatar
Paul committed
614
615

#endif