Re: [Libmesh-devel] EpetraMatrix sparsity pattern

2008-08-22 Thread Benjamin Kirk
> I know there's still a lot of work (and cleanup) left to do... but it > is technically working. One thing I know for sure is that it's > probably leaking memory like a bitch ;-) I'll deal with that next > week though ;-) ex4 on two processors doesn't work yet! ;-) > Thanks for all the help N

Re: [Libmesh-devel] EpetraMatrix sparsity pattern

2008-08-22 Thread Derek Gaston
Ok - I just committed the sparsity pattern stuff... and now you can use Trilinos to solve Example3! I know there's still a lot of work (and cleanup) left to do... but it is technically working. One thing I know for sure is that it's probably leaking memory like a bitch ;-) I'll deal with t

Re: [Libmesh-devel] EpetraMatrix sparsity pattern

2008-08-22 Thread Derek Gaston
Ok - one thing I'm confused about is what should go into update_sparsity_pattern and what should go into init() Will init() ever be called without calling update_sparsity_pattern() first? If not then I'm going to put both the Epetra_Map and Epetra_CrsGraph creation into update_sparsity_

Re: [Libmesh-devel] EpetraMatrix sparsity pattern

2008-08-22 Thread Derek Gaston
Great, thanks for the info... I am actually pretty close to having this done... Derek On Aug 22, 2008, at 2:09 PM, Benjamin Kirk wrote: > Should the trilinos list get queried? Is that a bug?? > > >> What do you think Ben? Shouldn't be too hard right? After a bit >> of looking >> around it

Re: [Libmesh-devel] EpetraMatrix sparsity pattern

2008-08-22 Thread Benjamin Kirk
Should the trilinos list get queried? Is that a bug?? > 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. Right.

Re: [Libmesh-devel] EpetraMatrix sparsity pattern

2008-08-22 Thread Derek Gaston
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_ int

Re: [Libmesh-devel] EpetraMatrix sparsity pattern

2008-08-22 Thread Derek Gaston
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,