Re: [Libmesh-users] Total system energy

2009-06-18 Thread Roy Stogner
On Thu, 18 Jun 2009, bkraczek wrote: On Thu, 18 Jun 2009 16:54:05 -0500, John Peterson wrote: On Thu, Jun 18, 2009 at 4:42 PM, Roy Stogner wrote: Doing matrix-based norms in LinearImplicitSystem rather than in SparseMatrix doesn't make sense to me.  Every LinearImplicitSystem user can use a

Re: [Libmesh-users] Total system energy

2009-06-18 Thread bkraczek
On Thu, 18 Jun 2009 16:54:05 -0500, John Peterson wrote: > On Thu, Jun 18, 2009 at 4:42 PM, Roy Stogner > wrote: >> >> On Thu, 18 Jun 2009, Brent Kraczek wrote: >> >>> Please let me know if you want anything changed. >> >> We can ignore my criticisms if other developers disagree, but: >> >> >>

Re: [Libmesh-users] Total system energy

2009-06-18 Thread John Peterson
On Thu, Jun 18, 2009 at 4:42 PM, Roy Stogner wrote: > > On Thu, 18 Jun 2009, Brent Kraczek wrote: > >> Please let me know if you want anything changed. > > We can ignore my criticisms if other developers disagree, but: > > > Doing matrix-based norms in LinearImplicitSystem rather than in > SparseMa

Re: [Libmesh-users] Total system energy

2009-06-18 Thread Roy Stogner
On Thu, 18 Jun 2009, Brent Kraczek wrote: Attached are modified versions of linear_implicit_system.{h,C} and ex3.C, which implement a system energy and an energy norm for a Poisson system (it made the most sense to me to add these to LinearImplicitSystem for what I am doing). The code passes

Re: [Libmesh-users] Total system energy

2009-06-18 Thread Brent Kraczek
Attached are modified versions of linear_implicit_system.{h,C} and ex3.C, which implement a system energy and an energy norm for a Poisson system (it made the most sense to me to add these to LinearImplicitSystem for what I am doing). The code passes the tests I have run. These additions were

Re: [Libmesh-users] Total system energy

2009-06-15 Thread Brent Kraczek
Sorry, it looks like I mis-read it. -Brent John Peterson wrote: > On Mon, Jun 15, 2009 at 11:00 AM, Brent Kraczek > wrote: > >> Roy, >> >> I am working toward having a patch to do this. >> >> But there is something which I don't understand here. It appears to me >> that SparseMatrix::vector_mu

Re: [Libmesh-users] Total system energy

2009-06-15 Thread John Peterson
On Mon, Jun 15, 2009 at 11:00 AM, Brent Kraczek wrote: > Roy, > > I am working toward having a patch to do this. > > But there is something which I don't understand here. It appears to me > that SparseMatrix::vector_mult doesn't call the appropriate PETSc > routine. It's not clear to me why SparseM

Re: [Libmesh-users] Total system energy

2009-06-15 Thread Brent Kraczek
Roy, I am working toward having a patch to do this. But there is something which I don't understand here. It appears to me that SparseMatrix::vector_mult doesn't call the appropriate PETSc routine. It's not clear to me why SparseMatrix::vector_mult or at least SparseMatrix::vector_mult_add has

Re: [Libmesh-users] Total system energy

2009-06-09 Thread Roy Stogner
On Tue, 9 Jun 2009, Brent Kraczek wrote: > I'm trying to calculate the total system energy for my Poisson solution, > which would be > >(solution)^T . matrix . solution > > Has this been implemented somewhere, I don't think so; you'll probably need a temporary vector. If you'd like to add t