Thanks Wolfgang for the reply. I created a small scalar example by modifying Step-6 and to my surprise, it works very well (can send if someone needs it). So, I have not figured out the error yet. I am still trying to narrow down the problem.
By the way, I am constructing the cell matrices using B'DB, where B is the strain-displacement matrix and D is the constitutive relation. Then I add it to the system matrix using distribute_local_to_gobal function. Since B and D are the same for global and adaptive case, I do not expect any difference in the cell matrices. The system_rhs is correct since I have only a point load and I can check it by printing the vector. Thus, the possible issue can be only in the assembly (for the local refinement one). My problem is vector-values. In case this information can lead to any other suggestion from someone, it would be great. Meanwhile I am going 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 <[email protected]> wrote: > On 10/05/2016 09:03 AM, Deepak Gupta wrote: > >> >> >> Since it works for the combination of p=1 and 2, I do not think hanging >> nodes >> are a problem. I looked at the hp-paper as well and do not see something >> that >> might be going wrong related to the hanging nodes. I am still trying to >> narrow >> down the problem and come up with a simplified example, however, I still >> thought of raising this question. Based on the above scenario, can someone >> figure out what might be going wrong? >> > > I am now aware of anything that would cause this, but it certainly sounds > like a bug. A small testcase would definitely be useful! > > In the past, I have debugged problems such as this by interpolating a > known and simple function onto the mesh, and then applying the constraints > via ConstraintMatrix::distribute() to the vector with the interpolant. If > the constraints are wrong, you will see this if you just output this as a > field (with sufficiently high n_subdivisions passed to > DataOut::build_patches()). > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: [email protected] > www: http://www.math.colostate.edu/~bangerth/ > > > -- > The deal.II project is located at http://www.dealii.org/ > For mailing list/forum options, see https://groups.google.com/d/fo > rum/dealii?hl=en > --- You received this message because you are subscribed to the Google > Groups "deal.II User Group" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to [email protected]. > For more options, visit https://groups.google.com/d/optout. > -- 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 Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. For more options, visit https://groups.google.com/d/optout.
