Exact diagonalization of Hubbard clusters at finite temperatures