[Libmesh-devel] Sparsity patterns and post-facto matrices

2012-08-21 Thread Roy Stogner
We currently have a bug in the interaction between post-initialization sparse matrix creation (as used by our reduced_basis code and triggered by that ex5) and sparse matrix subclasses that need a full sparsity pattern to examine (i.e. Laspack and Trilinos). Currently, the DofMap creates a sparsi

Re: [Libmesh-devel] Sparsity patterns and post-facto matrices

2012-08-21 Thread Derek Gaston
For now I'll vote for #2. One thing that you could add to #2 is a method that we could call to clear/delete the sparsity pattern in the case where we _know_ we won't need it again and we really do need the memory... Derek On Tue, Aug 21, 2012 at 9:49 AM, Roy Stogner wrote: > > We currently have

Re: [Libmesh-devel] Sparsity patterns and post-facto matrices

2012-08-21 Thread Roy Stogner
On Tue, 21 Aug 2012, Derek Gaston wrote: For now I'll vote for #2.  One thing that you could add to #2 is a method that we could call to clear/delete the sparsity pattern in the case where we _know_ we won't need it again and we really do need the memory... Nice idea - if we go with #2 then I

Re: [Libmesh-devel] Periodic BCs patch

2012-08-21 Thread John Peterson
On Tue, Aug 21, 2012 at 9:52 AM, John Peterson wrote: > On Mon, Aug 20, 2012 at 5:11 PM, John Peterson > wrote: >> On Mon, Aug 20, 2012 at 4:46 PM, John Peterson >> wrote: >>> On Mon, Aug 20, 2012 at 4:30 PM, Roy Stogner >>> wrote: This is against the SVN head. It seems to work fine

Re: [Libmesh-devel] Periodic BCs patch

2012-08-21 Thread Roy Stogner
On Tue, 21 Aug 2012, John Peterson wrote: > So there are small differences in several of the tests which I am > investigating One case in which small differences would be expected: if you have anything where the initial conditions are not periodic-consistent, then it's possible that the way in w

Re: [Libmesh-devel] Periodic BCs patch

2012-08-21 Thread Roy Stogner
On Tue, 21 Aug 2012, John Peterson wrote: > Note that the periodicity of v is set up with two sets of periodic BCs > (one for the left and one for the right half of the domain) so it > could have to do with a pair of periodic boundaries meeting at a > point? That's a tricky case, but one I speci

Re: [Libmesh-devel] Sparsity patterns and post-facto matrices

2012-08-21 Thread David Knezevic
On 08/21/2012 12:02 PM, Roy Stogner wrote: > > On Tue, 21 Aug 2012, Derek Gaston wrote: > >> For now I'll vote for #2. One thing that you could add to #2 is a >> method >> that we could call to clear/delete the sparsity pattern in the case >> where we >> _know_ we won't need it again and we rea

Re: [Libmesh-devel] Sparsity patterns and post-facto matrices

2012-08-21 Thread Roy Stogner
On Tue, 21 Aug 2012, David Knezevic wrote: > Another question: does this issue explain why the reduced basis examples > never worked with trilinos? Well, it's certainly one contributing factor. But it's never safe to assume that under any given bug is pure working code; see the periodic BCs dis

Re: [Libmesh-devel] Sparsity patterns and post-facto matrices

2012-08-21 Thread Kirk, Benjamin (JSC-EG311)
I'm glad Derek suggested this - it was my preferred approach too! -Ben On Aug 21, 2012, at 5:55 AM, "Derek Gaston" mailto:fried...@gmail.com>> wrote: For now I'll vote for #2. One thing that you could add to #2 is a method that we could call to clear/delete the sparsity pattern in the case