On Thu, Jul 9, 2020 at 11:06 AM Caleb M Phillips <calebphill...@utexas.edu>
wrote:

> John,
>
> I wasn't too sure how to accomplish that, but here is the main file. This
> is a module in a larger code, but the function I'm looking at is
> "assemble_pressure".
> https://github.com/CalebPhillips5/ABM_Ang/blob/master/pressure_ABM . In
> this instance I'm just trying to set the top boundary to 1 and the bottom
> boundary to 0 with 2nd order elements.
>

Looking quickly at the assembly code, it seems like it should be fine,
although I could be missing something subtle.

Are you hitting the "build_cube()" code path for the Mesh? If so, you will
need to change the element type to TRI6 in order to use SECOND-order FEs.
Similarly for any meshes you read in... one way to achieve this is by
generating/reading in a mesh with linear elements and then calling
all_second_order() to convert it.

If you aren't already, I'd also strongly suggest compiling and running in
dbg mode, as I believe an error like this should be caught by an assert in
dbg mode.

--
John

_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to