https://github.com/halide/Halide
Tip revision: fe8f71a70ad9370fae6c9448c250cb7142158f44 authored by Steven Johnson on 18 March 2019, 21:01:39 UTC
Merge branch 'master' into pr/3387
Merge branch 'master' into pr/3387
Tip revision: fe8f71a
CodeGen_C.cpp
#include <iostream>
#include <limits>
#include "CodeGen_C.h"
#include "CodeGen_Internal.h"
#include "Deinterleave.h"
#include "IROperator.h"
#include "Lerp.h"
#include "Param.h"
#include "Simplify.h"
#include "Substitute.h"
#include "Var.h"
namespace Halide {
namespace Internal {
using std::endl;
using std::map;
using std::ostream;
using std::ostringstream;
using std::string;
using std::vector;
extern "C" unsigned char halide_internal_initmod_inlined_c[];
extern "C" unsigned char halide_internal_runtime_header_HalideRuntime_h[];
extern "C" unsigned char halide_internal_runtime_header_HalideRuntimeCuda_h[];
extern "C" unsigned char halide_internal_runtime_header_HalideRuntimeHexagonHost_h[];
extern "C" unsigned char halide_internal_runtime_header_HalideRuntimeMetal_h[];
extern "C" unsigned char halide_internal_runtime_header_HalideRuntimeOpenCL_h[];
extern "C" unsigned char halide_internal_runtime_header_HalideRuntimeOpenGLCompute_h[];
extern "C" unsigned char halide_internal_runtime_header_HalideRuntimeOpenGL_h[];
extern "C" unsigned char halide_internal_runtime_header_HalideRuntimeQurt_h[];
extern "C" unsigned char halide_internal_runtime_header_HalideRuntimeD3D12Compute_h[];
namespace {
const string headers =
"#include <iostream>\n"
"#include <math.h>\n"
"#include <float.h>\n"
"#include <assert.h>\n"
"#include <string.h>\n"
"#include <stdio.h>\n"
"#include <stdint.h>\n";
// We now add definitions of things in the runtime which are
// intended to be inlined into every module but are only expressed
// in .ll. The redundancy is regrettable (FIXME).
const string globals = R"INLINE_CODE(
extern "C" {
int64_t halide_current_time_ns(void *ctx);
void halide_profiler_pipeline_end(void *, void *);
}
#ifdef _WIN32
__declspec(dllimport) float __cdecl roundf(float);
__declspec(dllimport) double __cdecl round(double);
#else
inline float asinh_f32(float x) {return asinhf(x);}
inline float acosh_f32(float x) {return acoshf(x);}
inline float atanh_f32(float x) {return atanhf(x);}
inline double asinh_f64(double x) {return asinh(x);}
inline double acosh_f64(double x) {return acosh(x);}
inline double atanh_f64(double x) {return atanh(x);}
#endif
inline float sqrt_f32(float x) {return sqrtf(x);}
inline float sin_f32(float x) {return sinf(x);}
inline float asin_f32(float x) {return asinf(x);}
inline float cos_f32(float x) {return cosf(x);}
inline float acos_f32(float x) {return acosf(x);}
inline float tan_f32(float x) {return tanf(x);}
inline float atan_f32(float x) {return atanf(x);}
inline float sinh_f32(float x) {return sinhf(x);}
inline float cosh_f32(float x) {return coshf(x);}
inline float tanh_f32(float x) {return tanhf(x);}
inline float hypot_f32(float x, float y) {return hypotf(x, y);}
inline float exp_f32(float x) {return expf(x);}
inline float log_f32(float x) {return logf(x);}
inline float pow_f32(float x, float y) {return powf(x, y);}
inline float floor_f32(float x) {return floorf(x);}
inline float ceil_f32(float x) {return ceilf(x);}
inline float round_f32(float x) {return roundf(x);}
inline double sqrt_f64(double x) {return sqrt(x);}
inline double sin_f64(double x) {return sin(x);}
inline double asin_f64(double x) {return asin(x);}
inline double cos_f64(double x) {return cos(x);}
inline double acos_f64(double x) {return acos(x);}
inline double tan_f64(double x) {return tan(x);}
inline double atan_f64(double x) {return atan(x);}
inline double sinh_f64(double x) {return sinh(x);}
inline double cosh_f64(double x) {return cosh(x);}
inline double tanh_f64(double x) {return tanh(x);}
inline double hypot_f64(double x, double y) {return hypot(x, y);}
inline double exp_f64(double x) {return exp(x);}
inline double log_f64(double x) {return log(x);}
inline double pow_f64(double x, double y) {return pow(x, y);}
inline double floor_f64(double x) {return floor(x);}
inline double ceil_f64(double x) {return ceil(x);}
inline double round_f64(double x) {return round(x);}
inline float nan_f32() {return NAN;}
inline float neg_inf_f32() {return -INFINITY;}
inline float inf_f32() {return INFINITY;}
inline bool is_nan_f32(float x) {return x != x;}
inline bool is_nan_f64(double x) {return x != x;}
template<typename A, typename B>
inline A reinterpret(const B &b) {
#if __cplusplus >= 201103L
static_assert(sizeof(A) == sizeof(B), "type size mismatch");
#endif
A a;
memcpy(&a, &b, sizeof(a));
return a;
}
inline float float_from_bits(uint32_t bits) {
return reinterpret<float, uint32_t>(bits);
}
template<typename T>
inline int halide_popcount(T a) {
int bits_set = 0;
while (a != 0) {
bits_set += a & 1;
a >>= 1;
}
return bits_set;
}
template<typename T>
inline int halide_count_leading_zeros(T a) {
int leading_zeros = 0;
int bit = sizeof(a) * 8 - 1;
while (bit >= 0 && (a & (((T)1) << bit)) == 0) {
leading_zeros++;
bit--;
}
return leading_zeros;
}
template<typename T>
inline int halide_count_trailing_zeros(T a) {
int trailing_zeros = 0;
constexpr int bits = sizeof(a) * 8;
int bit = 0;
while (bit < bits && (a & (((T)1) << bit)) == 0) {
trailing_zeros++;
bit++;
}
return trailing_zeros;
}
template<typename T>
inline T halide_cpp_max(const T &a, const T &b) {return (a > b) ? a : b;}
template<typename T>
inline T halide_cpp_min(const T &a, const T &b) {return (a < b) ? a : b;}
template<typename A, typename B>
const B &return_second(const A &a, const B &b) {
(void) a;
return b;
}
template<typename A, typename B>
inline auto quiet_div(const A &a, const B &b) -> decltype(a / b) {
return b == 0 ? static_cast<decltype(a / b)>(0) : (a / b);
}
template<typename A, typename B>
inline auto quiet_mod(const A &a, const B &b) -> decltype(a % b) {
return b == 0 ? static_cast<decltype(a % b)>(0) : (a % b);
}
namespace {
class HalideFreeHelper {
typedef void (*FreeFunction)(void *user_context, void *p);
void * user_context;
void *p;
FreeFunction free_function;
public:
HalideFreeHelper(void *user_context, void *p, FreeFunction free_function)
: user_context(user_context), p(p), free_function(free_function) {}
~HalideFreeHelper() { free(); }
void free() {
if (p) {
// TODO: do all free_functions guarantee to ignore a nullptr?
free_function(user_context, p);
p = nullptr;
}
}
};
} // namespace
)INLINE_CODE";
} // namespace
class TypeInfoGatherer : public IRGraphVisitor {
public:
std::set<ForType> for_types_used;
std::set<Type> vector_types_used;
using IRGraphVisitor::include;
using IRGraphVisitor::visit;
void include(const Expr &e) override {
if (e.type().is_vector()) {
if (e.type().is_bool()) {
// bool vectors are always emitted as uint8 in the C++ backend
// TODO: on some architectures, we could do better by choosing
// a bitwidth that matches the other vectors in use; EliminateBoolVectors
// could be used for this with a bit of work.
vector_types_used.insert(UInt(8).with_lanes(e.type().lanes()));
} else if (!e.type().is_handle()) {
// Vector-handle types can be seen when processing (e.g.)
// require() statements that are vectorized, but they
// will all be scalarized away prior to use, so don't emit
// them.
vector_types_used.insert(e.type());
}
}
IRGraphVisitor::include(e);
}
// GCC's __builtin_shuffle takes an integer vector of
// the size of its input vector. Make sure this type exists.
void visit(const Shuffle *op) override {
vector_types_used.insert(Int(32, op->vectors[0].type().lanes()));
IRGraphVisitor::visit(op);
}
void visit(const For *op) override {
for_types_used.insert(op->for_type);
IRGraphVisitor::visit(op);
}
};
CodeGen_C::CodeGen_C(ostream &s, Target t, OutputKind output_kind, const std::string &guard) :
IRPrinter(s), id("$$ BAD ID $$"), target(t), output_kind(output_kind), extern_c_open(false) {
if (is_header()) {
// If it's a header, emit an include guard.
stream << "#ifndef HALIDE_" << print_name(guard) << '\n'
<< "#define HALIDE_" << print_name(guard) << '\n'
<< "#include <stdint.h>\n"
<< "\n"
<< "// Forward declarations of the types used in the interface\n"
<< "// to the Halide pipeline.\n"
<< "//\n";
if (target.has_feature(Target::NoRuntime)) {
stream << "// For the definitions of these structs, include HalideRuntime.h\n";
} else {
stream << "// Definitions for these structs are below.\n";
}
stream << "\n"
<< "// Halide's representation of a multi-dimensional array.\n"
<< "// Halide::Runtime::Buffer is a more user-friendly wrapper\n"
<< "// around this. Its declaration is in HalideBuffer.h\n"
<< "struct halide_buffer_t;\n"
<< "\n"
<< "// Metadata describing the arguments to the generated function.\n"
<< "// Used to construct calls to the _argv version of the function.\n"
<< "struct halide_filter_metadata_t;\n"
<< "\n";
// We just forward declared the following types:
forward_declared.insert(type_of<halide_buffer_t *>().handle_type);
forward_declared.insert(type_of<halide_filter_metadata_t *>().handle_type);
if (t.has_feature(Target::LegacyBufferWrappers)) {
stream << "// The legacy buffer type. Do not use in new code.\n"
<< "struct buffer_t;\n"
<< "\n";
forward_declared.insert(type_of<buffer_t *>().handle_type);
}
} else {
// Include declarations of everything generated C source might want
stream
<< headers
<< globals
<< halide_internal_runtime_header_HalideRuntime_h << '\n'
<< halide_internal_initmod_inlined_c << '\n';
}
// Throw in a default (empty) definition of HALIDE_FUNCTION_ATTRS
// (some hosts may define this to e.g. __attribute__((warn_unused_result)))
stream << "#ifndef HALIDE_FUNCTION_ATTRS\n";
stream << "#define HALIDE_FUNCTION_ATTRS\n";
stream << "#endif\n";
}
CodeGen_C::~CodeGen_C() {
set_name_mangling_mode(NameMangling::Default);
if (is_header()) {
if (!target.has_feature(Target::NoRuntime)) {
stream << "\n"
<< "// The generated object file that goes with this header\n"
<< "// includes a full copy of the Halide runtime so that it\n"
<< "// can be used standalone. Declarations for the functions\n"
<< "// in the Halide runtime are below.\n";
if (target.os == Target::Windows) {
stream
<< "//\n"
<< "// The inclusion of this runtime means that it is not legal\n"
<< "// to link multiple Halide-generated object files together.\n"
<< "// This problem is Windows-specific. On other platforms, we\n"
<< "// use weak linkage.\n";
} else {
stream
<< "//\n"
<< "// The runtime is defined using weak linkage, so it is legal\n"
<< "// to link multiple Halide-generated object files together,\n"
<< "// or to clobber any of these functions with your own\n"
<< "// definition.\n";
}
stream << "//\n"
<< "// To generate an object file without a full copy of the\n"
<< "// runtime, use the -no_runtime target flag. To generate a\n"
<< "// standalone Halide runtime to use with such object files\n"
<< "// use the -r flag with any Halide generator binary, e.g.:\n"
<< "// $ ./my_generator -r halide_runtime -o . target=host\n"
<< "\n"
<< halide_internal_runtime_header_HalideRuntime_h << '\n';
if (target.has_feature(Target::CUDA)) {
stream << halide_internal_runtime_header_HalideRuntimeCuda_h << '\n';
}
if (target.has_feature(Target::HVX_128) ||
target.has_feature(Target::HVX_64)) {
stream << halide_internal_runtime_header_HalideRuntimeHexagonHost_h << '\n';
}
if (target.has_feature(Target::Metal)) {
stream << halide_internal_runtime_header_HalideRuntimeMetal_h << '\n';
}
if (target.has_feature(Target::OpenCL)) {
stream << halide_internal_runtime_header_HalideRuntimeOpenCL_h << '\n';
}
if (target.has_feature(Target::OpenGLCompute)) {
stream << halide_internal_runtime_header_HalideRuntimeOpenGLCompute_h << '\n';
}
if (target.has_feature(Target::OpenGL)) {
stream << halide_internal_runtime_header_HalideRuntimeOpenGL_h << '\n';
}
if (target.has_feature(Target::D3D12Compute)) {
stream << halide_internal_runtime_header_HalideRuntimeD3D12Compute_h << '\n';
}
}
stream << "#endif\n";
}
}
namespace {
string type_to_c_type(Type type, bool include_space, bool c_plus_plus = true) {
bool needs_space = true;
ostringstream oss;
if (type.is_float()) {
if (type.bits() == 32) {
oss << "float";
} else if (type.bits() == 64) {
oss << "double";
} else {
user_error << "Can't represent a float with this many bits in C: " << type << "\n";
}
if (type.is_vector()) {
oss << type.lanes();
}
} else if (type.is_handle()) {
needs_space = false;
// If there is no type info or is generating C (not C++) and
// the type is a class or in an inner scope, just use void *.
if (type.handle_type == NULL ||
(!c_plus_plus &&
(!type.handle_type->namespaces.empty() ||
!type.handle_type->enclosing_types.empty() ||
type.handle_type->inner_name.cpp_type_type == halide_cplusplus_type_name::Class))) {
oss << "void *";
} else {
if (type.handle_type->inner_name.cpp_type_type ==
halide_cplusplus_type_name::Struct) {
oss << "struct ";
}
if (!type.handle_type->namespaces.empty() ||
!type.handle_type->enclosing_types.empty()) {
oss << "::";
for (size_t i = 0; i < type.handle_type->namespaces.size(); i++) {
oss << type.handle_type->namespaces[i] << "::";
}
for (size_t i = 0; i < type.handle_type->enclosing_types.size(); i++) {
oss << type.handle_type->enclosing_types[i].name << "::";
}
}
oss << type.handle_type->inner_name.name;
if (type.handle_type->reference_type == halide_handle_cplusplus_type::LValueReference) {
oss << " &";
} else if (type.handle_type->reference_type == halide_handle_cplusplus_type::LValueReference) {
oss << " &&";
}
for (auto modifier : type.handle_type->cpp_type_modifiers) {
if (modifier & halide_handle_cplusplus_type::Const) {
oss << " const";
}
if (modifier & halide_handle_cplusplus_type::Volatile) {
oss << " volatile";
}
if (modifier & halide_handle_cplusplus_type::Restrict) {
oss << " restrict";
}
if (modifier & halide_handle_cplusplus_type::Pointer) {
oss << " *";
}
}
}
} else {
// This ends up using different type names than OpenCL does
// for the integer vector types. E.g. uint16x8_t rather than
// OpenCL's short8. Should be fine as CodeGen_C introduces
// typedefs for them and codegen always goes through this
// routine or its override in CodeGen_OpenCL to make the
// names. This may be the better bet as the typedefs are less
// likely to collide with built-in types (e.g. the OpenCL
// ones for a C compiler that decides to compile OpenCL).
// This code also supports arbitrary vector sizes where the
// OpenCL ones must be one of 2, 3, 4, 8, 16, which is too
// restrictive for already existing architectures.
switch (type.bits()) {
case 1:
// bool vectors are always emitted as uint8 in the C++ backend
if (type.is_vector()) {
oss << "uint8x" << type.lanes() << "_t";
} else {
oss << "bool";
}
break;
case 8: case 16: case 32: case 64:
if (type.is_uint()) {
oss << 'u';
}
oss << "int" << type.bits();
if (type.is_vector()) {
oss << "x" << type.lanes();
}
oss << "_t";
break;
default:
user_error << "Can't represent an integer with this many bits in C: " << type << "\n";
}
}
if (include_space && needs_space)
oss << " ";
return oss.str();
}
} // namespace
void CodeGen_C::add_vector_typedefs(const std::set<Type> &vector_types) {
if (!vector_types.empty()) {
// MSVC has a limit of ~16k for string literals, so split
// up these declarations accordingly
const char *cpp_vector_decl = R"INLINE_CODE(
#if !defined(__has_attribute)
#define __has_attribute(x) 0
#endif
#if !defined(__has_builtin)
#define __has_builtin(x) 0
#endif
template <typename ElementType_, size_t Lanes_>
class CppVector {
public:
typedef ElementType_ ElementType;
static const size_t Lanes = Lanes_;
typedef CppVector<ElementType, Lanes> Vec;
typedef CppVector<uint8_t, Lanes> Mask;
CppVector &operator=(const Vec &src) {
if (this != &src) {
for (size_t i = 0; i < Lanes; i++) {
elements[i] = src[i];
}
}
return *this;
}
/* not-explicit */ CppVector(const Vec &src) {
for (size_t i = 0; i < Lanes; i++) {
elements[i] = src[i];
}
}
CppVector() {
for (size_t i = 0; i < Lanes; i++) {
elements[i] = 0;
}
}
static Vec broadcast(const ElementType &v) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = v;
}
return r;
}
static Vec ramp(const ElementType &base, const ElementType &stride) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = base + stride * i;
}
return r;
}
static Vec load(const void *base, int32_t offset) {
Vec r(empty);
memcpy(&r.elements[0], ((const ElementType*)base + offset), sizeof(r.elements));
return r;
}
// gather
static Vec load(const void *base, const CppVector<int32_t, Lanes> &offset) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = ((const ElementType*)base)[offset[i]];
}
return r;
}
void store(void *base, int32_t offset) const {
memcpy(((ElementType*)base + offset), &this->elements[0], sizeof(this->elements));
}
// scatter
void store(void *base, const CppVector<int32_t, Lanes> &offset) const {
for (size_t i = 0; i < Lanes; i++) {
((ElementType*)base)[offset[i]] = elements[i];
}
}
static Vec shuffle(const Vec &a, const int32_t indices[Lanes]) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
if (indices[i] < 0) {
continue;
}
r.elements[i] = a[indices[i]];
}
return r;
}
template<size_t InputLanes>
static Vec concat(size_t count, const CppVector<ElementType, InputLanes> vecs[]) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = vecs[i / InputLanes][i % InputLanes];
}
return r;
}
Vec replace(size_t i, const ElementType &b) const {
Vec r = *this;
r.elements[i] = b;
return r;
}
ElementType operator[](size_t i) const {
return elements[i];
}
Vec operator~() const {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = ~elements[i];
}
return r;
}
friend Vec operator+(const Vec &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] + b[i];
}
return r;
}
friend Vec operator-(const Vec &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] - b[i];
}
return r;
}
friend Vec operator*(const Vec &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] * b[i];
}
return r;
}
friend Vec operator/(const Vec &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] / b[i];
}
return r;
}
friend Vec operator%(const Vec &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] % b[i];
}
return r;
}
friend Vec operator<<(const Vec &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] << b[i];
}
return r;
}
friend Vec operator>>(const Vec &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] >> b[i];
}
return r;
}
friend Vec operator&(const Vec &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] & b[i];
}
return r;
}
friend Vec operator|(const Vec &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] | b[i];
}
return r;
}
friend Vec operator+(const Vec &a, const ElementType &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] + b;
}
return r;
}
friend Vec operator-(const Vec &a, const ElementType &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] - b;
}
return r;
}
friend Vec operator*(const Vec &a, const ElementType &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] * b;
}
return r;
}
friend Vec operator/(const Vec &a, const ElementType &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] / b;
}
return r;
}
friend Vec operator%(const Vec &a, const ElementType &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] % b;
}
return r;
}
friend Vec operator>>(const Vec &a, const ElementType &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] >> b;
}
return r;
}
friend Vec operator<<(const Vec &a, const ElementType &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] << b;
}
return r;
}
friend Vec operator&(const Vec &a, const ElementType &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] & b;
}
return r;
}
friend Vec operator|(const Vec &a, const ElementType &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] | b;
}
return r;
}
friend Vec operator+(const ElementType &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a + b[i];
}
return r;
}
friend Vec operator-(const ElementType &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a - b[i];
}
return r;
}
friend Vec operator*(const ElementType &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a * b[i];
}
return r;
}
friend Vec operator/(const ElementType &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a / b[i];
}
return r;
}
friend Vec operator%(const ElementType &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a % b[i];
}
return r;
}
friend Vec operator>>(const ElementType &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a >> b[i];
}
return r;
}
friend Vec operator<<(const ElementType &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a << b[i];
}
return r;
}
friend Vec operator&(const ElementType &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a & b[i];
}
return r;
}
friend Vec operator|(const ElementType &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a | b[i];
}
return r;
}
friend Mask operator<(const Vec &a, const Vec &b) {
Mask r;
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] < b[i] ? 0xff : 0x00;
}
return r;
}
friend Mask operator<=(const Vec &a, const Vec &b) {
Mask r;
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] <= b[i] ? 0xff : 0x00;
}
return r;
}
friend Mask operator>(const Vec &a, const Vec &b) {
Mask r;
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] > b[i] ? 0xff : 0x00;
}
return r;
}
friend Mask operator>=(const Vec &a, const Vec &b) {
Mask r;
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] >= b[i] ? 0xff : 0x00;
}
return r;
}
friend Mask operator==(const Vec &a, const Vec &b) {
Mask r;
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] == b[i] ? 0xff : 0x00;
}
return r;
}
friend Mask operator!=(const Vec &a, const Vec &b) {
Mask r;
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = a[i] != b[i] ? 0xff : 0x00;
}
return r;
}
static Vec select(const Mask &cond, const Vec &true_value, const Vec &false_value) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = cond[i] ? true_value[i] : false_value[i];
}
return r;
}
template <typename OtherVec>
static Vec convert_from(const OtherVec &src) {
#if __cplusplus >= 201103L
static_assert(Vec::Lanes == OtherVec::Lanes, "Lanes mismatch");
#endif
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = static_cast<typename Vec::ElementType>(src[i]);
}
return r;
}
static Vec max(const Vec &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = ::halide_cpp_max(a[i], b[i]);
}
return r;
}
static Vec min(const Vec &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.elements[i] = ::halide_cpp_min(a[i], b[i]);
}
return r;
}
private:
template <typename, size_t> friend class CppVector;
ElementType elements[Lanes];
// Leave vector uninitialized for cases where we overwrite every entry
enum Empty { empty };
CppVector(Empty) {}
};
)INLINE_CODE";
const char *native_vector_decl = R"INLINE_CODE(
#if __has_attribute(ext_vector_type) || __has_attribute(vector_size)
template <typename ElementType_, size_t Lanes_>
class NativeVector {
public:
typedef ElementType_ ElementType;
static const size_t Lanes = Lanes_;
typedef NativeVector<ElementType, Lanes> Vec;
typedef NativeVector<uint8_t, Lanes> Mask;
#if __has_attribute(ext_vector_type)
typedef ElementType_ NativeVectorType __attribute__((ext_vector_type(Lanes), aligned(sizeof(ElementType))));
#elif __has_attribute(vector_size) || __GNUC__
typedef ElementType_ NativeVectorType __attribute__((vector_size(Lanes * sizeof(ElementType)), aligned(sizeof(ElementType))));
#endif
NativeVector &operator=(const Vec &src) {
if (this != &src) {
native_vector = src.native_vector;
}
return *this;
}
/* not-explicit */ NativeVector(const Vec &src) {
native_vector = src.native_vector;
}
NativeVector() {
native_vector = (NativeVectorType){};
}
static Vec broadcast(const ElementType &v) {
Vec zero; // Zero-initialized native vector.
return zero + v;
}
// TODO: this should be improved by taking advantage of native operator support.
static Vec ramp(const ElementType &base, const ElementType &stride) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.native_vector[i] = base + stride * i;
}
return r;
}
// TODO: could this be improved by taking advantage of native operator support?
static Vec load(const void *base, int32_t offset) {
Vec r(empty);
memcpy(&r.native_vector, ((const ElementType*)base + offset), sizeof(NativeVectorType));
return r;
}
// gather
// TODO: could this be improved by taking advantage of native operator support?
static Vec load(const void *base, const NativeVector<int32_t, Lanes> &offset) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.native_vector[i] = ((const ElementType*)base)[offset[i]];
}
return r;
}
// TODO: could this be improved by taking advantage of native operator support?
void store(void *base, int32_t offset) const {
memcpy(((ElementType*)base + offset), &native_vector, sizeof(NativeVectorType));
}
// scatter
// TODO: could this be improved by taking advantage of native operator support?
void store(void *base, const NativeVector<int32_t, Lanes> &offset) const {
for (size_t i = 0; i < Lanes; i++) {
((ElementType*)base)[offset[i]] = native_vector[i];
}
}
// TODO: this should be improved by taking advantage of native operator support.
static Vec shuffle(const Vec &a, const int32_t indices[Lanes]) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
if (indices[i] < 0) {
continue;
}
r.native_vector[i] = a[indices[i]];
}
return r;
}
// TODO: this should be improved by taking advantage of native operator support.
template<size_t InputLanes>
static Vec concat(size_t count, const NativeVector<ElementType, InputLanes> vecs[]) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.native_vector[i] = vecs[i / InputLanes][i % InputLanes];
}
return r;
}
// TODO: this should be improved by taking advantage of native operator support.
Vec replace(size_t i, const ElementType &b) const {
Vec r = *this;
r.native_vector[i] = b;
return r;
}
ElementType operator[](size_t i) const {
return native_vector[i];
}
Vec operator~() const {
return Vec(from_native_vector, ~native_vector);
}
friend Vec operator+(const Vec &a, const Vec &b) {
return Vec(from_native_vector, a.native_vector + b.native_vector);
}
friend Vec operator-(const Vec &a, const Vec &b) {
return Vec(from_native_vector, a.native_vector - b.native_vector);
}
friend Vec operator*(const Vec &a, const Vec &b) {
return Vec(from_native_vector, a.native_vector * b.native_vector);
}
friend Vec operator/(const Vec &a, const Vec &b) {
return Vec(from_native_vector, a.native_vector / b.native_vector);
}
friend Vec operator%(const Vec &a, const Vec &b) {
return Vec(from_native_vector, a.native_vector % b.native_vector);
}
friend Vec operator<<(const Vec &a, const Vec &b) {
return Vec(from_native_vector, a.native_vector << b.native_vector);
}
friend Vec operator>>(const Vec &a, const Vec &b) {
return Vec(from_native_vector, a.native_vector >> b.native_vector);
}
friend Vec operator&(const Vec &a, const Vec &b) {
return Vec(from_native_vector, a.native_vector & b.native_vector);
}
friend Vec operator|(const Vec &a, const Vec &b) {
return Vec(from_native_vector, a.native_vector | b.native_vector);
}
friend Vec operator+(const Vec &a, const ElementType &b) {
return Vec(from_native_vector, a.native_vector + b);
}
friend Vec operator-(const Vec &a, const ElementType &b) {
return Vec(from_native_vector, a.native_vector - b);
}
friend Vec operator*(const Vec &a, const ElementType &b) {
return Vec(from_native_vector, a.native_vector * b);
}
friend Vec operator/(const Vec &a, const ElementType &b) {
return Vec(from_native_vector, a.native_vector / b);
}
friend Vec operator%(const Vec &a, const ElementType &b) {
return Vec(from_native_vector, a.native_vector % b);
}
friend Vec operator<<(const Vec &a, const ElementType &b) {
return Vec(from_native_vector, a.native_vector << b);
}
friend Vec operator>>(const Vec &a, const ElementType &b) {
return Vec(from_native_vector, a.native_vector >> b);
}
friend Vec operator&(const Vec &a, const ElementType &b) {
return Vec(from_native_vector, a.native_vector & b);
}
friend Vec operator|(const Vec &a, const ElementType &b) {
return Vec(from_native_vector, a.native_vector | b);
}
friend Vec operator+(const ElementType &a, const Vec &b) {
return Vec(from_native_vector, a + b.native_vector);
}
friend Vec operator-(const ElementType &a, const Vec &b) {
return Vec(from_native_vector, a - b.native_vector);
}
friend Vec operator*(const ElementType &a, const Vec &b) {
return Vec(from_native_vector, a * b.native_vector);
}
friend Vec operator/(const ElementType &a, const Vec &b) {
return Vec(from_native_vector, a / b.native_vector);
}
friend Vec operator%(const ElementType &a, const Vec &b) {
return Vec(from_native_vector, a % b.native_vector);
}
friend Vec operator<<(const ElementType &a, const Vec &b) {
return Vec(from_native_vector, a << b.native_vector);
}
friend Vec operator>>(const ElementType &a, const Vec &b) {
return Vec(from_native_vector, a >> b.native_vector);
}
friend Vec operator&(const ElementType &a, const Vec &b) {
return Vec(from_native_vector, a & b.native_vector);
}
friend Vec operator|(const ElementType &a, const Vec &b) {
return Vec(from_native_vector, a | b.native_vector);
}
// TODO: this should be improved by taking advantage of native operator support.
friend Mask operator<(const Vec &a, const Vec &b) {
Mask r;
for (size_t i = 0; i < Lanes; i++) {
r.native_vector[i] = a[i] < b[i] ? 0xff : 0x00;
}
return r;
}
// TODO: this should be improved by taking advantage of native operator support.
friend Mask operator<=(const Vec &a, const Vec &b) {
Mask r;
for (size_t i = 0; i < Lanes; i++) {
r.native_vector[i] = a[i] <= b[i] ? 0xff : 0x00;
}
return r;
}
// TODO: this should be improved by taking advantage of native operator support.
friend Mask operator>(const Vec &a, const Vec &b) {
Mask r;
for (size_t i = 0; i < Lanes; i++) {
r.native_vector[i] = a[i] > b[i] ? 0xff : 0x00;
}
return r;
}
// TODO: this should be improved by taking advantage of native operator support.
friend Mask operator>=(const Vec &a, const Vec &b) {
Mask r;
for (size_t i = 0; i < Lanes; i++) {
r.native_vector[i] = a[i] >= b[i] ? 0xff : 0x00;
}
return r;
}
// TODO: this should be improved by taking advantage of native operator support.
friend Mask operator==(const Vec &a, const Vec &b) {
Mask r;
for (size_t i = 0; i < Lanes; i++) {
r.native_vector[i] = a[i] == b[i] ? 0xff : 0x00;
}
return r;
}
// TODO: this should be improved by taking advantage of native operator support.
friend Mask operator!=(const Vec &a, const Vec &b) {
Mask r;
for (size_t i = 0; i < Lanes; i++) {
r.native_vector[i] = a[i] != b[i] ? 0xff : 0x00;
}
return r;
}
// TODO: this should be improved by taking advantage of native operator support.
static Vec select(const Mask &cond, const Vec &true_value, const Vec &false_value) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.native_vector[i] = cond[i] ? true_value[i] : false_value[i];
}
return r;
}
template <typename OtherVec>
static Vec convert_from(const OtherVec &src) {
#if __cplusplus >= 201103L
static_assert(Vec::Lanes == OtherVec::Lanes, "Lanes mismatch");
#endif
#if 0 // __has_builtin(__builtin_convertvector)
// Disabled (for now) because __builtin_convertvector appears to have
// different float->int rounding behavior in at least some situations;
// for now we'll use the much-slower-but-correct explicit C++ code.
// (https://github.com/halide/Halide/issues/2080)
return Vec(from_native_vector, __builtin_convertvector(src.native_vector, NativeVectorType));
#else
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.native_vector[i] = static_cast<typename Vec::ElementType>(src.native_vector[i]);
}
return r;
#endif
}
// TODO: this should be improved by taking advantage of native operator support.
static Vec max(const Vec &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.native_vector[i] = ::halide_cpp_max(a[i], b[i]);
}
return r;
}
// TODO: this should be improved by taking advantage of native operator support.
static Vec min(const Vec &a, const Vec &b) {
Vec r(empty);
for (size_t i = 0; i < Lanes; i++) {
r.native_vector[i] = ::halide_cpp_min(a[i], b[i]);
}
return r;
}
private:
template<typename, size_t> friend class NativeVector;
NativeVectorType native_vector;
// Leave vector uninitialized for cases where we overwrite every entry
enum Empty { empty };
inline NativeVector(Empty) {}
// Syntactic sugar to avoid ctor overloading issues
enum FromNativeVector { from_native_vector };
inline NativeVector(FromNativeVector, const NativeVectorType &src) {
native_vector = src;
}
};
#endif // __has_attribute(ext_vector_type) || __has_attribute(vector_size)
)INLINE_CODE";
const char *vector_selection_decl = R"INLINE_CODE(
// Dec. 1, 2018: Apparently emscripten compilation runs with the __has_attribute true,
// then fails to handle the vector intrinsics later.
#if !defined(__EMSCRIPTEN__) && (__has_attribute(ext_vector_type) || __has_attribute(vector_size))
#if __GNUC__ && !__clang__
// GCC only allows powers-of-two; fall back to CppVector for other widths
#define halide_cpp_use_native_vector(type, lanes) ((lanes & (lanes - 1)) == 0)
#else
#define halide_cpp_use_native_vector(type, lanes) (true)
#endif
#else
// No NativeVector available
#define halide_cpp_use_native_vector(type, lanes) (false)
#endif // __has_attribute(ext_vector_type) || __has_attribute(vector_size)
// Failsafe to allow forcing non-native vectors in case of unruly compilers
#if HALIDE_CPP_ALWAYS_USE_CPP_VECTORS
#undef halide_cpp_use_native_vector
#define halide_cpp_use_native_vector(type, lanes) (false)
#endif
)INLINE_CODE";
// Vodoo fix: on at least one config (our arm32 buildbot running gcc 5.4),
// emitting this long text string was regularly garbled in a predictable pattern;
// flushing the stream before or after heals it. Since C++ codegen is rarely
// on a compilation critical path, we'll just band-aid it in this way.
stream << std::flush;
stream << cpp_vector_decl << native_vector_decl << vector_selection_decl;
stream << std::flush;
for (const auto &t : vector_types) {
string name = type_to_c_type(t, false, false);
string scalar_name = type_to_c_type(t.element_of(), false, false);
stream << "#if halide_cpp_use_native_vector(" << scalar_name << ", " << t.lanes() << ")\n";
stream << "typedef NativeVector<" << scalar_name << ", " << t.lanes() << "> " << name << ";\n";
// Useful for debugging which Vector implementation is being selected
// stream << "#pragma message \"using NativeVector for " << t << "\"\n";
stream << "#else\n";
stream << "typedef CppVector<" << scalar_name << ", " << t.lanes() << "> " << name << ";\n";
// Useful for debugging which Vector implementation is being selected
// stream << "#pragma message \"using CppVector for " << t << "\"\n";
stream << "#endif\n";
}
}
}
void CodeGen_C::set_name_mangling_mode(NameMangling mode) {
if (extern_c_open && mode != NameMangling::C) {
stream << "\n#ifdef __cplusplus\n";
stream << "} // extern \"C\"\n";
stream << "#endif\n\n";
extern_c_open = false;
} else if (!extern_c_open && mode == NameMangling::C) {
stream << "\n#ifdef __cplusplus\n";
stream << "extern \"C\" {\n";
stream << "#endif\n\n";
extern_c_open = true;
}
}
string CodeGen_C::print_type(Type type, AppendSpaceIfNeeded space_option) {
return type_to_c_type(type, space_option == AppendSpace);
}
string CodeGen_C::print_reinterpret(Type type, Expr e) {
ostringstream oss;
if (type.is_handle() || e.type().is_handle()) {
// Use a c-style cast if either src or dest is a handle --
// note that although Halide declares a "Handle" to always be 64 bits,
// the source "handle" might actually be a 32-bit pointer (from
// a function parameter), so calling reinterpret<> (which just memcpy's)
// would be garbage-producing.
oss << "(" << print_type(type) << ")";
} else {
oss << "reinterpret<" << print_type(type) << ">";
}
oss << "(" << print_expr(e) << ")";
return oss.str();
}
string CodeGen_C::print_name(const string &name) {
ostringstream oss;
// Prefix an underscore to avoid reserved words (e.g. a variable named "while")
if (isalpha(name[0])) {
oss << '_';
}
for (size_t i = 0; i < name.size(); i++) {
if (name[i] == '.') {
oss << '_';
} else if (name[i] == '$') {
oss << "__";
} else if (name[i] != '_' && !isalnum(name[i])) {
oss << "___";
}
else oss << name[i];
}
return oss.str();
}
namespace {
class ExternCallPrototypes : public IRGraphVisitor {
struct NamespaceOrCall {
const Call *call; // nullptr if this is a subnamespace
std::map<string, NamespaceOrCall> names;
NamespaceOrCall(const Call *call = nullptr) : call(call) { }
};
std::map<string, NamespaceOrCall> c_plus_plus_externs;
std::map<string, const Call *> c_externs;
std::set<std::string> processed;
std::set<std::string> internal_linkage;
std::set<std::string> destructors;
using IRGraphVisitor::visit;
void visit(const Call *op) override {
IRGraphVisitor::visit(op);
if (!processed.count(op->name)) {
if (op->call_type == Call::Extern || op->call_type == Call::PureExtern) {
c_externs.insert({op->name, op});
} else if (op->call_type == Call::ExternCPlusPlus) {
std::vector<std::string> namespaces;
std::string name = extract_namespaces(op->name, namespaces);
std::map<string, NamespaceOrCall> *namespace_map = &c_plus_plus_externs;
for (const auto &ns : namespaces) {
auto insertion = namespace_map->insert({ns, NamespaceOrCall()});
namespace_map = &insertion.first->second.names;
}
namespace_map->insert({name, NamespaceOrCall(op)});
}
processed.insert(op->name);
}
if (op->is_intrinsic(Call::register_destructor)) {
internal_assert(op->args.size() == 2);
const StringImm *fn = op->args[0].as<StringImm>();
internal_assert(fn);
destructors.insert(fn->value);
}
}
void visit(const Allocate *op) override {
IRGraphVisitor::visit(op);
if (!op->free_function.empty()) {
destructors.insert(op->free_function);
}
}
void emit_function_decl(ostream &stream, const Call *op, const std::string &name) const {
// op->name (rather than the name arg) since we need the fully-qualified C++ name
if (internal_linkage.count(op->name)) {
stream << "static ";
}
stream << type_to_c_type(op->type, /* append_space */ true) << name << "(";
if (function_takes_user_context(name)) {
stream << "void *";
if (!op->args.empty()) {
stream << ", ";
}
}
for (size_t i = 0; i < op->args.size(); i++) {
if (i > 0) {
stream << ", ";
}
if (op->args[i].as<StringImm>()) {
stream << "const char *";
} else {
stream << type_to_c_type(op->args[i].type(), true);
}
}
stream << ");\n";
}
void emit_namespace_or_call(ostream &stream, const NamespaceOrCall &ns_or_call, const std::string &name) const {
if (ns_or_call.call == nullptr) {
stream << "namespace " << name << " {\n";
for (const auto &ns_or_call_inner : ns_or_call.names) {
emit_namespace_or_call(stream, ns_or_call_inner.second, ns_or_call_inner.first);
}
stream << "} // namespace " << name << "\n";
} else {
emit_function_decl(stream, ns_or_call.call, name);
}
}
public:
ExternCallPrototypes() {
// Make sure we don't catch calls that are already in the global declarations
const char *strs[] = {globals.c_str(),
(const char *)halide_internal_runtime_header_HalideRuntime_h,
(const char *)halide_internal_initmod_inlined_c};
for (const char *str : strs) {
size_t j = 0;
for (size_t i = 0; str[i]; i++) {
char c = str[i];
if (c == '(' && i > j+1) {
// Could be the end of a function_name.
string name(str + j + 1, i-j-1);
processed.insert(name);
}
if (('A' <= c && c <= 'Z') ||
('a' <= c && c <= 'z') ||
c == '_' ||
('0' <= c && c <= '9')) {
// Could be part of a function name.
} else {
j = i;
}
}
}
}
void set_internal_linkage(const std::string &name) {
internal_linkage.insert(name);
}
bool has_c_declarations() const {
return !c_externs.empty();
}
bool has_c_plus_plus_declarations() const {
return !c_plus_plus_externs.empty();
}
void emit_c_declarations(ostream &stream) const {
for (const auto &call : c_externs) {
emit_function_decl(stream, call.second, call.first);
}
for (const auto &d : destructors) {
stream << "void " << d << "(void *, void *);\n";
}
stream << "\n";
}
void emit_c_plus_plus_declarations(ostream &stream) const {
for (const auto &ns_or_call : c_plus_plus_externs) {
emit_namespace_or_call(stream, ns_or_call.second, ns_or_call.first);
}
stream << "\n";
}
};
} // namespace
void CodeGen_C::forward_declare_type_if_needed(const Type &t) {
if (!t.handle_type ||
forward_declared.count(t.handle_type) ||
t.handle_type->inner_name.cpp_type_type == halide_cplusplus_type_name::Simple) {
return;
}
for (auto &ns : t.handle_type->namespaces) {
stream << "namespace " << ns << " { ";
}
switch (t.handle_type->inner_name.cpp_type_type) {
case halide_cplusplus_type_name::Simple:
// nothing
break;
case halide_cplusplus_type_name::Struct:
stream << "struct " << t.handle_type->inner_name.name << ";";
break;
case halide_cplusplus_type_name::Class:
stream << "class " << t.handle_type->inner_name.name << ";";
break;
case halide_cplusplus_type_name::Union:
stream << "union " << t.handle_type->inner_name.name << ";";
break;
case halide_cplusplus_type_name::Enum:
internal_error << "Passing pointers to enums is unsupported\n";
break;
}
for (auto &ns : t.handle_type->namespaces) {
(void) ns;
stream << " }";
}
stream << "\n";
forward_declared.insert(t.handle_type);
}
void CodeGen_C::compile(const Module &input) {
TypeInfoGatherer type_info;
for (const auto &f : input.functions()) {
if (f.body.defined()) {
f.body.accept(&type_info);
}
}
uses_gpu_for_loops = (type_info.for_types_used.count(ForType::GPUBlock) ||
type_info.for_types_used.count(ForType::GPUThread) ||
type_info.for_types_used.count(ForType::GPULane));
// Forward-declare all the types we need; this needs to happen before
// we emit function prototypes, since those may need the types.
stream << "\n";
for (const auto &f : input.functions()) {
for (auto &arg : f.args) {
forward_declare_type_if_needed(arg.type);
}
}
stream << "\n";
if (!is_header()) {
// Emit any external-code blobs that are C++.
for (const ExternalCode &code_blob : input.external_code()) {
if (code_blob.is_c_plus_plus_source()) {
stream << "\n";
stream << "// Begin External Code: " << code_blob.name() << "\n";
stream.write((const char *) code_blob.contents().data(), code_blob.contents().size());
stream << "\n";
stream << "// End External Code: " << code_blob.name() << "\n";
stream << "\n";
}
}
add_vector_typedefs(type_info.vector_types_used);
// Emit prototypes for all external and internal-only functions.
// Gather them up and do them all up front, to reduce duplicates,
// and to make it simpler to get internal-linkage functions correct.
ExternCallPrototypes e;
for (const auto &f : input.functions()) {
f.body.accept(&e);
if (f.linkage == LinkageType::Internal) {
// We can't tell at the call site if a LoweredFunc is intended to be internal
// or not, so mark them explicitly.
e.set_internal_linkage(f.name);
}
}
if (e.has_c_plus_plus_declarations()) {
set_name_mangling_mode(NameMangling::CPlusPlus);
e.emit_c_plus_plus_declarations(stream);
}
if (e.has_c_declarations()) {
set_name_mangling_mode(NameMangling::C);
e.emit_c_declarations(stream);
}
}
for (const auto &b : input.buffers()) {
compile(b);
}
for (const auto &f : input.functions()) {
compile(f);
}
}
void CodeGen_C::compile(const LoweredFunc &f) {
// Don't put non-external function declarations in headers.
if (is_header() && f.linkage == LinkageType::Internal) {
return;
}
const std::vector<LoweredArgument> &args = f.args;
have_user_context = false;
for (size_t i = 0; i < args.size(); i++) {
// TODO: check that its type is void *?
have_user_context |= (args[i].name == "__user_context");
}
NameMangling name_mangling = f.name_mangling;
if (name_mangling == NameMangling::Default) {
name_mangling = (target.has_feature(Target::CPlusPlusMangling) ?
NameMangling::CPlusPlus : NameMangling::C);
}
set_name_mangling_mode(name_mangling);
std::vector<std::string> namespaces;
std::string simple_name = extract_namespaces(f.name, namespaces);
if (!is_c_plus_plus_interface()) {
user_assert(namespaces.empty()) <<
"Namespace qualifiers not allowed on function name if not compiling with Target::CPlusPlusNameMangling.\n";
}
if (!namespaces.empty()) {
for (const auto &ns : namespaces) {
stream << "namespace " << ns << " {\n";
}
stream << "\n";
}
// Emit the function prototype
if (f.linkage == LinkageType::Internal) {
// If the function isn't public, mark it static.
stream << "static ";
}
stream << "int " << simple_name << "(";
for (size_t i = 0; i < args.size(); i++) {
if (args[i].is_buffer()) {
stream << "struct halide_buffer_t *"
<< print_name(args[i].name)
<< "_buffer";
} else {
stream << print_type(args[i].type, AppendSpace)
<< print_name(args[i].name);
}
if (i < args.size()-1) stream << ", ";
}
if (is_header()) {
stream << ") HALIDE_FUNCTION_ATTRS;\n";
} else {
stream << ") HALIDE_FUNCTION_ATTRS {\n";
indent += 1;
if (uses_gpu_for_loops) {
do_indent();
stream << "halide_error("
<< (have_user_context ? "__user_context_" : "nullptr")
<< ", \"C++ Backend does not support gpu_blocks() or gpu_threads() yet, "
<< "this function will always fail at runtime\");\n";
do_indent();
stream << "return halide_error_code_device_malloc_failed;\n";
} else {
// Emit a local user_context we can pass in all cases, either
// aliasing __user_context or nullptr.
if (!is_header()) {
do_indent();
stream << "void * const _ucon = "
<< (have_user_context ? "const_cast<void *>(__user_context)" : "nullptr")
<< ";\n";
}
// Emit the body
print(f.body);
// Return success.
do_indent();
stream << "return 0;\n";
}
indent -= 1;
stream << "}\n";
}
if (is_header() && f.linkage == LinkageType::ExternalPlusMetadata) {
// Emit the argv version
stream << "int " << simple_name << "_argv(void **args) HALIDE_FUNCTION_ATTRS;\n";
// And also the metadata.
stream << "const struct halide_filter_metadata_t *" << simple_name << "_metadata() HALIDE_FUNCTION_ATTRS;\n";
}
if (!namespaces.empty()) {
stream << "\n";
for (size_t i = namespaces.size(); i > 0; i--) {
stream << "} // namespace " << namespaces[i-1] << "\n";
}
stream << "\n";
}
}
void CodeGen_C::compile(const Buffer<> &buffer) {
// Don't define buffers in headers.
if (is_header()) {
return;
}
string name = print_name(buffer.name());
halide_buffer_t b = *(buffer.raw_buffer());
user_assert(b.host) << "Can't embed image: " << buffer.name() << " because it has a null host pointer\n";
user_assert(!b.device_dirty()) << "Can't embed image: " << buffer.name() << "because it has a dirty device pointer\n";
// Figure out the offset of the last pixel.
size_t num_elems = 1;
for (int d = 0; b.dim[d].extent; d++) {
num_elems += b.dim[d].stride * (b.dim[d].extent - 1);
}
// For now, we assume buffers that aren't scalar are constant,
// while scalars can be mutated. This accommodates all our existing
// use cases, which is that all buffers are constant, except those
// used to store stateful module information in offloading runtimes.
bool is_constant = buffer.dimensions() != 0;
// Emit the data
stream << "static " << (is_constant ? "const" : "") << " uint8_t " << name << "_data[] HALIDE_ATTRIBUTE_ALIGN(32) = {\n";
do_indent();
for (size_t i = 0; i < num_elems * b.type.bytes(); i++) {
if (i > 0) {
stream << ",";
if (i % 16 == 0) {
stream << "\n";
do_indent();
} else {
stream << " ";
}
}
stream << (int)(b.host[i]);
}
stream << "\n};\n";
// Emit the shape (constant even for scalar buffers)
stream << "static const halide_dimension_t " << name << "_buffer_shape[] = {";
for (int i = 0; i < buffer.dimensions(); i++) {
stream << "halide_dimension_t(" << buffer.dim(i).min()
<< ", " << buffer.dim(i).extent()
<< ", " << buffer.dim(i).stride() << ")";
if (i < buffer.dimensions() - 1) {
stream << ", ";
}
}
stream << "};\n";
Type t = buffer.type();
// Emit the buffer struct. Note that although our shape and (usually) our host
// data is const, the buffer itself isn't: embedded buffers in one pipeline
// can be passed to another pipeline (e.g. for an extern stage), in which
// case the buffer objects need to be non-const, because the constness
// (from the POV of the extern stage) is a runtime property.
stream << "static halide_buffer_t " << name << "_buffer_ = {"
<< "0, " // device
<< "nullptr, " // device_interface
<< "const_cast<uint8_t*>(&" << name << "_data[0]), " // host
<< "0, " // flags
<< "halide_type_t((halide_type_code_t)(" << (int)t.code() << "), " << t.bits() << ", " << t.lanes() << "), "
<< buffer.dimensions() << ", "
<< "const_cast<halide_dimension_t*>(" << name << "_buffer_shape)};\n";
// Make a global pointer to it.
stream << "static halide_buffer_t * const " << name << "_buffer = &" << name << "_buffer_;\n";
}
string CodeGen_C::print_expr(Expr e) {
id = "$$ BAD ID $$";
e.accept(this);
return id;
}
string CodeGen_C::print_cast_expr(const Type &t, Expr e) {
string value = print_expr(e);
string type = print_type(t);
if (t.is_vector() &&
t.lanes() == e.type().lanes() &&
t != e.type()) {
return print_assignment(t, type + "::convert_from<" + print_type(e.type()) + ">(" + value + ")");
} else {
return print_assignment(t, "(" + type + ")(" + value + ")");
}
}
void CodeGen_C::print_stmt(Stmt s) {
s.accept(this);
}
string CodeGen_C::print_assignment(Type t, const std::string &rhs) {
auto cached = cache.find(rhs);
if (cached == cache.end()) {
id = unique_name('_');
do_indent();
stream << print_type(t, AppendSpace) << id << " = " << rhs << ";\n";
cache[rhs] = id;
} else {
id = cached->second;
}
return id;
}
void CodeGen_C::open_scope() {
cache.clear();
do_indent();
indent++;
stream << "{\n";
}
void CodeGen_C::close_scope(const std::string &comment) {
cache.clear();
indent--;
do_indent();
if (!comment.empty()) {
stream << "} // " << comment << "\n";
} else {
stream << "}\n";
}
}
void CodeGen_C::visit(const Variable *op) {
id = print_name(op->name);
}
void CodeGen_C::visit(const Cast *op) {
id = print_cast_expr(op->type, op->value);
}
void CodeGen_C::visit_binop(Type t, Expr a, Expr b, const char * op) {
string sa = print_expr(a);
string sb = print_expr(b);
print_assignment(t, sa + " " + op + " " + sb);
}
void CodeGen_C::visit(const Add *op) {
visit_binop(op->type, op->a, op->b, "+");
}
void CodeGen_C::visit(const Sub *op) {
visit_binop(op->type, op->a, op->b, "-");
}
void CodeGen_C::visit(const Mul *op) {
visit_binop(op->type, op->a, op->b, "*");
}
void CodeGen_C::visit(const Div *op) {
int bits;
if (is_const_power_of_two_integer(op->b, &bits)) {
visit_binop(op->type, op->a, make_const(op->a.type(), bits), ">>");
} else if (op->type.is_int()) {
print_expr(lower_euclidean_div(op->a, op->b));
} else {
visit_binop(op->type, op->a, op->b, "/");
}
}
void CodeGen_C::visit(const Mod *op) {
int bits;
if (is_const_power_of_two_integer(op->b, &bits)) {
visit_binop(op->type, op->a, make_const(op->a.type(), (1 << bits)-1), "&");
} else if (op->type.is_int()) {
print_expr(lower_euclidean_mod(op->a, op->b));
} else {
visit_binop(op->type, op->a, op->b, "%");
}
}
void CodeGen_C::visit(const Max *op) {
// clang doesn't support the ternary operator on OpenCL style vectors.
// See: https://bugs.llvm.org/show_bug.cgi?id=33103
if (op->type.is_scalar()) {
print_expr(Call::make(op->type, "::halide_cpp_max", {op->a, op->b}, Call::Extern));
} else {
ostringstream rhs;
rhs << print_type(op->type) << "::max(" << print_expr(op->a) << ", " << print_expr(op->b) << ")";
print_assignment(op->type, rhs.str());
}
}
void CodeGen_C::visit(const Min *op) {
// clang doesn't support the ternary operator on OpenCL style vectors.
// See: https://bugs.llvm.org/show_bug.cgi?id=33103
if (op->type.is_scalar()) {
print_expr(Call::make(op->type, "::halide_cpp_min", {op->a, op->b}, Call::Extern));
} else {
ostringstream rhs;
rhs << print_type(op->type) << "::min(" << print_expr(op->a) << ", " << print_expr(op->b) << ")";
print_assignment(op->type, rhs.str());
}
}
void CodeGen_C::visit(const EQ *op) {
visit_binop(op->type, op->a, op->b, "==");
}
void CodeGen_C::visit(const NE *op) {
visit_binop(op->type, op->a, op->b, "!=");
}
void CodeGen_C::visit(const LT *op) {
visit_binop(op->type, op->a, op->b, "<");
}
void CodeGen_C::visit(const LE *op) {
visit_binop(op->type, op->a, op->b, "<=");
}
void CodeGen_C::visit(const GT *op) {
visit_binop(op->type, op->a, op->b, ">");
}
void CodeGen_C::visit(const GE *op) {
visit_binop(op->type, op->a, op->b, ">=");
}
void CodeGen_C::visit(const Or *op) {
visit_binop(op->type, op->a, op->b, "||");
}
void CodeGen_C::visit(const And *op) {
visit_binop(op->type, op->a, op->b, "&&");
}
void CodeGen_C::visit(const Not *op) {
print_assignment(op->type, "!(" + print_expr(op->a) + ")");
}
void CodeGen_C::visit(const IntImm *op) {
if (op->type == Int(32)) {
id = std::to_string(op->value);
} else {
print_assignment(op->type, "(" + print_type(op->type) + ")(" + std::to_string(op->value) + ")");
}
}
void CodeGen_C::visit(const UIntImm *op) {
print_assignment(op->type, "(" + print_type(op->type) + ")(" + std::to_string(op->value) + ")");
}
void CodeGen_C::visit(const StringImm *op) {
ostringstream oss;
oss << Expr(op);
id = oss.str();
}
// NaN is the only float/double for which this is true... and
// surprisingly, there doesn't seem to be a portable isnan function
// (dsharlet).
template <typename T>
static bool isnan(T x) { return x != x; }
template <typename T>
static bool isinf(T x)
{
return std::numeric_limits<T>::has_infinity && (
x == std::numeric_limits<T>::infinity() ||
x == -std::numeric_limits<T>::infinity());
}
void CodeGen_C::visit(const FloatImm *op) {
if (isnan(op->value)) {
id = "nan_f32()";
} else if (isinf(op->value)) {
if (op->value > 0) {
id = "inf_f32()";
} else {
id = "neg_inf_f32()";
}
} else {
// Write the constant as reinterpreted uint to avoid any bits lost in conversion.
union {
uint32_t as_uint;
float as_float;
} u;
u.as_float = op->value;
ostringstream oss;
if (op->type.bits() == 64) {
oss << "(double) ";
}
oss << "float_from_bits(" << u.as_uint << " /* " << u.as_float << " */)";
print_assignment(op->type, oss.str());
}
}
void CodeGen_C::visit(const Call *op) {
internal_assert(op->is_extern() || op->is_intrinsic())
<< "Can only codegen extern calls and intrinsics\n";
ostringstream rhs;
// Handle intrinsics first
if (op->is_intrinsic(Call::debug_to_file)) {
internal_assert(op->args.size() == 3);
const StringImm *string_imm = op->args[0].as<StringImm>();
internal_assert(string_imm);
string filename = string_imm->value;
string typecode = print_expr(op->args[1]);
string buffer = print_name(print_expr(op->args[2]));
rhs << "halide_debug_to_file(_ucon, "
<< "\"" << filename << "\", "
<< typecode
<< ", (struct halide_buffer_t *)" << buffer << ")";
} else if (op->is_intrinsic(Call::bitwise_and)) {
internal_assert(op->args.size() == 2);
string a0 = print_expr(op->args[0]);
string a1 = print_expr(op->args[1]);
rhs << a0 << " & " << a1;
} else if (op->is_intrinsic(Call::bitwise_xor)) {
internal_assert(op->args.size() == 2);
string a0 = print_expr(op->args[0]);
string a1 = print_expr(op->args[1]);
rhs << a0 << " ^ " << a1;
} else if (op->is_intrinsic(Call::bitwise_or)) {
internal_assert(op->args.size() == 2);
string a0 = print_expr(op->args[0]);
string a1 = print_expr(op->args[1]);
rhs << a0 << " | " << a1;
} else if (op->is_intrinsic(Call::bitwise_not)) {
internal_assert(op->args.size() == 1);
rhs << "~" << print_expr(op->args[0]);
} else if (op->is_intrinsic(Call::reinterpret)) {
internal_assert(op->args.size() == 1);
rhs << print_reinterpret(op->type, op->args[0]);
} else if (op->is_intrinsic(Call::shift_left)) {
internal_assert(op->args.size() == 2);
string a0 = print_expr(op->args[0]);
string a1 = print_expr(op->args[1]);
rhs << a0 << " << " << a1;
} else if (op->is_intrinsic(Call::shift_right)) {
internal_assert(op->args.size() == 2);
string a0 = print_expr(op->args[0]);
string a1 = print_expr(op->args[1]);
rhs << a0 << " >> " << a1;
} else if (op->is_intrinsic(Call::count_leading_zeros) ||
op->is_intrinsic(Call::count_trailing_zeros) ||
op->is_intrinsic(Call::popcount)) {
internal_assert(op->args.size() == 1);
if (op->args[0].type().is_vector()) {
rhs << print_scalarized_expr(op);
} else {
string a0 = print_expr(op->args[0]);
rhs << "halide_" << op->name << "(" << a0 << ")";
}
} else if (op->is_intrinsic(Call::lerp)) {
internal_assert(op->args.size() == 3);
Expr e = lower_lerp(op->args[0], op->args[1], op->args[2]);
rhs << print_expr(e);
} else if (op->is_intrinsic(Call::absd)) {
internal_assert(op->args.size() == 2);
Expr a = op->args[0];
Expr b = op->args[1];
Type t = op->type.with_code(op->type.is_int() ? Type::UInt : op->type.code());
Expr e = cast(t, select(a < b, b - a, a - b));
rhs << print_expr(e);
} else if (op->is_intrinsic(Call::return_second)) {
internal_assert(op->args.size() == 2);
string arg0 = print_expr(op->args[0]);
string arg1 = print_expr(op->args[1]);
rhs << "return_second(" << arg0 << ", " << arg1 << ")";
} else if (op->is_intrinsic(Call::if_then_else)) {
internal_assert(op->args.size() == 3);
string result_id = unique_name('_');
do_indent();
stream << print_type(op->args[1].type(), AppendSpace)
<< result_id << ";\n";
string cond_id = print_expr(op->args[0]);
do_indent();
stream << "if (" << cond_id << ")\n";
open_scope();
string true_case = print_expr(op->args[1]);
do_indent();
stream << result_id << " = " << true_case << ";\n";
close_scope("if " + cond_id);
do_indent();
stream << "else\n";
open_scope();
string false_case = print_expr(op->args[2]);
do_indent();
stream << result_id << " = " << false_case << ";\n";
close_scope("if " + cond_id + " else");
rhs << result_id;
} else if (op->is_intrinsic(Call::require)) {
internal_assert(op->args.size() == 3);
if (op->args[0].type().is_vector()) {
rhs << print_scalarized_expr(op);
} else {
create_assertion(op->args[0], op->args[2]);
rhs << print_expr(op->args[1]);
}
} else if (op->is_intrinsic(Call::abs)) {
internal_assert(op->args.size() == 1);
Expr a0 = op->args[0];
rhs << print_expr(cast(op->type, select(a0 > 0, a0, -a0)));
} else if (op->is_intrinsic(Call::memoize_expr)) {
internal_assert(op->args.size() >= 1);
string arg = print_expr(op->args[0]);
rhs << "(" << arg << ")";
} else if (op->is_intrinsic(Call::alloca)) {
internal_assert(op->args.size() == 1);
internal_assert(op->type.is_handle());
const Call *call = op->args[0].as<Call>();
if (op->type == type_of<struct halide_buffer_t *>() &&
call && call->is_intrinsic(Call::size_of_halide_buffer_t)) {
do_indent();
string buf_name = unique_name('b');
stream << "halide_buffer_t " << buf_name << ";\n";
rhs << "&" << buf_name;
} else {
// Make a stack of uint64_ts
string size = print_expr(simplify((op->args[0] + 7)/8));
do_indent();
string array_name = unique_name('a');
stream << "uint64_t " << array_name << "[" << size << "];";
rhs << "(" << print_type(op->type) << ")(&" << array_name << ")";
}
} else if (op->is_intrinsic(Call::make_struct)) {
if (op->args.empty()) {
internal_assert(op->type.handle_type);
// Add explicit cast so that different structs can't cache to the same value
rhs << "(" << print_type(op->type) << ")(NULL)";
} else {
// Emit a declaration like:
// struct {const int f_0, const char f_1, const int f_2} foo = {3, 'c', 4};
// Get the args
vector<string> values;
for (size_t i = 0; i < op->args.size(); i++) {
values.push_back(print_expr(op->args[i]));
}
do_indent();
stream << "struct {\n";
// List the types.
indent++;
for (size_t i = 0; i < op->args.size(); i++) {
do_indent();
stream << "const " << print_type(op->args[i].type()) << " f_" << i << ";\n";
}
indent--;
string struct_name = unique_name('s');
do_indent();
stream << "} " << struct_name << " = {\n";
// List the values.
indent++;
for (size_t i = 0; i < op->args.size(); i++) {
do_indent();
stream << values[i];
if (i < op->args.size() - 1) stream << ",";
stream << "\n";
}
indent--;
do_indent();
stream << "};\n";
// Return a pointer to it of the appropriate type
if (op->type.handle_type) {
rhs << "(" << print_type(op->type) << ")";
}
rhs << "(&" << struct_name << ")";
}
} else if (op->is_intrinsic(Call::stringify)) {
// Rewrite to an snprintf
vector<string> printf_args;
string format_string = "";
for (size_t i = 0; i < op->args.size(); i++) {
Type t = op->args[i].type();
printf_args.push_back(print_expr(op->args[i]));
if (t.is_int()) {
format_string += "%lld";
printf_args[i] = "(long long)(" + printf_args[i] + ")";
} else if (t.is_uint()) {
format_string += "%llu";
printf_args[i] = "(long long unsigned)(" + printf_args[i] + ")";
} else if (t.is_float()) {
if (t.bits() == 32) {
format_string += "%f";
} else {
format_string += "%e";
}
} else if (op->args[i].as<StringImm>()) {
format_string += "%s";
} else {
internal_assert(t.is_handle());
format_string += "%p";
}
}
string buf_name = unique_name('b');
do_indent();
stream << "char " << buf_name << "[1024];\n";
do_indent();
stream << "snprintf(" << buf_name << ", 1024, \"" << format_string << "\", " << with_commas(printf_args) << ");\n";
rhs << buf_name;
} else if (op->is_intrinsic(Call::register_destructor)) {
internal_assert(op->args.size() == 2);
const StringImm *fn = op->args[0].as<StringImm>();
internal_assert(fn);
string arg = print_expr(op->args[1]);
do_indent();
// Make a struct on the stack that calls the given function as a destructor
string struct_name = unique_name('s');
string instance_name = unique_name('d');
stream << "struct " << struct_name << " { "
<< "void * const ucon; "
<< "void * const arg; "
<< "" << struct_name << "(void *ucon, void *a) : ucon(ucon), arg((void *)a) {} "
<< "~" << struct_name << "() { " << fn->value + "(ucon, arg); } "
<< "} " << instance_name << "(_ucon, " << arg << ");\n";
rhs << print_expr(0);
} else if (op->is_intrinsic(Call::div_round_to_zero)) {
rhs << print_expr(op->args[0]) << " / " << print_expr(op->args[1]);
} else if (op->is_intrinsic(Call::mod_round_to_zero)) {
rhs << print_expr(op->args[0]) << " % " << print_expr(op->args[1]);
} else if (op->is_intrinsic(Call::signed_integer_overflow)) {
user_error << "Signed integer overflow occurred during constant-folding. Signed"
" integer overflow for int32 and int64 is undefined behavior in"
" Halide.\n";
} else if (op->is_intrinsic(Call::quiet_div)) {
internal_assert(op->args.size() == 2);
// Don't bother checking for zero denominator here; the quiet_div
// implementation will always do a runtime check and return zero
// (rather than failing at runtime).
string a = print_expr(op->args[0]);
string b = print_expr(op->args[1]);
rhs << "::quiet_div(" << a << ", " << b << ")";
} else if (op->is_intrinsic(Call::quiet_mod)) {
internal_assert(op->args.size() == 2);
// Don't bother checking for zero denominator here; the quiet_mod
// implementation will always do a runtime check and return zero
// (rather than failing at runtime).
string a = print_expr(op->args[0]);
string b = print_expr(op->args[1]);
rhs << "::quiet_mod(" << a << ", " << b << ")";
} else if (op->is_intrinsic(Call::prefetch)) {
user_assert((op->args.size() == 4) && is_one(op->args[2]))
<< "Only prefetch of 1 cache line is supported in C backend.\n";
const Variable *base = op->args[0].as<Variable>();
internal_assert(base && base->type.is_handle());
rhs << "__builtin_prefetch("
<< "((" << print_type(op->type) << " *)" << print_name(base->name)
<< " + " << print_expr(op->args[1]) << "), 1)";
} else if (op->is_intrinsic(Call::indeterminate_expression)) {
user_error << "Indeterminate expression occurred during constant-folding.\n";
} else if (op->is_intrinsic(Call::size_of_halide_buffer_t)) {
rhs << "(sizeof(halide_buffer_t))";
} else if (op->is_intrinsic(Call::strict_float)) {
internal_assert(op->args.size() == 1);
string arg0 = print_expr(op->args[0]);
rhs << "(" << arg0 << ")";
} else if (op->is_intrinsic()) {
// TODO: other intrinsics
internal_error << "Unhandled intrinsic in C backend: " << op->name << '\n';
} else {
// Generic extern calls
rhs << print_extern_call(op);
}
print_assignment(op->type, rhs.str());
}
string CodeGen_C::print_scalarized_expr(Expr e) {
Type t = e.type();
internal_assert(t.is_vector());
string v = unique_name('_');
do_indent();
stream << print_type(t, AppendSpace) << v << ";\n";
for (int lane = 0; lane < t.lanes(); lane++) {
Expr e2 = extract_lane(e, lane);
string elem = print_expr(e2);
ostringstream rhs;
rhs << v << ".replace(" << lane << ", " << elem << ")";
v = print_assignment(t, rhs.str());
}
return v;
}
string CodeGen_C::print_extern_call(const Call *op) {
if (op->type.is_vector()) {
// Need to split into multiple scalar calls.
return print_scalarized_expr(op);
}
ostringstream rhs;
vector<string> args(op->args.size());
for (size_t i = 0; i < op->args.size(); i++) {
args[i] = print_expr(op->args[i]);
// This substitution ensures const correctness for all calls
if (args[i] == "__user_context") {
args[i] = "_ucon";
}
}
if (function_takes_user_context(op->name)) {
args.insert(args.begin(), "_ucon");
}
rhs << op->name << "(" << with_commas(args) << ")";
return rhs.str();
}
void CodeGen_C::visit(const Load *op) {
user_assert(is_one(op->predicate)) << "Predicated load is not supported by C backend.\n";
// TODO: We could replicate the logic in the llvm codegen which decides whether
// the vector access can be aligned. Doing so would also require introducing
// aligned type equivalents for all the vector types.
ostringstream rhs;
Type t = op->type;
string name = print_name(op->name);
// If we're loading a contiguous ramp into a vector, just load the vector
Expr dense_ramp_base = strided_ramp_base(op->index, 1);
if (dense_ramp_base.defined()) {
internal_assert(t.is_vector());
string id_ramp_base = print_expr(dense_ramp_base);
rhs << print_type(t) + "::load(" << name << ", " << id_ramp_base << ")";
} else if (op->index.type().is_vector()) {
// If index is a vector, gather vector elements.
internal_assert(t.is_vector());
string id_index = print_expr(op->index);
rhs << print_type(t) + "::load(" << name << ", " << id_index << ")";
} else {
string id_index = print_expr(op->index);
bool type_cast_needed = !(allocations.contains(op->name) &&
allocations.get(op->name).type.element_of() == t.element_of());
if (type_cast_needed) {
rhs << "((const " << print_type(t.element_of()) << " *)" << name << ")";
} else {
rhs << name;
}
rhs << "[" << id_index << "]";
}
print_assignment(t, rhs.str());
}
void CodeGen_C::visit(const Store *op) {
user_assert(is_one(op->predicate)) << "Predicated store is not supported by C backend.\n";
Type t = op->value.type();
string id_value = print_expr(op->value);
string name = print_name(op->name);
// TODO: We could replicate the logic in the llvm codegen which decides whether
// the vector access can be aligned. Doing so would also require introducing
// aligned type equivalents for all the vector types.
// If we're writing a contiguous ramp, just store the vector.
Expr dense_ramp_base = strided_ramp_base(op->index, 1);
if (dense_ramp_base.defined()) {
internal_assert(op->value.type().is_vector());
string id_ramp_base = print_expr(dense_ramp_base);
do_indent();
stream << id_value + ".store(" << name << ", " << id_ramp_base << ");\n";
} else if (op->index.type().is_vector()) {
// If index is a vector, scatter vector elements.
internal_assert(t.is_vector());
string id_index = print_expr(op->index);
do_indent();
stream << id_value + ".store(" << name << ", " << id_index << ");\n";
} else {
bool type_cast_needed =
t.is_handle() ||
!allocations.contains(op->name) ||
allocations.get(op->name).type != t;
string id_index = print_expr(op->index);
do_indent();
if (type_cast_needed) {
stream << "((" << print_type(t) << " *)" << name << ")";
} else {
stream << name;
}
stream << "[" << id_index << "] = " << id_value << ";\n";
}
cache.clear();
}
void CodeGen_C::visit(const Let *op) {
string id_value = print_expr(op->value);
Expr body = op->body;
if (op->value.type().is_handle()) {
// The body might contain a Load that references this directly
// by name, so we can't rewrite the name.
do_indent();
stream << print_type(op->value.type())
<< " " << print_name(op->name)
<< " = " << id_value << ";\n";
} else {
Expr new_var = Variable::make(op->value.type(), id_value);
body = substitute(op->name, new_var, body);
}
print_expr(body);
}
void CodeGen_C::visit(const Select *op) {
ostringstream rhs;
string type = print_type(op->type);
string true_val = print_expr(op->true_value);
string false_val = print_expr(op->false_value);
string cond = print_expr(op->condition);
// clang doesn't support the ternary operator on OpenCL style vectors.
// See: https://bugs.llvm.org/show_bug.cgi?id=33103
if (op->condition.type().is_scalar()) {
rhs << "(" << type << ")"
<< "(" << cond
<< " ? " << true_val
<< " : " << false_val
<< ")";
} else {
rhs << type << "::select(" << cond << ", " << true_val << ", " << false_val << ")";
}
print_assignment(op->type, rhs.str());
}
void CodeGen_C::visit(const LetStmt *op) {
string id_value = print_expr(op->value);
Stmt body = op->body;
if (op->value.type().is_handle()) {
// The body might contain a Load or Store that references this
// directly by name, so we can't rewrite the name.
do_indent();
stream << print_type(op->value.type())
<< " " << print_name(op->name)
<< " = " << id_value << ";\n";
} else {
Expr new_var = Variable::make(op->value.type(), id_value);
body = substitute(op->name, new_var, body);
}
body.accept(this);
}
// Halide asserts have different semantics to C asserts. They're
// supposed to clean up and make the containing function return
// -1, so we can't use the C version of assert. Instead we convert
// to an if statement.
void CodeGen_C::create_assertion(const string &id_cond, const string &id_msg) {
if (target.has_feature(Target::NoAsserts)) return;
do_indent();
stream << "if (!" << id_cond << ")\n";
open_scope();
do_indent();
stream << "return " << id_msg << ";\n";
close_scope("");
}
void CodeGen_C::create_assertion(const string &id_cond, Expr message) {
internal_assert(!message.defined() || message.type() == Int(32))
<< "Assertion result is not an int: " << message;
if (target.has_feature(Target::NoAsserts)) return;
// don't call the create_assertion(string, string) version because
// we don't want to force evaluation of 'message' unless the condition fails
do_indent();
stream << "if (!" << id_cond << ") ";
open_scope();
string id_msg = print_expr(message);
do_indent();
stream << "return " << id_msg << ";\n";
close_scope("");
}
void CodeGen_C::create_assertion(Expr cond, Expr message) {
create_assertion(print_expr(cond), message);
}
void CodeGen_C::visit(const AssertStmt *op) {
create_assertion(op->condition, op->message);
}
void CodeGen_C::visit(const ProducerConsumer *op) {
do_indent();
if (op->is_producer) {
stream << "// produce " << op->name << '\n';
} else {
stream << "// consume " << op->name << '\n';
}
print_stmt(op->body);
}
void CodeGen_C::visit(const Fork *op) {
// TODO: This doesn't actually work with nested tasks
do_indent();
stream << "#pragma omp parallel\n";
open_scope();
do_indent();
stream << "#pragma omp single\n";
open_scope();
do_indent();
stream << "#pragma omp task\n";
open_scope();
print_stmt(op->first);
close_scope("");
do_indent();
stream << "#pragma omp task\n";
open_scope();
print_stmt(op->rest);
close_scope("");
do_indent();
stream << "#pragma omp taskwait\n";
close_scope("");
close_scope("");
}
void CodeGen_C::visit(const Acquire *op) {
string id_sem = print_expr(op->semaphore);
string id_count = print_expr(op->count);
open_scope();
do_indent();
stream << "while (!halide_semaphore_try_acquire(" << id_sem << ", " << id_count << "))\n";
open_scope();
do_indent();
stream << "#pragma omp taskyield\n";
close_scope("");
op->body.accept(this);
close_scope("");
}
void CodeGen_C::visit(const For *op) {
string id_min = print_expr(op->min);
string id_extent = print_expr(op->extent);
if (op->for_type == ForType::Parallel) {
do_indent();
stream << "#pragma omp parallel for\n";
} else {
internal_assert(op->for_type == ForType::Serial)
<< "Can only emit serial or parallel for loops to C\n";
}
do_indent();
stream << "for (int "
<< print_name(op->name)
<< " = " << id_min
<< "; "
<< print_name(op->name)
<< " < " << id_min
<< " + " << id_extent
<< "; "
<< print_name(op->name)
<< "++)\n";
open_scope();
op->body.accept(this);
close_scope("for " + print_name(op->name));
}
void CodeGen_C::visit(const Ramp *op) {
Type vector_type = op->type.with_lanes(op->lanes);
string id_base = print_expr(op->base);
string id_stride = print_expr(op->stride);
print_assignment(vector_type, print_type(vector_type) + "::ramp(" + id_base + ", " + id_stride + ")");
}
void CodeGen_C::visit(const Broadcast *op) {
Type vector_type = op->type.with_lanes(op->lanes);
string id_value = print_expr(op->value);
string rhs;
if (op->lanes > 1) {
rhs = print_type(vector_type) + "::broadcast(" + id_value + ")";
} else {
rhs = id_value;
}
print_assignment(vector_type, rhs);
}
void CodeGen_C::visit(const Provide *op) {
internal_error << "Cannot emit Provide statements as C\n";
}
void CodeGen_C::visit(const Allocate *op) {
open_scope();
string op_name = print_name(op->name);
string op_type = print_type(op->type, AppendSpace);
// For sizes less than 8k, do a stack allocation
bool on_stack = false;
int32_t constant_size;
string size_id;
if (op->new_expr.defined()) {
Allocation alloc;
alloc.type = op->type;
allocations.push(op->name, alloc);
heap_allocations.push(op->name);
stream << op_type << "*" << op_name << " = (" << print_expr(op->new_expr) << ");\n";
} else {
constant_size = op->constant_allocation_size();
if (constant_size > 0) {
int64_t stack_bytes = constant_size * op->type.bytes();
if (stack_bytes > ((int64_t(1) << 31) - 1)) {
user_error << "Total size for allocation "
<< op->name << " is constant but exceeds 2^31 - 1.\n";
} else {
size_id = print_expr(Expr(static_cast<int32_t>(constant_size)));
if (op->memory_type == MemoryType::Stack ||
(op->memory_type == MemoryType::Auto &&
can_allocation_fit_on_stack(stack_bytes))) {
on_stack = true;
}
}
} else {
// Check that the allocation is not scalar (if it were scalar
// it would have constant size).
internal_assert(op->extents.size() > 0);
size_id = print_assignment(Int(64), print_expr(op->extents[0]));
for (size_t i = 1; i < op->extents.size(); i++) {
// Make the code a little less cluttered for two-dimensional case
string new_size_id_rhs;
string next_extent = print_expr(op->extents[i]);
if (i > 1) {
new_size_id_rhs = "(" + size_id + " > ((int64_t(1) << 31) - 1)) ? " + size_id + " : (" + size_id + " * " + next_extent + ")";
} else {
new_size_id_rhs = size_id + " * " + next_extent;
}
size_id = print_assignment(Int(64), new_size_id_rhs);
}
do_indent();
stream << "if ((" << size_id << " > ((int64_t(1) << 31) - 1)) || ((" << size_id <<
" * sizeof(" << op_type << ")) > ((int64_t(1) << 31) - 1)))\n";
open_scope();
do_indent();
// TODO: call halide_error_buffer_allocation_too_large() here instead
// TODO: call create_assertion() so that NoAssertions works
stream << "halide_error(_ucon, "
<< "\"32-bit signed overflow computing size of allocation " << op->name << "\\n\");\n";
do_indent();
stream << "return -1;\n";
close_scope("overflow test " + op->name);
}
// Check the condition to see if this allocation should actually be created.
// If the allocation is on the stack, the only condition we can respect is
// unconditional false (otherwise a non-constant-sized array declaration
// will be generated).
if (!on_stack || is_zero(op->condition)) {
Expr conditional_size = Select::make(op->condition,
Var(size_id),
Expr(static_cast<int32_t>(0)));
conditional_size = simplify(conditional_size);
size_id = print_assignment(Int(64), print_expr(conditional_size));
}
Allocation alloc;
alloc.type = op->type;
allocations.push(op->name, alloc);
do_indent();
stream << op_type;
if (on_stack) {
stream << op_name
<< "[" << size_id << "];\n";
} else {
stream << "*"
<< op_name
<< " = ("
<< op_type
<< " *)halide_malloc(_ucon, sizeof("
<< op_type
<< ")*" << size_id << ");\n";
heap_allocations.push(op->name);
}
}
if (!on_stack) {
create_assertion(op_name, "halide_error_out_of_memory(_ucon)");
do_indent();
string free_function = op->free_function.empty() ? "halide_free" : op->free_function;
stream << "HalideFreeHelper " << op_name << "_free(_ucon, "
<< op_name << ", " << free_function << ");\n";
}
op->body.accept(this);
// Should have been freed internally
internal_assert(!allocations.contains(op->name));
close_scope("alloc " + print_name(op->name));
}
void CodeGen_C::visit(const Free *op) {
if (heap_allocations.contains(op->name)) {
do_indent();
stream << print_name(op->name) << "_free.free();\n";
heap_allocations.pop(op->name);
}
allocations.pop(op->name);
}
void CodeGen_C::visit(const Realize *op) {
internal_error << "Cannot emit realize statements to C\n";
}
void CodeGen_C::visit(const Prefetch *op) {
internal_error << "Cannot emit prefetch statements to C\n";
}
void CodeGen_C::visit(const IfThenElse *op) {
string cond_id = print_expr(op->condition);
do_indent();
stream << "if (" << cond_id << ")\n";
open_scope();
op->then_case.accept(this);
close_scope("if " + cond_id);
if (op->else_case.defined()) {
do_indent();
stream << "else\n";
open_scope();
op->else_case.accept(this);
close_scope("if " + cond_id + " else");
}
}
void CodeGen_C::visit(const Evaluate *op) {
if (is_const(op->value)) return;
string id = print_expr(op->value);
do_indent();
stream << "(void)" << id << ";\n";
}
void CodeGen_C::visit(const Shuffle *op) {
internal_assert(op->vectors.size() >= 1);
internal_assert(op->vectors[0].type().is_vector());
for (size_t i = 1; i < op->vectors.size(); i++) {
internal_assert(op->vectors[0].type() == op->vectors[i].type());
}
internal_assert(op->type.lanes() == (int) op->indices.size());
const int max_index = (int) (op->vectors[0].type().lanes() * op->vectors.size());
for (int i : op->indices) {
internal_assert(i >= -1 && i < max_index);
}
std::vector<string> vecs;
for (Expr v : op->vectors) {
vecs.push_back(print_expr(v));
}
string src = vecs[0];
if (op->vectors.size() > 1) {
ostringstream rhs;
string storage_name = unique_name('_');
do_indent();
stream << "const " << print_type(op->vectors[0].type()) << " " << storage_name << "[] = { " << with_commas(vecs) << " };\n";
rhs << print_type(op->type) << "::concat(" << op->vectors.size() << ", " << storage_name << ")";
src = print_assignment(op->type, rhs.str());
}
ostringstream rhs;
if (op->type.is_scalar()) {
rhs << src << "[" << op->indices[0] << "]";
} else {
string indices_name = unique_name('_');
do_indent();
stream << "const int32_t " << indices_name << "[" << op->indices.size() << "] = { " << with_commas(op->indices) << " };\n";
rhs << print_type(op->type) << "::shuffle(" << src << ", " << indices_name << ")";
}
print_assignment(op->type, rhs.str());
}
void CodeGen_C::test() {
LoweredArgument buffer_arg("buf", Argument::OutputBuffer, Int(32), 3, ArgumentEstimates{});
LoweredArgument float_arg("alpha", Argument::InputScalar, Float(32), 0, ArgumentEstimates{});
LoweredArgument int_arg("beta", Argument::InputScalar, Int(32), 0, ArgumentEstimates{});
LoweredArgument user_context_arg("__user_context", Argument::InputScalar, type_of<const void*>(), 0, ArgumentEstimates{});
vector<LoweredArgument> args = { buffer_arg, float_arg, int_arg, user_context_arg };
Var x("x");
Param<float> alpha("alpha");
Param<int> beta("beta");
Expr e = Select::make(alpha > 4.0f, print_when(x < 1, 3), 2);
Stmt s = Store::make("buf", e, x, Parameter(), const_true(), ModulusRemainder());
s = LetStmt::make("x", beta+1, s);
s = Block::make(s, Free::make("tmp.stack"));
s = Allocate::make("tmp.stack", Int(32), MemoryType::Stack, {127}, const_true(), s);
s = Block::make(s, Free::make("tmp.heap"));
s = Allocate::make("tmp.heap", Int(32), MemoryType::Heap, {43, beta}, const_true(), s);
Expr buf = Variable::make(Handle(), "buf.buffer");
s = LetStmt::make("buf", Call::make(Handle(), Call::buffer_get_host, {buf}, Call::Extern), s);
Module m("", get_host_target());
m.append(LoweredFunc("test1", args, s, LinkageType::External));
ostringstream source;
{
CodeGen_C cg(source, Target("host"), CodeGen_C::CImplementation);
cg.compile(m);
}
string src = source.str();
string correct_source =
headers +
globals +
string((const char *)halide_internal_runtime_header_HalideRuntime_h) + '\n' +
string((const char *)halide_internal_initmod_inlined_c) + R"GOLDEN_CODE(
#ifndef HALIDE_FUNCTION_ATTRS
#define HALIDE_FUNCTION_ATTRS
#endif
#ifdef __cplusplus
extern "C" {
#endif
int test1(struct halide_buffer_t *_buf_buffer, float _alpha, int32_t _beta, void const *__user_context) HALIDE_FUNCTION_ATTRS {
void * const _ucon = const_cast<void *>(__user_context);
void *_0 = _halide_buffer_get_host(_buf_buffer);
void * _buf = _0;
{
int64_t _1 = 43;
int64_t _2 = _1 * _beta;
if ((_2 > ((int64_t(1) << 31) - 1)) || ((_2 * sizeof(int32_t )) > ((int64_t(1) << 31) - 1)))
{
halide_error(_ucon, "32-bit signed overflow computing size of allocation tmp.heap\n");
return -1;
} // overflow test tmp.heap
int64_t _3 = _2;
int32_t *_tmp_heap = (int32_t *)halide_malloc(_ucon, sizeof(int32_t )*_3);
if (!_tmp_heap)
{
return halide_error_out_of_memory(_ucon);
}
HalideFreeHelper _tmp_heap_free(_ucon, _tmp_heap, halide_free);
{
int32_t _tmp_stack[127];
int32_t _4 = _beta + 1;
int32_t _5;
bool _6 = _4 < 1;
if (_6)
{
char b0[1024];
snprintf(b0, 1024, "%lld%s", (long long)(3), "\n");
char const *_7 = b0;
int32_t _8 = halide_print(_ucon, _7);
int32_t _9 = return_second(_8, 3);
_5 = _9;
} // if _6
else
{
_5 = 3;
} // if _6 else
int32_t _10 = _5;
float _11 = float_from_bits(1082130432 /* 4 */);
bool _12 = _alpha > _11;
int32_t _13 = (int32_t)(_12 ? _10 : 2);
((int32_t *)_buf)[_4] = _13;
} // alloc _tmp_stack
_tmp_heap_free.free();
} // alloc _tmp_heap
return 0;
}
#ifdef __cplusplus
} // extern "C"
#endif
)GOLDEN_CODE";
if (src != correct_source) {
int diff = 0;
while (src[diff] == correct_source[diff]) diff++;
int diff_end = diff + 1;
while (diff > 0 && src[diff] != '\n') diff--;
while (diff_end < (int)src.size() && src[diff_end] != '\n') diff_end++;
internal_error
<< "Correct source code:\n" << correct_source
<< "Actual source code:\n" << src
<< "Difference starts at:\n"
<< "Correct: " << correct_source.substr(diff, diff_end - diff) << "\n"
<< "Actual: " << src.substr(diff, diff_end - diff) << "\n";
}
std::cout << "CodeGen_C test passed\n";
}
} // namespace Internal
} // namespace Halide