template <int dim>
inline
-void VectorBoundaryValues<dim>::vector_value (const Point<dim> &p,
+void VectorBoundaryValues<dim>::vector_value (const Point<dim> &,
Vector<double> &values) const
{
Assert (values.size() == dim+1,
ExcDimensionMismatch (values.size(), dim+1));
- for (unsigned int i=0; i<dim; ++i)
- values(i) = 0;
- values(dim) = p(0)*p(0);
+ values = 13;
}
{
if (fixed_dofs[i] == true)
{
- Assert (dirichlet_dofs[i] == 0, ExcInternalError());
+ Assert (dirichlet_dofs[i] == 13, ExcInternalError());
}
else
{
--- /dev/null
+
+DEAL::Number of active cells: 4
+DEAL::Total number of cells: 5
+DEAL::Number of degrees of freedom: 21
+DEAL::0 13.0000
+DEAL::1 13.0000
+DEAL::2 13.0000
+DEAL::8 13.0000
+DEAL::9 13.0000
+DEAL::13 13.0000
+DEAL::14 13.0000
+DEAL::18 13.0000
+DEAL::Number of active cells: 8
+DEAL::Total number of cells: 9
+DEAL::Number of degrees of freedom: 63
+DEAL::0 13.0000
+DEAL::1 13.0000
+DEAL::2 13.0000
+DEAL::3 13.0000
+DEAL::4 13.0000
+DEAL::5 13.0000
+DEAL::6 13.0000
+DEAL::14 13.0000
+DEAL::15 13.0000
+DEAL::16 13.0000
+DEAL::17 13.0000
+DEAL::23 13.0000
+DEAL::24 13.0000
+DEAL::25 13.0000
+DEAL::26 13.0000
+DEAL::32 13.0000
+DEAL::33 13.0000
+DEAL::38 13.0000
+DEAL::39 13.0000
+DEAL::40 13.0000
+DEAL::41 13.0000
+DEAL::47 13.0000
+DEAL::48 13.0000
+DEAL::53 13.0000
+DEAL::54 13.0000
+DEAL::59 13.0000
In the beginning the Universe was created. This has made a lot of
people very angry and has been widely regarded as a bad move.
Douglas Adams