I am not sure you need to flatten (you are not using any internal structure of bc yet in your code), preprocess is probably needed (to do unaliasing). Now for 1:length(r) is probably safe, but a normal thing would be to use axes(bcf2, 1) (but I think it should be the same in this case as you are not accepting non standard indexing in bcf2 anyway - right?).
2 Likes