I think the Functional works like this if you have one functional (number is double or int depending on your functional)
MeshWorker::Assembler::Functional<double> assembler;assembler.initialize(1); ... call integration_loop then assembler (0) should give you the result. praveen On Mon, Feb 7, 2011 at 7:26 AM, mihai alexe <[email protected]> wrote: > Hi all, > > I am trying to evaluate an integral-type functional, say \int_\Omega > -(\Delta q) * delta_q dx (+ some other boundary/volume terms) where \Delta > is the Laplacian. I'm looking into the MeshWorker framework for this... > which class would be the best choice to put together this integral given > that both q and delta_q are known (I have them as Vector<double> and > Function<dim>, respectively). I have tried doing things with > MeshWorker::Assembler::ResidualSimple, but I cannot initialize it with a > dummy Vector<double>. Instead, it wants a NamedData object whose purpose I > don't quite understand. > > Upon inspection of the MeshWorker framework I found the Functional class > that looks (at least by name) to be more suitable to what I'm trying to do. > I'm not quite sure how to initialize that. It asks for an integer n. What is > that used for? > > I could use assembler::MatrixSimple but I figure that may be overkill in my > case. Also, I am trying to avoid writing the assembly loop by hand . > > Thanks! > -- Mihai > > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii > >
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
