On solving large systems of polynomial equations appearing in discrete differential geometry