for i in range(len(r)):
        energy_residual[i] = norm(r[i], "fro")