loxodrome_direct stability near asymptote

Thanks to @noritada for finding and making tests for loxodrome_direct issues. Always appreciate hearing where the test suite isn't covering a broken case.

used numpy.broadcast_arrays to simplify code.