[dune-pdelab] Instationary heat equation
conf86 at web.de
conf86 at web.de
Sat Jul 11 21:36:00 CEST 2015
Hello,
I want to solve the instationary equation
\partial_t u + a * \laplace u = f
Where and how can I define this factor a?
When I try this in operator.hh:
RF a = 1.0;
RF factor = it->weight()*eg.geometry().integrationElement(it->position());
for (size_type i=0; i<lfsu.size(); ++i)
r.accumulate(lfsu, i, (a * gradu*gradphi[i] - f * phi[i]) * factor);
it doesn't work (sorry, I'm currently not at my office, so I don't remember the exact error message).
Thank you in advance.
More information about the dune-pdelab
mailing list