Remove error messages from lagrange shape computations in 1D to allow their inlining

refs #4272

Author: GiudGiud

include/fe/fe_lagrange_shape_1D.h

@@ -39,11 +39,8 @@ Real fe_lagrange_1D_linear_shape(const unsigned int i,
     case 0:
       return .5*(1. - xi);
 
-    case 1:
-      return .5*(1. + xi);
-
     default:
-      libmesh_error_msg("Invalid shape function index i = " << i);
+      return .5*(1. + xi);
     }
 }
 
@@ -63,11 +60,8 @@ Real fe_lagrange_1D_quadratic_shape(const unsigned int i,
     case 1:
       return .5*xi*(xi + 1);
 
-    case 2:
-      return (1. - xi*xi);
-
     default:
-      libmesh_error_msg("Invalid shape function index i = " << i);
+      return (1. - xi*xi);
     }
 }
 
@@ -90,11 +84,8 @@ Real fe_lagrange_1D_cubic_shape(const unsigned int i,
     case 2:
       return 27./16.*(1.-xi*xi)*(1./3.-xi);
 
-    case 3:
-      return 27./16.*(1.-xi*xi)*(1./3.+xi);
-
     default:
-      libmesh_error_msg("Invalid shape function index i = " << i);
+      return 27./16.*(1.-xi*xi)*(1./3.+xi);
     }
 }
 
@@ -116,11 +107,8 @@ Real fe_lagrange_1D_shape(const Order order,
       return fe_lagrange_1D_quadratic_shape(i, xi);
 
       // Lagrange cubics
-    case THIRD:
-      return fe_lagrange_1D_cubic_shape(i, xi);
-
     default:
-      libmesh_error_msg("ERROR: Unsupported polynomial order = " << order);
+      return fe_lagrange_1D_cubic_shape(i, xi);
     }
 }
 
@@ -141,11 +129,8 @@ Real fe_lagrange_1D_linear_shape_deriv(const unsigned int i,
     case 0:
       return -.5;
 
-    case 1:
-      return .5;
-
     default:
-      libmesh_error_msg("Invalid shape function index i = " << i);
+      return .5;
     }
 }
 
@@ -168,11 +153,8 @@ Real fe_lagrange_1D_quadratic_shape_deriv(const unsigned int i,
     case 1:
       return xi+.5;
 
-    case 2:
-      return -2.*xi;
-
     default:
-      libmesh_error_msg("Invalid shape function index i = " << i);
+      return -2.*xi;
     }
 }
 
@@ -198,11 +180,8 @@ Real fe_lagrange_1D_cubic_shape_deriv(const unsigned int i,
     case 2:
       return 27./16.*(3.*xi*xi-2./3.*xi-1.);
 
-    case 3:
-      return 27./16.*(-3.*xi*xi-2./3.*xi+1.);
-
     default:
-      libmesh_error_msg("Invalid shape function index i = " << i);
+      return 27./16.*(-3.*xi*xi-2./3.*xi+1.);
     }
 }
 
@@ -222,11 +201,8 @@ Real fe_lagrange_1D_shape_deriv(const Order order,
     case SECOND:
       return fe_lagrange_1D_quadratic_shape_deriv(i, j, xi);
 
-    case THIRD:
-      return fe_lagrange_1D_cubic_shape_deriv(i, j, xi);
-
     default:
-      libmesh_error_msg("ERROR: Unsupported polynomial order = " << order);
+      return fe_lagrange_1D_cubic_shape_deriv(i, j, xi);
     }
 }
 
@@ -253,11 +229,8 @@ Real fe_lagrange_1D_quadratic_shape_second_deriv(const unsigned int i,
     case 1:
       return 1.;
 
-    case 2:
-      return -2.;
-
     default:
-      libmesh_error_msg("Invalid shape function index i = " << i);
+      return -2.;
     }
 }
 
@@ -282,11 +255,8 @@ Real fe_lagrange_1D_cubic_shape_second_deriv(const unsigned int i,
     case 2:
       return 27./16.*(6*xi-2./3.);
 
-    case 3:
-      return 27./16.*(-6*xi-2./3.);
-
     default:
-      libmesh_error_msg("Invalid shape function index i = " << i);
+      return 27./16.*(-6*xi-2./3.);
     }
 }
 
@@ -307,11 +277,8 @@ Real fe_lagrange_1D_shape_second_deriv(const Order order,
     case SECOND:
       return fe_lagrange_1D_quadratic_shape_second_deriv(i, j, xi);
 
-    case THIRD:
-      return fe_lagrange_1D_cubic_shape_second_deriv(i, j, xi);
-
     default:
-      libmesh_error_msg("ERROR: Unsupported polynomial order = " << order);
+      return fe_lagrange_1D_cubic_shape_second_deriv(i, j, xi);
     } // end switch (order)
 }