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
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
