src/panfrost/midgard/midgard_nir_algebraic.py

   1 #
   2 # Copyright (C) 2018 Alyssa Rosenzweig
   3 #
   4 # Copyright (C) 2016 Intel Corporation
   5 #
   6 # Permission is hereby granted, free of charge, to any person obtaining a
   7 # copy of this software and associated documentation files (the "Software"),
   8 # to deal in the Software without restriction, including without limitation
   9 # the rights to use, copy, modify, merge, publish, distribute, sublicense,
  10 # and/or sell copies of the Software, and to permit persons to whom the
  11 # Software is furnished to do so, subject to the following conditions:
  12 #
  13 # The above copyright notice and this permission notice (including the next
  14 # paragraph) shall be included in all copies or substantial portions of the
  15 # Software.
  16 #
  17 # THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
  18 # IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
  19 # FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.  IN NO EVENT SHALL
  20 # THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
  21 # LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
  22 # FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS
  23 # IN THE SOFTWARE.
  24
  25 import argparse
  26 import sys
  27 import math
  28
  29 a = 'a'
  30 b = 'b'
  31 c = 'c'
  32
  33 algebraic = [
  34    (('pack_unorm_4x8', a), ('pack_32_4x8', ('f2u8', ('fround_even', ('fmul', ('fsat', a), 255.0)))))
  35 ]
  36
  37 algebraic_late = [
  38     # ineg must be lowered late, but only for integers; floats will try to
  39     # have modifiers attached... hence why this has to be here rather than
  40     # a more standard lower_negate approach
  41
  42     (('ineg', a), ('isub', 0, a)),
  43
  44     # Likewise we want fsub lowered but not isub
  45     (('fsub', a, b), ('fadd', a, ('fneg', b))),
  46
  47     # These two special-cases save space/an op than the actual csel op +
  48     # scheduler flexibility
  49
  50     (('b32csel', a, 'b@32', 0), ('iand', a, b)),
  51     (('b32csel', a, 0, 'b@32'), ('iand', ('inot', a), b)),
  52
  53     # Fuse sat_signed. This should probably be shared with Bifrost
  54     (('~fmin', ('fmax', a, -1.0), 1.0), ('fsat_signed', a)),
  55     (('~fmax', ('fmin', a, 1.0), -1.0), ('fsat_signed', a)),
  56
  57     # Fuse clamp_positive. This should probably be shared with Utgard/bifrost
  58     (('fmax', a, 0.0), ('fclamp_pos', a)),
  59
  60     (('ishl', 'a@16', b), ('u2u16', ('ishl', ('u2u32', a), b))),
  61     (('ishr', 'a@16', b), ('i2i16', ('ishr', ('i2i32', a), b))),
  62     (('ushr', 'a@16', b), ('u2u16', ('ushr', ('u2u32', a), b))),
  63
  64     (('ishl', 'a@8', b), ('u2u8', ('u2u16', ('ishl', ('u2u32', ('u2u16', a)), b)))),
  65     (('ishr', 'a@8', b), ('i2i8', ('i2i16', ('ishr', ('i2i32', ('i2i16', a)), b)))),
  66     (('ushr', 'a@8', b), ('u2u8', ('u2u16', ('ushr', ('u2u32', ('u2u16', a)), b)))),
  67
  68     # Canonical form. The scheduler will convert back if it makes sense.
  69     (('fmul', a, 2.0), ('fadd', a, a))
  70 ]
  71
  72
  73 # Midgard is able to type convert down by only one "step" per instruction; if
  74 # NIR wants more than one step, we need to break up into multiple instructions.
  75 # Nevertheless, we can do both a size step and a floating/int step at once.
  76
  77 converts = []
  78
  79 for op in ('u2u', 'i2i', 'f2f', 'i2f', 'u2f', 'f2i', 'f2u'):
  80     srcsz_max = 64
  81     dstsz_max = 64
  82     # 8 bit float doesn't exist
  83     srcsz_min = 8 if op[0] != 'f' else 16
  84     dstsz_min = 8 if op[2] != 'f' else 16
  85     dstsz = dstsz_min
  86     # Iterate over all possible destination and source sizes
  87     while dstsz <= dstsz_max:
  88         srcsz = srcsz_min
  89         while srcsz <= srcsz_max:
  90             # Size converter lowering is only needed if src and dst sizes are
  91             # spaced by a factor > 2.
  92             if srcsz != dstsz and (srcsz * 2 != dstsz and srcsz != dstsz * 2):
  93                 cursz = srcsz
  94                 rule = a
  95                 # When converting down we first do the type conversion followed
  96                 # by one or more size conversions. When converting up, we do
  97                 # the type conversion at the end. This way we don't have to
  98                 # deal with the fact that f2f8 doesn't exists.
  99                 sizeconvop = op[0] + '2' + op[0] if srcsz < dstsz else op[2] + '2' + op[2]
 100                 if srcsz > dstsz and op[0] != op[2]:
 101                     rule = (op + str(int(cursz)), rule)
 102                 while cursz != dstsz:
 103                     cursz = cursz / 2 if dstsz < srcsz else cursz * 2
 104                     rule = (sizeconvop + str(int(cursz)), rule)
 105                 if srcsz < dstsz and op[0] != op[2]:
 106                     rule = (op + str(int(cursz)), rule)
 107                 converts += [((op + str(int(dstsz)), 'a@' + str(int(srcsz))), rule)]
 108             srcsz *= 2
 109         dstsz *= 2
 110
 111 # Try to force constants to the right
 112 constant_switch = [
 113         # fge gets flipped to fle, so we invert to keep the order
 114         (('fge', 'a', '#b'), (('inot', ('flt', a, b)))),
 115         (('fge32', 'a', '#b'), (('inot', ('flt32', a, b)))),
 116         (('ige32', 'a', '#b'), (('inot', ('ilt32', a, b)))),
 117         (('uge32', 'a', '#b'), (('inot', ('ult32', a, b)))),
 118
 119         # fge gets mapped to fle with a flip
 120         (('flt32', '#a', 'b'), ('inot', ('fge32', a, b))),
 121         (('ilt32', '#a', 'b'), ('inot', ('ige32', a, b))),
 122         (('ult32', '#a', 'b'), ('inot', ('uge32', a, b)))
 123 ]
 124
 125 # ..since the above switching happens after algebraic stuff is done
 126 cancel_inot = [
 127         (('inot', ('inot', a)), a)
 128 ]
 129
 130 # Midgard scales fsin/fcos arguments by pi.
 131 # Pass must be run only once, after the main loop
 132
 133 scale_trig = [
 134         (('fsin', a), ('fsin', ('fdiv', a, math.pi))),
 135         (('fcos', a), ('fcos', ('fdiv', a, math.pi))),
 136 ]
 137
 138 def main():
 139     parser = argparse.ArgumentParser()
 140     parser.add_argument('-p', '--import-path', required=True)
 141     args = parser.parse_args()
 142     sys.path.insert(0, args.import_path)
 143     run()
 144
 145
 146 def run():
 147     import nir_algebraic  # pylint: disable=import-error
 148
 149     print('#include "midgard_nir.h"')
 150
 151     print(nir_algebraic.AlgebraicPass("midgard_nir_lower_algebraic_early",
 152                                       algebraic).render())
 153
 154     print(nir_algebraic.AlgebraicPass("midgard_nir_lower_algebraic_late",
 155                                       algebraic_late + converts + constant_switch).render())
 156
 157     print(nir_algebraic.AlgebraicPass("midgard_nir_scale_trig",
 158                                       scale_trig).render())
 159
 160     print(nir_algebraic.AlgebraicPass("midgard_nir_cancel_inot",
 161                                       cancel_inot).render())
 162
 163
 164 if __name__ == '__main__':
 165     main()