26     CeedScalar wdetJb, norm[3];  in Slip()  local
27     QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, NULL, norm);  in Slip()
31     const CeedScalar vel_normal = Dot3(s.Y.velocity, norm);  in Slip()
32     for (CeedInt j = 0; j < 3; j++) vel_reflect[j] = s.Y.velocity[j] - 2. * norm[j] * vel_normal;  in Slip()
36     StateConservative flux = RiemannFlux_HLLC(newt_ctx, s, s_reflect, norm);  in Slip()
73     CeedScalar wdetJb, norm[3];  in Slip_Jacobian()  local
74     QdataBoundaryUnpack_3D(Q, i, q_data_sur, &wdetJb, NULL, norm);  in Slip_Jacobian()
84     const CeedScalar vel_normal = Dot3(s.Y.velocity, norm);  in Slip_Jacobian()
85     for (CeedInt j = 0; j < 3; j++) vel_reflect[j] = s.Y.velocity[j] - 2. * norm[j] * vel_normal;  in Slip_Jacobian()
90     const CeedScalar dvel_normal = Dot3(ds.Y.velocity, norm);  in Slip_Jacobian()
91     for (CeedInt j = 0; j < 3; j++) dvel_reflect[j] = ds.Y.velocity[j] - 2. * norm[j] * dvel_normal;  in Slip_Jacobian()
95     StateConservative dflux = RiemannFlux_HLLC_fwd(newt_ctx, s, ds, s_reflect, ds_reflect, norm);  in Slip_Jacobian()