[dune-pdelab] Time-dependent right hand side
conf86 at web.de
conf86 at web.de
Fri Aug 28 12:57:56 CEST 2015
> I have the heat equation:
> \partial_t u - \Delta u = f(t) = (e.g.) t
>
> So the Integral in *LocalOperator is (example02_operator.hh):
>
> RF factor = it->weight()*eg.geometry().integrationElement(it->position());
> for (size_type i=0; i<lfsu.size(); ++i)
> r.accumulate(lfsu, i, (gradu*gradphi[i] - f*phi[i]) * factor);
>
> My questions are: Is *LocalOperator actually the right place to define the time-dependent f? (If not, where does it go?)
> And how can I get the time variable for RF f = time?
Hi,
I hope I'm not getting annoying, but I cant't figure out how to implement this.
Is there some obvious example where this is done which I missed or is it a difficult/unclear question?
Thank you in advance!
Best regards
More information about the dune-pdelab
mailing list