Parallellize distance ansatz calculation

[ravioli1369]

This implements a parallelized root solver function to solve for the distance ansatz for all pixels at once.

However, this implementation is ~2x slower than the original scipy implementation.