BFieldZone(int id, double zmin, double zmax, double rmin, double rmax, double phimin, double phimax, double scale)
void addBiotSavart(const double *ATH_RESTRICT xyz, double *ATH_RESTRICT B, double *ATH_RESTRICT deriv=nullptr) const
compute Biot-Savart magnetic field and add to B[3]