[Trilinos-Users] experiences with Trilinos and OpenMP

Jonas.Thies at dlr.de Jonas.Thies at dlr.de
Wed Mar 5 07:30:14 MST 2014


Hi,

we're getting a bit frustrated trying to combine Epetra and Teuchos with
OpenMP. A few things just seem to be wrong, for instance CrsMatrix::SumIntoGlobalValues
has a line

#ifdef EPETRA_HAVE_OMP_NONASSOCIATIVE
#pragma omp atomic
#endif

I think this should always be done atomic because FECrsMatrix uses this function
and several threads thus might sum into the same row simultaneously. At least some
synchronization is needed here.

Also it seems that Teuchos::RCP does not work correctly if used in an OpenMP section,
but I can't give a simple example of that right now. Is this a known problem?

Does anyone have similar experience with OpenMP and Epetra?

thanks,
Jonas

--------------------------
Deutsches Zentrum für Luft- und Raumfahrt e.V. (DLR)
Simulations- und Softwaretechnik | Linder Höhe | 51147 Köln

Jonas Thies
Telefon +49-2203-601-4145 | Telefax +49-2203-601-3070 | jonas.thies at dlr.de<mailto:jonas.thies at dlr.de>
www.DLR.de/sc<http://www.dlr.de/sc>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://software.sandia.gov/pipermail/trilinos-users/attachments/20140305/57c00221/attachment.html 


More information about the Trilinos-Users mailing list