Re: [deal.II] Does MappingFEField apply to all cells ?

2016-08-22 Thread Wolfgang Bangerth
Praveen, I am using MappingFEField to get high order meshes which I construct by solving Winslow equations. I set initial condition and save it to file with 3 patches. As the figure shows, it seems to me that the MappingFEField is used only for the cells adjacent to the boundary. Is this what

Re: [deal.II] Does MappingFEField apply to all cells ?

2016-08-22 Thread Praveen C
Thanks Wolfgang. I should have checked the built_patches function. It works fine if I give the third argument. But now I have to find the real reason why my code is crashing with MappingFEField :-( Best praveen On Mon, Aug 22, 2016 at 6:36 PM, Wolfgang Bangerth wrote: > > Praveen, > > I am usin

Re: [deal.II] Does MappingFEField apply to all cells ?

2016-08-22 Thread Wolfgang Bangerth
On 08/22/2016 07:30 AM, Praveen C wrote: Thanks Wolfgang. I should have checked the built_patches function. It works fine if I give the third argument. But now I have to find the real reason why my code is crashing with MappingFEField :-( There is always one more bug. Always. Best W. --

[deal.II] Are you funded to write applications based on deal.II?

2016-08-22 Thread Wolfgang Bangerth
All, Timo and I are writing a proposal that will (hopefully) fund further development of deal.II in the coming years. It would be helpful if we had an overview of what other projects there are where people are funded to write applications based on deal.II (whether or not they will be made avai

Re: [deal.II] weakly enforcing an identity - vector mean curvature as the surface laplacian of the identity on a codim-1 manifold

2016-08-22 Thread Wolfgang Bangerth
On 08/19/2016 01:15 PM, thomas stephens wrote: Wolfgang, I appreciate your assistance so far. I'm getting closer, but I still need some help. I would like to look at this problem as it appears in [1], Equation 27 on p. 11, as the bilinear forms are written as integrals, making things more clear

[deal.II] Re: Vector-valued gradient of solution vector

2016-08-22 Thread Joel Davidsson
Dear Daniel, Thank you for a very good answer, adding the fe_values_scalar and cell_scalar fixed the problem. I have a follow-up question about the solution I get out, how is the data organized in the solution vector? Say for example I want to loop over all the x components, how would I do tha

[deal.II] Re: Vector-valued gradient of solution vector

2016-08-22 Thread Daniel Arndt
Joel, You should not rely on any particular sorting of the dofs in the solution vector. Instead you can ask the FiniteElement on each cell for the component of a local dof by const unsigned int component = fe.system_to_component_index(i).first; Apart from that the DataOut object knows about co

[deal.II] Some bug with deal.II from git

2016-08-22 Thread Praveen C
Dear all I have convergence problem with one of my codes using deal.II from git. Using the Aug 4 version (commit ed1bdc6fc1931a664bc30920f66f054907a49cce) the code converges to 1.0e-10 in 20 iterations. Using the Aug 16 version (commit 2ca845c83a1bceb8d1d09df2ee729b4ec6d5547b) does not converge,

Re: [deal.II] Some bug with deal.II from git

2016-08-22 Thread Matthias Maier
> Using the Aug 4 version (commit ed1bdc6fc1931a664bc30920f66f054907a49cce) > the code converges to 1.0e-10 in 20 iterations. > > Using the Aug 16 version (commit 2ca845c83a1bceb8d1d09df2ee729b4ec6d5547b) > does not converge, residual is 1.8 after 20 iterations, which is a huge > difference. You c