Thankyou!
that works :)
On Saturday, 18 March 2017 12:03:55 UTC+5:30, Sumedh Yadav wrote:
>
> Hello,
> I need to do some post-processing on computed solution data (say
> Vector Solution(dof_handler.n_dofs()) is the solution vector). For
> post-processing field (say Vector Curvature(triangulatio
Dr. Bangerth,
Thank you very much for your kind reply.
> What is the content of 'constraints' in your case?
>
> DoFTools::make_hanging_node_constraints (dof_handler,
constraints);
VectorTools::interpolate_boundary_values (dof_handler,
I need to do some post-processing on computed solution data (say
Vector Solution(dof_handler.n_dofs()) is the solution vector). For
post-processing field (say Vector Curvature(triangulation.n_dofs()) is
the scalar field to be computed) I decided to have one entry per cell. Once I
calculate this
On 03/17/2017 08:03 AM, Jaekwang Kim wrote:
After the first cycle, I am doing
template
voidnonlinear::refine_grid ()
{
Vector estimated_error_per_cell (triangulation.n_active_cells());
KellyErrorEstimator::estimate (dof_handler,
QGauss(3),
Kyusik,
I'm studying Inhomogenous Dirichlet boundary conditions watching Lecture 21.65.
To practice this boundary condition, I'm trying to reproduce the result of
step-4(in which Inhomogenous Dirichlet boundary condition is used) in step-6.
As I learned in the lecture...
First, I assembled A
I assign such a boundary condition at points by giving boundary values
to std::map boundary_values directly, and it looks
working well.
Yes, this is the correct approach.
However, for handling conflict between hanging node and
boundary condition with ConstraintMatrix, just doing something l
