How do feel about making this a method on the CSD model? It can take a gradient table and FOD coefficients as arguments, and it'll get the response function (R matrix) from the model.