Add SVE nontemporal stores and scatters, including masked variants

@holzer pointed out that SVE has nontemporal stores, which I overlooked when I implemented SVE support three years ago. So there are actually eight different kinds of stores we have to support (store, stream, mask-store, mask-stream, scatter, stream-scatter, mask-scatter, mask-stream-scatter).

(Mask-)stream-scatter requires SVE2; adding support for that introduced changes to a number of unrelated files. SVE2 is a superset of SVE, so the automatic detection won‘t return SVE if SVE2 is also available.

The added test coverage revealed a few bugs on master that I also fixed:

• We used an incorrect argument order for maskStoreS on RISC-V-V. This meant that maskStoreS just wouldn't compile on RISC-V-V.
• Our emulation of maskStore (via blendv) on ARM Neon and POWER VSX was zeroing the masked-out elements when nontemporal mode was selected (i.e. which due to the lack of real nontemporal stores maps to an emulation via cacheline zeroing). This lead to incorrect results.
• Our emulation of maskStore on POWER VSX was ignoring the mask for the last vector of each cacheline. This lead to even more incorrect results when nontemporal mode was selected.

src/pystencils/backends/cbackend.py

@@ -280,9 +280,11 @@ class CBackend:
             if type(lhs_type) is VectorType and isinstance(node.lhs, CastFunc):
                 arg, data_type, aligned, nontemporal, mask, stride = node.lhs.args
                 instr = 'storeU'
-                if aligned:
-                    instr = 'stream' if nontemporal and 'stream' in self._vector_instruction_set else \
-                        'storeA' if 'storeA' in self._vector_instruction_set else 'storeU'
+                if nontemporal and 'storeA' not in self._vector_instruction_set and \
+                        'stream' in self._vector_instruction_set:
+                    instr = 'stream'
+                elif aligned:
+                    instr = 'stream' if nontemporal and 'stream' in self._vector_instruction_set else 'storeA'
                 if mask != True:  # NOQA
                     instr = 'maskStream' if nontemporal and 'maskStream' in self._vector_instruction_set else \
                             'maskStoreA' if aligned else 'maskStoreU'
@@ -319,12 +321,14 @@ class CBackend:
 
                 if stride != 1:
                     instr = ('maskStreamS' if nontemporal and 'maskStreamS' in self._vector_instruction_set else
-                             'maskStoreS') if mask != True else ('streamS' if nontemporal else 'storeS') # NOQA
+                             'maskStoreS') if mask != True else \
+                            ('streamS' if nontemporal and 'streamS' in self._vector_instruction_set else 'storeS')  # NOQA
                     return self._vector_instruction_set[instr].format(ptr, self.sympy_printer.doprint(rhs),
                                                                       stride, printed_mask, **self._kwargs) + ';'
 
                 pre_code = ''
-                if nontemporal and 'cachelineZero' in self._vector_instruction_set:
+                blended = self._vector_instruction_set['blendv'].split('{')[0] in self._vector_instruction_set[instr]
+                if nontemporal and 'cachelineZero' in self._vector_instruction_set and not blended:
                     first_cond = f"((uintptr_t) {ptr} & {CachelineSize.mask_symbol}) == 0"
                     offset = sp.Add(*[sp.Symbol(LoopOverCoordinate.get_loop_counter_name(i))
                                       * node.lhs.args[0].field.spatial_strides[i] for i in
@@ -701,8 +705,7 @@ class VectorizedCustomSympyPrinter(CustomSympyPrinter):
             arg, data_type, aligned, _, mask, stride = expr.args
             if stride != 1:
                 return self.instruction_set['loadS'].format(f"& {self._print(arg)}", stride, **self._kwargs)
-            instruction = self.instruction_set['loadA'] if aligned and 'loadA' in self.instruction_set \
-                else self.instruction_set['loadU']
+            instruction = self.instruction_set['loadA'] if aligned else self.instruction_set['loadU']
             return instruction.format(f"& {self._print(arg)}", **self._kwargs)
         elif expr.func == DivFunc:
             result = self._scalarFallback('_print_Function', expr)