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(); > >for (; cell!=endc; +

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: > Dear all, > > if I hav

Re: [deal.II] Changing geometry

Hi, I think (a guess) you forgot to change the Coefficient function. For [-1, 1], the radius of 0.5 is around (0, 0). For [0, 1], I guess it needs to be adjusted accordingly. Best, Deepak On Wed, Nov 2, 2016 at 4:24 PM, A Banerjee wrote: > Dear All, > > Let me tell you that I am a beginner in

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

together (one after the other) in the whole dof_list. However, this was not true beyond p=2 and that is where my code failed. Best regards Deepak On Fri, Oct 7, 2016 at 2:22 PM, Deepak Gupta wrote: > Thanks. > > Just some add on, I figured out that the replacement for the lines > comment

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, Deepak Gupta wrote: >&

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

ive quadrature formula from FEValues. > - Run this in debug mode to check that there are no other mistakes that we > can't pick up by just looking at the code. > > If you do all of this and still have no success, then post the updated > code and I'll fiddle with it la

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

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 I >>

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

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 &g

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

n/deal.II/classDoFAccessor.html#a5560151b5407e4851d5c1009c7753764> > function > to retrieve the DoF associated with a vertex. There are a number of posts > <https://groups.google.com/forum/#!searchin/dealii/vertex_dof_index%7Csort:relevance> > that discuss the use of this func

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

ng 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 wrote: > On 10/05/2016 09:03 AM, Deepak Gupta wrote: > >> >> &g

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

Dear All, I am trying to solve an elasticity problem and am dealing with p-refinement. With global refinement, my solution is as expected (correct). I am using an initial mesh of 4X2 elements. If I keep p=1 for all the elements, and make p = 2 for the 7th element, the solution is still correct

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 wrote: > On 09/27/2016 08:52 AM, Deepak Gupta wrote: > >> >> I created a simple example but it worked corre

Re: [deal.II] Hanging node constraints

Thanks, after this explanation, it sounds extremely silly on my side :) Best, Deepak On Tue, Sep 27, 2016 at 6:08 PM, Wolfgang Bangerth wrote: > On 09/27/2016 10:05 AM, Deepak Gupta wrote: > >> >> >> If I remove the else part of the if-condition, it works, else the

Re: [deal.II] Hanging node constraints

: > On 09/27/2016 06:44 AM, Deepak Gupta wrote: > >> 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 out

[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] Strange issue with KellyErrorEstimator functionality

that this might be problem. Do I need to pass this information to the estimate() function, if yes, could you help me on how to do this? Best Deepak On Fri, Aug 26, 2016 at 4:04 PM, Wolfgang Bangerth wrote: > On 08/26/2016 05:23 AM, Deepak GUpta wrote: > >> >> I am trying to use K

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 re

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, so

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,

[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: Reference article for modified Kelly Estimator

Thanks Wolfgang, I had already seen these videos quickly and they are quite helpful. May be I will take another deep look to get more insight. Best Deepak On 25-08-16 00:58, Wolfgang Bangerth wrote: Deepak, I am using the KellyEstimator as an 'indicator' for simple hp-adaptive p

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

acent cells being on the same refinement level (i.e. there are no actual hanging _nodes_). Perhaps someone else could clarify this? J-P On Tuesday, August 23, 2016 at 9:45:41 AM UTC+2, Deepak Gupta wrote: Dear Daniel and others, Below is the part of my code that might explain better:

[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] 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

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 ou

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

Dear All, I had been using *create_point_source_vector()* function with* DoFHandler* and it had been working fine. Recently, I modified my code for hp-version and am now using *hp::DoFHandler*. With p=1, the code works fine. However, with p>1, the program breaks during runtime at *create_point_sou

Re: [deal.II] Re: Decoupling FECollection and QCollection

Thanks Simon and Jean-Paul, Your answers clarified by doubt and thanks Wolfgang for updating it already. Best regards Deepak On Sun, Jul 24, 2016 at 5:39 AM, Wolfgang Bangerth wrote: > > Deepak & J-P, > > Its possible that the initial implementation of the functionality has >> changed, >> and

Re: [deal.II] Re: Decoupling FECollection and QCollection

Dear Simon Thanks for the reply. I was guessing the same, but then this sentence from the documentation seems a bit ambiguous. Isn't it? "/The quadrature rules have to be added in the same order as for the FECollection

[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] 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 suppo

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

Dear All, Since I didn't receive any response to my question, I thought of raising it up as a reminder if someone can help. Best Deepak On Friday, July 15, 2016 at 2:02:43 PM UTC+2, Deepak Gupta wrote: > > Dear All, > > I am currently dealing with hanging support points in

[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 support