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