Superconduting wire current-phase relation¶

We want to calculate the current-phase relation of a superconducting wire. The relation is multivalued, and so we use pseudo-arclength continuation to generate the $(\varphi, I)$ trace. Since it is an equilibrium calculation, we do it in imaginary time.