Then, calling DoFTools::map_dofs_to_support_points(MappingQ1<dim>(),
dof_handler) gives an error in *debug mode*, but *no *errors in *release mode*!
There are ~14,000 checks in the library that only run in debug mode. They
typically verify that input arguments of functions are reasonable. You ran
into one of them. But just because you don't get an error in release mode
(because the check is not run in release mode) does not mean that the result
you get is reasonable.
The error message is as follows:
--------------------------------------------------------
Calling MPI_Abort now.
To break execution in a GDB session, execute 'break MPI_Abort' before running.
You can also put the following into your ~/.gdbinit:
set breakpoint pending on
break MPI_Abort
set breakpoint pending auto
--------------------------------------------------------
An error occurred in line <2791> of file
</proj/petsc/dealii-9.5.2/source/fe/fe_values.cc> in function
dealii::FEValuesBase<dim, spacedim>::FEValuesBase(unsigned int, unsigned
int, dealii::UpdateFlags, const dealii::Mapping<dim, spacedim>&, const
dealii::FiniteElement<dim, spacedim>&) [with int dim = 3; int spacedim = 3]
The violated condition was:
n_q_points > 0
Additional information:
There is nothing useful you can do with an FEValues object when using
a quadrature formula with zero quadrature points!
[...]
I'd consider this a bug. I opened a github issue here:
https://github.com/dealii/dealii/issues/18279
I don't know when someone will get around to fixing it, but at least the issue
is recorded.
Minimum code to reproduce the issue:
Great little test, thanks for coming up with a minimum piece of code that
shows the issue!
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
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/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 dealii+unsubscr...@googlegroups.com.
To view this discussion visit
https://groups.google.com/d/msgid/dealii/a3d5d5b3-390c-4121-99fe-fc03fd831aff%40colostate.edu.
