1 // Copyright (c) 2017-2026, Lawrence Livermore National Security, LLC and other CEED contributors.
2 // All Rights Reserved. See the top-level LICENSE and NOTICE files for details.
4 // SPDX-License-Identifier: BSD-2-Clause
20   CeedScalar(*jac_data_sur)               = newt_ctx->is_implicit ? out[1] : NULL;  in Slip()
24     State            s     = StateFromQ(newt_ctx, qi, state_var);  in Slip()  local
28     wdetJb *= newt_ctx->is_implicit ? -1. : 1.;  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()
33 …const CeedScalar Y_reflect[5] = {s.Y.pressure, vel_reflect[0], vel_reflect[1], vel_reflect[2], s.Y…  in Slip()
36     StateConservative flux = RiemannFlux_HLLC(newt_ctx, s, s_reflect, norm);  in Slip()
40     for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * Flux[j];  in Slip()
42     if (newt_ctx->is_implicit) {  in Slip()
75     wdetJb *= newt_ctx->is_implicit ? -1. : 1.;  in Slip_Jacobian()
80     State s  = StateFromQ(newt_ctx, qi, state_var);  in Slip_Jacobian()  local
81     State ds = StateFromQ_fwd(newt_ctx, s, dqi, state_var);  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()
86 …const CeedScalar Y_reflect[5] = {s.Y.pressure, vel_reflect[0], vel_reflect[1], vel_reflect[2], s.Y…  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()
99     for (CeedInt j = 0; j < 5; j++) v[j][i] = -wdetJb * dFlux[j];  in Slip_Jacobian()