def METHOD_NAME(ph_nacl_nofcsym: Phonopy):
    """Test of total DOS with tetrahedron method."""
    phonon = ph_nacl_nofcsym
    phonon.run_mesh([5, 5, 5])
    phonon.run_total_dos(freq_pitch=1, use_tetrahedron_method=True)
    dos = phonon.total_dos.dos
    freqs = phonon.total_dos.frequency_points
    data_ref = np.reshape([float(x) for x in tdos_thm_str.split()], (-1, 2))
    np.testing.assert_allclose(data_ref, np.c_[freqs, dos], atol=1e-5)
    # for f, d in zip(freqs, dos):
    #     print("%f %f" % (f, d))