Sambit,
just for confirmation:
the first problem with your test is that calling the constructor FE_Q with
a (one-dimensional) Quadrature formula requires the first point to be 0 and
the last one 1 as described in the documentation.
<https://www.dealii.org/8.5.1/doxygen/deal.II/classFE__Q.html#a32ccedfe518cfe3fca02cdddf9e664f6>
Apart from that you probably get the error message:
An error occurred in line <1510> of file
<../include/deal.II/lac/constraint_matrix.h> in function
void dealii::ConstraintMatrix::add_line(const
dealii::ConstraintMatrix::size_type)
The violated condition was:
sorted==false
Additional information:
(none)
telling you that ConstraintMatrix::close() was called and you are not
allowed to add any constraints anymore as described here
<https://www.dealii.org/8.5.1/doxygen/deal.II/classConstraintMatrix.html#ada29d922494c2b645e418cc311ec872d>
.
In particular, removing the call to DoFTools::make_periodicity_constraints
makes the test run fine even when calling ConstraintMatrix::close()
a second time.
Can you confirm this?
Best,
Daniel
--
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.
