def METHOD_NAME(filename, dm_name):
    """Get trace of spin up and down 1RDM as function tau. Useful for debugging.
    Parameters
    ----------
    filename : string
        QMCPACK output containing density matrix (*.h5 file).
    dm_name: string
        Estimator type. Mixed or BackPropagated. Default: Mixed.
    Returns
    -------
    dmu : :class:`numpy.ndarray`
        Trace of spin up QMC 1RDM as function of sample block.
    dmd : :class:`numpy.ndarray`
        Trace of spin down QMC 1RDM as function of sample block.
    """
    (rdm, weights) = extract_qmc_dm(filename, dm_name)
    dmu = []
    dmd = []
    for dm in rdm:
        dmu.append(dm[0].trace())
        dmd.append(dm[1].trace())
    return (dmu, dmd)