34#ifndef __FLUID_EULER_RES_KERNEL__
35#define __FLUID_EULER_RES_KERNEL__
48 for (
int i = idx;
i < n;
i +=
str) {
67 for (
int i = idx;
i < n;
i +=
str) {
68 f_x[
i] = m_x[
i] * m_x[
i] / rho_field[
i] + p[
i];
69 f_y[
i] = m_x[
i] * m_y[
i] / rho_field[
i];
70 f_z[
i] = m_x[
i] * m_z[
i] / rho_field[
i];
88 for (
int i = idx;
i < n;
i +=
str) {
89 f_x[
i] = m_y[
i] * m_x[
i] / rho_field[
i];
90 f_y[
i] = m_y[
i] * m_y[
i] / rho_field[
i] + p[
i];
91 f_z[
i] = m_y[
i] * m_z[
i] / rho_field[
i];
109 for (
int i = idx;
i < n;
i +=
str) {
110 f_x[
i] = m_z[
i] * m_x[
i] / rho_field[
i];
111 f_y[
i] = m_z[
i] * m_y[
i] / rho_field[
i];
112 f_z[
i] = m_z[
i] * m_z[
i] / rho_field[
i] + p[
i];
116template<
typename T >
131 for (
int i = idx;
i < n;
i +=
str) {
132 f_x[
i] = (
E[
i] + p[
i]) * m_x[
i] / rho_field[
i];
133 f_y[
i] = (
E[
i] + p[
i]) * m_y[
i] / rho_field[
i];
134 f_z[
i] = (
E[
i] + p[
i]) * m_z[
i] / rho_field[
i];
138template<
typename T >
150 for (
int i = idx;
i < n;
i +=
str) {
159template<
typename T >
177 for (
int i = idx;
i < n;
i +=
str) {
__global__ void const T *__restrict__ const T *__restrict__ const T *__restrict__ const T *__restrict__ dt
__global__ void dirichlet_apply_scalar_kernel(const int *__restrict__ msk, T *__restrict__ x, const T g, const int m)
__global__ void euler_res_part_visc_kernel(T *__restrict__ rhs, const T *__restrict__ Binv, const T *__restrict__ lap_sol, const T *__restrict__ h, const T c_avisc, const int n)
__global__ void euler_res_part_rk_sum_kernel(T *__restrict__ rho, T *__restrict__ m_x, T *__restrict__ m_y, T *__restrict__ m_z, T *__restrict__ E, const T *__restrict__ k_rho_i, const T *__restrict__ k_m_x_i, const T *__restrict__ k_m_y_i, const T *__restrict__ k_m_z_i, const T *__restrict__ k_E_i, const T dt, const T c, const int n)
__global__ void euler_res_part_my_flux_kernel(T *__restrict__ f_x, T *__restrict__ f_y, T *__restrict__ f_z, const T *__restrict__ m_x, const T *__restrict__ m_y, const T *__restrict__ m_z, const T *__restrict__ rho_field, const T *__restrict__ p, const int n)
__global__ void euler_res_part_E_flux_kernel(T *__restrict__ f_x, T *__restrict__ f_y, T *__restrict__ f_z, const T *__restrict__ m_x, const T *__restrict__ m_y, const T *__restrict__ m_z, const T *__restrict__ rho_field, const T *__restrict__ p, const T *__restrict__ E, const int n)
__global__ void euler_res_part_mz_flux_kernel(T *__restrict__ f_x, T *__restrict__ f_y, T *__restrict__ f_z, const T *__restrict__ m_x, const T *__restrict__ m_y, const T *__restrict__ m_z, const T *__restrict__ rho_field, const T *__restrict__ p, const int n)
__global__ void euler_res_part_coef_mult_kernel(T *__restrict__ rhs_rho, T *__restrict__ rhs_m_x, T *__restrict__ rhs_m_y, T *__restrict__ rhs_m_z, T *__restrict__ rhs_E, const T *__restrict__ mult, const int n)
__global__ void euler_res_part_mx_flux_kernel(T *__restrict__ f_x, T *__restrict__ f_y, T *__restrict__ f_z, const T *__restrict__ m_x, const T *__restrict__ m_y, const T *__restrict__ m_z, const T *__restrict__ rho_field, const T *__restrict__ p, const int n)