Bug in generate_pack_info_from_kernel
generate_pack_info_from_kernel throws the error "Slice dimension does not match sizes" sometimes. Basically this happens when in the forcing terms also field accesses to other fields are done. I recreated the situation by applying the gradient of one field as forcing terms to an lattice Boltzmann update rule
Important note:
The error can only be reproduced when cuda is set to true in the code generation context. This means with the standard ManualCodeGenerationContext it does not work. I just set cuda True in the ManualCodeGenerationContext to reproduce the error. Then the error can be seen in the following notebook. Bug_in_generate_pack_info_from_kernel.ipynb
Side note: I tried to run the whole notebook in loop to be able to reproduce the error more often. Apparently the error appears either in the first loop iteration or never (100 iterations). So restart and run all was the only change for me to see the error in about one of three runs.