ENH: interpolate in cylindrical coordinates: See #2145

This intercepts all vector/tensor AMI interpolations and does the interpolation in cylindrical coordinates.

Use (component) 'coupled' linear solver to enable this in the linear solver sweeps as well (however this is not very stable since the vector solution is still in Cartesian coordinates - only the AMI interpolation is in cylindrical)