add_in_contribs_to_energy Subroutine

    subroutine add_in_contribs_to_energy()

        integer :: i, j
        real(dp) :: beta

        if (iProcIndex /= root) return

        do i = 1, n_lanc_vecs_ftlm
            beta = 0.0_dp
            do j = 1, nbeta_ftlm + 1
                ftlm_trace(j) = ftlm_trace(j) + (ftlm_hamil(1, i)**2) * exp(-beta * ftlm_h_eigv(i))
                ftlm_e_num(j) = ftlm_e_num(j) + (ftlm_hamil(1, i)**2) * ftlm_h_eigv(i) * exp(-beta * ftlm_h_eigv(i))
                beta = beta + delta_beta_ftlm
            end do
        end do

    end subroutine add_in_contribs_to_energy