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

OK, I see. Anyway, as always, Thanks for educating me free of charge :-) Kyusik. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Googl

[deal.II] Re: More questions about FE_enriched & PUM FEM

Hi Edith, On Monday, November 28, 2016 at 7:19:56 PM UTC+1, Edith Sotelo wrote: > > Hi all, > I have more questions regarding Fe_enriched & PUM > > 1) Are enriched functions shifted, so that the enriched part evaluates to > zero at the mesh nodes? > no, nothing is shifted. Enrichment function

[deal.II] Re: Ordering of entries in FEValuesExtractors::SymmetricTensor

Dear Arjun, The ordering of the extracted data should align with how the components are internally unrolled. Both the Tensor and SymmetricTensor classes have two functions to assist with this type of operation, namely unrolled_to_component_indices

[deal.II] Ordering of entries in FEValuesExtractors::SymmetricTensor

Hello everyone, I am trying to follow Wolfgang's suggestion here (https://groups.google.com/forum/#!topic/dealii/KrSp6k0FSx8) to implement a mixed formulation for elasticity. I'm experimenting with the Hu-Washizu 3-field variational form. The most straightforward approach appears to be to decl

Re: [deal.II] Re: Problem with get_back_interpolation_matrix when dim=2, spacedim=3

On 10/19/2016 07:59 AM, Praveen C wrote: Well, thinking about it, it seems that spacedim=3 is not relevant here. I can do constFE_DGQArbitraryNodes fe1(QGauss<1>(fe.degree+1)); constFE_DGQArbitraryNodes fe2(QGauss<1>(fe.degree)); FETools::get_back_interpolation_matrix(fe1, fe

[deal.II] problem in periodic boundary condition

Dear all, I am trying to implement the periodic boundary condition for a nonlinear parabolic equation (scalar). I have considered a *parallelepiped *as shown in the attached pdf. I want to impose periodic BC at x=0 and x=w surfaces. I have used an offset vector as appearing in the piece of code

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

Kyusik, As you said above, I only touched the points on the boundary as follows... (I attached the file for detail) GridTools::copy_boundary_to_manifold_id(triangulation); const SphericalManifold<2> boundary_description(Point<2>(0,0)); triangulation.set_manifold (1, boundary_description);

[deal.II] More questions about FE_enriched & PUM FEM

Hi all, I have more questions regarding Fe_enriched & PUM 1) Are enriched functions shifted, so that the enriched part evaluates to zero at the mesh nodes? 2) About strong boundary conditions: I am using the function interpolate_boundary_values,I think this function is also interpolating va

Re: [deal.II] Issue on Windows x64 Debug linkage

I used both methods that i said earlier and both worked. But to avoid having to set a variable inside the visual studio command prompt and opening the project from there, at the end, i prefer to use the CMake solution in https://github.com/preillyme/hiphop-php/blob/master/CMake/VisualStudioTool

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

Daniel, thank you very much for the quick reply. Now I know how to approach the hanging nodes correctly at once in one single system! Best regards, Heinrich -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii