Hi all,
I am trying to implement robin boundary condition as the folowing:
du/dy=u;
I implement the exact solution as
template <int dim>
double Solution<dim>::value (const Point<dim> &p,
const unsigned int) const
{
double return_value = 0.5*(p(0)*p(0)+p(1)*p(1));
return return_value;
}
I implement robin values in the suitable edges like the folowing
for (unsigned int q_point=0; q_point<n_q_points; ++q_point)
const double Robin_value
=
(exact_solution.value(fe_face_values.quadrature_point(q_point)));
//fe_face_values.normal_vector(q_point));
for (unsigned int i=0; i<dofs_per_cell; ++i)
{
for (unsigned int j=0; j<dofs_per_cell; ++j)
cell_matrix(i,j) += (Robin_value*
fe_values.shape_value(i,q_point) *
fe_values.shape_value(j,q_point) *
fe_values.JxW(q_point));
}
but it gave me the folowing errror:
"Robin_value’ was not declared in this scope"
could you please help me to correct this error?
Best Regards,
Ahmad
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
