I feel like this linearity check is very messy, and might fail for some odd corner cases, I feel like you'd be better served doing matching with `Wild`, similar to the matching that is done in `classify_ode`, to get `func_coeff` and `linearity`. When I get home tonight, I will try to sit down with this tonight and work through this, to figure out what exactly this is doing.