def METHOD_NAME(self):     
    order = 25
    ppp = 10
    curves, currents, ma = get_ncsx_data(Nt_coils=order, ppp=ppp)
    with ScratchDir("."):
        coils_to_makegrid("coils.file_to_load", curves, currents, nfp=1)
        loaded_coils = load_coils_from_makegrid_file("coils.file_to_load", order, ppp)
    gamma = [curve.gamma() for curve in curves]
    loaded_gamma = [coil.curve.gamma() for coil in loaded_coils]
    loaded_currents = [coil.current for coil in loaded_coils]
    coils = [Coil(curve, current) for curve, current in zip(curves, currents)]
    for j_coil in range(len(coils)):
        np.testing.assert_allclose(
            currents[j_coil].get_value(),
            loaded_currents[j_coil].get_value()
        )
        np.testing.assert_allclose(curves[j_coil].x, loaded_coils[j_coil].curve.x)
    np.random.seed(1)
    bs = BiotSavart(coils)
    loaded_bs = BiotSavart(loaded_coils)
    points = np.asarray(17 * [[0.9, 0.4, -0.85]])
    points += 0.01 * (np.random.rand(*points.shape) - 0.5)
    bs.set_points(points)
    loaded_bs.set_points(points)
    B = bs.B()
    loaded_B = loaded_bs.B()
    np.testing.assert_allclose(B, loaded_B)
    np.testing.assert_allclose(gamma, loaded_gamma)