## 4 Discussion
### Concluding Remarks  
The source-growth algorithm-based inversion scheme herein presented has shown to be an efficient method to invert the geometry of geobodies in highly complex zones such as the deepwater Gulf of Mexico. Part of the success is the incorporation of constrictions to mass density or magnetic susceptibility (pre-fixed) contrasts, which optimizes the inversion of large-scale bodies. On the other hand, the incorporation of FTG data provided an excellent approximation to the geometry of salt structures and their surroundings within the study region. While the velocities models of the zone under study obtained by seismics presented severe limitations in their interpretation due to the existence of a highly complex geology with intense allochthonous salt and clay tectonics--in particular because the lateral and vertical anisotropic changes are very abrupt--the geometric elements provided by potential data inversion allowed to reconstruct the interest horizons with a good approximation. The selection of an adequate \(\lambda\) parameter and contrasts for the inversion of data coming from the geological acquaintance of the region under study are fundamental to estimate accurate results. In these cases, this information aids in establishing adequate and more concise search spaces which highly reduces the non-unicity problems or the appearance of fictitious sources. The significant reduction of the space search allows the generation of non-prohibitive solutions from the computational burden point of view. This, because the initial model was built with appropriate distributions of mass densities/magnetic susceptibilities according to the current geology of the zone under study, which provided an initial model which is close to the global optimum of the objective function. In this sense, the pre-fixed mass density and magnetic susceptibility for the inversion algorithm were of 2200 kg/m\({}^{3}\) and \(-\)0.00036 [1], respectively. In all cases, the total number of parameters for the inversion algorithm is 72,269 prisms; the observations mesh was located 100 m above the origin, and the misfit between computed and observed values was evaluated through the \(L_{2}\) norm. In the case of FTG data, the joined inversion of the \(G_{uv}\), \(G_{xy}\) and \(G_{zz}\) components of the gradiometric tensor allowed to diminish the memory allocation for the Jacobian used in the simultaneous inversion, which in turn also speeds up the convergence of the process. As gradiometric data represent a better sensitivity to characterize the limits of bodies than gravimetric or magnetic data, their inversion provides results of great aid for the interpreter of reservoir engineering, due to the importance of possessing clear information to define geometries of bodies which volutmetries are later to be estimated if they are of economic interest. Although simultaneous inversion of gravity, magnetic and FTG were not performed due to the capacity of the available workstation, it is important to point out that the proposed method for the simultaneous inversion of FTG components can be easily extended to provide a simultaneous inversion of all potential field data in order to yield to a univocal and better determined imaging of 3D geological structures. With respect to the parallelization strategy, as OpenMP could be considered a standard for MT programming on SMP systems, sustained by a combination of functions and compiler directives, it reduces the lines of code and complexity with respect to low-level coding (like with the usage of pthreads), easing the migration of serial codes to a parallel shared memory environment. Additionally, the performance is really acceptable; in this regard, for the performance results in this work, we can observe an almost perfect speed-up until 12 threads; beyond this limit, the performance no longer scales on logical cores, this means that the parallel OpenMP-based implementation was well developed.