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.
>
>
I have used MeshWorker to compute the right hand side like this
// Attach rhs vector to assembler MeshWorker::Assembler::ResidualSimple<
Vector<double> >
assembler; NamedData< Vector<double>* > rhs;
Vector<double>* data = &right_hand_side;
rhs.add (data, "RHS");
assembler.initialize (rhs);
All steps above are necessary, though I too dont understand their purpose.
praveen
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
