Refactor gpu indexing

To map an iteration space to GPU Threads indexing classes are used. These indexing classes receive a field and iterations slice to determine the iteration space. This MR refactors the indexing classes to directly receive an iteration space. With this, the indexing classes are more general and not dependent on pystencils Fields.

Further improvements/fixes:

• Line indexing works now with iteration slices. This did not work at all before
• Both indexing schemes calculate a correct block and grid size for iteration slices. This means if for example if only every second element is touched (due to a given iteration slice) the number of threads will be half. This removes modulo calculation that was needed before
• Both indexing schemes now support up to 4 dimensions

pystencils/gpu/indexing.py

@@ -143,7 +143,7 @@ class BlockIndexing(AbstractIndexing):
                  maximum_block_size=(1024, 1024, 64), device_number=None):
         super(BlockIndexing, self).__init__(iteration_space, data_layout)
 
-        if self._dim < 4:
+        if self._dim > 4:
             raise NotImplementedError("This indexing scheme supports at most 4 spatial dimensions")
 
         if permute_block_size_dependent_on_layout and self._dim < 4: