Re: [deal.II] Help with triangulation::face_iterator

2016-11-27 Thread Wolfgang Bangerth
On 11/27/2016 01:48 PM, RAJAT ARORA wrote: That worked but I am not sure why a reference didn't work but a const. reference did. C++ says that you can convert an argument from one type to another if a function takes a const-reference, but not if the function takes a non-const reference. It's

Re: [deal.II] Help with triangulation::face_iterator

2016-11-27 Thread RAJAT ARORA
Thank you Professor. That worked but I am not sure why a reference didn't work but a const. reference did. On Tuesday, November 22, 2016 at 7:31:47 PM UTC-5, Wolfgang Bangerth wrote: > > > Rajat, > I don't know the exact cause of the problem, but... > > > Can we pass Triangulation::face_iter

[deal.II] Re: Problems with hanging nodes and "distribute_local_to_global"

2016-11-27 Thread Daniel Arndt
Heinrich, The way constraints are resolved in deal.II is as follows (related to *M. S. Shephard: Linear multipoint constraints applied via transformation as part of a direct stiffness assembly process. Int. J. Numer. Meth. Engrg., vol. 20 (1984), pp. 2107-2112.*): You have a linear system Ax=b

Re: [deal.II] GridTools::laplace_transform and interior points

2016-11-27 Thread hanks0227
Dr. Bangerth, Thanks for your reply again. > Yes, but that's because you have a small bug in your code: you fix *all* > vertices before you call laplace_transform. You should only fix the ones > that > are are at the boundary of the domain. > > As you said above, I only touched the points on