[Trilinos-Users] [EXTERNAL] Re: Find enclosing element in EXODUS II tetrahedral mesh.

Pawlowski, Roger P rppawlo at sandia.gov
Thu Jan 16 09:10:17 MST 2014

Hi Michael,

There is a Trilinos package called DataTransferKit (DTK) that does this exact computation (locating the global element index for an arbitrary point) as part of a data transfer in multiphysics code-to-code coupling.  It does not have the exact interface to return the data to you in the form you expect, so you might have to modify the interfaces (or pull the data out of the source evaluation function). 

DTK is a trilinos package but it is housed in a separate repository.  You can download from:


You can build this as part of trilinos by putting this repository in the top level of trilinos and telling trilnos about the extra repository at configure time with the configure flag:

-D Trilinos_EXTRA_REPOSITORIES="DataTransferKit"

The way to do this is to build the data transfer kit and use the STK mesh interface to load the exodus mesh.  Once in STK form, pull the element gids and coordinates out of STK and pass to DTK.  Then you can try to do a data transfer requesting an arbitrary list of points.  In the source evaluation, the gids for each point are passed into that function.


From: trilinos-users-bounces at software.sandia.gov [trilinos-users-bounces at software.sandia.gov] on behalf of Sjaardema, Gregory D [gdsjaar at sandia.gov]
Sent: Thursday, January 16, 2014 10:35 AM
To: Heroux, Mike; Michael Afanasiev; trilinos-users at software.sandia.gov
Subject: Re: [Trilinos-Users] [EXTERNAL] Re: Find enclosing element in EXODUS II tetrahedral mesh.

The exodus library does not have this functionality either.

On 1/16/14 8:26 AM, "Heroux, Mike" <MHeroux at csbsju.edu> wrote:

>I don't know if you got a private response or not, but I am unaware of
>functionality in Trilinos that would address your need.
>Best regards,
>From: Michael Afanasiev
><mafanasiev at gmail.com<mailto:mafanasiev at gmail.com>>
>Date: Monday, January 13, 2014 1:26 PM
>"trilinos-users at software.sandia.gov<mailto:trilinos-users at software.sandia.
><trilinos-users at software.sandia.gov<mailto:trilinos-users at software.sandia.
>Subject: [Trilinos-Users] Find enclosing element in EXODUS II tetrahedral
>Hi everyone,
>My name is Michael Afanasiev, and I am a PhD student at ETH Zürich, in
>I¹m working with a large 3-D tetrahedral finite element mesh, which is
>stored as an EXODUS II file. I need to find the smallest simplex that
>encloses an arbitrary 3-D point. Alternatively, the enclosing element
>would do just fine. The purpose is for interpolation over the element.
>I¹m currently accomplishing this by recursively searching a KDTree around
>my requested point, and finding the smallest enclosing simplex. This
>works pretty well in most cases. I was wondering if there was any
>functionality built into Trilinos or EXODUS that would accomplish this
>instead, or if anyone has an alternative idea.
>Cheers, and all the best. Apologies if this comes through as a duplicate.
>Michael Afanasiev
>Ph.D. Candidate
>Computational Seismology
>Institut für Geophysik
>ETH Zürich
>Sonneggstrasse 5, NO H 52.2
>CH-8092 Zürich
>michael.afanasiev at erdw.ethz.ch<mailto:michael.afanasiev at erdw.ethz.ch>
>Trilinos-Users mailing list
>Trilinos-Users at software.sandia.gov

Trilinos-Users mailing list
Trilinos-Users at software.sandia.gov

More information about the Trilinos-Users mailing list