Hi All, I wanted to update DoFTools::make_hanging_node_constraints is not the issue above. Even I comment out that line and call ConstraintMatrix.close() prior to calling make_periodicity_constraints, I still get the same error. This means calling ConstraintMatrix.close() on a constraintMatrix with no constraints is leading to the error.
Thanks, Sambit On Sunday, December 10, 2017 at 4:30:52 PM UTC-5, Sambit Das wrote: > > Hi All, > > I am facing the following issue: I am setting periodic boundary conditions > in all directions on a dofHandler object attached to a single element > triangulation. When I print the ConstraintMatrix I observe that a trivial > call to DoFTools::make_hanging_node_constraints(..) which I call prior to > calling DoFTools::make_periodicity_constraints gives an erroneous > ConstraintMatrix as shown below for the 8 degrees of freedom corresponding > to the 8 nodes. > > 0 1: 1 > 2 3: 1 > 4 5: 1 > 6 7: 1 > 1 3: 1 > 5 7: 1 > 3 7: 1 > > However, if the don't call ConstraintMatrix.close() (line 102 in the > minimal example) after the trivial call to > DoFTools::make_hanging_node_constraints(..), I get the correct > ConstraintMatrix. Likewise, if I don't use > make_hanging_node_constraints(..) the ConstraintMatrix is correct too > 0 7: 1 > 1 7: 1 > 2 7: 1 > 3 7: 1 > 4 7: 1 > 5 7: 1 > 6 7: 1 > > I also checked that if I printed the ConstraintMatrix after just setting > the hanging_node_constraints, it doesn't print anything. > I wonder if I am making a mistake here. I have provided a minimal working > example file which reproduces this error. > > Thanks, > Sambit > -- 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.
