def METHOD_NAME(ph_nacl_gruneisen):
    """Test of mode Grueneisen parameter calculation along band paths."""
    paths = get_band_qpoints([[[0.05, 0.05, 0.05], [0.5, 0.5, 0.5]]], npoints=10)
    ph0, ph_minus, ph_plus = ph_nacl_gruneisen
    phg = PhonopyGruneisen(ph0, ph_minus, ph_plus)
    phg.set_band_structure(paths)
    qpoints, distances, freqs, _, gammas = phg.get_band_structure()
    # print(", ".join(["%.5f" % v for v in distances[0]]))
    np.testing.assert_allclose(distances[0], band_distances, atol=1e-5)
    np.testing.assert_allclose(qpoints[0], band_qpoints, atol=1e-5)
    # for line in freqs[0].reshape(-1, 6):
    #     print(", ".join(["%.5f" % v for v in line]))
    for freqs_q, band_freqs_q in zip(freqs[0], band_freqs):
        np.testing.assert_allclose(np.sort(freqs_q), np.sort(band_freqs_q), atol=1e-5)
    # for line in gammas[0].reshape(-1, 6):
    #     print(", ".join(["%.5f" % v for v in line]))
    for gammas_q, band_gammas_q in zip(gammas[0], band_gammas):
        np.testing.assert_allclose(np.sort(gammas_q), np.sort(band_gammas_q), atol=1e-5)