Re: [deal.II] Dealing with conflicting constraints without using VectorTools::interpolate_boundary_values()

2017-03-20 Thread Sogo Shiozawa
Thank you for your quick replay. It seems like the way I assigned the condition in if ("certain condition") was not appropriate, and neither do the boundary conditions. That was the reason the matrix was invertible. After fixing that, it looks working well now. Best. Sogo On Sunday, March 19,

Re: [deal.II] Dealing with conflicting constraints without using VectorTools::interpolate_boundary_values()

2017-03-19 Thread Wolfgang Bangerth
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

[deal.II] Dealing with conflicting constraints without using VectorTools::interpolate_boundary_values()

2017-03-18 Thread Sogo Shiozawa
Hello, I am relatively new to deal.ii. I wonder if I can deal with conflicting constraints without using VectorTools::interpolate_boundary_values() . The explanation here (https://www.dea

Re: [deal.II] Dealing with conflicting constraints

2016-07-05 Thread Ce Qin
Dear Wolfgang, I think the first one does the correct thing whereas the second may not. I think so. But in my case, when following the instructions in the first one, the linear solver does not converge at all. The second one does. Since @Markus implemented the preject_boundary_values_curl_confor

Re: [deal.II] Dealing with conflicting constraints

2016-07-05 Thread Wolfgang Bangerth
On 07/04/2016 08:52 AM, Ce Qin wrote: Dear all, I am now solving a problem using adaptive finite element method. The boundary condition is Dirichlet condition. When the cells on boundary are locally refined, the hanging node constraints are conflicting with inhomogeneous constraints. In Const

[deal.II] Dealing with conflicting constraints

2016-07-04 Thread Ce Qin
Dear all, I am now solving a problem using adaptive finite element method. The boundary condition is Dirichlet condition. When the cells on boundary are locally refined, the hanging node constraints are conflicting with inhomogeneous constraints. In Constraints on degrees of freedom