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
