Skip to content

Bug in element coordinate replacement if cubes loop pattern is chosen

The function _compute_micro_element_coordinates computes the element coordinates from the loop counters. To enable vectorization, the replacement is currently done in a custom code node.

The issue arises during loop cutting for the cubes loop, since loops with end - start are automatically replaced with their loop bodies by pystencils. However, the custom code node requires the presence of the loop counter, which has been removed through the loop cutting.

A possible solution would be to not remove these short loops during cutting.