[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