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

Reply via email to