Three-dimensional Cauchy-condition surface method to identify the shape of the last closed magnetic surface in the Large Helical Device