Hello Yuliang,
the modification is still not very complicated. You just have to add
additionally the line
VectorTools::interpolate_boundary_values
<http://www.dealii.org/developer/doxygen/deal.II/classVectorTools.html#ae81c0e62e5bdfae7f3805c0109542c0c>
(dof_handler,
1,
Solution<dim>(),
boundary_values);
in assemble_system after the almost identical line for boundary indicator 0.
Best Regards,
Markus
Am 01.03.11 00:50, schrieb Yuliang Wang:
Thanks for such a quick reply.
So the example program cannot be easily modified to include the 1D
case. This is a little bit strange. I hope the authors will implement
that in later releases.
On Mon, Feb 28, 2011 at 6:17 PM, Markus Bürg <[email protected]
<mailto:[email protected]>> wrote:
Dear Yuliang,
in 1D you cannot set boundary indicators. They are fixed to 0 on
the left and to 1 on the right boundary as you observed.
Best Regards,
Markus
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
