Mercurial > vec
diff gen/gendouble.c @ 46:31cee67540b5
f32/f64: add floor, ceil, round, and trunc operations
we also need saturated add sub mul etc
author | Paper <paper@tflc.us> |
---|---|
date | Sat, 09 Aug 2025 15:55:59 -0400 |
parents | 7955bed1d169 |
children |
line wrap: on
line diff
--- a/gen/gendouble.c Wed Apr 30 18:36:38 2025 -0400 +++ b/gen/gendouble.c Sat Aug 09 15:55:59 2025 -0400 @@ -24,123 +24,6 @@ #include "genlib.h" -#define ARRAY_SIZE(x) (sizeof(x)/sizeof((x)[0])) - -/* XXX: would it be faster to unroll literally everything instead of defining everything, - * and then unpacking it all? */ -static const char *header = - "/**\n" - " * vec - a tiny SIMD vector library in C99\n" - " * \n" - " * Copyright (c) 2024-2025 Paper\n" - " * \n" - " * Permission is hereby granted, free of charge, to any person obtaining a copy\n" - " * of this software and associated documentation files (the \"Software\"), to deal\n" - " * in the Software without restriction, including without limitation the rights\n" - " * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell\n" - " * copies of the Software, and to permit persons to whom the Software is\n" - " * furnished to do so, subject to the following conditions:\n" - " * \n" - " * The above copyright notice and this permission notice shall be included in all\n" - " * copies or substantial portions of the Software.\n" - " * \n" - " * THE SOFTWARE IS PROVIDED \"AS IS\", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR\n" - " * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,\n" - " * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE\n" - " * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER\n" - " * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,\n" - " * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE\n" - " * SOFTWARE.\n" - "**/\n" - "\n" - "/* This file is automatically generated! Do not edit it directly!\n" - " * Edit the code that generates it in utils/gendbl.c --paper */\n" - "\n" - "#ifndef VEC_IMPL_DOUBLE_H_\n" - "#define VEC_IMPL_DOUBLE_H_\n" - "\n" - "#define VEC_DOUBLE_SPLAT(sign, bits, size, halfsize) \\\n" - " VEC_FUNC_IMPL v##sign##int##bits##x##size v##sign##int##bits##x##size##_splat(vec_##sign##int##bits x) \\\n" - " { \\\n" - " v##sign##int##bits##x##size vec; \\\n" - " \\\n" - " vec.dbl[0] = v##sign##int##bits##x##halfsize##_splat(x); \\\n" - " vec.dbl[1] = v##sign##int##bits##x##halfsize##_splat(x); \\\n" - " \\\n" - " return vec; \\\n" - " }\n" - "\n" - "#define VEC_DOUBLE_LOAD_EX(name, sign, bits, size, halfsize) \\\n" - " VEC_FUNC_IMPL v##sign##int##bits##x##size v##sign##int##bits##x##size##_##name(const vec_##sign##int##bits x[size]) \\\n" - " { \\\n" - " v##sign##int##bits##x##size vec; \\\n" - " \\\n" - " vec.dbl[0] = v##sign##int##bits##x##halfsize##_##name(x); \\\n" - " vec.dbl[1] = v##sign##int##bits##x##halfsize##_##name(x + halfsize); \\\n" - " \\\n" - " return vec; \\\n" - " }\n" - "\n" - "#define VEC_DOUBLE_LOAD(sign, bits, size, halfsize) VEC_DOUBLE_LOAD_EX(load, sign, bits, size, halfsize)\n" - "#define VEC_DOUBLE_LOAD_ALIGNED(sign, bits, size, halfsize) VEC_DOUBLE_LOAD_EX(load_aligned, sign, bits, size, halfsize)\n" - "\n" - "#define VEC_DOUBLE_STORE_EX(name, sign, bits, size, halfsize) \\\n" - " VEC_FUNC_IMPL void v##sign##int##bits##x##size##_##name(v##sign##int##bits##x##size vec, vec_##sign##int##bits x[size]) \\\n" - " { \\\n" - " v##sign##int##bits##x##halfsize##_##name(vec.dbl[0], x); \\\n" - " v##sign##int##bits##x##halfsize##_##name(vec.dbl[1], x + halfsize); \\\n" - " }\n" - "\n" - "#define VEC_DOUBLE_STORE(sign, bits, size, halfsize) VEC_DOUBLE_STORE_EX(store, sign, bits, size, halfsize)\n" - "#define VEC_DOUBLE_STORE_ALIGNED(sign, bits, size, halfsize) VEC_DOUBLE_STORE_EX(store_aligned, sign, bits, size, halfsize)\n" - "\n" - "#define VEC_DOUBLE_OP(name, sign, bits, size, halfsize, secondsign) \\\n" - " VEC_FUNC_IMPL v##sign##int##bits##x##size v##sign##int##bits##x##size##_##name(v##sign##int##bits##x##size vec1, v##secondsign##int##bits##x##size vec2) \\\n" - " { \\\n" - " vec1.dbl[0] = v##sign##int##bits##x##halfsize##_##name(vec1.dbl[0], vec2.dbl[0]); \\\n" - " vec1.dbl[1] = v##sign##int##bits##x##halfsize##_##name(vec1.dbl[1], vec2.dbl[1]); \\\n" - " \\\n" - " return vec1; \\\n" - " }\n" - "\n" - "#define VEC_DOUBLE_ADD(sign, bits, size, halfsize) VEC_DOUBLE_OP(add, sign, bits, size, halfsize, sign)\n" - "#define VEC_DOUBLE_SUB(sign, bits, size, halfsize) VEC_DOUBLE_OP(sub, sign, bits, size, halfsize, sign)\n" - "#define VEC_DOUBLE_MUL(sign, bits, size, halfsize) VEC_DOUBLE_OP(mul, sign, bits, size, halfsize, sign)\n" - "#define VEC_DOUBLE_DIV(sign, bits, size, halfsize) VEC_DOUBLE_OP(div, sign, bits, size, halfsize, sign)\n" - "#define VEC_DOUBLE_MOD(sign, bits, size, halfsize) VEC_DOUBLE_OP(mod, sign, bits, size, halfsize, sign)\n" - "#define VEC_DOUBLE_AVG(sign, bits, size, halfsize) VEC_DOUBLE_OP(avg, sign, bits, size, halfsize, sign)\n" - "#define VEC_DOUBLE_LSHIFT(sign, bits, size, halfsize) VEC_DOUBLE_OP(lshift, sign, bits, size, halfsize, u)\n" - "#define VEC_DOUBLE_RSHIFT(sign, bits, size, halfsize) VEC_DOUBLE_OP(rshift, sign, bits, size, halfsize, u)\n" - "#define VEC_DOUBLE_LRSHIFT(sign, bits, size, halfsize) VEC_DOUBLE_OP(lrshift, sign, bits, size, halfsize, u)\n" - "#define VEC_DOUBLE_AND(sign, bits, size, halfsize) VEC_DOUBLE_OP(and, sign, bits, size, halfsize, sign)\n" - "#define VEC_DOUBLE_OR(sign, bits, size, halfsize) VEC_DOUBLE_OP(or, sign, bits, size, halfsize, sign)\n" - "#define VEC_DOUBLE_XOR(sign, bits, size, halfsize) VEC_DOUBLE_OP(xor, sign, bits, size, halfsize, sign)\n" - "#define VEC_DOUBLE_MIN(sign, bits, size, halfsize) VEC_DOUBLE_OP(min, sign, bits, size, halfsize, sign)\n" - "#define VEC_DOUBLE_MAX(sign, bits, size, halfsize) VEC_DOUBLE_OP(max, sign, bits, size, halfsize, sign)\n" - "#define VEC_DOUBLE_CMPLT(sign, bits, size, halfsize) VEC_DOUBLE_OP(cmplt, sign, bits, size, halfsize, sign)\n" - "#define VEC_DOUBLE_CMPLE(sign, bits, size, halfsize) VEC_DOUBLE_OP(cmple, sign, bits, size, halfsize, sign)\n" - "#define VEC_DOUBLE_CMPEQ(sign, bits, size, halfsize) VEC_DOUBLE_OP(cmpeq, sign, bits, size, halfsize, sign)\n" - "#define VEC_DOUBLE_CMPGE(sign, bits, size, halfsize) VEC_DOUBLE_OP(cmpge, sign, bits, size, halfsize, sign)\n" - "#define VEC_DOUBLE_CMPGT(sign, bits, size, halfsize) VEC_DOUBLE_OP(cmpgt, sign, bits, size, halfsize, sign)\n" - "\n" - "#define VEC_DOUBLE_NOT(sign, bits, size, halfsize) \\\n" - " VEC_FUNC_IMPL v##sign##int##bits##x##size v##sign##int##bits##x##size##_not(v##sign##int##bits##x##size vec) \\\n" - " { \\\n" - " vec.dbl[0] = v##sign##int##bits##x##halfsize##_not(vec.dbl[0]); \\\n" - " vec.dbl[1] = v##sign##int##bits##x##halfsize##_not(vec.dbl[1]); \\\n" - " \\\n" - " return vec; \\\n" - " }\n" - "\n" - "#endif /* VEC_IMPL_DOUBLE_H_ */ \n" - "\n" - "/* ------------------------------------------------------------------------ */\n" - "/* PREPROCESSOR HELL INCOMING */\n" - ""; - -static const char *footer = - "" /* nothing */; - /* ------------------------------------------------------------------------ */ static void op_print_pp_halfsize(int op, int type, int bits, int size) @@ -173,7 +56,7 @@ gen_print_vtype(type, bits, size / 2); printf("_%s(vec.dbl[0]);\n", op_info->l); - printf("\tvec1.dbl[1] = "); + printf("\tvec.dbl[1] = "); gen_print_vtype(type, bits, size / 2); printf("_%s(vec.dbl[1]);\n", op_info->l); @@ -267,6 +150,11 @@ [OP_CMPEQ] = {NULL, op_print_pp_halfsize, op_print_twoop}, [OP_CMPGE] = {NULL, op_print_pp_halfsize, op_print_twoop}, [OP_CMPGT] = {NULL, op_print_pp_halfsize, op_print_twoop}, + + [OP_FLOOR] = {NULL, op_print_pp_halfsize, op_print_unoop}, + [OP_CEIL] = {NULL, op_print_pp_halfsize, op_print_unoop}, + [OP_ROUND] = {NULL, op_print_pp_halfsize, op_print_unoop}, + [OP_TRUNC] = {NULL, op_print_pp_halfsize, op_print_unoop}, }; int main(void)