Seems like this was added in “65cde6a” on 5/17/14. Was this intended?
-Manav > On Oct 18, 2016, at 9:40 PM, Manav Bhatia <[email protected]> wrote: > > Hi, > > I am using SZABAB shape functions with QUAD9 element, and am tripping the > error in the "default:” case for the enum in lines 208-228 in > fe_szabab_shape_2D.C (shown below). > > I had successfully used this combination about 3 years ago, so I am not > sure what, if anything, has changed in libMesh for this. > > Could someone please comment on the intent of the error message? > > Thanks, > Manav > > > > > switch(i) > { > case 5: // edge 0 points > if (elem->point(0) > elem->point(1))f = -1.; > break; > case 7: // edge 1 points > if (elem->point(1) > elem->point(2))f = -1.; > break; > case 9: // edge 2 points > if (elem->point(3) > elem->point(2))f = -1.; > break; > case 11: // edge 3 points > if (elem->point(0) > elem->point(3))f = -1.; > break; > > default: > libmesh_error_msg("Invalid shape function index i = " << i); > } > > return f*(FE<1,SZABAB>::shape(EDGE3, totalorder, i0[i], xi)* > FE<1,SZABAB>::shape(EDGE3, totalorder, i1[i], eta)); > ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, SlashDot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
