On 06/24/2016 11:57 PM, Praveen C wrote:
MappingFEField<dim> map(dh_euler, euler_vector);
When “map" is destroyed, I get some error
An error occurred in line <128> of file
</Users/praveen/Applications/deal.II/git/source/base/subscriptor.cc
<https://urldefense.proofpoint.com/v2/url?u=http-3A__subscriptor.cc&d=CwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=bn1clvQAiDQrfAC4yKbN0PlSr7RWRs-U3WJ0zRBB2qM&m=Vf-2N0p6VmAXKK4t0RHa8AEgald8F3wLtzHtbbHTHDk&s=Q5qaqH8k_jfyOzP2rprL-TgG2uQe5NpA4zddILDSvuc&e=>>
in function
void dealii::Subscriptor::check_no_subscribers() const
The violated condition was:
counter == 0
The name and call sequence of the exception was:
ExcInUse (counter, object_info->name(), infostring)
Additional Information:
Object of class N6dealii6VectorIdEE is still used by 1 other objects.
What a coincidence -- I've explained this very error just 10 minutes to
someone sitting here next to me at the ASPECT hackathon :-)
The issue is that MappingFEField takes a dealii::Vector as second argument.
What you pass is a vector of a different class, so the constructor creates a
temporary vector to pass to the constructor, and then destroys it again after
the call. But the MappingFEField class has already taken its address, and so
you get the error above.
What you need to do is use the second template argument of the MappingFEField
class and set it to TrilinosWrappers::MPI::Vector. This should make it all work!
Best
Wolfgang
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@math.tamu.edu
www: http://www.math.tamu.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.
For more options, visit https://groups.google.com/d/optout.
