diff --git a/laghos_assembly.cpp b/laghos_assembly.cpp
index 60ee07f3..930fe8d5 100644
--- a/laghos_assembly.cpp
+++ b/laghos_assembly.cpp
@@ -538,6 +538,7 @@ static void ForceMult(const int DIM, const int D1D, const int Q1D,
       {0x234,&ForceMult2D<2,3,4,2>},
       {0x246,&ForceMult2D<2,4,6,3>},
       {0x258,&ForceMult2D<2,5,8,4>},
+      {0x26A,&ForceMult2D<2,6,10,5>},
       // 3D
       {0x334,&ForceMult3D<3,3,4,2>},
       {0x346,&ForceMult3D<3,4,6,3>},
@@ -944,6 +945,7 @@ static void ForceMultTranspose(const int DIM, const int D1D, const int Q1D,
       {0x234,&ForceMultTranspose2D<2,3,4,2>},
       {0x246,&ForceMultTranspose2D<2,4,6,3>},
       {0x258,&ForceMultTranspose2D<2,5,8,4>},
+      {0x26A,&ForceMultTranspose2D<2,6,10,5>},
       {0x334,&ForceMultTranspose3D<3,3,4,2>},
       {0x346,&ForceMultTranspose3D<3,4,6,3>},
       {0x358,&ForceMultTranspose3D<3,5,8,4>}
diff --git a/laghos_solver.cpp b/laghos_solver.cpp
index fe9f9acd..18cde356 100644
--- a/laghos_solver.cpp
+++ b/laghos_solver.cpp
@@ -1318,7 +1318,10 @@ void QUpdate::UpdateQuadratureData(const Vector &S, QuadratureData &qdata)
                             Vector &dt_est, DenseTensor &stressJinvT);
    static std::unordered_map<int, fQKernel> qupdate =
    {
-      {0x24,&QKernel<2,4>}, {0x26,&QKernel<2,6>}, {0x28,&QKernel<2,8>},
+      // 2D.
+      {0x24,&QKernel<2,4>}, {0x26,&QKernel<2,6>},
+      {0x28,&QKernel<2,8>}, {0x2A,&QKernel<2,10>},
+      // 3D.
       {0x34,&QKernel<3,4>}, {0x36,&QKernel<3,6>}, {0x38,&QKernel<3,8>}
    };
    if (!qupdate[id])