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

Hello 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 is hap

Re: [deal.II] Efficient implementation of varying Dirichlet BCs when everything else in the system is constant

On 08/19/2016 05:04 PM, David F wrote: Thank you, I'll take a look in detail at those examples. However, those examples don't use PETSc and the problem is that I can't find any method in PETScWrappers::SparseMatrix to copy a matrix. Which way would you suggest? Oh, I see, that is a problem. I t

Re: [deal.II] Efficient implementation of varying Dirichlet BCs when everything else in the system is constant

Thank you, I'll take a look in detail at those examples. However, those examples don't use PETSc and the problem is that I can't find any method in PETScWrappers::SparseMatrix to copy a matrix. Which way would you suggest? David. On Saturday, 20 August 2016 00:41:05 UTC+2, Wolfgang Bangerth w

Re: [deal.II] Efficient implementation of varying Dirichlet BCs when everything else in the system is constant

On 08/19/2016 04:23 PM, David F wrote: P.S.: I tried to simply repeat the process to re-apply the BCs with a new value using the same matrix, but I get always the first BCs I apply (I think because once I apply it, the matrix knows it has been already condensed and ignores the next calls). Ye

[deal.II] Efficient implementation of varying Dirichlet BCs when everything else in the system is constant

Hello everyone, I looked in the forum, in the documentation and in the tutorials but I couldn't find a way to solve this problem in the way I need. I would appreciate a lot any help you could provide me. I would like to have Dirichlet BCs applied to faces with some boundary id. This id is not

[deal.II] Efficient implementation of varying Dirichlet BCs when everything else in the system is constant

Hello everyone, I looked in the forum, in the documentation and in the tutorials but I couldn't find a way to solve this problem in the way I need. I would appreciate a lot any help you could provide me. I would like to have Dirichlet BCs applied to faces with some boundary id. This id is not

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

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. Below I have set up what I think to be the

Re: [deal.II] Re: Getting number of hanging support points

Deepak, I am aware of the fact that ConstraintMatrix.n_constraints() gives the > number of hanging nodes for h-refinement, but if I only use p-refinement, > can it give the number of hanging support points occurring due to the > difference in the order of bases between two adjacent elements? I