ENH: Add fitting of current sheets

The `fit()` function only takes observations of B right now. I don't see why we can't do a similar method that we have in `predict()` and have both a `fit_B` and `fit_J` function. I think this would just require the T matrix being the J transform instead.