Rory Conlin

could do something like: ```python for coil in coils: errorif(not isinstance(coil, _Coil), TypeError, f"Expected object of type Coil, got {type(coil)}", ) ``` that way it actually says what the offending...

actually im now doubting myself. Can plasma current produce toroidal flux? i thought that wasn't possible (or maybe it just cant produce "enough" flux?)

div/grad/curl in jax (assumes cartesian coordinates, will need some extra stuff to work in cylindrical) ```python def grad(fun, argnum): return jax.jacfwd(fun, argnum) def div(fun, argnum): def divf(*args, **kwargs): return jnp.trace(jax.jacfwd(fun,...

not really sure what you mean by > `factorize_linear_constraints` normalizes the optimization variables `x_reduced` to order unity,

I think this might be as simple as redefining things as ``` x = xp + Z @ D @ x_reduced x_reduced = D^-1 @ Z.T @ (x - xp)...

In general the SVD is unique up to a permutation of the columns of $U$ and $V$. If we enforce that the singular values are ordered in decreasing magnitude, this...

Might be more useful as a general IOAble family for arbitrary things, not just coils. And could make EquilibriaFamily a subclass

I don't think we should get rid of it, since that would likely break every existing saved `EquilibriaFamily`

Another thought: how useful are circular cross sections in practice? Do people use them? Would just rectangular be enough? (I would assume rectangular is closer to reality?)

Yeah, I think the order would be something like ```python Coil(MagneticField): # base class FramedCoil(Coil): # coil plus rotating reference frame, but no width/height - still filamentary FiniteBuildCoil(FramedCoil): # as...