Two new methods for solving large scale least squares in geodetic surveying computations