Okay, so I have started my teaching duties for the summer, I have also started working with the original code to add inverse-Compton scatter--This may be helpful for finally finishing up the model fit of GRB 070125: The original reason I started playing with boxfit . Both of these have slowed progess on the CUDA port, so it is behind where I wanted it to be in the original roadmap. I do bear some good news, though. I think I understand why the static is causing such a big problem with CUDA. Originally, I was confused because I had overloaded the new and delete operators to automatically move any declared variables to device memory using CUDAMallocManaged() (I know this is not an efficient means of allocated device memory, but it is useful for getting a prototype up quickly.) The problem is that static members of a class are identical across all instances of the class, so there is no way to seperately call the array fo...