The structural inverse gravity problem in a multilayer medium is one of the most important geophysics problem. Until recently, the problem was reduced to the separation of gravitational fields and the restoration ofunknown layers independently. Now the methods are in demand that allow find unknown layers simultaneously. For solving Urysohn integral equation of the first kind describing the problem regularized algorithmsLevenberg-Marquardt type with weight factors are investigated. A new Levenberg-Marquardt type methodbased on Levenberg-Marquardt scheme is proposed. A regularized Levenberg-Marquardt type method comparedwith classic Levenberg-Marquardt method. For classic Levenberg-Marquardt method some computationaloptimizations are offered. The numerical experiments using model gravitational data allow to compareconvergence rates, relative errors and program execution times of classic Levenberg-Marquardt algorithm andLevenberg-Marquardt method. The parallel programs implementing the algorithms are developed using CUDAand OpenMP technologies.