The process of converting between ECI and LLA is generally known to be complex, as it involves a intermediate conversion to ECEF (Earth-Centered, Earth-Fixed). This, in turn, involves the consideration of precession, nutation and polar motion to transform the ECI coordinate frame to an ECEF one.
However, I have seen in a few commonly used software libraries the use of a CelesTrak article's algorithm for a direct ECI to LLA operation, whereby they calculate the Longitude, Latitude and Altitude of a satellite using values from an ECI frame (without going through ECEF) by instead calculating the LLA of the sub-satellite point - that being, the point on the Earth's surface directly below the satellite. It also involves the usage of the Julian Date in order to determine the rotation of the Earth beneath the satellite.
How accurate is the use of this method in calculating the LLA of a satellite's position from its ECI values? This includes both how it compares to the ECEF method and the legitimacy of this usage in general (it looks like it came from a book made in 1996, before the developments of the IAU SOFA library.) Moreover, would it be valid to use the CelesTrak method if its only use was to plot a satellite's position above the map of an Earth?