Proceedings of: Second International Workshop on Sustainable Ultrascale Computing Systems (NESUS 2015). Krakow (Poland), September 10-11, 2015. In this paper results of recent updates of a simple algorithm for the inversion of gravity anomalies for 3D geosections in parallel computer systems are presented. A relaxation iterative principle was used updating step by step the geosection distribution of mass density. Selection of updates was done on basis of least squares error match of the update effect with the observed anomaly. Locally weighted least squares combined with the linear trend were used to obtain good inversion results for two-body geosections.