[deal.II] Re: How to change material ids during iteration.

Thanks for your reply. you are right, I shouldn't have been in a hurry to show some results to my boss... I should have asked, after looking through several examples especially step49 that includes GridTools::laplace_transform. I guess, using GridTools::laplace_transform, I can solve my proble

Re: [deal.II] Interpolating p::d::Vector

Hello Martin I am patching many charts to define a manifold to cover the sphere. I use MappingManifold. I suspect at the interface of two charts, which also coincides with the mesh partition boundary, there may be some round off error. Best praveen On Wed, Nov 23, 2016 at 7:12 PM, Martin Kronbic

Re: [deal.II] Interpolating p::d::Vector

Dear Praveen, Called compress(VectorOperation::insert), but the element received from a remote processor, value 3.343515457745695e-15, does not match with the value -3.343515457745695e-15 on the owner processor 0 I dont understand why this error is triggered. The error says that two processo

[deal.II] Interpolating p::d::Vector

Dear all I have an FE_Q space and a ghosted p::d::Vector. When I call VectorTools::interpolate on this, I get this error An error occurred in line <679> of file in function void dealii::LinearAlgebra::distributed::Vector::compress_finish(::dealii::VectorOperation::values) The violated cond

[deal.II] Calculate Stresses in Step-14

Hello, I would like to modify step-14 and calculate stress at nodes, but I found step-14 complicated and also unfortunately my c++ is not good :( I would like to add stress_strain_nodal_values method to classes Base, PrimalSolver and WeightedResidual ( similar to output_solution() ) then insid

Re: [deal.II] Relation of FE_Q support points and quadrature points

Wonderful. It solved my problem. Best praveen On Wed, Nov 23, 2016 at 3:42 PM, luca.heltai wrote: > Start with a brute force approach: > > FE_Q fe(degree); > QGaussLobatto quad(degree+1); > > const std::vector > &q = fe.get_unit_support_points() > > std::vector w(q.size()); > > tol = 1e-10; > >

Re: [deal.II] Relation of FE_Q support points and quadrature points

Start with a brute force approach: FE_Q fe(degree); QGaussLobatto quad(degree+1); const std::vector > &q = fe.get_unit_support_points() std::vector w(q.size()); tol = 1e-10; // check everything is fine AssertDimension(w.size(), quad.size()); for(unsigned int i=0; i new_quad(q,w); Now you ar

Re: [deal.II] Relation of FE_Q support points and quadrature points

Thanks Luca. I needed this property to efficiently implement a spectral element method, where support points are used for quadrature also. Now I have to figure out some other way to do this efficiently. I want to avoid use of get_function_values. Best praveen On Wed, Nov 23, 2016 at 3:27 PM, lu

Re: [deal.II] Relation of FE_Q support points and quadrature points

Not necessarily. In fact the ordering is reshuffled so that vertex support points comes first, then edge support points (3D), face support points, and lastly cell support points. To see the actual ordering, take a look at fe.get_unit_support_points() https://www.dealii.org/8.4.1/doxygen/deal.

[deal.II] Relation of FE_Q support points and quadrature points

Dear all Suppose I have FE_Q fe(degree); QGaussLobatto quad(degree+1); Within a cell, I get FEValues fe_values(quad, update_quadrature_points); fe_values.reinit(cell); cell->get_dof_indices(local_dof_indices); Then can I assume that, in the current cell, solution(local_dof_indices[q]) is the