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 CAFFE2_OPERATORS_MATMUL_OP_H_
#define CAFFE2_OPERATORS_MATMUL_OP_H_
#include "caffe2/core/context.h"
#include "caffe2/core/operator.h"
#include "caffe2/utils/math.h"
namespace caffe2 {
template <typename T, class Context, class Engine = DefaultEngine>
class MatMulOp final : public Operator<Context> {
public:
USE_OPERATOR_CONTEXT_FUNCTIONS;
MatMulOp(const OperatorDef& operator_def, Workspace* ws)
: Operator<Context>(operator_def, ws),
axis_a_(this->template GetSingleArgument<int>("axis_a", 1)),
axis_b_(this->template GetSingleArgument<int>("axis_b", 1)),
trans_a_(this->template GetSingleArgument<int>("trans_a", 0)),
trans_b_(this->template GetSingleArgument<int>("trans_b", 0)) {}
~MatMulOp() {}
bool RunOnDevice() override {
const auto& A = Input(0);
const auto& B = Input(1);
const auto canonical_axis_a = A.canonical_axis_index(axis_a_);
const auto canonical_axis_b = B.canonical_axis_index(axis_b_);
int A_dim0 = A.size_to_dim(canonical_axis_a);
int A_dim1 = A.size_from_dim(canonical_axis_a);
int B_dim0 = B.size_to_dim(canonical_axis_b);
int B_dim1 = B.size_from_dim(canonical_axis_b);
int a_dim0, a_dim1, b_dim0, b_dim1;
if (trans_a_) {
a_dim0 = A_dim1;
a_dim1 = A_dim0;
} else {
a_dim0 = A_dim0;
a_dim1 = A_dim1;
}
if (trans_b_) {
b_dim0 = B_dim1;
b_dim1 = B_dim0;
} else {
b_dim0 = B_dim0;
b_dim1 = B_dim1;
}
auto dimErrorString = [&]() {
return c10::str(
"Dimension mismatch: ",
trans_a_ ? "trans(A): " : "A: ",
a_dim0,
" ",
a_dim1,
trans_b_ ? ", trans(B): " : ", B: ",
b_dim0,
" ",
b_dim1);
};
// Error checking
CAFFE_ENFORCE(a_dim1 == b_dim0, dimErrorString());
Y_shape_cache_[0] = a_dim0;
Y_shape_cache_[1] = b_dim1;
auto* Y = Output(0, Y_shape_cache_, at::dtype<T>());
CAFFE_ENFORCE(a_dim0 * b_dim1 == Y->numel(), dimErrorString());
// Y = A * B
math::Gemm<T, Context, Engine>(
trans_a_ ? CblasTrans : CblasNoTrans,
trans_b_ ? CblasTrans : CblasNoTrans,
a_dim0,
b_dim1,
a_dim1,
1,
A.template data<T>(),
B.template data<T>(),
0,
Y->template mutable_data<T>(),
&context_);
if (InputSize() == 3) {
// In gradient op, resize to input
Y->ResizeLike(Input(2));
}
return true;
}
protected:
// A local vector to cache the output shape so we don't need to recreate
// a vector object every time we run Run().
vector<int64_t> Y_shape_cache_{0, 0};
int axis_a_{1};
int axis_b_{1};
bool trans_a_;
bool trans_b_;
};
} // namespace caffe2
#endif // CAFFE2_OPERATORS_MATMUL_OP_H_
|
