Hello and thanks a lot for what is possibly my favorite Julia package!
Unless I am mistaken, there was a small bug introduced with the pointer access. The last element of Dkβ was stored outside the actual array, since size(Dkp1,1) == size(Dk,1) - 1
On a related note, is this pointer business still necessary for high performance? Or could one of the newer language constructs (such as SubArrays) also achieve the same result?
Hello and thanks a lot for what is possibly my favorite Julia package!
Unless I am mistaken, there was a small bug introduced with the pointer access. The last element of
Dkβ
was stored outside the actual array, sincesize(Dkp1,1) == size(Dk,1) - 1
On a related note, is this pointer business still necessary for high performance? Or could one of the newer language constructs (such as SubArrays) also achieve the same result?