With the commit I just made.... Trilinos support for solving linear problems
is basically there.... there is just one piece missing....
The piece that is missing has to do with the way Epetra_FECrsMatrix works
with SumIntoGlobalValues. It turns out that if there hasn't been a value
_INSERTED_ into that position before you try to sum into it.... it will
ignore the summed value! I ran into this about a year ago when I was doing
my own application.
What I did back then was just do an InsertGlobalValues() of all zeroes for
everywhere that I needed to sum in. In reality, I think what needs to be
done is an Epetra_CrsGraph() needs to get filled out with the exact non-zero
sparsity pattern.
What do you think Ben? Shouldn't be too hard right? After a bit of looking
around it appears that the dof_map can generate the full sparsity pattern
which shouldn't be too hard to translate into the CrsGraph that Epetra
wants.
I'm going to work on this a bit for the next hour or so, but if you can get
it done quicker than I am... feel free!
Derek
On Fri, Aug 22, 2008 at 7:35 AM, Derek Gaston <[EMAIL PROTECTED]> wrote:
> Sweet! I'll take a look.
>
> I had to do presentations all day yesterday... but I'm back on working on
> this today.
>
> Derek
>
>
> On Aug 21, 2008, at 10:58 AM, Benjamin Kirk wrote:
>
> Derek,
>>
>> I just checked in a change which should construct the map properly for a
>> sparse matrix in trilinos. Now, I remember thinking this was gonna be
>> harder than what I just did, so maybe I missed something?
>>
>> In any case, check it out...
>>
>> -Ben
>>
>>
>> -------------------------------------------------------------------------
>> This SF.Net email is sponsored by the Moblin Your Move Developer's
>> challenge
>> Build the coolest Linux based applications with Moblin SDK & win great
>> prizes
>> Grand prize is a trip for two to an Open Source event anywhere in the
>> world
>> http://moblin-contest.org/redirect.php?banner_id=100&url=/
>> _______________________________________________
>> Libmesh-devel mailing list
>> [email protected]
>> https://lists.sourceforge.net/lists/listinfo/libmesh-devel
>>
>
>
-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Libmesh-devel mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-devel
