This will be tricky to implement on GPU but can be done on C++/Python side, i.e. simulated BOLD is calculated on GPU, returned to CPU [C++/Python] to be filtered, and then FC and FCD are again calculated on GPU. To do this the mean_bold and ssd_bold calculations should not be done within the bnm kernel (which is the current behavior), but rather should be done in a separate kernel after filtering.
This will be tricky to implement on GPU but can be done on C++/Python side, i.e. simulated BOLD is calculated on GPU, returned to CPU [C++/Python] to be filtered, and then FC and FCD are again calculated on GPU. To do this the
mean_boldandssd_boldcalculations should not be done within thebnmkernel (which is the current behavior), but rather should be done in a separate kernel after filtering.