I'd like to be able to impose boundary conditions on an edge. We added
support for this in libMesh, but the only way I could come up with for
creating edge BCs in a mesh generator is to:
1. Create a mesh with extra 1D elements along the desired edge
2. Convert those 1D elements into edge BCs and delete the 1D elements
I've attached code that does this, with an associated test mesh. But after
the edge to BC conversion is complete, the code hits an assert in dbg mode
when I make a FEMContext (stack trace is below), and in opt mode it reports
a memory leak.
If anyone has any thoughts on what might be wrong here, I'd be interested
to hear.
David
----------------
Stack trace:
Assertion `q' failed.
Stack frames: 7
0: libMesh::print_trace(std::ostream&)
1: libMesh::MacroFunctions::report_error(char const*, int, char const*,
char const*)
2: libMesh::FE<1u,
(libMesh::FEFamily)0>::attach_quadrature_rule(libMesh::QBase*)
3: libMesh::FEMContext::FEMContext(libMesh::System const&)
4: ./example-dbg() [0x4416ac]
5: __libc_start_main
6: ./example-dbg() [0x43ff89]
[0] /home/dknez/software/libmesh-src/src/fe/fe.C, line 46, compiled Apr 22
2015 at 15:28:47
------------------------------------------------------------------------------
One dashboard for servers and applications across Physical-Virtual-Cloud
Widest out-of-the-box monitoring support with 50+ applications
Performance metrics, stats and reports that give you Actionable Insights
Deep dive visibility with transaction tracing using APM Insight.
http://ad.doubleclick.net/ddm/clk/290420510;117567292;y
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
