Re: [deal.II] Re: Using create_point_source_vector() function with hp::DoFHandler

Dear Daniel, I will try to come up with a minimal example. In the meantime, below is the piece of code and I will try to explain using this. I tried running in Debug mode as well but not able to figure out the error. For the code below in combination with the rest of the code, the last line of

[deal.II] Decoupling FECollection and QCollection

I tried finding functionalities such as active_quad_index() similar to active_fe_index(). But, I found that they are coupled when inserted using push_back(). Currently, I am not able to see a way out of it. If someone knows how to do this or has some trick in mind, please let me know. Best -- *Deepak Gupta*

[deal.II] Getting number of hanging support points

Dear All, I am currently dealing with hanging support points in my problem and I have a couple of questions. These are: 1. For my case, different elements of the mesh have different polynomial orders of the basis. Due to support points at the element edges, I believe there will be hanging

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

Dear All, Thanks for your replies. Indeed I had thought of the same as suggested by Jean and am implementing it. Thanks for helping out. Best regards Deepak On Friday, July 15, 2016 at 2:02:43 PM UTC+2, Deepak Gupta wrote: > > Dear All, > > I am currently dealing with hanging su

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

Dear All, 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 am

[deal.II] Reference article for modified Kelly Estimator

Dear All, I am using the KellyEstimator as an 'indicator' for simple hp-adaptive problems. I read the documentation that for hp cases, CF = hF/2PF is considered. I would like to know if this has been published somewhere which I can refer. In addition, I am looking for a reference where it is

Re: [deal.II] Using p-refinement with high order elements

to extend the simple example for a basic version of the elastic problem I am trying to solve. Hope then I can figure out the error. Best regards Deepak On Wed, Oct 5, 2016 at 6:03 PM, Wolfgang Bangerth <bange...@colostate.edu> wrote: > On 10/05/2016 09:03 AM, Deepak Gupta wrote: > &g

Re: [deal.II] Using p-refinement with high order elements

it later. > J-P > > > On Friday, October 7, 2016 at 12:19:39 PM UTC+2, Deepak Gupta wrote: >> >> Thanks JP for the elaborate clarification. I quickly tried to test a >> *distributed >> source term* commenting out the point load function and adding the >>

Re: [deal.II] Using p-refinement with high order elements

> VectorTools::interpolate_boundary_values (dof_handler, >> 0, >> ZeroFunction(), >> constraints); > > > On Friday, October 7, 2016 at 2:02:18 PM UTC+2,

Re: [deal.II] Using p-refinement with high order elements

s of using it in this context are (i.e. the interaction of > hp-constraints). > > Regards, > J-P > > On Friday, October 7, 2016 at 8:41:42 AM UTC+2, Deepak Gupta wrote: >> >> Dear JP, >> >> Thanks for looking at the example. The code itself is not a simple &

Re: [deal.II] Using p-refinement with high order elements

system RHS vector. > > I hope that this explains better what I think you might need to do to > correct the issue. > > J-P > > On Friday, October 7, 2016 at 11:26:18 AM UTC+2, Deepak Gupta wrote: >> >> Thanks JP. Indeed a concise example should help and that is where

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

is correct). However, now I have to figure out why set_boundary_indicator affects the number of constraints. Anyways have a nice weekend. Best Deepak On 27-08-16 03:57, Wolfgang Bangerth wrote: On 08/26/2016 03:49 PM, Deepak GUpta wrote: Indeed Wolfgang is right. Can't help but respond

[deal.II] Strange issue with KellyErrorEstimator functionality

Dear All, I am trying to use KellyErrorEstimator of deal.II and I saw various tutorials which use it. I ran those examples and they work well. Recently, I tried it with my work and when I print the estimated_error_per_cell vector, all of the values in that vector are zero. I have been trying

Re: [deal.II] Re: Strange issue with KellyErrorEstimator functionality

. Have you tried increasing the number of face quadrature points (I notice that qrule starts at 1, but without knowing which type and degree FE you are using, its hard to judge if the way you choose the quadrature order is correct). Regards, J-P On Friday, August 26, 2016 at 1:23:03 PM UTC+2,

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

Indeed Wolfgang is right. With the stated change, the example code of JP works now. So I need to relate the error in mine. Thanks Wolfgang and Jean-Paul for the help. Best Deepak On 26-08-16 23:08, Wolfgang Bangerth wrote: On 08/25/2016 02:34 AM, Jean-Paul Pelteret wrote: Hi Deepak, Ok,

[deal.II] Hanging node constraints

Dear All, Below is a simple piece of code where one finite element has a different p-order compared to the rest. Thus, I expect certain hanging_node_constraints (which is 4) in the output. What I cannot figure out is why I get *zero *constraints when I check the hanging constraints a second time?

Re: [deal.II] Hanging node constraints

e, but I would like to know why it was a problem (what goes on in the background as per the implementation?) Best Deepak On Tue, Sep 27, 2016 at 5:02 PM, Wolfgang Bangerth <bange...@colostate.edu> wrote: > On 09/27/2016 08:52 AM, Deepak Gupta wrote: > >> >> I created a simple

Re: [deal.II] coordinates of dof (Q1)

Dear Julian, This question has been asked in the past and there are already some threads which might be of help to you: https://groups.google.com/forum/#!searchin/dealii/dof$20coordinates%7Csort:relevance Best, Deepak On Thu, Nov 10, 2016 at 2:41 PM, Julian Dorn wrote:

Re: [deal.II] problem with user_index()

Dear Sudharshan, Probably you are looking for cell->index(). Am I right? best Deepak On Sun, Nov 13, 2016 at 7:03 PM, Sudarshan Kumar wrote: > > typename DoFHandler::active_cell_ > iterator > cell = dof_handler.begin_active(), > endc = dof_handler.end(); > >