Hello,
I am solving a non-linear equation with DG.
As per step 15, we set the boundary values initially using the following
function:
template <int dim>
void MinimalSurfaceProblem<dim>::set_boundary_values()
{
std::map<types::global_dof_index, double> boundary_values;
VectorTools::interpolate_boundary_values
<https://www.dealii.org/current/doxygen/deal.II/namespaceVectorTools.html#ab2562d41bb26f362043f9719a8cd9b87>
(dof_handler, 0, BoundaryValues<dim(), boundary_values);
for (auto &boundary_value : boundary_values)
current_solution(boundary_value.first) = boundary_value.second;
hanging_node_constraints.distribute(current_solution);
}
Something like this doesn't work with DG.
Loosely speaking, my initial approach would be to solve a system with just
the BC terms but I am not sure. Do I have other options?
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/9e893d7a-2629-4378-aa8e-ce457f0c0e51n%40googlegroups.com.
