I'm working through things related to #203 and #204 and I believe I have identified the principle problem, which is that there is some behavioral difference in the initialization of center_to_polar_bound_vectors_ and center_to_azimuthal_bound_vectors_ when num_polar_sections == num_azimuthal_sections and when the converse. The more correct (appearing) behavior is when they are equal, but I'm not sure I can completely convince myself of the correctness of that. The thing that gives me the most pause is that in the case that the count is equal, I don't see why it would get the correct answer for the values for both polar and azimuthal, since they span different values.
Anyway, I'm working on it, and have a partially completed pull request in the works.