@misc {Simula.simula.2004,
	title = {Developing and Implementing Parallel Algorithms for Fold Simulations in 3D},
	howpublished = {Presented during a workshop at Simula: Using the GPU for Scientific Computing: coding concepts and a look at applications in structural geology http://simula.no/intranet/news/using-the-gpu-for-scie\ ntific-computing-coding-concepts-and-a-look-at-appl},
	year = {2013},
	note = {Presented during a workshop at Simula: Using the GPU for Scientific Computing: coding concepts and a look at applications in structural geology http://simula.no/intranet/news/using-the-gpu-for-scientific-computing-coding-concepts-and-a-look-at-applications-in-structural-geology},
	month = {February},
	abstract = {The Compound Earth Simulator is an exploration software developed by Kalkulo for Statoil. In the Compound methodology, large folded geological structures are simulated repeatedly as geological models are investigated. Folded layers are modeled as positions of a front moving through a geological, and are computed using first-arrival-time equations. The numerical simulation times must be low for the software to function at interactive times, which is a challenge in 3D as the number of unknown nodes increases rapidly. This talk presents a new set of algorithms for such simulations, specifically developed for massively parallel compute architectures such as GPUs. A couple of optimization approaches and their effects are presented for the GPU implementation. Finally, some early result from a multicore CPU and a GPU implementation will be given.},
	keywords = {Workshop},
	author = {Gillberg, Tor}
}
@article {Simula.simula.622,
	title = {Accuracy and Efficiency of Stencils for the Eikonal Equation in Earth Modelling},
	journal = {Computational Geosciences},
	volume = {16},
	number = {4},
	year = {2012},
	month = {September},
	pages = {933-952},
	abstract = {Motivated by the needs for creating fast and accurate models of complex geological scenarios, accuracy and efficiency of three stencils for the isotropic eikonal equation on rectangular grids are evaluated using a Fast Marching implementation. The stencils are derived by direct modelling of the wave front, resulting in new and valuable insight in terms of improved upwind and causality conditions. After intro- ducing a method for generalising first order upwind stencils to higher order, a new second order Diagonal stencil is presented. Similarly to the Multistencils Fast Marching ap- proach, the Diagonal stencil makes use of nodes in the diagonal directions, whereas the traditional Godunov stencil uses solely edge-connected neighbours. The Diagonal stencil uses nodes close to each other, reaching upwind, to get a more accurate estimate of the angle of incidence of the arriving wave front. Although the stencils are evaluated in a Fast Marching setting, they can be adapted to other efficient eikonal solvers. All first and second order stencils are evaluated in a range of tests. The first test case models a folded structure from the Zagros fold belt in Iran. The other test cases are constructed to investigate specific properties of the examined stencils. The numerical investigation considers convergence rates and CPU times for non-constant and constant speed first-arrival computations. In conclusion, the Diagonal stencil is the most efficient and accurate of the three alternatives.},
	keywords = {Workshop},
	doi = {10.1007/s10596-012-9296-0},
	author = {Gillberg, Tor and Hjelle, {\O}yvind and Are Magnus {Bruaset}}
}
@misc {Simula.simula.1126,
	title = {Massively Parallel Front Tracking},
	howpublished = {Computational Geoscience Workshop, 18-20 January 2012 in Geilo},
	year = {2012},
	month = {January 2012},
	abstract = {We present a novel method for three dimensional anisotropic front propagation, and apply it to the simulation of geological folding. The iterative algorithm has a simple structure, and a high level of fine grained parallelism. The source-code translator Mint is used to generate Cuda code executable on Graphic Processor Units (GPUs) from our C implementation. Our numerical section present results of such an approach, and that of a multicore implementation. The implementations have been tested and compared on several examples of three dimensional geological folding.},
	keywords = {Workshop},
	author = {Gillberg, Tor}
}
