[CustomOp] Support complex dtype in custom op (#31657)
* support custom complex op * fix detail error * add inference support * fix setup windows failedtest_benchmark_ci
Chen Weihang
4 years ago
committed by
GitHub
parent
fe241fd02f
commit
87852616aa
@ -0,0 +1,94 @@
|
|||||||
|
// Copyright (c) 2021 PaddlePaddle Authors. All Rights Reserved.
|
||||||
|
//
|
||||||
|
// Licensed under the Apache License, Version 2.0 (the "License");
|
||||||
|
// you may not use this file except in compliance with the License.
|
||||||
|
// You may obtain a copy of the License at
|
||||||
|
//
|
||||||
|
// http://www.apache.org/licenses/LICENSE-2.0
|
||||||
|
//
|
||||||
|
// Unless required by applicable law or agreed to in writing, software
|
||||||
|
// distributed under the License is distributed on an "AS IS" BASIS,
|
||||||
|
// WIdata_tHOUdata_t WARRANdata_tIES OR CONDIdata_tIONS OF ANY KIND, either
|
||||||
|
// express or implied.
|
||||||
|
// See the License for the specific language governing permissions and
|
||||||
|
// limitations under the License.
|
||||||
|
|
||||||
|
#include <iostream>
|
||||||
|
#include <vector>
|
||||||
|
|
||||||
|
#include "paddle/extension.h"
|
||||||
|
|
||||||
|
#define CHECK_INPUT(x) \
|
||||||
|
PD_CHECK(x.place() == paddle::PlaceType::kCPU, #x " must be a CPU Tensor.")
|
||||||
|
|
||||||
|
template <typename data_t>
|
||||||
|
using EnableComplex = typename std::enable_if<
|
||||||
|
std::is_same<data_t, paddle::complex64>::value ||
|
||||||
|
std::is_same<data_t, paddle::complex128>::value>::type;
|
||||||
|
|
||||||
|
template <typename data_t>
|
||||||
|
using DisableComplex = typename std::enable_if<
|
||||||
|
!std::is_same<data_t, paddle::complex64>::value &&
|
||||||
|
!std::is_same<data_t, paddle::complex128>::value>::type;
|
||||||
|
|
||||||
|
template <typename data_t, typename Enable = void>
|
||||||
|
struct ConjFunctor;
|
||||||
|
|
||||||
|
template <typename data_t>
|
||||||
|
struct ConjFunctor<data_t, EnableComplex<data_t>> {
|
||||||
|
ConjFunctor(const data_t* input, int64_t numel, data_t* output)
|
||||||
|
: input_(input), numel_(numel), output_(output) {}
|
||||||
|
|
||||||
|
void operator()(size_t idx) const {
|
||||||
|
output_[idx] = data_t(input_[idx].real, -input_[idx].imag);
|
||||||
|
}
|
||||||
|
|
||||||
|
const data_t* input_;
|
||||||
|
int64_t numel_;
|
||||||
|
data_t* output_;
|
||||||
|
};
|
||||||
|
|
||||||
|
template <typename data_t>
|
||||||
|
struct ConjFunctor<data_t, DisableComplex<data_t>> {
|
||||||
|
ConjFunctor(const data_t* input, int64_t numel, data_t* output)
|
||||||
|
: input_(input), numel_(numel), output_(output) {}
|
||||||
|
|
||||||
|
void operator()(size_t idx) const { output_[idx] = input_[idx]; }
|
||||||
|
|
||||||
|
const data_t* input_;
|
||||||
|
int64_t numel_;
|
||||||
|
data_t* output_;
|
||||||
|
};
|
||||||
|
|
||||||
|
template <typename data_t>
|
||||||
|
void ConjCPUKernel(const data_t* x_data, int64_t numel, data_t* out_data) {
|
||||||
|
ConjFunctor<data_t> conj(x_data, numel, out_data);
|
||||||
|
for (int64_t i = 0; i < numel; ++i) {
|
||||||
|
conj(i);
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
std::vector<paddle::Tensor> ConjFunction(const paddle::Tensor& x) {
|
||||||
|
CHECK_INPUT(x);
|
||||||
|
|
||||||
|
paddle::Tensor out(x.place());
|
||||||
|
out.reshape(x.shape());
|
||||||
|
|
||||||
|
PD_DISPATCH_FLOATING_AND_COMPLEX_TYPES(
|
||||||
|
x.type(), "ConjCPUKernel", ([&] {
|
||||||
|
ConjCPUKernel<data_t>(
|
||||||
|
x.data<data_t>(), x.size(), out.mutable_data<data_t>());
|
||||||
|
}));
|
||||||
|
|
||||||
|
return {out};
|
||||||
|
}
|
||||||
|
|
||||||
|
PD_BUILD_OP(custom_conj)
|
||||||
|
.Inputs({"X"})
|
||||||
|
.Outputs({"Out"})
|
||||||
|
.SetKernelFn(PD_KERNEL(ConjFunction));
|
||||||
|
|
||||||
|
PD_BUILD_GRAD_OP(custom_conj)
|
||||||
|
.Inputs({paddle::Grad("Out")})
|
||||||
|
.Outputs({paddle::Grad("X")})
|
||||||
|
.SetKernelFn(PD_KERNEL(ConjFunction));
|
||||||
@ -0,0 +1,136 @@
|
|||||||
|
# Copyright (c) 2021 PaddlePaddle Authors. All Rights Reserved.
|
||||||
|
#
|
||||||
|
# Licensed under the Apache License, Version 2.0 (the "License");
|
||||||
|
# you may not use this file except in compliance with the License.
|
||||||
|
# You may obtain a copy of the License at
|
||||||
|
#
|
||||||
|
# http://www.apache.org/licenses/LICENSE-2.0
|
||||||
|
#
|
||||||
|
# Unless required by applicable law or agreed to in writing, software
|
||||||
|
# distributed under the License is distributed on an "AS IS" BASIS,
|
||||||
|
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
||||||
|
# See the License for the specific language governing permissions and
|
||||||
|
# limitations under the License.
|
||||||
|
|
||||||
|
import os
|
||||||
|
import unittest
|
||||||
|
import numpy as np
|
||||||
|
|
||||||
|
import paddle
|
||||||
|
import paddle.static as static
|
||||||
|
from paddle.utils.cpp_extension import load, get_build_directory
|
||||||
|
from paddle.utils.cpp_extension.extension_utils import run_cmd
|
||||||
|
from utils import paddle_includes, extra_cc_args, extra_nvcc_args
|
||||||
|
|
||||||
|
# Because Windows don't use docker, the shared lib already exists in the
|
||||||
|
# cache dir, it will not be compiled again unless the shared lib is removed.
|
||||||
|
file = '{}\\custom_relu_module_jit\\custom_relu_module_jit.pyd'.format(
|
||||||
|
get_build_directory())
|
||||||
|
if os.name == 'nt' and os.path.isfile(file):
|
||||||
|
cmd = 'del {}'.format(file)
|
||||||
|
run_cmd(cmd, True)
|
||||||
|
|
||||||
|
custom_ops = load(
|
||||||
|
name='custom_conj_jit',
|
||||||
|
sources=['custom_conj_op.cc'],
|
||||||
|
extra_include_paths=paddle_includes, # add for Coverage CI
|
||||||
|
extra_cxx_cflags=extra_cc_args, # test for cc flags
|
||||||
|
extra_cuda_cflags=extra_nvcc_args, # test for nvcc flags
|
||||||
|
verbose=True)
|
||||||
|
|
||||||
|
|
||||||
|
def is_complex(dtype):
|
||||||
|
return dtype == paddle.fluid.core.VarDesc.VarType.COMPLEX64 or \
|
||||||
|
dtype == paddle.fluid.core.VarDesc.VarType.COMPLEX128
|
||||||
|
|
||||||
|
|
||||||
|
def to_complex(dtype):
|
||||||
|
if dtype == "float32":
|
||||||
|
return np.complex64
|
||||||
|
elif dtype == "float64":
|
||||||
|
return np.complex128
|
||||||
|
else:
|
||||||
|
return dtype
|
||||||
|
|
||||||
|
|
||||||
|
def conj_dynamic(func, dtype, np_input):
|
||||||
|
paddle.set_device("cpu")
|
||||||
|
x = paddle.to_tensor(np_input)
|
||||||
|
out = func(x)
|
||||||
|
out.stop_gradient = False
|
||||||
|
sum_out = paddle.sum(out)
|
||||||
|
if is_complex(sum_out.dtype):
|
||||||
|
sum_out.real().backward()
|
||||||
|
else:
|
||||||
|
sum_out.backward()
|
||||||
|
return out.numpy(), x.grad
|
||||||
|
|
||||||
|
|
||||||
|
def conj_static(func, shape, dtype, np_input):
|
||||||
|
paddle.enable_static()
|
||||||
|
paddle.set_device("cpu")
|
||||||
|
with static.scope_guard(static.Scope()):
|
||||||
|
with static.program_guard(static.Program()):
|
||||||
|
x = static.data(name="x", shape=shape, dtype=dtype)
|
||||||
|
x.stop_gradient = False
|
||||||
|
out = func(x)
|
||||||
|
sum_out = paddle.sum(out)
|
||||||
|
static.append_backward(sum_out)
|
||||||
|
|
||||||
|
exe = static.Executor()
|
||||||
|
exe.run(static.default_startup_program())
|
||||||
|
|
||||||
|
out_v, x_grad_v = exe.run(static.default_main_program(),
|
||||||
|
feed={"x": np_input},
|
||||||
|
fetch_list=[out.name, x.name + "@GRAD"])
|
||||||
|
paddle.disable_static()
|
||||||
|
return out_v, x_grad_v
|
||||||
|
|
||||||
|
|
||||||
|
class TestCustomConjJit(unittest.TestCase):
|
||||||
|
def setUp(self):
|
||||||
|
self.dtypes = ['float32', 'float64']
|
||||||
|
self.shape = [2, 20, 2, 3]
|
||||||
|
|
||||||
|
def check_output(self, out, pd_out, name):
|
||||||
|
self.assertTrue(
|
||||||
|
np.array_equal(out, pd_out),
|
||||||
|
"custom op {}: {},\n paddle api {}: {}".format(name, out, name,
|
||||||
|
pd_out))
|
||||||
|
|
||||||
|
def run_dynamic(self, dtype, np_input):
|
||||||
|
out, x_grad = conj_dynamic(custom_ops.custom_conj, dtype, np_input)
|
||||||
|
pd_out, pd_x_grad = conj_dynamic(paddle.conj, dtype, np_input)
|
||||||
|
|
||||||
|
self.check_output(out, pd_out, "out")
|
||||||
|
self.check_output(x_grad, pd_x_grad, "x's grad")
|
||||||
|
|
||||||
|
def run_static(self, dtype, np_input):
|
||||||
|
out, x_grad = conj_static(custom_ops.custom_conj, self.shape, dtype,
|
||||||
|
np_input)
|
||||||
|
pd_out, pd_x_grad = conj_static(paddle.conj, self.shape, dtype,
|
||||||
|
np_input)
|
||||||
|
|
||||||
|
self.check_output(out, pd_out, "out")
|
||||||
|
self.check_output(x_grad, pd_x_grad, "x's grad")
|
||||||
|
|
||||||
|
def test_dynamic(self):
|
||||||
|
for dtype in self.dtypes:
|
||||||
|
np_input = np.random.random(self.shape).astype(dtype)
|
||||||
|
self.run_dynamic(dtype, np_input)
|
||||||
|
|
||||||
|
def test_static(self):
|
||||||
|
for dtype in self.dtypes:
|
||||||
|
np_input = np.random.random(self.shape).astype(dtype)
|
||||||
|
self.run_static(dtype, np_input)
|
||||||
|
|
||||||
|
# complex only used in dynamic mode now
|
||||||
|
def test_complex_dynamic(self):
|
||||||
|
for dtype in self.dtypes:
|
||||||
|
np_input = np.random.random(self.shape).astype(
|
||||||
|
dtype) + 1j * np.random.random(self.shape).astype(dtype)
|
||||||
|
self.run_dynamic(to_complex(dtype), np_input)
|
||||||
|
|
||||||
|
|
||||||
|
if __name__ == "__main__":
|
||||||
|
unittest.main()
|
||||||
Loading…
Reference in new issue