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