diff --git a/README.md b/README.md index 8c3a850..ed40b38 100644 --- a/README.md +++ b/README.md @@ -49,7 +49,7 @@ and runtime measurements of the presented code. The implementation is based on the paper [Tsoulis, D., 2012. Analytical computation of the full gravity tensor of a homogeneous arbitrarily shaped polyhedral source using line integrals. Geophysics, 77(2), pp.F1-F11.](http://dx.doi.org/10.1190/geo2010-0334.1) -and its corresponding [implementation in FORTRAN](https://software.seg.org/2012/0001/index.html). +and its corresponding implementation in FORTRAN. Supplementary details can be found in the more recent paper [TSOULIS, Dimitrios; GAVRIILIDOU, Georgia. A computational review of the line integral analytical formulation of the polyhedral gravity signal. Geophysical Prospecting, 2021, 69. Jg., Nr. 8-9, S. 1745-1760.](https://doi.org/10.1111/1365-2478.13134) diff --git a/docs/background/approach.rst b/docs/background/approach.rst index 9a3005d..f405195 100644 --- a/docs/background/approach.rst +++ b/docs/background/approach.rst @@ -12,11 +12,11 @@ and runtime measurements of the presented code. The implementation is based on the paper `Tsoulis, D., 2012. Analytical computation of the full gravity tensor of a homogeneous arbitrarily shaped polyhedral source using line integrals. Geophysics, 77(2), pp.F1-F11. `__ -and its corresponding implementation in FORTRAN_. +and its corresponding implementation in FORTRAN_ (last accessed: 12.09.2022). Supplementary details can be found in the more recent paper `TSOULIS, Dimitrios; GAVRIILIDOU, Georgia. A computational review of the line integral analytical formulation of the polyhedral gravity signal. Geophysical Prospecting, 2021, 69. Jg., Nr. 8-9, S. 1745-1760 `__ -and its corresponding implementation in MATLAB_, +and its corresponding implementation in MATLAB_ (last accessed: 28.03.2024), which is strongly based on the former implementation in FORTRAN. diff --git a/paper/paper.md b/paper/paper.md index 00c6960..235c717 100644 --- a/paper/paper.md +++ b/paper/paper.md @@ -50,7 +50,7 @@ The complex gravitational fields of irregular bodies, such as asteroids and come In contrast, polyhedral gravity models provide an analytic solution for the computation of the gravitational potential, acceleration (and second derivative) given a mesh of the body [@tsoulis2012analytical;@tsoulis2021computational] with the only assumption of homogeneous density. The computation of the gravitational potential and acceleration is a computationally expensive task, especially for large meshes, which can however benefit from parallelization either over computed target points for which we seek potential and acceleration or over the mesh. Thus, a high-performance implementation of a polyhedral gravity model is desirable. -While some research code for these models exists, they are not focused on usability and are limited to FORTRAN\footnote{\url{https://software.seg.org/2012/0001/index.html}} and proprietary software like MATLAB\footnote{\url{https://github.com/Gavriilidou/GPolyhedron}}. There is a lack of well-documented, actively maintained open-source implementations, particularly in modern programming languages, and with a focus on scalability and performance. +While some research code for these models exists, they are not focused on usability and are limited to FORTRAN\footnote{\url{https://software.seg.org/2012/0001/index.html}, last accessed: 12.09.2022} and proprietary software like MATLAB\footnote{\url{https://github.com/Gavriilidou/GPolyhedron}, last accessed: 28.03.2024}. There is a lack of well-documented, actively maintained open-source implementations, particularly in modern programming languages, and with a focus on scalability and performance. This circumstance and the fact that polyhedral models are often used in studying gravitational fields, e.g., for Eros [@zhang2010modeling], or as a reference for creating new neural models [@martin2023physics] make an easy-to-install implementation necessary. diff --git a/test/resources/mod.py b/test/resources/mod.py deleted file mode 100644 index fc268f3..0000000 --- a/test/resources/mod.py +++ /dev/null @@ -1,22 +0,0 @@ -FILE_1 = "analytic_cube_solution_density1_old.txt" -FILE_1_OUT = "analytic_cube_solution_density1.txt" -FILE_42 = "analytic_cube_solution_density42_old.txt" -FILE_42_OUT = "analytic_cube_solution_density42.txt" - -def mod(input_name, output_name): - with open(input_name, 'r') as read_file, open(output_name, 'w') as write_file: - lines = read_file.readlines()[1:] # Skip the first line - - for line in lines: - elements = line.split() # Splits the line into multiple elements - - # Changes the sign of the last three elements - elements[-1] = str(float(elements[-1]) * -1) - elements[-2] = str(float(elements[-2]) * -1) - elements[-3] = str(float(elements[-3]) * -1) - - write_file.write(' '.join(elements) + '\n') - - -if __name__ == '__main__': - mod(FILE_42, FILE_42_OUT)