[Trilinos-Users] Epetra_FECrsGraph question

kent-and at simula.no kent-and at simula.no
Wed Aug 13 10:28:58 MDT 2008


Dear Alan,

The :InsertGlobalIndices call is performed one time element so it should
be the same number of calls in all examples (except that the 2D and 3D
mesh are different).
I have attached the Dolfin code for creating the sparsity pattern.


Thanks,
Kent



>
>
>> -----Original Message-----
>> From: trilinos-users-bounces at software.sandia.gov
>> [mailto:trilinos-users-bounces at software.sandia.gov] On Behalf
>> Of Kent-Andre Mardal
>> Hello,
>>
>> But what I don't understand is that there is such a great
>> difference in
>> the case where the finite element method contains vector type elements
>> in 3D. While in the scalar case, Epetra is comparable to the other
>> backends.
>
> I see what you mean, the timings are substantially different for those
> last two cases.
> The pseudo-code you sent indicates you're calling
> CrsGraph::InsertGlobalIndices, which takes a single row. If you have
> vector fields it just means that you call InsertGlobalIndices more times.
> Is it possible that you could send more details? Perhaps more of your
> code, or some profiling data.
>
> Alan
>
>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: EpetraSparsityPattern.cpp
Url: http://software.sandia.gov/mailman/private/trilinos-users/attachments/20080813/fdd5738f/attachment-0001.pl 


More information about the Trilinos-Users mailing list