Problem Solved! Thank you!!
2016년 11월 13일 일요일 오후 5시 25분 20초 UTC-6, Daniel Arndt 님의 말:
>
> Jaekwang,
>
> It seems that you are missing
> constraints.distribute(solution);
> after solving the linear system. The constrained dofs are set to zero
> during assembly and
> this call is necessary to
Jaekwang,
It seems that you are missing
constraints.distribute(solution);
after solving the linear system. The constrained dofs are set to zero
during assembly and
this call is necessary to get the correct values for these in the solution
as well. You might want to have a look at
Thank you for being kindness!
I followed your suggestions but I couldn't resolve the problem yet.
Or did I make other mistake in another places?
>From a slight modification from STEP-4 tutorial , I have been trying solve
Poisson Equation.
Actually, I have manufactured solution to see my
Yes, well this is going to give strange results because you're setting the
boundary ID of internal faces (x=0.5 is the centreline of your geometry).
Before you set any boundary ID's, you should always first check to see that
you're actually on a boundary. The first few tutorials do go through
Thank you for quick reply!!!
I tried your suggestions but I couldn't resolve the problem.
The more weird thing is if I commands as..
GridGenerator::hyper_cube (triangulation, 0, 1);
triangulation.refine_global ();
for (typename Triangulation::active_cell_iterator
Hi Jaekwang,
With what you've shown I'm going to guess that your problem might simply be
related to numerical / round-off error associated with floating point
calculations.
if ( cell->face(f)->center()[0]==1 ) // The LHS value is a computed double,
> and the RHS value is an integer which
