Last active
October 13, 2022 23:41
-
-
Save kaushikcfd/13dcb39130c7e8f5195c224b751a8b90 to your computer and use it in GitHub Desktop.
2D DG-FEM Grad Operator over 16 ranks
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs(double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_all, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_all, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm4_0; | |
| double acc_idof_ensm4_1; | |
| double acc_iface_ensm3_idof_ensm3; | |
| double acc_iface_ensm3_idof_ensm3_0; | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| for (int iel_ensm1 = 0; iel_ensm1 <= 1; ++iel_ensm1) | |
| cse_0[idof_ensm1 + 4 * iel_ensm1] = 0.5 * (_pt_part_ph_id_3[idof_ensm1 + 4 * iel_ensm1] + _pt_dist_id_0[4 * from_el_indices_3[iel_ensm1] + _pt_data_7[idof_ensm1 + 4 * iel_ensm1]]); | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| for (int iel_ensm2 = 0; iel_ensm2 <= 4; ++iel_ensm2) | |
| cse[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_2[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id[4 * from_el_indices_1[iel_ensm2] + _pt_data_5[idof_ensm2 + 4 * iel_ensm2]]); | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| for (int iel_ensm0 = 0; iel_ensm0 <= 31; ++iel_ensm0) | |
| for (int iface_ensm0 = 0; iface_ensm0 <= 2; ++iface_ensm0) | |
| { | |
| _pt_temp_1[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0] = (from_el_present[iel_ensm0 + 32 * iface_ensm0] ? normal_1_b_all[from_el_indices[iel_ensm0 + 32 * iface_ensm0]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]]) : 0.0) + (from_el_present_2[iel_ensm0 + 32 * iface_ensm0] ? cse_0[4 * from_el_indices_4[iel_ensm0 + 32 * iface_ensm0] + _pt_data_8[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_0[iel_ensm0 + 32 * iface_ensm0] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm0 + 32 * iface_ensm0] + _pt_data_4[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_face_restr_interior[from_el_indices_0[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_1[iel_ensm0 + 32 * iface_ensm0] ? cse[4 * from_el_indices_2[iel_ensm0 + 32 * iface_ensm0] + _pt_data_6[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION[from_el_indices_2[iel_ensm0 + 32 * iface_ensm0]] : 0.0); | |
| _pt_temp[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0] = (from_el_present[iel_ensm0 + 32 * iface_ensm0] ? normal_2_b_all[from_el_indices[iel_ensm0 + 32 * iface_ensm0]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]]) : 0.0) + (from_el_present_2[iel_ensm0 + 32 * iface_ensm0] ? cse_0[4 * from_el_indices_4[iel_ensm0 + 32 * iface_ensm0] + _pt_data_8[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_0[iel_ensm0 + 32 * iface_ensm0] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm0 + 32 * iface_ensm0] + _pt_data_4[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_face_restr_interior[from_el_indices_0[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_1[iel_ensm0 + 32 * iface_ensm0] ? cse[4 * from_el_indices_2[iel_ensm0 + 32 * iface_ensm0] + _pt_data_6[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION[from_el_indices_2[iel_ensm0 + 32 * iface_ensm0]] : 0.0); | |
| } | |
| for (int idof_ensm3_0 = 0; idof_ensm3_0 <= 9; ++idof_ensm3_0) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 31; ++iel_ensm3) | |
| { | |
| acc_iface_ensm3_idof_ensm3_0 = 0.0; | |
| acc_iface_ensm3_idof_ensm3 = 0.0; | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 3; ++idof_ensm3) | |
| for (int iface_ensm3 = 0; iface_ensm3 <= 2; ++iface_ensm3) | |
| { | |
| acc_iface_ensm3_idof_ensm3_0 = acc_iface_ensm3_idof_ensm3_0 + _pt_data_1[idof_ensm3 + 4 * iface_ensm3 + 12 * idof_ensm3_0] * _pt_data_2[32 * iface_ensm3 + iel_ensm3] * _pt_temp_1[idof_ensm3 + 128 * iface_ensm3 + 4 * iel_ensm3]; | |
| acc_iface_ensm3_idof_ensm3 = acc_iface_ensm3_idof_ensm3 + _pt_data_1[idof_ensm3 + 4 * iface_ensm3 + 12 * idof_ensm3_0] * _pt_data_2[32 * iface_ensm3 + iel_ensm3] * _pt_temp[idof_ensm3 + 128 * iface_ensm3 + 4 * iel_ensm3]; | |
| } | |
| face_mass_0[idof_ensm3_0 + 10 * iel_ensm3] = acc_iface_ensm3_idof_ensm3_0; | |
| face_mass[idof_ensm3_0 + 10 * iel_ensm3] = acc_iface_ensm3_idof_ensm3; | |
| } | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 9; ++idof_ensm5) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| { | |
| mass_inv_inp_0[idof_ensm5 + 10 * iel_ensm5] = _pt_part_ph_id[idof_ensm5 + 10 * iel_ensm5] + -1.0 * face_mass_0[idof_ensm5 + 10 * iel_ensm5]; | |
| mass_inv_inp[idof_ensm5 + 10 * iel_ensm5] = _pt_part_ph_id[320 + idof_ensm5 + 10 * iel_ensm5] + -1.0 * face_mass[idof_ensm5 + 10 * iel_ensm5]; | |
| } | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 9; ++idof_ensm4) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 31; ++iel_ensm4) | |
| { | |
| acc_idof_ensm4_0 = 0.0; | |
| acc_idof_ensm4_1 = 0.0; | |
| for (int idof_ensm4_0 = 0; idof_ensm4_0 <= 9; ++idof_ensm4_0) | |
| { | |
| acc_idof_ensm4_0 = acc_idof_ensm4_0 + _pt_data[iel_ensm4] * _pt_data_0[idof_ensm4_0 + 10 * idof_ensm4] * mass_inv_inp[idof_ensm4_0 + 10 * iel_ensm4]; | |
| acc_idof_ensm4_1 = acc_idof_ensm4_1 + _pt_data[iel_ensm4] * _pt_data_0[idof_ensm4_0 + 10 * idof_ensm4] * mass_inv_inp_0[idof_ensm4_0 + 10 * iel_ensm4]; | |
| } | |
| _pt_temp_0[idof_ensm4 + 10 * iel_ensm4] = acc_idof_ensm4_0; | |
| _pt_temp_2[idof_ensm4 + 10 * iel_ensm4] = acc_idof_ensm4_1; | |
| } | |
| for (int idof_ensm6 = 0; idof_ensm6 <= 9; ++idof_ensm6) | |
| for (int iel_ensm6 = 0; iel_ensm6 <= 31; ++iel_ensm6) | |
| { | |
| _pt_out_0_0[idof_ensm6 + 10 * iel_ensm6] = -1.0 * _pt_temp_2[idof_ensm6 + 10 * iel_ensm6]; | |
| _pt_out_1_0[idof_ensm6 + 10 * iel_ensm6] = -1.0 * _pt_temp_0[idof_ensm6 + 10 * iel_ensm6]; | |
| } | |
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs(double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_all, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, char const *__restrict__ from_el_present_3, double const *__restrict__ _pt_part_ph_id_4, double const *__restrict__ _pt_dist_id_1, int const *__restrict__ from_el_indices_5, int const *__restrict__ _pt_data_9, double const *__restrict__ normal_2_b_BTAG_PARTITION_1, int const *__restrict__ from_el_indices_6, int const *__restrict__ _pt_data_10, char const *__restrict__ from_el_present_4, double const *__restrict__ _pt_part_ph_id_5, double const *__restrict__ _pt_dist_id_2, int const *__restrict__ from_el_indices_7, int const *__restrict__ _pt_data_11, double const *__restrict__ normal_2_b_BTAG_PARTITION_2, int const *__restrict__ from_el_indices_8, int const *__restrict__ _pt_data_12, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_all, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double const *__restrict__ normal_1_b_BTAG_PARTITION_1, double const *__restrict__ normal_1_b_BTAG_PARTITION_2, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ cse_1, double *__restrict__ cse_2, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm6_0; | |
| double acc_idof_ensm6_1; | |
| double acc_iface_ensm4_idof_ensm4; | |
| double acc_iface_ensm4_idof_ensm4_0; | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| for (int iel_ensm1 = 0; iel_ensm1 <= 1; ++iel_ensm1) | |
| { | |
| cse_1[idof_ensm1 + 4 * iel_ensm1] = 0.5 * (_pt_part_ph_id_4[idof_ensm1 + 4 * iel_ensm1] + _pt_dist_id_1[4 * from_el_indices_5[iel_ensm1] + _pt_data_9[idof_ensm1 + 4 * iel_ensm1]]); | |
| cse_2[idof_ensm1 + 4 * iel_ensm1] = 0.5 * (_pt_part_ph_id_5[idof_ensm1 + 4 * iel_ensm1] + _pt_dist_id_2[4 * from_el_indices_7[iel_ensm1] + _pt_data_11[idof_ensm1 + 4 * iel_ensm1]]); | |
| } | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| for (int iel_ensm2 = 0; iel_ensm2 <= 2; ++iel_ensm2) | |
| cse_0[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_3[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id_0[4 * from_el_indices_3[iel_ensm2] + _pt_data_7[idof_ensm2 + 4 * iel_ensm2]]); | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 3; ++idof_ensm3) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 4; ++iel_ensm3) | |
| cse[idof_ensm3 + 4 * iel_ensm3] = 0.5 * (_pt_part_ph_id_2[idof_ensm3 + 4 * iel_ensm3] + _pt_dist_id[4 * from_el_indices_1[iel_ensm3] + _pt_data_5[idof_ensm3 + 4 * iel_ensm3]]); | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| for (int iel_ensm0 = 0; iel_ensm0 <= 31; ++iel_ensm0) | |
| for (int iface_ensm0 = 0; iface_ensm0 <= 2; ++iface_ensm0) | |
| { | |
| _pt_temp_1[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0] = (from_el_present[iel_ensm0 + 32 * iface_ensm0] ? normal_1_b_all[from_el_indices[iel_ensm0 + 32 * iface_ensm0]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]]) : 0.0) + (from_el_present_4[iel_ensm0 + 32 * iface_ensm0] ? cse_2[4 * from_el_indices_8[iel_ensm0 + 32 * iface_ensm0] + _pt_data_12[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION_2[from_el_indices_8[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_3[iel_ensm0 + 32 * iface_ensm0] ? cse_1[4 * from_el_indices_6[iel_ensm0 + 32 * iface_ensm0] + _pt_data_10[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_2[iel_ensm0 + 32 * iface_ensm0] ? cse_0[4 * from_el_indices_4[iel_ensm0 + 32 * iface_ensm0] + _pt_data_8[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_0[iel_ensm0 + 32 * iface_ensm0] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm0 + 32 * iface_ensm0] + _pt_data_4[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_face_restr_interior[from_el_indices_0[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_1[iel_ensm0 + 32 * iface_ensm0] ? cse[4 * from_el_indices_2[iel_ensm0 + 32 * iface_ensm0] + _pt_data_6[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION[from_el_indices_2[iel_ensm0 + 32 * iface_ensm0]] : 0.0); | |
| _pt_temp[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0] = (from_el_present[iel_ensm0 + 32 * iface_ensm0] ? normal_2_b_all[from_el_indices[iel_ensm0 + 32 * iface_ensm0]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]]) : 0.0) + (from_el_present_4[iel_ensm0 + 32 * iface_ensm0] ? cse_2[4 * from_el_indices_8[iel_ensm0 + 32 * iface_ensm0] + _pt_data_12[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION_2[from_el_indices_8[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_3[iel_ensm0 + 32 * iface_ensm0] ? cse_1[4 * from_el_indices_6[iel_ensm0 + 32 * iface_ensm0] + _pt_data_10[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_2[iel_ensm0 + 32 * iface_ensm0] ? cse_0[4 * from_el_indices_4[iel_ensm0 + 32 * iface_ensm0] + _pt_data_8[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_0[iel_ensm0 + 32 * iface_ensm0] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm0 + 32 * iface_ensm0] + _pt_data_4[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_face_restr_interior[from_el_indices_0[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_1[iel_ensm0 + 32 * iface_ensm0] ? cse[4 * from_el_indices_2[iel_ensm0 + 32 * iface_ensm0] + _pt_data_6[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION[from_el_indices_2[iel_ensm0 + 32 * iface_ensm0]] : 0.0); | |
| } | |
| for (int idof_ensm4_0 = 0; idof_ensm4_0 <= 9; ++idof_ensm4_0) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 31; ++iel_ensm4) | |
| { | |
| acc_iface_ensm4_idof_ensm4_0 = 0.0; | |
| acc_iface_ensm4_idof_ensm4 = 0.0; | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 3; ++idof_ensm4) | |
| for (int iface_ensm4 = 0; iface_ensm4 <= 2; ++iface_ensm4) | |
| { | |
| acc_iface_ensm4_idof_ensm4_0 = acc_iface_ensm4_idof_ensm4_0 + _pt_data_1[idof_ensm4 + 4 * iface_ensm4 + 12 * idof_ensm4_0] * _pt_data_2[32 * iface_ensm4 + iel_ensm4] * _pt_temp_1[idof_ensm4 + 128 * iface_ensm4 + 4 * iel_ensm4]; | |
| acc_iface_ensm4_idof_ensm4 = acc_iface_ensm4_idof_ensm4 + _pt_data_1[idof_ensm4 + 4 * iface_ensm4 + 12 * idof_ensm4_0] * _pt_data_2[32 * iface_ensm4 + iel_ensm4] * _pt_temp[idof_ensm4 + 128 * iface_ensm4 + 4 * iel_ensm4]; | |
| } | |
| face_mass_0[idof_ensm4_0 + 10 * iel_ensm4] = acc_iface_ensm4_idof_ensm4_0; | |
| face_mass[idof_ensm4_0 + 10 * iel_ensm4] = acc_iface_ensm4_idof_ensm4; | |
| } | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 9; ++idof_ensm5) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| { | |
| mass_inv_inp_0[idof_ensm5 + 10 * iel_ensm5] = _pt_part_ph_id[idof_ensm5 + 10 * iel_ensm5] + -1.0 * face_mass_0[idof_ensm5 + 10 * iel_ensm5]; | |
| mass_inv_inp[idof_ensm5 + 10 * iel_ensm5] = _pt_part_ph_id[320 + idof_ensm5 + 10 * iel_ensm5] + -1.0 * face_mass[idof_ensm5 + 10 * iel_ensm5]; | |
| } | |
| for (int idof_ensm6 = 0; idof_ensm6 <= 9; ++idof_ensm6) | |
| for (int iel_ensm6 = 0; iel_ensm6 <= 31; ++iel_ensm6) | |
| { | |
| acc_idof_ensm6_0 = 0.0; | |
| acc_idof_ensm6_1 = 0.0; | |
| for (int idof_ensm6_0 = 0; idof_ensm6_0 <= 9; ++idof_ensm6_0) | |
| { | |
| acc_idof_ensm6_0 = acc_idof_ensm6_0 + _pt_data[iel_ensm6] * _pt_data_0[idof_ensm6_0 + 10 * idof_ensm6] * mass_inv_inp[idof_ensm6_0 + 10 * iel_ensm6]; | |
| acc_idof_ensm6_1 = acc_idof_ensm6_1 + _pt_data[iel_ensm6] * _pt_data_0[idof_ensm6_0 + 10 * idof_ensm6] * mass_inv_inp_0[idof_ensm6_0 + 10 * iel_ensm6]; | |
| } | |
| _pt_temp_0[idof_ensm6 + 10 * iel_ensm6] = acc_idof_ensm6_0; | |
| _pt_temp_2[idof_ensm6 + 10 * iel_ensm6] = acc_idof_ensm6_1; | |
| } | |
| for (int idof_ensm7 = 0; idof_ensm7 <= 9; ++idof_ensm7) | |
| for (int iel_ensm7 = 0; iel_ensm7 <= 31; ++iel_ensm7) | |
| { | |
| _pt_out_0_0[idof_ensm7 + 10 * iel_ensm7] = -1.0 * _pt_temp_2[idof_ensm7 + 10 * iel_ensm7]; | |
| _pt_out_1_0[idof_ensm7 + 10 * iel_ensm7] = -1.0 * _pt_temp_0[idof_ensm7 + 10 * iel_ensm7]; | |
| } | |
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs(double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_all, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, char const *__restrict__ from_el_present_3, double const *__restrict__ _pt_part_ph_id_4, double const *__restrict__ _pt_dist_id_1, int const *__restrict__ from_el_indices_5, int const *__restrict__ _pt_data_9, double const *__restrict__ normal_2_b_BTAG_PARTITION_1, int const *__restrict__ from_el_indices_6, int const *__restrict__ _pt_data_10, char const *__restrict__ from_el_present_4, double const *__restrict__ _pt_part_ph_id_5, double const *__restrict__ _pt_dist_id_2, int const *__restrict__ from_el_indices_7, int const *__restrict__ _pt_data_11, double const *__restrict__ normal_2_b_BTAG_PARTITION_2, int const *__restrict__ from_el_indices_8, int const *__restrict__ _pt_data_12, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_all, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double const *__restrict__ normal_1_b_BTAG_PARTITION_1, double const *__restrict__ normal_1_b_BTAG_PARTITION_2, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ cse_1, double *__restrict__ cse_2, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm8_0; | |
| double acc_idof_ensm8_1; | |
| double acc_iface_ensm5_idof_ensm5; | |
| double acc_iface_ensm5_idof_ensm5_0; | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| for (int iel_ensm1 = 0; iel_ensm1 <= 6; ++iel_ensm1) | |
| cse_0[idof_ensm1 + 4 * iel_ensm1] = 0.5 * (_pt_part_ph_id_3[idof_ensm1 + 4 * iel_ensm1] + _pt_dist_id_0[4 * from_el_indices_3[iel_ensm1] + _pt_data_7[idof_ensm1 + 4 * iel_ensm1]]); | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| { | |
| int const iel_ensm2 = 0; | |
| cse_1[idof_ensm2] = 0.5 * (_pt_part_ph_id_4[idof_ensm2] + _pt_dist_id_1[4 * from_el_indices_5[0] + _pt_data_9[idof_ensm2]]); | |
| } | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 3; ++idof_ensm3) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 3; ++iel_ensm3) | |
| cse_2[idof_ensm3 + 4 * iel_ensm3] = 0.5 * (_pt_part_ph_id_5[idof_ensm3 + 4 * iel_ensm3] + _pt_dist_id_2[4 * from_el_indices_7[iel_ensm3] + _pt_data_11[idof_ensm3 + 4 * iel_ensm3]]); | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 3; ++idof_ensm4) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 5; ++iel_ensm4) | |
| cse[idof_ensm4 + 4 * iel_ensm4] = 0.5 * (_pt_part_ph_id_2[idof_ensm4 + 4 * iel_ensm4] + _pt_dist_id[4 * from_el_indices_1[iel_ensm4] + _pt_data_5[idof_ensm4 + 4 * iel_ensm4]]); | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| for (int iel_ensm0 = 0; iel_ensm0 <= 31; ++iel_ensm0) | |
| for (int iface_ensm0 = 0; iface_ensm0 <= 2; ++iface_ensm0) | |
| { | |
| _pt_temp_1[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0] = (from_el_present[iel_ensm0 + 32 * iface_ensm0] ? normal_1_b_all[from_el_indices[iel_ensm0 + 32 * iface_ensm0]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]]) : 0.0) + (from_el_present_4[iel_ensm0 + 32 * iface_ensm0] ? cse_2[4 * from_el_indices_8[iel_ensm0 + 32 * iface_ensm0] + _pt_data_12[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION_2[from_el_indices_8[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_3[iel_ensm0 + 32 * iface_ensm0] ? cse_1[4 * from_el_indices_6[iel_ensm0 + 32 * iface_ensm0] + _pt_data_10[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_2[iel_ensm0 + 32 * iface_ensm0] ? cse_0[4 * from_el_indices_4[iel_ensm0 + 32 * iface_ensm0] + _pt_data_8[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_0[iel_ensm0 + 32 * iface_ensm0] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm0 + 32 * iface_ensm0] + _pt_data_4[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_face_restr_interior[from_el_indices_0[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_1[iel_ensm0 + 32 * iface_ensm0] ? cse[4 * from_el_indices_2[iel_ensm0 + 32 * iface_ensm0] + _pt_data_6[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION[from_el_indices_2[iel_ensm0 + 32 * iface_ensm0]] : 0.0); | |
| _pt_temp[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0] = (from_el_present[iel_ensm0 + 32 * iface_ensm0] ? normal_2_b_all[from_el_indices[iel_ensm0 + 32 * iface_ensm0]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]]) : 0.0) + (from_el_present_4[iel_ensm0 + 32 * iface_ensm0] ? cse_2[4 * from_el_indices_8[iel_ensm0 + 32 * iface_ensm0] + _pt_data_12[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION_2[from_el_indices_8[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_3[iel_ensm0 + 32 * iface_ensm0] ? cse_1[4 * from_el_indices_6[iel_ensm0 + 32 * iface_ensm0] + _pt_data_10[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_2[iel_ensm0 + 32 * iface_ensm0] ? cse_0[4 * from_el_indices_4[iel_ensm0 + 32 * iface_ensm0] + _pt_data_8[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_0[iel_ensm0 + 32 * iface_ensm0] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm0 + 32 * iface_ensm0] + _pt_data_4[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_face_restr_interior[from_el_indices_0[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_1[iel_ensm0 + 32 * iface_ensm0] ? cse[4 * from_el_indices_2[iel_ensm0 + 32 * iface_ensm0] + _pt_data_6[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION[from_el_indices_2[iel_ensm0 + 32 * iface_ensm0]] : 0.0); | |
| } | |
| for (int idof_ensm5_0 = 0; idof_ensm5_0 <= 9; ++idof_ensm5_0) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| { | |
| acc_iface_ensm5_idof_ensm5_0 = 0.0; | |
| acc_iface_ensm5_idof_ensm5 = 0.0; | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 3; ++idof_ensm5) | |
| for (int iface_ensm5 = 0; iface_ensm5 <= 2; ++iface_ensm5) | |
| { | |
| acc_iface_ensm5_idof_ensm5_0 = acc_iface_ensm5_idof_ensm5_0 + _pt_data_1[idof_ensm5 + 4 * iface_ensm5 + 12 * idof_ensm5_0] * _pt_data_2[32 * iface_ensm5 + iel_ensm5] * _pt_temp_1[idof_ensm5 + 128 * iface_ensm5 + 4 * iel_ensm5]; | |
| acc_iface_ensm5_idof_ensm5 = acc_iface_ensm5_idof_ensm5 + _pt_data_1[idof_ensm5 + 4 * iface_ensm5 + 12 * idof_ensm5_0] * _pt_data_2[32 * iface_ensm5 + iel_ensm5] * _pt_temp[idof_ensm5 + 128 * iface_ensm5 + 4 * iel_ensm5]; | |
| } | |
| face_mass_0[idof_ensm5_0 + 10 * iel_ensm5] = acc_iface_ensm5_idof_ensm5_0; | |
| face_mass[idof_ensm5_0 + 10 * iel_ensm5] = acc_iface_ensm5_idof_ensm5; | |
| } | |
| for (int idof_ensm7 = 0; idof_ensm7 <= 9; ++idof_ensm7) | |
| for (int iel_ensm7 = 0; iel_ensm7 <= 31; ++iel_ensm7) | |
| { | |
| mass_inv_inp_0[idof_ensm7 + 10 * iel_ensm7] = _pt_part_ph_id[idof_ensm7 + 10 * iel_ensm7] + -1.0 * face_mass_0[idof_ensm7 + 10 * iel_ensm7]; | |
| mass_inv_inp[idof_ensm7 + 10 * iel_ensm7] = _pt_part_ph_id[320 + idof_ensm7 + 10 * iel_ensm7] + -1.0 * face_mass[idof_ensm7 + 10 * iel_ensm7]; | |
| } | |
| for (int idof_ensm8 = 0; idof_ensm8 <= 9; ++idof_ensm8) | |
| for (int iel_ensm8 = 0; iel_ensm8 <= 31; ++iel_ensm8) | |
| { | |
| acc_idof_ensm8_0 = 0.0; | |
| acc_idof_ensm8_1 = 0.0; | |
| for (int idof_ensm8_0 = 0; idof_ensm8_0 <= 9; ++idof_ensm8_0) | |
| { | |
| acc_idof_ensm8_0 = acc_idof_ensm8_0 + _pt_data[iel_ensm8] * _pt_data_0[idof_ensm8_0 + 10 * idof_ensm8] * mass_inv_inp[idof_ensm8_0 + 10 * iel_ensm8]; | |
| acc_idof_ensm8_1 = acc_idof_ensm8_1 + _pt_data[iel_ensm8] * _pt_data_0[idof_ensm8_0 + 10 * idof_ensm8] * mass_inv_inp_0[idof_ensm8_0 + 10 * iel_ensm8]; | |
| } | |
| _pt_temp_0[idof_ensm8 + 10 * iel_ensm8] = acc_idof_ensm8_0; | |
| _pt_temp_2[idof_ensm8 + 10 * iel_ensm8] = acc_idof_ensm8_1; | |
| } | |
| for (int idof_ensm6 = 0; idof_ensm6 <= 9; ++idof_ensm6) | |
| for (int iel_ensm6 = 0; iel_ensm6 <= 31; ++iel_ensm6) | |
| { | |
| _pt_out_0_0[idof_ensm6 + 10 * iel_ensm6] = -1.0 * _pt_temp_2[idof_ensm6 + 10 * iel_ensm6]; | |
| _pt_out_1_0[idof_ensm6 + 10 * iel_ensm6] = -1.0 * _pt_temp_0[idof_ensm6 + 10 * iel_ensm6]; | |
| } | |
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs( double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_all, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, char const *__restrict__ from_el_present_3, double const *__restrict__ _pt_part_ph_id_4, double const *__restrict__ _pt_dist_id_1, int const *__restrict__ from_el_indices_5, int const *__restrict__ _pt_data_9, double const *__restrict__ normal_2_b_BTAG_PARTITION_1, int const *__restrict__ from_el_indices_6, int const *__restrict__ _pt_data_10, char const *__restrict__ from_el_present_4, double const *__restrict__ _pt_part_ph_id_5, double const *__restrict__ _pt_dist_id_2, int const *__restrict__ from_el_indices_7, int const *__restrict__ _pt_data_11, double const *__restrict__ normal_2_b_BTAG_PARTITION_2, int const *__restrict__ from_el_indices_8, int const *__restrict__ _pt_data_12, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_all, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double const *__restrict__ normal_1_b_BTAG_PARTITION_1, double const *__restrict__ normal_1_b_BTAG_PARTITION_2, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ cse_1, double *__restrict__ cse_2, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm5_0; | |
| double acc_idof_ensm5_1; | |
| double acc_iface_ensm4_idof_ensm4; | |
| double acc_iface_ensm4_idof_ensm4_0; | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| for (int iel_ensm0 = 0; iel_ensm0 <= 2; ++iel_ensm0) | |
| cse_2[idof_ensm0 + 4 * iel_ensm0] = 0.5 * (_pt_part_ph_id_5[idof_ensm0 + 4 * iel_ensm0] + _pt_dist_id_2[4 * from_el_indices_7[iel_ensm0] + _pt_data_11[idof_ensm0 + 4 * iel_ensm0]]); | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| for (int iel_ensm1 = 0; iel_ensm1 <= 4; ++iel_ensm1) | |
| cse[idof_ensm1 + 4 * iel_ensm1] = 0.5 * (_pt_part_ph_id_2[idof_ensm1 + 4 * iel_ensm1] + _pt_dist_id[4 * from_el_indices_1[iel_ensm1] + _pt_data_5[idof_ensm1 + 4 * iel_ensm1]]); | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| for (int iel_ensm2 = 0; iel_ensm2 <= 6; ++iel_ensm2) | |
| { | |
| cse_0[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_3[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id_0[4 * from_el_indices_3[iel_ensm2] + _pt_data_7[idof_ensm2 + 4 * iel_ensm2]]); | |
| cse_1[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_4[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id_1[4 * from_el_indices_5[iel_ensm2] + _pt_data_9[idof_ensm2 + 4 * iel_ensm2]]); | |
| } | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 3; ++idof_ensm3) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 31; ++iel_ensm3) | |
| for (int iface_ensm3 = 0; iface_ensm3 <= 2; ++iface_ensm3) | |
| { | |
| _pt_temp_1[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3] = (from_el_present[iel_ensm3 + 32 * iface_ensm3] ? normal_1_b_all[from_el_indices[iel_ensm3 + 32 * iface_ensm3]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm3 + 32 * iface_ensm3] + _pt_data_3[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm3 + 32 * iface_ensm3] + _pt_data_3[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]]) : 0.0) + (from_el_present_4[iel_ensm3 + 32 * iface_ensm3] ? cse_2[4 * from_el_indices_8[iel_ensm3 + 32 * iface_ensm3] + _pt_data_12[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION_2[from_el_indices_8[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_3[iel_ensm3 + 32 * iface_ensm3] ? cse_1[4 * from_el_indices_6[iel_ensm3 + 32 * iface_ensm3] + _pt_data_10[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_2[iel_ensm3 + 32 * iface_ensm3] ? cse_0[4 * from_el_indices_4[iel_ensm3 + 32 * iface_ensm3] + _pt_data_8[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_0[iel_ensm3 + 32 * iface_ensm3] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm3 + 32 * iface_ensm3] + _pt_data_4[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_face_restr_interior[from_el_indices_0[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_1[iel_ensm3 + 32 * iface_ensm3] ? cse[4 * from_el_indices_2[iel_ensm3 + 32 * iface_ensm3] + _pt_data_6[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION[from_el_indices_2[iel_ensm3 + 32 * iface_ensm3]] : 0.0); | |
| _pt_temp[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3] = (from_el_present[iel_ensm3 + 32 * iface_ensm3] ? normal_2_b_all[from_el_indices[iel_ensm3 + 32 * iface_ensm3]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm3 + 32 * iface_ensm3] + _pt_data_3[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm3 + 32 * iface_ensm3] + _pt_data_3[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]]) : 0.0) + (from_el_present_4[iel_ensm3 + 32 * iface_ensm3] ? cse_2[4 * from_el_indices_8[iel_ensm3 + 32 * iface_ensm3] + _pt_data_12[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION_2[from_el_indices_8[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_3[iel_ensm3 + 32 * iface_ensm3] ? cse_1[4 * from_el_indices_6[iel_ensm3 + 32 * iface_ensm3] + _pt_data_10[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_2[iel_ensm3 + 32 * iface_ensm3] ? cse_0[4 * from_el_indices_4[iel_ensm3 + 32 * iface_ensm3] + _pt_data_8[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_0[iel_ensm3 + 32 * iface_ensm3] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm3 + 32 * iface_ensm3] + _pt_data_4[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_face_restr_interior[from_el_indices_0[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_1[iel_ensm3 + 32 * iface_ensm3] ? cse[4 * from_el_indices_2[iel_ensm3 + 32 * iface_ensm3] + _pt_data_6[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION[from_el_indices_2[iel_ensm3 + 32 * iface_ensm3]] : 0.0); | |
| } | |
| for (int idof_ensm4_0 = 0; idof_ensm4_0 <= 9; ++idof_ensm4_0) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 31; ++iel_ensm4) | |
| { | |
| acc_iface_ensm4_idof_ensm4_0 = 0.0; | |
| acc_iface_ensm4_idof_ensm4 = 0.0; | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 3; ++idof_ensm4) | |
| for (int iface_ensm4 = 0; iface_ensm4 <= 2; ++iface_ensm4) | |
| { | |
| acc_iface_ensm4_idof_ensm4_0 = acc_iface_ensm4_idof_ensm4_0 + _pt_data_1[idof_ensm4 + 4 * iface_ensm4 + 12 * idof_ensm4_0] * _pt_data_2[32 * iface_ensm4 + iel_ensm4] * _pt_temp_1[idof_ensm4 + 128 * iface_ensm4 + 4 * iel_ensm4]; | |
| acc_iface_ensm4_idof_ensm4 = acc_iface_ensm4_idof_ensm4 + _pt_data_1[idof_ensm4 + 4 * iface_ensm4 + 12 * idof_ensm4_0] * _pt_data_2[32 * iface_ensm4 + iel_ensm4] * _pt_temp[idof_ensm4 + 128 * iface_ensm4 + 4 * iel_ensm4]; | |
| } | |
| face_mass_0[idof_ensm4_0 + 10 * iel_ensm4] = acc_iface_ensm4_idof_ensm4_0; | |
| face_mass[idof_ensm4_0 + 10 * iel_ensm4] = acc_iface_ensm4_idof_ensm4; | |
| } | |
| for (int idof_ensm6 = 0; idof_ensm6 <= 9; ++idof_ensm6) | |
| for (int iel_ensm6 = 0; iel_ensm6 <= 31; ++iel_ensm6) | |
| { | |
| mass_inv_inp_0[idof_ensm6 + 10 * iel_ensm6] = _pt_part_ph_id[idof_ensm6 + 10 * iel_ensm6] + -1.0 * face_mass_0[idof_ensm6 + 10 * iel_ensm6]; | |
| mass_inv_inp[idof_ensm6 + 10 * iel_ensm6] = _pt_part_ph_id[320 + idof_ensm6 + 10 * iel_ensm6] + -1.0 * face_mass[idof_ensm6 + 10 * iel_ensm6]; | |
| } | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 9; ++idof_ensm5) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| { | |
| acc_idof_ensm5_0 = 0.0; | |
| acc_idof_ensm5_1 = 0.0; | |
| for (int idof_ensm5_0 = 0; idof_ensm5_0 <= 9; ++idof_ensm5_0) | |
| { | |
| acc_idof_ensm5_0 = acc_idof_ensm5_0 + _pt_data[iel_ensm5] * _pt_data_0[idof_ensm5_0 + 10 * idof_ensm5] * mass_inv_inp[idof_ensm5_0 + 10 * iel_ensm5]; | |
| acc_idof_ensm5_1 = acc_idof_ensm5_1 + _pt_data[iel_ensm5] * _pt_data_0[idof_ensm5_0 + 10 * idof_ensm5] * mass_inv_inp_0[idof_ensm5_0 + 10 * iel_ensm5]; | |
| } | |
| _pt_temp_0[idof_ensm5 + 10 * iel_ensm5] = acc_idof_ensm5_0; | |
| _pt_temp_2[idof_ensm5 + 10 * iel_ensm5] = acc_idof_ensm5_1; | |
| } | |
| for (int idof_ensm7 = 0; idof_ensm7 <= 9; ++idof_ensm7) | |
| for (int iel_ensm7 = 0; iel_ensm7 <= 31; ++iel_ensm7) | |
| { | |
| _pt_out_0_0[idof_ensm7 + 10 * iel_ensm7] = -1.0 * _pt_temp_2[idof_ensm7 + 10 * iel_ensm7]; | |
| _pt_out_1_0[idof_ensm7 + 10 * iel_ensm7] = -1.0 * _pt_temp_0[idof_ensm7 + 10 * iel_ensm7]; | |
| } | |
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs(double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_1, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, double const *__restrict__ normal_2_b_BTAG_PARTITION_1, int const *__restrict__ from_el_indices_5, int const *__restrict__ _pt_data_9, char const *__restrict__ from_el_present_3, double const *__restrict__ _pt_part_ph_id_4, double const *__restrict__ _pt_dist_id_2, int const *__restrict__ from_el_indices_6, int const *__restrict__ _pt_data_10, double const *__restrict__ normal_2_b_BTAG_PARTITION_2, int const *__restrict__ from_el_indices_7, int const *__restrict__ _pt_data_11, char const *__restrict__ from_el_present_4, double const *__restrict__ _pt_part_ph_id_5, double const *__restrict__ _pt_dist_id_3, int const *__restrict__ from_el_indices_8, int const *__restrict__ _pt_data_12, double const *__restrict__ normal_2_b_BTAG_PARTITION_3, int const *__restrict__ from_el_indices_9, int const *__restrict__ _pt_data_13, char const *__restrict__ from_el_present_5, double const *__restrict__ _pt_part_ph_id_6, double const *__restrict__ _pt_dist_id_4, int const *__restrict__ from_el_indices_10, int const *__restrict__ _pt_data_14, double const *__restrict__ normal_2_b_BTAG_PARTITION_4, int const *__restrict__ from_el_indices_11, int const *__restrict__ _pt_data_15, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double const *__restrict__ normal_1_b_BTAG_PARTITION_1, double const *__restrict__ normal_1_b_BTAG_PARTITION_2, double const *__restrict__ normal_1_b_BTAG_PARTITION_3, double const *__restrict__ normal_1_b_BTAG_PARTITION_4, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ cse_1, double *__restrict__ cse_2, double *__restrict__ cse_3, double *__restrict__ cse_4, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm8_0; | |
| double acc_idof_ensm8_1; | |
| double acc_iface_ensm5_idof_ensm5; | |
| double acc_iface_ensm5_idof_ensm5_0; | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| { | |
| int const iel_ensm0 = 0; | |
| cse_0[idof_ensm0] = 0.5 * (_pt_part_ph_id_2[idof_ensm0] + _pt_dist_id_0[4 * from_el_indices_2[0] + _pt_data_6[idof_ensm0]]); | |
| cse_1[idof_ensm0] = 0.5 * (_pt_part_ph_id_3[idof_ensm0] + _pt_dist_id_1[4 * from_el_indices_4[0] + _pt_data_8[idof_ensm0]]); | |
| } | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| for (int iel_ensm1 = 0; iel_ensm1 <= 1; ++iel_ensm1) | |
| { | |
| cse_4[idof_ensm1 + 4 * iel_ensm1] = 0.5 * (_pt_part_ph_id_6[idof_ensm1 + 4 * iel_ensm1] + _pt_dist_id_4[4 * from_el_indices_10[iel_ensm1] + _pt_data_14[idof_ensm1 + 4 * iel_ensm1]]); | |
| cse[idof_ensm1 + 4 * iel_ensm1] = 0.5 * (_pt_part_ph_id_1[idof_ensm1 + 4 * iel_ensm1] + _pt_dist_id[4 * from_el_indices_0[iel_ensm1] + _pt_data_4[idof_ensm1 + 4 * iel_ensm1]]); | |
| } | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 3; ++idof_ensm3) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 7; ++iel_ensm3) | |
| cse_2[idof_ensm3 + 4 * iel_ensm3] = 0.5 * (_pt_part_ph_id_4[idof_ensm3 + 4 * iel_ensm3] + _pt_dist_id_2[4 * from_el_indices_6[iel_ensm3] + _pt_data_10[idof_ensm3 + 4 * iel_ensm3]]); | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 3; ++idof_ensm4) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 3; ++iel_ensm4) | |
| cse_3[idof_ensm4 + 4 * iel_ensm4] = 0.5 * (_pt_part_ph_id_5[idof_ensm4 + 4 * iel_ensm4] + _pt_dist_id_3[4 * from_el_indices_8[iel_ensm4] + _pt_data_12[idof_ensm4 + 4 * iel_ensm4]]); | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| for (int iel_ensm2 = 0; iel_ensm2 <= 31; ++iel_ensm2) | |
| for (int iface_ensm2 = 0; iface_ensm2 <= 2; ++iface_ensm2) | |
| { | |
| _pt_temp_1[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2] = (from_el_present_5[iel_ensm2 + 32 * iface_ensm2] ? cse_4[4 * from_el_indices_11[iel_ensm2 + 32 * iface_ensm2] + _pt_data_15[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2]] * normal_1_b_BTAG_PARTITION_4[from_el_indices_11[iel_ensm2 + 32 * iface_ensm2]] : 0.0) + (from_el_present_4[iel_ensm2 + 32 * iface_ensm2] ? cse_3[4 * from_el_indices_9[iel_ensm2 + 32 * iface_ensm2] + _pt_data_13[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2]] * normal_1_b_BTAG_PARTITION_3[from_el_indices_9[iel_ensm2 + 32 * iface_ensm2]] : 0.0) + (from_el_present_3[iel_ensm2 + 32 * iface_ensm2] ? cse_2[4 * from_el_indices_7[iel_ensm2 + 32 * iface_ensm2] + _pt_data_11[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2]] * normal_1_b_BTAG_PARTITION_2[from_el_indices_7[iel_ensm2 + 32 * iface_ensm2]] : 0.0) + (from_el_present_2[iel_ensm2 + 32 * iface_ensm2] ? cse_1[4 * from_el_indices_5[iel_ensm2 + 32 * iface_ensm2] + _pt_data_9[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2]] * normal_1_b_BTAG_PARTITION_1[from_el_indices_5[iel_ensm2 + 32 * iface_ensm2]] : 0.0) + (from_el_present_1[iel_ensm2 + 32 * iface_ensm2] ? cse_0[4 * from_el_indices_3[iel_ensm2 + 32 * iface_ensm2] + _pt_data_7[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_3[iel_ensm2 + 32 * iface_ensm2]] : 0.0) + (from_el_present[iel_ensm2 + 32 * iface_ensm2] ? _pt_part_ph_id_0[4 * from_el_indices[iel_ensm2 + 32 * iface_ensm2] + _pt_data_3[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2]] * normal_1_b_face_restr_interior[from_el_indices[iel_ensm2 + 32 * iface_ensm2]] : 0.0) + (from_el_present_0[iel_ensm2 + 32 * iface_ensm2] ? cse[4 * from_el_indices_1[iel_ensm2 + 32 * iface_ensm2] + _pt_data_5[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2]] * normal_1_b_BTAG_PARTITION[from_el_indices_1[iel_ensm2 + 32 * iface_ensm2]] : 0.0); | |
| _pt_temp[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2] = (from_el_present_5[iel_ensm2 + 32 * iface_ensm2] ? cse_4[4 * from_el_indices_11[iel_ensm2 + 32 * iface_ensm2] + _pt_data_15[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2]] * normal_2_b_BTAG_PARTITION_4[from_el_indices_11[iel_ensm2 + 32 * iface_ensm2]] : 0.0) + (from_el_present_4[iel_ensm2 + 32 * iface_ensm2] ? cse_3[4 * from_el_indices_9[iel_ensm2 + 32 * iface_ensm2] + _pt_data_13[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2]] * normal_2_b_BTAG_PARTITION_3[from_el_indices_9[iel_ensm2 + 32 * iface_ensm2]] : 0.0) + (from_el_present_3[iel_ensm2 + 32 * iface_ensm2] ? cse_2[4 * from_el_indices_7[iel_ensm2 + 32 * iface_ensm2] + _pt_data_11[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2]] * normal_2_b_BTAG_PARTITION_2[from_el_indices_7[iel_ensm2 + 32 * iface_ensm2]] : 0.0) + (from_el_present_2[iel_ensm2 + 32 * iface_ensm2] ? cse_1[4 * from_el_indices_5[iel_ensm2 + 32 * iface_ensm2] + _pt_data_9[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2]] * normal_2_b_BTAG_PARTITION_1[from_el_indices_5[iel_ensm2 + 32 * iface_ensm2]] : 0.0) + (from_el_present_1[iel_ensm2 + 32 * iface_ensm2] ? cse_0[4 * from_el_indices_3[iel_ensm2 + 32 * iface_ensm2] + _pt_data_7[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_3[iel_ensm2 + 32 * iface_ensm2]] : 0.0) + (from_el_present[iel_ensm2 + 32 * iface_ensm2] ? _pt_part_ph_id_0[4 * from_el_indices[iel_ensm2 + 32 * iface_ensm2] + _pt_data_3[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2]] * normal_2_b_face_restr_interior[from_el_indices[iel_ensm2 + 32 * iface_ensm2]] : 0.0) + (from_el_present_0[iel_ensm2 + 32 * iface_ensm2] ? cse[4 * from_el_indices_1[iel_ensm2 + 32 * iface_ensm2] + _pt_data_5[idof_ensm2 + 4 * iel_ensm2 + 128 * iface_ensm2]] * normal_2_b_BTAG_PARTITION[from_el_indices_1[iel_ensm2 + 32 * iface_ensm2]] : 0.0); | |
| } | |
| for (int idof_ensm5_0 = 0; idof_ensm5_0 <= 9; ++idof_ensm5_0) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| { | |
| acc_iface_ensm5_idof_ensm5_0 = 0.0; | |
| acc_iface_ensm5_idof_ensm5 = 0.0; | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 3; ++idof_ensm5) | |
| for (int iface_ensm5 = 0; iface_ensm5 <= 2; ++iface_ensm5) | |
| { | |
| acc_iface_ensm5_idof_ensm5_0 = acc_iface_ensm5_idof_ensm5_0 + _pt_data_1[idof_ensm5 + 4 * iface_ensm5 + 12 * idof_ensm5_0] * _pt_data_2[32 * iface_ensm5 + iel_ensm5] * _pt_temp_1[idof_ensm5 + 128 * iface_ensm5 + 4 * iel_ensm5]; | |
| acc_iface_ensm5_idof_ensm5 = acc_iface_ensm5_idof_ensm5 + _pt_data_1[idof_ensm5 + 4 * iface_ensm5 + 12 * idof_ensm5_0] * _pt_data_2[32 * iface_ensm5 + iel_ensm5] * _pt_temp[idof_ensm5 + 128 * iface_ensm5 + 4 * iel_ensm5]; | |
| } | |
| face_mass_0[idof_ensm5_0 + 10 * iel_ensm5] = acc_iface_ensm5_idof_ensm5_0; | |
| face_mass[idof_ensm5_0 + 10 * iel_ensm5] = acc_iface_ensm5_idof_ensm5; | |
| } | |
| for (int idof_ensm7 = 0; idof_ensm7 <= 9; ++idof_ensm7) | |
| for (int iel_ensm7 = 0; iel_ensm7 <= 31; ++iel_ensm7) | |
| { | |
| mass_inv_inp_0[idof_ensm7 + 10 * iel_ensm7] = _pt_part_ph_id[idof_ensm7 + 10 * iel_ensm7] + -1.0 * face_mass_0[idof_ensm7 + 10 * iel_ensm7]; | |
| mass_inv_inp[idof_ensm7 + 10 * iel_ensm7] = _pt_part_ph_id[320 + idof_ensm7 + 10 * iel_ensm7] + -1.0 * face_mass[idof_ensm7 + 10 * iel_ensm7]; | |
| } | |
| for (int idof_ensm8 = 0; idof_ensm8 <= 9; ++idof_ensm8) | |
| for (int iel_ensm8 = 0; iel_ensm8 <= 31; ++iel_ensm8) | |
| { | |
| acc_idof_ensm8_0 = 0.0; | |
| acc_idof_ensm8_1 = 0.0; | |
| for (int idof_ensm8_0 = 0; idof_ensm8_0 <= 9; ++idof_ensm8_0) | |
| { | |
| acc_idof_ensm8_0 = acc_idof_ensm8_0 + _pt_data[iel_ensm8] * _pt_data_0[idof_ensm8_0 + 10 * idof_ensm8] * mass_inv_inp[idof_ensm8_0 + 10 * iel_ensm8]; | |
| acc_idof_ensm8_1 = acc_idof_ensm8_1 + _pt_data[iel_ensm8] * _pt_data_0[idof_ensm8_0 + 10 * idof_ensm8] * mass_inv_inp_0[idof_ensm8_0 + 10 * iel_ensm8]; | |
| } | |
| _pt_temp_0[idof_ensm8 + 10 * iel_ensm8] = acc_idof_ensm8_0; | |
| _pt_temp_2[idof_ensm8 + 10 * iel_ensm8] = acc_idof_ensm8_1; | |
| } | |
| for (int idof_ensm6 = 0; idof_ensm6 <= 9; ++idof_ensm6) | |
| for (int iel_ensm6 = 0; iel_ensm6 <= 31; ++iel_ensm6) | |
| { | |
| _pt_out_0_0[idof_ensm6 + 10 * iel_ensm6] = -1.0 * _pt_temp_2[idof_ensm6 + 10 * iel_ensm6]; | |
| _pt_out_1_0[idof_ensm6 + 10 * iel_ensm6] = -1.0 * _pt_temp_0[idof_ensm6 + 10 * iel_ensm6]; | |
| } | |
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs(double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_1, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, double const *__restrict__ normal_2_b_BTAG_PARTITION_1, int const *__restrict__ from_el_indices_5, int const *__restrict__ _pt_data_9, char const *__restrict__ from_el_present_3, double const *__restrict__ _pt_part_ph_id_4, double const *__restrict__ _pt_dist_id_2, int const *__restrict__ from_el_indices_6, int const *__restrict__ _pt_data_10, double const *__restrict__ normal_2_b_BTAG_PARTITION_2, int const *__restrict__ from_el_indices_7, int const *__restrict__ _pt_data_11, char const *__restrict__ from_el_present_4, double const *__restrict__ _pt_part_ph_id_5, double const *__restrict__ _pt_dist_id_3, int const *__restrict__ from_el_indices_8, int const *__restrict__ _pt_data_12, double const *__restrict__ normal_2_b_BTAG_PARTITION_3, int const *__restrict__ from_el_indices_9, int const *__restrict__ _pt_data_13, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double const *__restrict__ normal_1_b_BTAG_PARTITION_1, double const *__restrict__ normal_1_b_BTAG_PARTITION_2, double const *__restrict__ normal_1_b_BTAG_PARTITION_3, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ cse_1, double *__restrict__ cse_2, double *__restrict__ cse_3, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm7_0; | |
| double acc_idof_ensm7_1; | |
| double acc_iface_ensm6_idof_ensm6; | |
| double acc_iface_ensm6_idof_ensm6_0; | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| { | |
| int const iel_ensm0 = 0; | |
| cse_0[idof_ensm0] = 0.5 * (_pt_part_ph_id_2[idof_ensm0] + _pt_dist_id_0[4 * from_el_indices_2[0] + _pt_data_6[idof_ensm0]]); | |
| } | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| for (int iel_ensm1 = 0; iel_ensm1 <= 3; ++iel_ensm1) | |
| cse_1[idof_ensm1 + 4 * iel_ensm1] = 0.5 * (_pt_part_ph_id_3[idof_ensm1 + 4 * iel_ensm1] + _pt_dist_id_1[4 * from_el_indices_4[iel_ensm1] + _pt_data_8[idof_ensm1 + 4 * iel_ensm1]]); | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| for (int iel_ensm2 = 0; iel_ensm2 <= 2; ++iel_ensm2) | |
| cse_2[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_4[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id_2[4 * from_el_indices_6[iel_ensm2] + _pt_data_10[idof_ensm2 + 4 * iel_ensm2]]); | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 3; ++idof_ensm3) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 7; ++iel_ensm3) | |
| cse_3[idof_ensm3 + 4 * iel_ensm3] = 0.5 * (_pt_part_ph_id_5[idof_ensm3 + 4 * iel_ensm3] + _pt_dist_id_3[4 * from_el_indices_8[iel_ensm3] + _pt_data_12[idof_ensm3 + 4 * iel_ensm3]]); | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 3; ++idof_ensm4) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 5; ++iel_ensm4) | |
| cse[idof_ensm4 + 4 * iel_ensm4] = 0.5 * (_pt_part_ph_id_1[idof_ensm4 + 4 * iel_ensm4] + _pt_dist_id[4 * from_el_indices_0[iel_ensm4] + _pt_data_4[idof_ensm4 + 4 * iel_ensm4]]); | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 3; ++idof_ensm5) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| for (int iface_ensm5 = 0; iface_ensm5 <= 2; ++iface_ensm5) | |
| { | |
| _pt_temp_1[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5] = | |
| (from_el_present_4[iel_ensm5 + 32 * iface_ensm5] ? | |
| cse_3[4 * from_el_indices_9[iel_ensm5 + 32 * iface_ensm5] + _pt_data_13[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_1_b_BTAG_PARTITION_3[from_el_indices_9[iel_ensm5 + 32 * iface_ensm5]] | |
| : 0.0) | |
| + (from_el_present_3[iel_ensm5 + 32 * iface_ensm5] ? | |
| cse_2[4 * from_el_indices_7[iel_ensm5 + 32 * iface_ensm5] + _pt_data_11[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_1_b_BTAG_PARTITION_2[from_el_indices_7[iel_ensm5 + 32 * iface_ensm5]] | |
| : 0.0) | |
| + (from_el_present_2[iel_ensm5 + 32 * iface_ensm5] ? | |
| cse_1[4 * from_el_indices_5[iel_ensm5 + 32 * iface_ensm5] + _pt_data_9[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_1_b_BTAG_PARTITION_1[from_el_indices_5[iel_ensm5 + 32 * iface_ensm5]] | |
| : 0.0) | |
| + (from_el_present_1[iel_ensm5 + 32 * iface_ensm5] ? | |
| cse_0[4 * from_el_indices_3[iel_ensm5 + 32 * iface_ensm5] + _pt_data_7[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_3[iel_ensm5 + 32 * iface_ensm5]] | |
| : 0.0) | |
| + (from_el_present[iel_ensm5 + 32 * iface_ensm5] ? | |
| _pt_part_ph_id_0[4 * from_el_indices[iel_ensm5 + 32 * iface_ensm5] + _pt_data_3[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_1_b_face_restr_interior[from_el_indices[iel_ensm5 + 32 * iface_ensm5]] | |
| : 0.0) | |
| + (from_el_present_0[iel_ensm5 + 32 * iface_ensm5] ? | |
| cse[4 * from_el_indices_1[iel_ensm5 + 32 * iface_ensm5] + _pt_data_5[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_1_b_BTAG_PARTITION[from_el_indices_1[iel_ensm5 + 32 * iface_ensm5]] | |
| : 0.0); | |
| _pt_temp[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5] = (from_el_present_4[iel_ensm5 + 32 * iface_ensm5] ? cse_3[4 * from_el_indices_9[iel_ensm5 + 32 * iface_ensm5] + _pt_data_13[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_2_b_BTAG_PARTITION_3[from_el_indices_9[iel_ensm5 + 32 * iface_ensm5]] : 0.0) + (from_el_present_3[iel_ensm5 + 32 * iface_ensm5] ? cse_2[4 * from_el_indices_7[iel_ensm5 + 32 * iface_ensm5] + _pt_data_11[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_2_b_BTAG_PARTITION_2[from_el_indices_7[iel_ensm5 + 32 * iface_ensm5]] : 0.0) + (from_el_present_2[iel_ensm5 + 32 * iface_ensm5] ? cse_1[4 * from_el_indices_5[iel_ensm5 + 32 * iface_ensm5] + _pt_data_9[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_2_b_BTAG_PARTITION_1[from_el_indices_5[iel_ensm5 + 32 * iface_ensm5]] : 0.0) + (from_el_present_1[iel_ensm5 + 32 * iface_ensm5] ? cse_0[4 * from_el_indices_3[iel_ensm5 + 32 * iface_ensm5] + _pt_data_7[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_3[iel_ensm5 + 32 * iface_ensm5]] : 0.0) + (from_el_present[iel_ensm5 + 32 * iface_ensm5] ? _pt_part_ph_id_0[4 * from_el_indices[iel_ensm5 + 32 * iface_ensm5] + _pt_data_3[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_2_b_face_restr_interior[from_el_indices[iel_ensm5 + 32 * iface_ensm5]] : 0.0) + (from_el_present_0[iel_ensm5 + 32 * iface_ensm5] ? cse[4 * from_el_indices_1[iel_ensm5 + 32 * iface_ensm5] + _pt_data_5[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_2_b_BTAG_PARTITION[from_el_indices_1[iel_ensm5 + 32 * iface_ensm5]] : 0.0); | |
| } | |
| for (int idof_ensm6_0 = 0; idof_ensm6_0 <= 9; ++idof_ensm6_0) | |
| for (int iel_ensm6 = 0; iel_ensm6 <= 31; ++iel_ensm6) | |
| { | |
| acc_iface_ensm6_idof_ensm6_0 = 0.0; | |
| acc_iface_ensm6_idof_ensm6 = 0.0; | |
| for (int idof_ensm6 = 0; idof_ensm6 <= 3; ++idof_ensm6) | |
| for (int iface_ensm6 = 0; iface_ensm6 <= 2; ++iface_ensm6) | |
| { | |
| acc_iface_ensm6_idof_ensm6_0 = acc_iface_ensm6_idof_ensm6_0 + _pt_data_1[idof_ensm6 + 4 * iface_ensm6 + 12 * idof_ensm6_0] * _pt_data_2[32 * iface_ensm6 + iel_ensm6] * _pt_temp_1[idof_ensm6 + 128 * iface_ensm6 + 4 * iel_ensm6]; | |
| acc_iface_ensm6_idof_ensm6 = acc_iface_ensm6_idof_ensm6 + _pt_data_1[idof_ensm6 + 4 * iface_ensm6 + 12 * idof_ensm6_0] * _pt_data_2[32 * iface_ensm6 + iel_ensm6] * _pt_temp[idof_ensm6 + 128 * iface_ensm6 + 4 * iel_ensm6]; | |
| } | |
| face_mass_0[idof_ensm6_0 + 10 * iel_ensm6] = acc_iface_ensm6_idof_ensm6_0; | |
| face_mass[idof_ensm6_0 + 10 * iel_ensm6] = acc_iface_ensm6_idof_ensm6; | |
| } | |
| for (int idof_ensm8 = 0; idof_ensm8 <= 9; ++idof_ensm8) | |
| for (int iel_ensm8 = 0; iel_ensm8 <= 31; ++iel_ensm8) | |
| { | |
| mass_inv_inp_0[idof_ensm8 + 10 * iel_ensm8] = _pt_part_ph_id[idof_ensm8 + 10 * iel_ensm8] + -1.0 * face_mass_0[idof_ensm8 + 10 * iel_ensm8]; | |
| mass_inv_inp[idof_ensm8 + 10 * iel_ensm8] = _pt_part_ph_id[320 + idof_ensm8 + 10 * iel_ensm8] + -1.0 * face_mass[idof_ensm8 + 10 * iel_ensm8]; | |
| } | |
| for (int idof_ensm7 = 0; idof_ensm7 <= 9; ++idof_ensm7) | |
| for (int iel_ensm7 = 0; iel_ensm7 <= 31; ++iel_ensm7) | |
| { | |
| acc_idof_ensm7_0 = 0.0; | |
| acc_idof_ensm7_1 = 0.0; | |
| for (int idof_ensm7_0 = 0; idof_ensm7_0 <= 9; ++idof_ensm7_0) | |
| { | |
| acc_idof_ensm7_0 = acc_idof_ensm7_0 + _pt_data[iel_ensm7] * _pt_data_0[idof_ensm7_0 + 10 * idof_ensm7] * mass_inv_inp[idof_ensm7_0 + 10 * iel_ensm7]; | |
| acc_idof_ensm7_1 = acc_idof_ensm7_1 + _pt_data[iel_ensm7] * _pt_data_0[idof_ensm7_0 + 10 * idof_ensm7] * mass_inv_inp_0[idof_ensm7_0 + 10 * iel_ensm7]; | |
| } | |
| _pt_temp_0[idof_ensm7 + 10 * iel_ensm7] = acc_idof_ensm7_0; | |
| _pt_temp_2[idof_ensm7 + 10 * iel_ensm7] = acc_idof_ensm7_1; | |
| } | |
| for (int idof_ensm9 = 0; idof_ensm9 <= 9; ++idof_ensm9) | |
| for (int iel_ensm9 = 0; iel_ensm9 <= 31; ++iel_ensm9) | |
| { | |
| _pt_out_0_0[idof_ensm9 + 10 * iel_ensm9] = -1.0 * _pt_temp_2[idof_ensm9 + 10 * iel_ensm9]; | |
| _pt_out_1_0[idof_ensm9 + 10 * iel_ensm9] = -1.0 * _pt_temp_0[idof_ensm9 + 10 * iel_ensm9]; | |
| } | |
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs(double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_all, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, char const *__restrict__ from_el_present_3, double const *__restrict__ _pt_part_ph_id_4, double const *__restrict__ _pt_dist_id_1, int const *__restrict__ from_el_indices_5, int const *__restrict__ _pt_data_9, double const *__restrict__ normal_2_b_BTAG_PARTITION_1, int const *__restrict__ from_el_indices_6, int const *__restrict__ _pt_data_10, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_all, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double const *__restrict__ normal_1_b_BTAG_PARTITION_1, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ cse_1, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm6_0; | |
| double acc_idof_ensm6_1; | |
| double acc_iface_ensm3_idof_ensm3; | |
| double acc_iface_ensm3_idof_ensm3_0; | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| for (int iel_ensm1 = 0; iel_ensm1 <= 2; ++iel_ensm1) | |
| cse[idof_ensm1 + 4 * iel_ensm1] = 0.5 * (_pt_part_ph_id_2[idof_ensm1 + 4 * iel_ensm1] + _pt_dist_id[4 * from_el_indices_1[iel_ensm1] + _pt_data_5[idof_ensm1 + 4 * iel_ensm1]]); | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| for (int iel_ensm2 = 0; iel_ensm2 <= 3; ++iel_ensm2) | |
| { | |
| cse_0[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_3[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id_0[4 * from_el_indices_3[iel_ensm2] + _pt_data_7[idof_ensm2 + 4 * iel_ensm2]]); | |
| cse_1[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_4[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id_1[4 * from_el_indices_5[iel_ensm2] + _pt_data_9[idof_ensm2 + 4 * iel_ensm2]]); | |
| } | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| for (int iel_ensm0 = 0; iel_ensm0 <= 31; ++iel_ensm0) | |
| for (int iface_ensm0 = 0; iface_ensm0 <= 2; ++iface_ensm0) | |
| { | |
| _pt_temp_1[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0] = (from_el_present[iel_ensm0 + 32 * iface_ensm0] ? normal_1_b_all[from_el_indices[iel_ensm0 + 32 * iface_ensm0]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]]) : 0.0) + (from_el_present_3[iel_ensm0 + 32 * iface_ensm0] ? cse_1[4 * from_el_indices_6[iel_ensm0 + 32 * iface_ensm0] + _pt_data_10[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_2[iel_ensm0 + 32 * iface_ensm0] ? cse_0[4 * from_el_indices_4[iel_ensm0 + 32 * iface_ensm0] + _pt_data_8[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_0[iel_ensm0 + 32 * iface_ensm0] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm0 + 32 * iface_ensm0] + _pt_data_4[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_face_restr_interior[from_el_indices_0[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_1[iel_ensm0 + 32 * iface_ensm0] ? cse[4 * from_el_indices_2[iel_ensm0 + 32 * iface_ensm0] + _pt_data_6[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION[from_el_indices_2[iel_ensm0 + 32 * iface_ensm0]] : 0.0); | |
| _pt_temp[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0] = (from_el_present[iel_ensm0 + 32 * iface_ensm0] ? normal_2_b_all[from_el_indices[iel_ensm0 + 32 * iface_ensm0]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]]) : 0.0) + (from_el_present_3[iel_ensm0 + 32 * iface_ensm0] ? cse_1[4 * from_el_indices_6[iel_ensm0 + 32 * iface_ensm0] + _pt_data_10[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_2[iel_ensm0 + 32 * iface_ensm0] ? cse_0[4 * from_el_indices_4[iel_ensm0 + 32 * iface_ensm0] + _pt_data_8[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_0[iel_ensm0 + 32 * iface_ensm0] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm0 + 32 * iface_ensm0] + _pt_data_4[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_face_restr_interior[from_el_indices_0[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_1[iel_ensm0 + 32 * iface_ensm0] ? cse[4 * from_el_indices_2[iel_ensm0 + 32 * iface_ensm0] + _pt_data_6[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION[from_el_indices_2[iel_ensm0 + 32 * iface_ensm0]] : 0.0); | |
| } | |
| for (int idof_ensm3_0 = 0; idof_ensm3_0 <= 9; ++idof_ensm3_0) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 31; ++iel_ensm3) | |
| { | |
| acc_iface_ensm3_idof_ensm3_0 = 0.0; | |
| acc_iface_ensm3_idof_ensm3 = 0.0; | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 3; ++idof_ensm3) | |
| for (int iface_ensm3 = 0; iface_ensm3 <= 2; ++iface_ensm3) | |
| { | |
| acc_iface_ensm3_idof_ensm3_0 = acc_iface_ensm3_idof_ensm3_0 + _pt_data_1[idof_ensm3 + 4 * iface_ensm3 + 12 * idof_ensm3_0] * _pt_data_2[32 * iface_ensm3 + iel_ensm3] * _pt_temp_1[idof_ensm3 + 128 * iface_ensm3 + 4 * iel_ensm3]; | |
| acc_iface_ensm3_idof_ensm3 = acc_iface_ensm3_idof_ensm3 + _pt_data_1[idof_ensm3 + 4 * iface_ensm3 + 12 * idof_ensm3_0] * _pt_data_2[32 * iface_ensm3 + iel_ensm3] * _pt_temp[idof_ensm3 + 128 * iface_ensm3 + 4 * iel_ensm3]; | |
| } | |
| face_mass_0[idof_ensm3_0 + 10 * iel_ensm3] = acc_iface_ensm3_idof_ensm3_0; | |
| face_mass[idof_ensm3_0 + 10 * iel_ensm3] = acc_iface_ensm3_idof_ensm3; | |
| } | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 9; ++idof_ensm4) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 31; ++iel_ensm4) | |
| { | |
| mass_inv_inp_0[idof_ensm4 + 10 * iel_ensm4] = _pt_part_ph_id[idof_ensm4 + 10 * iel_ensm4] + -1.0 * face_mass_0[idof_ensm4 + 10 * iel_ensm4]; | |
| mass_inv_inp[idof_ensm4 + 10 * iel_ensm4] = _pt_part_ph_id[320 + idof_ensm4 + 10 * iel_ensm4] + -1.0 * face_mass[idof_ensm4 + 10 * iel_ensm4]; | |
| } | |
| for (int idof_ensm6 = 0; idof_ensm6 <= 9; ++idof_ensm6) | |
| for (int iel_ensm6 = 0; iel_ensm6 <= 31; ++iel_ensm6) | |
| { | |
| acc_idof_ensm6_0 = 0.0; | |
| acc_idof_ensm6_1 = 0.0; | |
| for (int idof_ensm6_0 = 0; idof_ensm6_0 <= 9; ++idof_ensm6_0) | |
| { | |
| acc_idof_ensm6_0 = acc_idof_ensm6_0 + _pt_data[iel_ensm6] * _pt_data_0[idof_ensm6_0 + 10 * idof_ensm6] * mass_inv_inp[idof_ensm6_0 + 10 * iel_ensm6]; | |
| acc_idof_ensm6_1 = acc_idof_ensm6_1 + _pt_data[iel_ensm6] * _pt_data_0[idof_ensm6_0 + 10 * idof_ensm6] * mass_inv_inp_0[idof_ensm6_0 + 10 * iel_ensm6]; | |
| } | |
| _pt_temp_0[idof_ensm6 + 10 * iel_ensm6] = acc_idof_ensm6_0; | |
| _pt_temp_2[idof_ensm6 + 10 * iel_ensm6] = acc_idof_ensm6_1; | |
| } | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 9; ++idof_ensm5) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| { | |
| _pt_out_0_0[idof_ensm5 + 10 * iel_ensm5] = -1.0 * _pt_temp_2[idof_ensm5 + 10 * iel_ensm5]; | |
| _pt_out_1_0[idof_ensm5 + 10 * iel_ensm5] = -1.0 * _pt_temp_0[idof_ensm5 + 10 * iel_ensm5]; | |
| } | |
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs(double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_all, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, char const *__restrict__ from_el_present_3, double const *__restrict__ _pt_part_ph_id_4, double const *__restrict__ _pt_dist_id_1, int const *__restrict__ from_el_indices_5, int const *__restrict__ _pt_data_9, double const *__restrict__ normal_2_b_BTAG_PARTITION_1, int const *__restrict__ from_el_indices_6, int const *__restrict__ _pt_data_10, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_all, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double const *__restrict__ normal_1_b_BTAG_PARTITION_1, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ cse_1, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm6_0; | |
| double acc_idof_ensm6_1; | |
| double acc_iface_ensm3_idof_ensm3; | |
| double acc_iface_ensm3_idof_ensm3_0; | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| for (int iel_ensm0 = 0; iel_ensm0 <= 1; ++iel_ensm0) | |
| cse_0[idof_ensm0 + 4 * iel_ensm0] = 0.5 * (_pt_part_ph_id_3[idof_ensm0 + 4 * iel_ensm0] + _pt_dist_id_0[4 * from_el_indices_3[iel_ensm0] + _pt_data_7[idof_ensm0 + 4 * iel_ensm0]]); | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| for (int iel_ensm2 = 0; iel_ensm2 <= 3; ++iel_ensm2) | |
| { | |
| cse_1[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_4[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id_1[4 * from_el_indices_5[iel_ensm2] + _pt_data_9[idof_ensm2 + 4 * iel_ensm2]]); | |
| cse[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_2[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id[4 * from_el_indices_1[iel_ensm2] + _pt_data_5[idof_ensm2 + 4 * iel_ensm2]]); | |
| } | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| for (int iel_ensm1 = 0; iel_ensm1 <= 31; ++iel_ensm1) | |
| for (int iface_ensm1 = 0; iface_ensm1 <= 2; ++iface_ensm1) | |
| { | |
| _pt_temp_1[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1] = (from_el_present[iel_ensm1 + 32 * iface_ensm1] ? normal_1_b_all[from_el_indices[iel_ensm1 + 32 * iface_ensm1]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm1 + 32 * iface_ensm1] + _pt_data_3[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm1 + 32 * iface_ensm1] + _pt_data_3[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]]) : 0.0) + (from_el_present_3[iel_ensm1 + 32 * iface_ensm1] ? cse_1[4 * from_el_indices_6[iel_ensm1 + 32 * iface_ensm1] + _pt_data_10[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_1_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm1 + 32 * iface_ensm1]] : 0.0) + (from_el_present_2[iel_ensm1 + 32 * iface_ensm1] ? cse_0[4 * from_el_indices_4[iel_ensm1 + 32 * iface_ensm1] + _pt_data_8[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm1 + 32 * iface_ensm1]] : 0.0) + (from_el_present_0[iel_ensm1 + 32 * iface_ensm1] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm1 + 32 * iface_ensm1] + _pt_data_4[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_1_b_face_restr_interior[from_el_indices_0[iel_ensm1 + 32 * iface_ensm1]] : 0.0) + (from_el_present_1[iel_ensm1 + 32 * iface_ensm1] ? cse[4 * from_el_indices_2[iel_ensm1 + 32 * iface_ensm1] + _pt_data_6[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_1_b_BTAG_PARTITION[from_el_indices_2[iel_ensm1 + 32 * iface_ensm1]] : 0.0); | |
| _pt_temp[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1] = (from_el_present[iel_ensm1 + 32 * iface_ensm1] ? normal_2_b_all[from_el_indices[iel_ensm1 + 32 * iface_ensm1]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm1 + 32 * iface_ensm1] + _pt_data_3[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm1 + 32 * iface_ensm1] + _pt_data_3[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]]) : 0.0) + (from_el_present_3[iel_ensm1 + 32 * iface_ensm1] ? cse_1[4 * from_el_indices_6[iel_ensm1 + 32 * iface_ensm1] + _pt_data_10[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_2_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm1 + 32 * iface_ensm1]] : 0.0) + (from_el_present_2[iel_ensm1 + 32 * iface_ensm1] ? cse_0[4 * from_el_indices_4[iel_ensm1 + 32 * iface_ensm1] + _pt_data_8[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm1 + 32 * iface_ensm1]] : 0.0) + (from_el_present_0[iel_ensm1 + 32 * iface_ensm1] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm1 + 32 * iface_ensm1] + _pt_data_4[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_2_b_face_restr_interior[from_el_indices_0[iel_ensm1 + 32 * iface_ensm1]] : 0.0) + (from_el_present_1[iel_ensm1 + 32 * iface_ensm1] ? cse[4 * from_el_indices_2[iel_ensm1 + 32 * iface_ensm1] + _pt_data_6[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_2_b_BTAG_PARTITION[from_el_indices_2[iel_ensm1 + 32 * iface_ensm1]] : 0.0); | |
| } | |
| for (int idof_ensm3_0 = 0; idof_ensm3_0 <= 9; ++idof_ensm3_0) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 31; ++iel_ensm3) | |
| { | |
| acc_iface_ensm3_idof_ensm3_0 = 0.0; | |
| acc_iface_ensm3_idof_ensm3 = 0.0; | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 3; ++idof_ensm3) | |
| for (int iface_ensm3 = 0; iface_ensm3 <= 2; ++iface_ensm3) | |
| { | |
| acc_iface_ensm3_idof_ensm3_0 = acc_iface_ensm3_idof_ensm3_0 + _pt_data_1[idof_ensm3 + 4 * iface_ensm3 + 12 * idof_ensm3_0] * _pt_data_2[32 * iface_ensm3 + iel_ensm3] * _pt_temp_1[idof_ensm3 + 128 * iface_ensm3 + 4 * iel_ensm3]; | |
| acc_iface_ensm3_idof_ensm3 = acc_iface_ensm3_idof_ensm3 + _pt_data_1[idof_ensm3 + 4 * iface_ensm3 + 12 * idof_ensm3_0] * _pt_data_2[32 * iface_ensm3 + iel_ensm3] * _pt_temp[idof_ensm3 + 128 * iface_ensm3 + 4 * iel_ensm3]; | |
| } | |
| face_mass_0[idof_ensm3_0 + 10 * iel_ensm3] = acc_iface_ensm3_idof_ensm3_0; | |
| face_mass[idof_ensm3_0 + 10 * iel_ensm3] = acc_iface_ensm3_idof_ensm3; | |
| } | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 9; ++idof_ensm5) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| { | |
| mass_inv_inp_0[idof_ensm5 + 10 * iel_ensm5] = _pt_part_ph_id[idof_ensm5 + 10 * iel_ensm5] + -1.0 * face_mass_0[idof_ensm5 + 10 * iel_ensm5]; | |
| mass_inv_inp[idof_ensm5 + 10 * iel_ensm5] = _pt_part_ph_id[320 + idof_ensm5 + 10 * iel_ensm5] + -1.0 * face_mass[idof_ensm5 + 10 * iel_ensm5]; | |
| } | |
| for (int idof_ensm6 = 0; idof_ensm6 <= 9; ++idof_ensm6) | |
| for (int iel_ensm6 = 0; iel_ensm6 <= 31; ++iel_ensm6) | |
| { | |
| acc_idof_ensm6_0 = 0.0; | |
| acc_idof_ensm6_1 = 0.0; | |
| for (int idof_ensm6_0 = 0; idof_ensm6_0 <= 9; ++idof_ensm6_0) | |
| { | |
| acc_idof_ensm6_0 = acc_idof_ensm6_0 + _pt_data[iel_ensm6] * _pt_data_0[idof_ensm6_0 + 10 * idof_ensm6] * mass_inv_inp[idof_ensm6_0 + 10 * iel_ensm6]; | |
| acc_idof_ensm6_1 = acc_idof_ensm6_1 + _pt_data[iel_ensm6] * _pt_data_0[idof_ensm6_0 + 10 * idof_ensm6] * mass_inv_inp_0[idof_ensm6_0 + 10 * iel_ensm6]; | |
| } | |
| _pt_temp_0[idof_ensm6 + 10 * iel_ensm6] = acc_idof_ensm6_0; | |
| _pt_temp_2[idof_ensm6 + 10 * iel_ensm6] = acc_idof_ensm6_1; | |
| } | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 9; ++idof_ensm4) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 31; ++iel_ensm4) | |
| { | |
| _pt_out_0_0[idof_ensm4 + 10 * iel_ensm4] = -1.0 * _pt_temp_2[idof_ensm4 + 10 * iel_ensm4]; | |
| _pt_out_1_0[idof_ensm4 + 10 * iel_ensm4] = -1.0 * _pt_temp_0[idof_ensm4 + 10 * iel_ensm4]; | |
| } | |
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs(double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_all, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, char const *__restrict__ from_el_present_3, double const *__restrict__ _pt_part_ph_id_4, double const *__restrict__ _pt_dist_id_1, int const *__restrict__ from_el_indices_5, int const *__restrict__ _pt_data_9, double const *__restrict__ normal_2_b_BTAG_PARTITION_1, int const *__restrict__ from_el_indices_6, int const *__restrict__ _pt_data_10, char const *__restrict__ from_el_present_4, double const *__restrict__ _pt_part_ph_id_5, double const *__restrict__ _pt_dist_id_2, int const *__restrict__ from_el_indices_7, int const *__restrict__ _pt_data_11, double const *__restrict__ normal_2_b_BTAG_PARTITION_2, int const *__restrict__ from_el_indices_8, int const *__restrict__ _pt_data_12, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_all, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double const *__restrict__ normal_1_b_BTAG_PARTITION_1, double const *__restrict__ normal_1_b_BTAG_PARTITION_2, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ cse_1, double *__restrict__ cse_2, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm7_0; | |
| double acc_idof_ensm7_1; | |
| double acc_iface_ensm4_idof_ensm4; | |
| double acc_iface_ensm4_idof_ensm4_0; | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| for (int iel_ensm0 = 0; iel_ensm0 <= 3; ++iel_ensm0) | |
| cse_1[idof_ensm0 + 4 * iel_ensm0] = 0.5 * (_pt_part_ph_id_4[idof_ensm0 + 4 * iel_ensm0] + _pt_dist_id_1[4 * from_el_indices_5[iel_ensm0] + _pt_data_9[idof_ensm0 + 4 * iel_ensm0]]); | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| for (int iel_ensm1 = 0; iel_ensm1 <= 1; ++iel_ensm1) | |
| cse[idof_ensm1 + 4 * iel_ensm1] = 0.5 * (_pt_part_ph_id_2[idof_ensm1 + 4 * iel_ensm1] + _pt_dist_id[4 * from_el_indices_1[iel_ensm1] + _pt_data_5[idof_ensm1 + 4 * iel_ensm1]]); | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| for (int iel_ensm2 = 0; iel_ensm2 <= 2; ++iel_ensm2) | |
| { | |
| cse_0[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_3[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id_0[4 * from_el_indices_3[iel_ensm2] + _pt_data_7[idof_ensm2 + 4 * iel_ensm2]]); | |
| cse_2[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_5[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id_2[4 * from_el_indices_7[iel_ensm2] + _pt_data_11[idof_ensm2 + 4 * iel_ensm2]]); | |
| } | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 3; ++idof_ensm3) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 31; ++iel_ensm3) | |
| for (int iface_ensm3 = 0; iface_ensm3 <= 2; ++iface_ensm3) | |
| { | |
| _pt_temp_1[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3] = (from_el_present[iel_ensm3 + 32 * iface_ensm3] ? normal_1_b_all[from_el_indices[iel_ensm3 + 32 * iface_ensm3]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm3 + 32 * iface_ensm3] + _pt_data_3[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm3 + 32 * iface_ensm3] + _pt_data_3[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]]) : 0.0) + (from_el_present_4[iel_ensm3 + 32 * iface_ensm3] ? cse_2[4 * from_el_indices_8[iel_ensm3 + 32 * iface_ensm3] + _pt_data_12[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION_2[from_el_indices_8[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_3[iel_ensm3 + 32 * iface_ensm3] ? cse_1[4 * from_el_indices_6[iel_ensm3 + 32 * iface_ensm3] + _pt_data_10[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_2[iel_ensm3 + 32 * iface_ensm3] ? cse_0[4 * from_el_indices_4[iel_ensm3 + 32 * iface_ensm3] + _pt_data_8[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_0[iel_ensm3 + 32 * iface_ensm3] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm3 + 32 * iface_ensm3] + _pt_data_4[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_face_restr_interior[from_el_indices_0[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_1[iel_ensm3 + 32 * iface_ensm3] ? cse[4 * from_el_indices_2[iel_ensm3 + 32 * iface_ensm3] + _pt_data_6[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION[from_el_indices_2[iel_ensm3 + 32 * iface_ensm3]] : 0.0); | |
| _pt_temp[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3] = (from_el_present[iel_ensm3 + 32 * iface_ensm3] ? normal_2_b_all[from_el_indices[iel_ensm3 + 32 * iface_ensm3]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm3 + 32 * iface_ensm3] + _pt_data_3[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm3 + 32 * iface_ensm3] + _pt_data_3[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]]) : 0.0) + (from_el_present_4[iel_ensm3 + 32 * iface_ensm3] ? cse_2[4 * from_el_indices_8[iel_ensm3 + 32 * iface_ensm3] + _pt_data_12[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION_2[from_el_indices_8[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_3[iel_ensm3 + 32 * iface_ensm3] ? cse_1[4 * from_el_indices_6[iel_ensm3 + 32 * iface_ensm3] + _pt_data_10[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_2[iel_ensm3 + 32 * iface_ensm3] ? cse_0[4 * from_el_indices_4[iel_ensm3 + 32 * iface_ensm3] + _pt_data_8[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_0[iel_ensm3 + 32 * iface_ensm3] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm3 + 32 * iface_ensm3] + _pt_data_4[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_face_restr_interior[from_el_indices_0[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_1[iel_ensm3 + 32 * iface_ensm3] ? cse[4 * from_el_indices_2[iel_ensm3 + 32 * iface_ensm3] + _pt_data_6[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION[from_el_indices_2[iel_ensm3 + 32 * iface_ensm3]] : 0.0); | |
| } | |
| for (int idof_ensm4_0 = 0; idof_ensm4_0 <= 9; ++idof_ensm4_0) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 31; ++iel_ensm4) | |
| { | |
| acc_iface_ensm4_idof_ensm4_0 = 0.0; | |
| acc_iface_ensm4_idof_ensm4 = 0.0; | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 3; ++idof_ensm4) | |
| for (int iface_ensm4 = 0; iface_ensm4 <= 2; ++iface_ensm4) | |
| { | |
| acc_iface_ensm4_idof_ensm4_0 = acc_iface_ensm4_idof_ensm4_0 + _pt_data_1[idof_ensm4 + 4 * iface_ensm4 + 12 * idof_ensm4_0] * _pt_data_2[32 * iface_ensm4 + iel_ensm4] * _pt_temp_1[idof_ensm4 + 128 * iface_ensm4 + 4 * iel_ensm4]; | |
| acc_iface_ensm4_idof_ensm4 = acc_iface_ensm4_idof_ensm4 + _pt_data_1[idof_ensm4 + 4 * iface_ensm4 + 12 * idof_ensm4_0] * _pt_data_2[32 * iface_ensm4 + iel_ensm4] * _pt_temp[idof_ensm4 + 128 * iface_ensm4 + 4 * iel_ensm4]; | |
| } | |
| face_mass_0[idof_ensm4_0 + 10 * iel_ensm4] = acc_iface_ensm4_idof_ensm4_0; | |
| face_mass[idof_ensm4_0 + 10 * iel_ensm4] = acc_iface_ensm4_idof_ensm4; | |
| } | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 9; ++idof_ensm5) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| { | |
| mass_inv_inp_0[idof_ensm5 + 10 * iel_ensm5] = _pt_part_ph_id[idof_ensm5 + 10 * iel_ensm5] + -1.0 * face_mass_0[idof_ensm5 + 10 * iel_ensm5]; | |
| mass_inv_inp[idof_ensm5 + 10 * iel_ensm5] = _pt_part_ph_id[320 + idof_ensm5 + 10 * iel_ensm5] + -1.0 * face_mass[idof_ensm5 + 10 * iel_ensm5]; | |
| } | |
| for (int idof_ensm7 = 0; idof_ensm7 <= 9; ++idof_ensm7) | |
| for (int iel_ensm7 = 0; iel_ensm7 <= 31; ++iel_ensm7) | |
| { | |
| acc_idof_ensm7_0 = 0.0; | |
| acc_idof_ensm7_1 = 0.0; | |
| for (int idof_ensm7_0 = 0; idof_ensm7_0 <= 9; ++idof_ensm7_0) | |
| { | |
| acc_idof_ensm7_0 = acc_idof_ensm7_0 + _pt_data[iel_ensm7] * _pt_data_0[idof_ensm7_0 + 10 * idof_ensm7] * mass_inv_inp[idof_ensm7_0 + 10 * iel_ensm7]; | |
| acc_idof_ensm7_1 = acc_idof_ensm7_1 + _pt_data[iel_ensm7] * _pt_data_0[idof_ensm7_0 + 10 * idof_ensm7] * mass_inv_inp_0[idof_ensm7_0 + 10 * iel_ensm7]; | |
| } | |
| _pt_temp_0[idof_ensm7 + 10 * iel_ensm7] = acc_idof_ensm7_0; | |
| _pt_temp_2[idof_ensm7 + 10 * iel_ensm7] = acc_idof_ensm7_1; | |
| } | |
| for (int idof_ensm6 = 0; idof_ensm6 <= 9; ++idof_ensm6) | |
| for (int iel_ensm6 = 0; iel_ensm6 <= 31; ++iel_ensm6) | |
| { | |
| _pt_out_0_0[idof_ensm6 + 10 * iel_ensm6] = -1.0 * _pt_temp_2[idof_ensm6 + 10 * iel_ensm6]; | |
| _pt_out_1_0[idof_ensm6 + 10 * iel_ensm6] = -1.0 * _pt_temp_0[idof_ensm6 + 10 * iel_ensm6]; | |
| } | |
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs(double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_all, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, char const *__restrict__ from_el_present_3, double const *__restrict__ _pt_part_ph_id_4, double const *__restrict__ _pt_dist_id_1, int const *__restrict__ from_el_indices_5, int const *__restrict__ _pt_data_9, double const *__restrict__ normal_2_b_BTAG_PARTITION_1, int const *__restrict__ from_el_indices_6, int const *__restrict__ _pt_data_10, char const *__restrict__ from_el_present_4, double const *__restrict__ _pt_part_ph_id_5, double const *__restrict__ _pt_dist_id_2, int const *__restrict__ from_el_indices_7, int const *__restrict__ _pt_data_11, double const *__restrict__ normal_2_b_BTAG_PARTITION_2, int const *__restrict__ from_el_indices_8, int const *__restrict__ _pt_data_12, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_all, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double const *__restrict__ normal_1_b_BTAG_PARTITION_1, double const *__restrict__ normal_1_b_BTAG_PARTITION_2, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ cse_1, double *__restrict__ cse_2, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm6_0; | |
| double acc_idof_ensm6_1; | |
| double acc_iface_ensm3_idof_ensm3; | |
| double acc_iface_ensm3_idof_ensm3_0; | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| for (int iel_ensm0 = 0; iel_ensm0 <= 3; ++iel_ensm0) | |
| { | |
| cse_1[idof_ensm0 + 4 * iel_ensm0] = 0.5 * (_pt_part_ph_id_4[idof_ensm0 + 4 * iel_ensm0] + _pt_dist_id_1[4 * from_el_indices_5[iel_ensm0] + _pt_data_9[idof_ensm0 + 4 * iel_ensm0]]); | |
| cse_2[idof_ensm0 + 4 * iel_ensm0] = 0.5 * (_pt_part_ph_id_5[idof_ensm0 + 4 * iel_ensm0] + _pt_dist_id_2[4 * from_el_indices_7[iel_ensm0] + _pt_data_11[idof_ensm0 + 4 * iel_ensm0]]); | |
| cse[idof_ensm0 + 4 * iel_ensm0] = 0.5 * (_pt_part_ph_id_2[idof_ensm0 + 4 * iel_ensm0] + _pt_dist_id[4 * from_el_indices_1[iel_ensm0] + _pt_data_5[idof_ensm0 + 4 * iel_ensm0]]); | |
| } | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| for (int iel_ensm2 = 0; iel_ensm2 <= 1; ++iel_ensm2) | |
| cse_0[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_3[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id_0[4 * from_el_indices_3[iel_ensm2] + _pt_data_7[idof_ensm2 + 4 * iel_ensm2]]); | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| for (int iel_ensm1 = 0; iel_ensm1 <= 31; ++iel_ensm1) | |
| for (int iface_ensm1 = 0; iface_ensm1 <= 2; ++iface_ensm1) | |
| { | |
| _pt_temp_1[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1] = (from_el_present[iel_ensm1 + 32 * iface_ensm1] ? normal_1_b_all[from_el_indices[iel_ensm1 + 32 * iface_ensm1]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm1 + 32 * iface_ensm1] + _pt_data_3[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm1 + 32 * iface_ensm1] + _pt_data_3[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]]) : 0.0) + (from_el_present_4[iel_ensm1 + 32 * iface_ensm1] ? cse_2[4 * from_el_indices_8[iel_ensm1 + 32 * iface_ensm1] + _pt_data_12[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_1_b_BTAG_PARTITION_2[from_el_indices_8[iel_ensm1 + 32 * iface_ensm1]] : 0.0) + (from_el_present_3[iel_ensm1 + 32 * iface_ensm1] ? cse_1[4 * from_el_indices_6[iel_ensm1 + 32 * iface_ensm1] + _pt_data_10[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_1_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm1 + 32 * iface_ensm1]] : 0.0) + (from_el_present_2[iel_ensm1 + 32 * iface_ensm1] ? cse_0[4 * from_el_indices_4[iel_ensm1 + 32 * iface_ensm1] + _pt_data_8[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm1 + 32 * iface_ensm1]] : 0.0) + (from_el_present_0[iel_ensm1 + 32 * iface_ensm1] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm1 + 32 * iface_ensm1] + _pt_data_4[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_1_b_face_restr_interior[from_el_indices_0[iel_ensm1 + 32 * iface_ensm1]] : 0.0) + (from_el_present_1[iel_ensm1 + 32 * iface_ensm1] ? cse[4 * from_el_indices_2[iel_ensm1 + 32 * iface_ensm1] + _pt_data_6[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_1_b_BTAG_PARTITION[from_el_indices_2[iel_ensm1 + 32 * iface_ensm1]] : 0.0); | |
| _pt_temp[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1] = (from_el_present[iel_ensm1 + 32 * iface_ensm1] ? normal_2_b_all[from_el_indices[iel_ensm1 + 32 * iface_ensm1]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm1 + 32 * iface_ensm1] + _pt_data_3[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm1 + 32 * iface_ensm1] + _pt_data_3[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]]) : 0.0) + (from_el_present_4[iel_ensm1 + 32 * iface_ensm1] ? cse_2[4 * from_el_indices_8[iel_ensm1 + 32 * iface_ensm1] + _pt_data_12[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_2_b_BTAG_PARTITION_2[from_el_indices_8[iel_ensm1 + 32 * iface_ensm1]] : 0.0) + (from_el_present_3[iel_ensm1 + 32 * iface_ensm1] ? cse_1[4 * from_el_indices_6[iel_ensm1 + 32 * iface_ensm1] + _pt_data_10[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_2_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm1 + 32 * iface_ensm1]] : 0.0) + (from_el_present_2[iel_ensm1 + 32 * iface_ensm1] ? cse_0[4 * from_el_indices_4[iel_ensm1 + 32 * iface_ensm1] + _pt_data_8[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm1 + 32 * iface_ensm1]] : 0.0) + (from_el_present_0[iel_ensm1 + 32 * iface_ensm1] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm1 + 32 * iface_ensm1] + _pt_data_4[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_2_b_face_restr_interior[from_el_indices_0[iel_ensm1 + 32 * iface_ensm1]] : 0.0) + (from_el_present_1[iel_ensm1 + 32 * iface_ensm1] ? cse[4 * from_el_indices_2[iel_ensm1 + 32 * iface_ensm1] + _pt_data_6[idof_ensm1 + 4 * iel_ensm1 + 128 * iface_ensm1]] * normal_2_b_BTAG_PARTITION[from_el_indices_2[iel_ensm1 + 32 * iface_ensm1]] : 0.0); | |
| } | |
| for (int idof_ensm3_0 = 0; idof_ensm3_0 <= 9; ++idof_ensm3_0) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 31; ++iel_ensm3) | |
| { | |
| acc_iface_ensm3_idof_ensm3_0 = 0.0; | |
| acc_iface_ensm3_idof_ensm3 = 0.0; | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 3; ++idof_ensm3) | |
| for (int iface_ensm3 = 0; iface_ensm3 <= 2; ++iface_ensm3) | |
| { | |
| acc_iface_ensm3_idof_ensm3_0 = acc_iface_ensm3_idof_ensm3_0 + _pt_data_1[idof_ensm3 + 4 * iface_ensm3 + 12 * idof_ensm3_0] * _pt_data_2[32 * iface_ensm3 + iel_ensm3] * _pt_temp_1[idof_ensm3 + 128 * iface_ensm3 + 4 * iel_ensm3]; | |
| acc_iface_ensm3_idof_ensm3 = acc_iface_ensm3_idof_ensm3 + _pt_data_1[idof_ensm3 + 4 * iface_ensm3 + 12 * idof_ensm3_0] * _pt_data_2[32 * iface_ensm3 + iel_ensm3] * _pt_temp[idof_ensm3 + 128 * iface_ensm3 + 4 * iel_ensm3]; | |
| } | |
| face_mass_0[idof_ensm3_0 + 10 * iel_ensm3] = acc_iface_ensm3_idof_ensm3_0; | |
| face_mass[idof_ensm3_0 + 10 * iel_ensm3] = acc_iface_ensm3_idof_ensm3; | |
| } | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 9; ++idof_ensm5) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| { | |
| mass_inv_inp_0[idof_ensm5 + 10 * iel_ensm5] = _pt_part_ph_id[idof_ensm5 + 10 * iel_ensm5] + -1.0 * face_mass_0[idof_ensm5 + 10 * iel_ensm5]; | |
| mass_inv_inp[idof_ensm5 + 10 * iel_ensm5] = _pt_part_ph_id[320 + idof_ensm5 + 10 * iel_ensm5] + -1.0 * face_mass[idof_ensm5 + 10 * iel_ensm5]; | |
| } | |
| for (int idof_ensm6 = 0; idof_ensm6 <= 9; ++idof_ensm6) | |
| for (int iel_ensm6 = 0; iel_ensm6 <= 31; ++iel_ensm6) | |
| { | |
| acc_idof_ensm6_0 = 0.0; | |
| acc_idof_ensm6_1 = 0.0; | |
| for (int idof_ensm6_0 = 0; idof_ensm6_0 <= 9; ++idof_ensm6_0) | |
| { | |
| acc_idof_ensm6_0 = acc_idof_ensm6_0 + _pt_data[iel_ensm6] * _pt_data_0[idof_ensm6_0 + 10 * idof_ensm6] * mass_inv_inp[idof_ensm6_0 + 10 * iel_ensm6]; | |
| acc_idof_ensm6_1 = acc_idof_ensm6_1 + _pt_data[iel_ensm6] * _pt_data_0[idof_ensm6_0 + 10 * idof_ensm6] * mass_inv_inp_0[idof_ensm6_0 + 10 * iel_ensm6]; | |
| } | |
| _pt_temp_0[idof_ensm6 + 10 * iel_ensm6] = acc_idof_ensm6_0; | |
| _pt_temp_2[idof_ensm6 + 10 * iel_ensm6] = acc_idof_ensm6_1; | |
| } | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 9; ++idof_ensm4) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 31; ++iel_ensm4) | |
| { | |
| _pt_out_0_0[idof_ensm4 + 10 * iel_ensm4] = -1.0 * _pt_temp_2[idof_ensm4 + 10 * iel_ensm4]; | |
| _pt_out_1_0[idof_ensm4 + 10 * iel_ensm4] = -1.0 * _pt_temp_0[idof_ensm4 + 10 * iel_ensm4]; | |
| } | |
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs(double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_1, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, double const *__restrict__ normal_2_b_BTAG_PARTITION_1, int const *__restrict__ from_el_indices_5, int const *__restrict__ _pt_data_9, char const *__restrict__ from_el_present_3, double const *__restrict__ _pt_part_ph_id_4, double const *__restrict__ _pt_dist_id_2, int const *__restrict__ from_el_indices_6, int const *__restrict__ _pt_data_10, double const *__restrict__ normal_2_b_BTAG_PARTITION_2, int const *__restrict__ from_el_indices_7, int const *__restrict__ _pt_data_11, char const *__restrict__ from_el_present_4, double const *__restrict__ _pt_part_ph_id_5, double const *__restrict__ _pt_dist_id_3, int const *__restrict__ from_el_indices_8, int const *__restrict__ _pt_data_12, double const *__restrict__ normal_2_b_BTAG_PARTITION_3, int const *__restrict__ from_el_indices_9, int const *__restrict__ _pt_data_13, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double const *__restrict__ normal_1_b_BTAG_PARTITION_1, double const *__restrict__ normal_1_b_BTAG_PARTITION_2, double const *__restrict__ normal_1_b_BTAG_PARTITION_3, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ cse_1, double *__restrict__ cse_2, double *__restrict__ cse_3, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm7_0; | |
| double acc_idof_ensm7_1; | |
| double acc_iface_ensm5_idof_ensm5; | |
| double acc_iface_ensm5_idof_ensm5_0; | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| for (int iel_ensm0 = 0; iel_ensm0 <= 2; ++iel_ensm0) | |
| cse_1[idof_ensm0 + 4 * iel_ensm0] = 0.5 * (_pt_part_ph_id_3[idof_ensm0 + 4 * iel_ensm0] + _pt_dist_id_1[4 * from_el_indices_4[iel_ensm0] + _pt_data_8[idof_ensm0 + 4 * iel_ensm0]]); | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| { | |
| int const iel_ensm1 = 0; | |
| cse_2[idof_ensm1] = 0.5 * (_pt_part_ph_id_4[idof_ensm1] + _pt_dist_id_2[4 * from_el_indices_6[0] + _pt_data_10[idof_ensm1]]); | |
| } | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| for (int iel_ensm2 = 0; iel_ensm2 <= 5; ++iel_ensm2) | |
| cse_3[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_5[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id_3[4 * from_el_indices_8[iel_ensm2] + _pt_data_12[idof_ensm2 + 4 * iel_ensm2]]); | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 3; ++idof_ensm4) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 4; ++iel_ensm4) | |
| { | |
| cse_0[idof_ensm4 + 4 * iel_ensm4] = 0.5 * (_pt_part_ph_id_2[idof_ensm4 + 4 * iel_ensm4] + _pt_dist_id_0[4 * from_el_indices_2[iel_ensm4] + _pt_data_6[idof_ensm4 + 4 * iel_ensm4]]); | |
| cse[idof_ensm4 + 4 * iel_ensm4] = 0.5 * (_pt_part_ph_id_1[idof_ensm4 + 4 * iel_ensm4] + _pt_dist_id[4 * from_el_indices_0[iel_ensm4] + _pt_data_4[idof_ensm4 + 4 * iel_ensm4]]); | |
| } | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 3; ++idof_ensm3) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 31; ++iel_ensm3) | |
| for (int iface_ensm3 = 0; iface_ensm3 <= 2; ++iface_ensm3) | |
| { | |
| _pt_temp_1[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3] = (from_el_present_4[iel_ensm3 + 32 * iface_ensm3] ? cse_3[4 * from_el_indices_9[iel_ensm3 + 32 * iface_ensm3] + _pt_data_13[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION_3[from_el_indices_9[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_3[iel_ensm3 + 32 * iface_ensm3] ? cse_2[4 * from_el_indices_7[iel_ensm3 + 32 * iface_ensm3] + _pt_data_11[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION_2[from_el_indices_7[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_2[iel_ensm3 + 32 * iface_ensm3] ? cse_1[4 * from_el_indices_5[iel_ensm3 + 32 * iface_ensm3] + _pt_data_9[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION_1[from_el_indices_5[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_1[iel_ensm3 + 32 * iface_ensm3] ? cse_0[4 * from_el_indices_3[iel_ensm3 + 32 * iface_ensm3] + _pt_data_7[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_3[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present[iel_ensm3 + 32 * iface_ensm3] ? _pt_part_ph_id_0[4 * from_el_indices[iel_ensm3 + 32 * iface_ensm3] + _pt_data_3[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_face_restr_interior[from_el_indices[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_0[iel_ensm3 + 32 * iface_ensm3] ? cse[4 * from_el_indices_1[iel_ensm3 + 32 * iface_ensm3] + _pt_data_5[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION[from_el_indices_1[iel_ensm3 + 32 * iface_ensm3]] : 0.0); | |
| _pt_temp[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3] = (from_el_present_4[iel_ensm3 + 32 * iface_ensm3] ? cse_3[4 * from_el_indices_9[iel_ensm3 + 32 * iface_ensm3] + _pt_data_13[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION_3[from_el_indices_9[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_3[iel_ensm3 + 32 * iface_ensm3] ? cse_2[4 * from_el_indices_7[iel_ensm3 + 32 * iface_ensm3] + _pt_data_11[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION_2[from_el_indices_7[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_2[iel_ensm3 + 32 * iface_ensm3] ? cse_1[4 * from_el_indices_5[iel_ensm3 + 32 * iface_ensm3] + _pt_data_9[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION_1[from_el_indices_5[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_1[iel_ensm3 + 32 * iface_ensm3] ? cse_0[4 * from_el_indices_3[iel_ensm3 + 32 * iface_ensm3] + _pt_data_7[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_3[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present[iel_ensm3 + 32 * iface_ensm3] ? _pt_part_ph_id_0[4 * from_el_indices[iel_ensm3 + 32 * iface_ensm3] + _pt_data_3[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_face_restr_interior[from_el_indices[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_0[iel_ensm3 + 32 * iface_ensm3] ? cse[4 * from_el_indices_1[iel_ensm3 + 32 * iface_ensm3] + _pt_data_5[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION[from_el_indices_1[iel_ensm3 + 32 * iface_ensm3]] : 0.0); | |
| } | |
| for (int idof_ensm5_0 = 0; idof_ensm5_0 <= 9; ++idof_ensm5_0) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| { | |
| acc_iface_ensm5_idof_ensm5_0 = 0.0; | |
| acc_iface_ensm5_idof_ensm5 = 0.0; | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 3; ++idof_ensm5) | |
| for (int iface_ensm5 = 0; iface_ensm5 <= 2; ++iface_ensm5) | |
| { | |
| acc_iface_ensm5_idof_ensm5_0 = acc_iface_ensm5_idof_ensm5_0 + _pt_data_1[idof_ensm5 + 4 * iface_ensm5 + 12 * idof_ensm5_0] * _pt_data_2[32 * iface_ensm5 + iel_ensm5] * _pt_temp_1[idof_ensm5 + 128 * iface_ensm5 + 4 * iel_ensm5]; | |
| acc_iface_ensm5_idof_ensm5 = acc_iface_ensm5_idof_ensm5 + _pt_data_1[idof_ensm5 + 4 * iface_ensm5 + 12 * idof_ensm5_0] * _pt_data_2[32 * iface_ensm5 + iel_ensm5] * _pt_temp[idof_ensm5 + 128 * iface_ensm5 + 4 * iel_ensm5]; | |
| } | |
| face_mass_0[idof_ensm5_0 + 10 * iel_ensm5] = acc_iface_ensm5_idof_ensm5_0; | |
| face_mass[idof_ensm5_0 + 10 * iel_ensm5] = acc_iface_ensm5_idof_ensm5; | |
| } | |
| for (int idof_ensm6 = 0; idof_ensm6 <= 9; ++idof_ensm6) | |
| for (int iel_ensm6 = 0; iel_ensm6 <= 31; ++iel_ensm6) | |
| { | |
| mass_inv_inp_0[idof_ensm6 + 10 * iel_ensm6] = _pt_part_ph_id[idof_ensm6 + 10 * iel_ensm6] + -1.0 * face_mass_0[idof_ensm6 + 10 * iel_ensm6]; | |
| mass_inv_inp[idof_ensm6 + 10 * iel_ensm6] = _pt_part_ph_id[320 + idof_ensm6 + 10 * iel_ensm6] + -1.0 * face_mass[idof_ensm6 + 10 * iel_ensm6]; | |
| } | |
| for (int idof_ensm7 = 0; idof_ensm7 <= 9; ++idof_ensm7) | |
| for (int iel_ensm7 = 0; iel_ensm7 <= 31; ++iel_ensm7) | |
| { | |
| acc_idof_ensm7_0 = 0.0; | |
| acc_idof_ensm7_1 = 0.0; | |
| for (int idof_ensm7_0 = 0; idof_ensm7_0 <= 9; ++idof_ensm7_0) | |
| { | |
| acc_idof_ensm7_0 = acc_idof_ensm7_0 + _pt_data[iel_ensm7] * _pt_data_0[idof_ensm7_0 + 10 * idof_ensm7] * mass_inv_inp[idof_ensm7_0 + 10 * iel_ensm7]; | |
| acc_idof_ensm7_1 = acc_idof_ensm7_1 + _pt_data[iel_ensm7] * _pt_data_0[idof_ensm7_0 + 10 * idof_ensm7] * mass_inv_inp_0[idof_ensm7_0 + 10 * iel_ensm7]; | |
| } | |
| _pt_temp_0[idof_ensm7 + 10 * iel_ensm7] = acc_idof_ensm7_0; | |
| _pt_temp_2[idof_ensm7 + 10 * iel_ensm7] = acc_idof_ensm7_1; | |
| } | |
| for (int idof_ensm8 = 0; idof_ensm8 <= 9; ++idof_ensm8) | |
| for (int iel_ensm8 = 0; iel_ensm8 <= 31; ++iel_ensm8) | |
| { | |
| _pt_out_0_0[idof_ensm8 + 10 * iel_ensm8] = -1.0 * _pt_temp_2[idof_ensm8 + 10 * iel_ensm8]; | |
| _pt_out_1_0[idof_ensm8 + 10 * iel_ensm8] = -1.0 * _pt_temp_0[idof_ensm8 + 10 * iel_ensm8]; | |
| } | |
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs(double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_all, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, char const *__restrict__ from_el_present_3, double const *__restrict__ _pt_part_ph_id_4, double const *__restrict__ _pt_dist_id_1, int const *__restrict__ from_el_indices_5, int const *__restrict__ _pt_data_9, double const *__restrict__ normal_2_b_BTAG_PARTITION_1, int const *__restrict__ from_el_indices_6, int const *__restrict__ _pt_data_10, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_all, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double const *__restrict__ normal_1_b_BTAG_PARTITION_1, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ cse_1, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm6_0; | |
| double acc_idof_ensm6_1; | |
| double acc_iface_ensm4_idof_ensm4; | |
| double acc_iface_ensm4_idof_ensm4_0; | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| for (int iel_ensm1 = 0; iel_ensm1 <= 4; ++iel_ensm1) | |
| cse_0[idof_ensm1 + 4 * iel_ensm1] = 0.5 * (_pt_part_ph_id_3[idof_ensm1 + 4 * iel_ensm1] + _pt_dist_id_0[4 * from_el_indices_3[iel_ensm1] + _pt_data_7[idof_ensm1 + 4 * iel_ensm1]]); | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| for (int iel_ensm2 = 0; iel_ensm2 <= 1; ++iel_ensm2) | |
| cse_1[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_4[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id_1[4 * from_el_indices_5[iel_ensm2] + _pt_data_9[idof_ensm2 + 4 * iel_ensm2]]); | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 3; ++idof_ensm3) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 5; ++iel_ensm3) | |
| cse[idof_ensm3 + 4 * iel_ensm3] = 0.5 * (_pt_part_ph_id_2[idof_ensm3 + 4 * iel_ensm3] + _pt_dist_id[4 * from_el_indices_1[iel_ensm3] + _pt_data_5[idof_ensm3 + 4 * iel_ensm3]]); | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| for (int iel_ensm0 = 0; iel_ensm0 <= 31; ++iel_ensm0) | |
| for (int iface_ensm0 = 0; iface_ensm0 <= 2; ++iface_ensm0) | |
| { | |
| _pt_temp_1[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0] = (from_el_present[iel_ensm0 + 32 * iface_ensm0] ? normal_1_b_all[from_el_indices[iel_ensm0 + 32 * iface_ensm0]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]]) : 0.0) + (from_el_present_3[iel_ensm0 + 32 * iface_ensm0] ? cse_1[4 * from_el_indices_6[iel_ensm0 + 32 * iface_ensm0] + _pt_data_10[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_2[iel_ensm0 + 32 * iface_ensm0] ? cse_0[4 * from_el_indices_4[iel_ensm0 + 32 * iface_ensm0] + _pt_data_8[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_0[iel_ensm0 + 32 * iface_ensm0] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm0 + 32 * iface_ensm0] + _pt_data_4[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_face_restr_interior[from_el_indices_0[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_1[iel_ensm0 + 32 * iface_ensm0] ? cse[4 * from_el_indices_2[iel_ensm0 + 32 * iface_ensm0] + _pt_data_6[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION[from_el_indices_2[iel_ensm0 + 32 * iface_ensm0]] : 0.0); | |
| _pt_temp[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0] = (from_el_present[iel_ensm0 + 32 * iface_ensm0] ? normal_2_b_all[from_el_indices[iel_ensm0 + 32 * iface_ensm0]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]]) : 0.0) + (from_el_present_3[iel_ensm0 + 32 * iface_ensm0] ? cse_1[4 * from_el_indices_6[iel_ensm0 + 32 * iface_ensm0] + _pt_data_10[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_2[iel_ensm0 + 32 * iface_ensm0] ? cse_0[4 * from_el_indices_4[iel_ensm0 + 32 * iface_ensm0] + _pt_data_8[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_0[iel_ensm0 + 32 * iface_ensm0] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm0 + 32 * iface_ensm0] + _pt_data_4[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_face_restr_interior[from_el_indices_0[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_1[iel_ensm0 + 32 * iface_ensm0] ? cse[4 * from_el_indices_2[iel_ensm0 + 32 * iface_ensm0] + _pt_data_6[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION[from_el_indices_2[iel_ensm0 + 32 * iface_ensm0]] : 0.0); | |
| } | |
| for (int idof_ensm4_0 = 0; idof_ensm4_0 <= 9; ++idof_ensm4_0) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 31; ++iel_ensm4) | |
| { | |
| acc_iface_ensm4_idof_ensm4_0 = 0.0; | |
| acc_iface_ensm4_idof_ensm4 = 0.0; | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 3; ++idof_ensm4) | |
| for (int iface_ensm4 = 0; iface_ensm4 <= 2; ++iface_ensm4) | |
| { | |
| acc_iface_ensm4_idof_ensm4_0 = acc_iface_ensm4_idof_ensm4_0 + _pt_data_1[idof_ensm4 + 4 * iface_ensm4 + 12 * idof_ensm4_0] * _pt_data_2[32 * iface_ensm4 + iel_ensm4] * _pt_temp_1[idof_ensm4 + 128 * iface_ensm4 + 4 * iel_ensm4]; | |
| acc_iface_ensm4_idof_ensm4 = acc_iface_ensm4_idof_ensm4 + _pt_data_1[idof_ensm4 + 4 * iface_ensm4 + 12 * idof_ensm4_0] * _pt_data_2[32 * iface_ensm4 + iel_ensm4] * _pt_temp[idof_ensm4 + 128 * iface_ensm4 + 4 * iel_ensm4]; | |
| } | |
| face_mass_0[idof_ensm4_0 + 10 * iel_ensm4] = acc_iface_ensm4_idof_ensm4_0; | |
| face_mass[idof_ensm4_0 + 10 * iel_ensm4] = acc_iface_ensm4_idof_ensm4; | |
| } | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 9; ++idof_ensm5) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| { | |
| mass_inv_inp_0[idof_ensm5 + 10 * iel_ensm5] = _pt_part_ph_id[idof_ensm5 + 10 * iel_ensm5] + -1.0 * face_mass_0[idof_ensm5 + 10 * iel_ensm5]; | |
| mass_inv_inp[idof_ensm5 + 10 * iel_ensm5] = _pt_part_ph_id[320 + idof_ensm5 + 10 * iel_ensm5] + -1.0 * face_mass[idof_ensm5 + 10 * iel_ensm5]; | |
| } | |
| for (int idof_ensm6 = 0; idof_ensm6 <= 9; ++idof_ensm6) | |
| for (int iel_ensm6 = 0; iel_ensm6 <= 31; ++iel_ensm6) | |
| { | |
| acc_idof_ensm6_0 = 0.0; | |
| acc_idof_ensm6_1 = 0.0; | |
| for (int idof_ensm6_0 = 0; idof_ensm6_0 <= 9; ++idof_ensm6_0) | |
| { | |
| acc_idof_ensm6_0 = acc_idof_ensm6_0 + _pt_data[iel_ensm6] * _pt_data_0[idof_ensm6_0 + 10 * idof_ensm6] * mass_inv_inp[idof_ensm6_0 + 10 * iel_ensm6]; | |
| acc_idof_ensm6_1 = acc_idof_ensm6_1 + _pt_data[iel_ensm6] * _pt_data_0[idof_ensm6_0 + 10 * idof_ensm6] * mass_inv_inp_0[idof_ensm6_0 + 10 * iel_ensm6]; | |
| } | |
| _pt_temp_0[idof_ensm6 + 10 * iel_ensm6] = acc_idof_ensm6_0; | |
| _pt_temp_2[idof_ensm6 + 10 * iel_ensm6] = acc_idof_ensm6_1; | |
| } | |
| for (int idof_ensm7 = 0; idof_ensm7 <= 9; ++idof_ensm7) | |
| for (int iel_ensm7 = 0; iel_ensm7 <= 31; ++iel_ensm7) | |
| { | |
| _pt_out_0_0[idof_ensm7 + 10 * iel_ensm7] = -1.0 * _pt_temp_2[idof_ensm7 + 10 * iel_ensm7]; | |
| _pt_out_1_0[idof_ensm7 + 10 * iel_ensm7] = -1.0 * _pt_temp_0[idof_ensm7 + 10 * iel_ensm7]; | |
| } | |
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs(double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_all, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, char const *__restrict__ from_el_present_3, double const *__restrict__ _pt_part_ph_id_4, double const *__restrict__ _pt_dist_id_1, int const *__restrict__ from_el_indices_5, int const *__restrict__ _pt_data_9, double const *__restrict__ normal_2_b_BTAG_PARTITION_1, int const *__restrict__ from_el_indices_6, int const *__restrict__ _pt_data_10, char const *__restrict__ from_el_present_4, double const *__restrict__ _pt_part_ph_id_5, double const *__restrict__ _pt_dist_id_2, int const *__restrict__ from_el_indices_7, int const *__restrict__ _pt_data_11, double const *__restrict__ normal_2_b_BTAG_PARTITION_2, int const *__restrict__ from_el_indices_8, int const *__restrict__ _pt_data_12, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_all, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double const *__restrict__ normal_1_b_BTAG_PARTITION_1, double const *__restrict__ normal_1_b_BTAG_PARTITION_2, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ cse_1, double *__restrict__ cse_2, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm6_0; | |
| double acc_idof_ensm6_1; | |
| double acc_iface_ensm4_idof_ensm4; | |
| double acc_iface_ensm4_idof_ensm4_0; | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| for (int iel_ensm0 = 0; iel_ensm0 <= 1; ++iel_ensm0) | |
| { | |
| cse_1[idof_ensm0 + 4 * iel_ensm0] = 0.5 * (_pt_part_ph_id_4[idof_ensm0 + 4 * iel_ensm0] + _pt_dist_id_1[4 * from_el_indices_5[iel_ensm0] + _pt_data_9[idof_ensm0 + 4 * iel_ensm0]]); | |
| cse[idof_ensm0 + 4 * iel_ensm0] = 0.5 * (_pt_part_ph_id_2[idof_ensm0 + 4 * iel_ensm0] + _pt_dist_id[4 * from_el_indices_1[iel_ensm0] + _pt_data_5[idof_ensm0 + 4 * iel_ensm0]]); | |
| } | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| for (int iel_ensm1 = 0; iel_ensm1 <= 4; ++iel_ensm1) | |
| cse_0[idof_ensm1 + 4 * iel_ensm1] = 0.5 * (_pt_part_ph_id_3[idof_ensm1 + 4 * iel_ensm1] + _pt_dist_id_0[4 * from_el_indices_3[iel_ensm1] + _pt_data_7[idof_ensm1 + 4 * iel_ensm1]]); | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| for (int iel_ensm2 = 0; iel_ensm2 <= 5; ++iel_ensm2) | |
| cse_2[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_5[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id_2[4 * from_el_indices_7[iel_ensm2] + _pt_data_11[idof_ensm2 + 4 * iel_ensm2]]); | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 3; ++idof_ensm3) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 31; ++iel_ensm3) | |
| for (int iface_ensm3 = 0; iface_ensm3 <= 2; ++iface_ensm3) | |
| { | |
| _pt_temp_1[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3] = | |
| (from_el_present[iel_ensm3 + 32 * iface_ensm3] ? | |
| normal_1_b_all[from_el_indices[iel_ensm3 + 32 * iface_ensm3]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm3 + 32 * iface_ensm3] + _pt_data_3[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm3 + 32 * iface_ensm3] + _pt_data_3[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]]) | |
| : 0.0) | |
| + (from_el_present_4[iel_ensm3 + 32 * iface_ensm3] ? | |
| cse_2[4 * from_el_indices_8[iel_ensm3 + 32 * iface_ensm3] + _pt_data_12[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION_2[from_el_indices_8[iel_ensm3 + 32 * iface_ensm3]] | |
| : 0.0) | |
| + (from_el_present_3[iel_ensm3 + 32 * iface_ensm3] ? | |
| cse_1[4 * from_el_indices_6[iel_ensm3 + 32 * iface_ensm3] + _pt_data_10[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm3 + 32 * iface_ensm3]] | |
| : 0.0) | |
| + (from_el_present_2[iel_ensm3 + 32 * iface_ensm3] ? | |
| cse_0[4 * from_el_indices_4[iel_ensm3 + 32 * iface_ensm3] + _pt_data_8[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm3 + 32 * iface_ensm3]] | |
| : 0.0) | |
| + (from_el_present_0[iel_ensm3 + 32 * iface_ensm3] ? | |
| _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm3 + 32 * iface_ensm3] + _pt_data_4[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_face_restr_interior[from_el_indices_0[iel_ensm3 + 32 * iface_ensm3]] | |
| : 0.0) | |
| + (from_el_present_1[iel_ensm3 + 32 * iface_ensm3] ? cse[4 * from_el_indices_2[iel_ensm3 + 32 * iface_ensm3] + _pt_data_6[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_1_b_BTAG_PARTITION[from_el_indices_2[iel_ensm3 + 32 * iface_ensm3]] | |
| : 0.0); | |
| _pt_temp[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3] = (from_el_present[iel_ensm3 + 32 * iface_ensm3] ? normal_2_b_all[from_el_indices[iel_ensm3 + 32 * iface_ensm3]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm3 + 32 * iface_ensm3] + _pt_data_3[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm3 + 32 * iface_ensm3] + _pt_data_3[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]]) : 0.0) + (from_el_present_4[iel_ensm3 + 32 * iface_ensm3] ? cse_2[4 * from_el_indices_8[iel_ensm3 + 32 * iface_ensm3] + _pt_data_12[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION_2[from_el_indices_8[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_3[iel_ensm3 + 32 * iface_ensm3] ? cse_1[4 * from_el_indices_6[iel_ensm3 + 32 * iface_ensm3] + _pt_data_10[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_2[iel_ensm3 + 32 * iface_ensm3] ? cse_0[4 * from_el_indices_4[iel_ensm3 + 32 * iface_ensm3] + _pt_data_8[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_0[iel_ensm3 + 32 * iface_ensm3] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm3 + 32 * iface_ensm3] + _pt_data_4[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_face_restr_interior[from_el_indices_0[iel_ensm3 + 32 * iface_ensm3]] : 0.0) + (from_el_present_1[iel_ensm3 + 32 * iface_ensm3] ? cse[4 * from_el_indices_2[iel_ensm3 + 32 * iface_ensm3] + _pt_data_6[idof_ensm3 + 4 * iel_ensm3 + 128 * iface_ensm3]] * normal_2_b_BTAG_PARTITION[from_el_indices_2[iel_ensm3 + 32 * iface_ensm3]] : 0.0); | |
| } | |
| for (int idof_ensm4_0 = 0; idof_ensm4_0 <= 9; ++idof_ensm4_0) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 31; ++iel_ensm4) | |
| { | |
| acc_iface_ensm4_idof_ensm4_0 = 0.0; | |
| acc_iface_ensm4_idof_ensm4 = 0.0; | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 3; ++idof_ensm4) | |
| for (int iface_ensm4 = 0; iface_ensm4 <= 2; ++iface_ensm4) | |
| { | |
| acc_iface_ensm4_idof_ensm4_0 = acc_iface_ensm4_idof_ensm4_0 + _pt_data_1[idof_ensm4 + 4 * iface_ensm4 + 12 * idof_ensm4_0] * _pt_data_2[32 * iface_ensm4 + iel_ensm4] * _pt_temp_1[idof_ensm4 + 128 * iface_ensm4 + 4 * iel_ensm4]; | |
| acc_iface_ensm4_idof_ensm4 = acc_iface_ensm4_idof_ensm4 + _pt_data_1[idof_ensm4 + 4 * iface_ensm4 + 12 * idof_ensm4_0] * _pt_data_2[32 * iface_ensm4 + iel_ensm4] * _pt_temp[idof_ensm4 + 128 * iface_ensm4 + 4 * iel_ensm4]; | |
| } | |
| face_mass_0[idof_ensm4_0 + 10 * iel_ensm4] = acc_iface_ensm4_idof_ensm4_0; | |
| face_mass[idof_ensm4_0 + 10 * iel_ensm4] = acc_iface_ensm4_idof_ensm4; | |
| } | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 9; ++idof_ensm5) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| { | |
| mass_inv_inp_0[idof_ensm5 + 10 * iel_ensm5] = _pt_part_ph_id[idof_ensm5 + 10 * iel_ensm5] + -1.0 * face_mass_0[idof_ensm5 + 10 * iel_ensm5]; | |
| mass_inv_inp[idof_ensm5 + 10 * iel_ensm5] = _pt_part_ph_id[320 + idof_ensm5 + 10 * iel_ensm5] + -1.0 * face_mass[idof_ensm5 + 10 * iel_ensm5]; | |
| } | |
| for (int idof_ensm6 = 0; idof_ensm6 <= 9; ++idof_ensm6) | |
| for (int iel_ensm6 = 0; iel_ensm6 <= 31; ++iel_ensm6) | |
| { | |
| acc_idof_ensm6_0 = 0.0; | |
| acc_idof_ensm6_1 = 0.0; | |
| for (int idof_ensm6_0 = 0; idof_ensm6_0 <= 9; ++idof_ensm6_0) | |
| { | |
| acc_idof_ensm6_0 = acc_idof_ensm6_0 + _pt_data[iel_ensm6] * _pt_data_0[idof_ensm6_0 + 10 * idof_ensm6] * mass_inv_inp[idof_ensm6_0 + 10 * iel_ensm6]; | |
| acc_idof_ensm6_1 = acc_idof_ensm6_1 + _pt_data[iel_ensm6] * _pt_data_0[idof_ensm6_0 + 10 * idof_ensm6] * mass_inv_inp_0[idof_ensm6_0 + 10 * iel_ensm6]; | |
| } | |
| _pt_temp_0[idof_ensm6 + 10 * iel_ensm6] = acc_idof_ensm6_0; | |
| _pt_temp_2[idof_ensm6 + 10 * iel_ensm6] = acc_idof_ensm6_1; | |
| } | |
| for (int idof_ensm7 = 0; idof_ensm7 <= 9; ++idof_ensm7) | |
| for (int iel_ensm7 = 0; iel_ensm7 <= 31; ++iel_ensm7) | |
| { | |
| _pt_out_0_0[idof_ensm7 + 10 * iel_ensm7] = -1.0 * _pt_temp_2[idof_ensm7 + 10 * iel_ensm7]; | |
| _pt_out_1_0[idof_ensm7 + 10 * iel_ensm7] = -1.0 * _pt_temp_0[idof_ensm7 + 10 * iel_ensm7]; | |
| } | |
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs(double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_1, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, double const *__restrict__ normal_2_b_BTAG_PARTITION_1, int const *__restrict__ from_el_indices_5, int const *__restrict__ _pt_data_9, char const *__restrict__ from_el_present_3, double const *__restrict__ _pt_part_ph_id_4, double const *__restrict__ _pt_dist_id_2, int const *__restrict__ from_el_indices_6, int const *__restrict__ _pt_data_10, double const *__restrict__ normal_2_b_BTAG_PARTITION_2, int const *__restrict__ from_el_indices_7, int const *__restrict__ _pt_data_11, char const *__restrict__ from_el_present_4, double const *__restrict__ _pt_part_ph_id_5, double const *__restrict__ _pt_dist_id_3, int const *__restrict__ from_el_indices_8, int const *__restrict__ _pt_data_12, double const *__restrict__ normal_2_b_BTAG_PARTITION_3, int const *__restrict__ from_el_indices_9, int const *__restrict__ _pt_data_13, char const *__restrict__ from_el_present_5, double const *__restrict__ _pt_part_ph_id_6, double const *__restrict__ _pt_dist_id_4, int const *__restrict__ from_el_indices_10, int const *__restrict__ _pt_data_14, double const *__restrict__ normal_2_b_BTAG_PARTITION_4, int const *__restrict__ from_el_indices_11, int const *__restrict__ _pt_data_15, char const *__restrict__ from_el_present_6, double const *__restrict__ _pt_part_ph_id_7, double const *__restrict__ _pt_dist_id_5, int const *__restrict__ from_el_indices_12, int const *__restrict__ _pt_data_16, double const *__restrict__ normal_2_b_BTAG_PARTITION_5, int const *__restrict__ from_el_indices_13, int const *__restrict__ _pt_data_17, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double const *__restrict__ normal_1_b_BTAG_PARTITION_1, double const *__restrict__ normal_1_b_BTAG_PARTITION_2, double const *__restrict__ normal_1_b_BTAG_PARTITION_3, double const *__restrict__ normal_1_b_BTAG_PARTITION_4, double const *__restrict__ normal_1_b_BTAG_PARTITION_5, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ cse_1, double *__restrict__ cse_2, double *__restrict__ cse_3, double *__restrict__ cse_4, double *__restrict__ cse_5, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm9_0; | |
| double acc_idof_ensm9_1; | |
| double acc_iface_ensm6_idof_ensm6; | |
| double acc_iface_ensm6_idof_ensm6_0; | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| { | |
| int const iel_ensm0 = 0; | |
| cse_4[idof_ensm0] = 0.5 * (_pt_part_ph_id_6[idof_ensm0] + _pt_dist_id_4[4 * from_el_indices_10[0] + _pt_data_14[idof_ensm0]]); | |
| cse_5[idof_ensm0] = 0.5 * (_pt_part_ph_id_7[idof_ensm0] + _pt_dist_id_5[4 * from_el_indices_12[0] + _pt_data_16[idof_ensm0]]); | |
| } | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| for (int iel_ensm1 = 0; iel_ensm1 <= 3; ++iel_ensm1) | |
| cse_1[idof_ensm1 + 4 * iel_ensm1] = 0.5 * (_pt_part_ph_id_3[idof_ensm1 + 4 * iel_ensm1] + _pt_dist_id_1[4 * from_el_indices_4[iel_ensm1] + _pt_data_8[idof_ensm1 + 4 * iel_ensm1]]); | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| for (int iel_ensm2 = 0; iel_ensm2 <= 5; ++iel_ensm2) | |
| cse_3[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_5[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id_3[4 * from_el_indices_8[iel_ensm2] + _pt_data_12[idof_ensm2 + 4 * iel_ensm2]]); | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 3; ++idof_ensm3) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 1; ++iel_ensm3) | |
| cse[idof_ensm3 + 4 * iel_ensm3] = 0.5 * (_pt_part_ph_id_1[idof_ensm3 + 4 * iel_ensm3] + _pt_dist_id[4 * from_el_indices_0[iel_ensm3] + _pt_data_4[idof_ensm3 + 4 * iel_ensm3]]); | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 3; ++idof_ensm4) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 2; ++iel_ensm4) | |
| { | |
| cse_0[idof_ensm4 + 4 * iel_ensm4] = 0.5 * (_pt_part_ph_id_2[idof_ensm4 + 4 * iel_ensm4] + _pt_dist_id_0[4 * from_el_indices_2[iel_ensm4] + _pt_data_6[idof_ensm4 + 4 * iel_ensm4]]); | |
| cse_2[idof_ensm4 + 4 * iel_ensm4] = 0.5 * (_pt_part_ph_id_4[idof_ensm4 + 4 * iel_ensm4] + _pt_dist_id_2[4 * from_el_indices_6[iel_ensm4] + _pt_data_10[idof_ensm4 + 4 * iel_ensm4]]); | |
| } | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 3; ++idof_ensm5) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| for (int iface_ensm5 = 0; iface_ensm5 <= 2; ++iface_ensm5) | |
| { | |
| _pt_temp_1[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5] = | |
| (from_el_present_6[iel_ensm5 + 32 * iface_ensm5] ? | |
| cse_5[4 * from_el_indices_13[iel_ensm5 + 32 * iface_ensm5] + _pt_data_17[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_1_b_BTAG_PARTITION_5[from_el_indices_13[iel_ensm5 + 32 * iface_ensm5]] | |
| : 0.0) | |
| + (from_el_present_5[iel_ensm5 + 32 * iface_ensm5] ? | |
| cse_4[4 * from_el_indices_11[iel_ensm5 + 32 * iface_ensm5] + _pt_data_15[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_1_b_BTAG_PARTITION_4[from_el_indices_11[iel_ensm5 + 32 * iface_ensm5]] | |
| : 0.0) | |
| + (from_el_present_4[iel_ensm5 + 32 * iface_ensm5] ? | |
| cse_3[4 * from_el_indices_9[iel_ensm5 + 32 * iface_ensm5] + _pt_data_13[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_1_b_BTAG_PARTITION_3[from_el_indices_9[iel_ensm5 + 32 * iface_ensm5]] | |
| : 0.0) | |
| + (from_el_present_3[iel_ensm5 + 32 * iface_ensm5] ? | |
| cse_2[4 * from_el_indices_7[iel_ensm5 + 32 * iface_ensm5] + _pt_data_11[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_1_b_BTAG_PARTITION_2[from_el_indices_7[iel_ensm5 + 32 * iface_ensm5]] | |
| : 0.0) | |
| + (from_el_present_2[iel_ensm5 + 32 * iface_ensm5] ? | |
| cse_1[4 * from_el_indices_5[iel_ensm5 + 32 * iface_ensm5] + _pt_data_9[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_1_b_BTAG_PARTITION_1[from_el_indices_5[iel_ensm5 + 32 * iface_ensm5]] | |
| : 0.0) | |
| + (from_el_present_1[iel_ensm5 + 32 * iface_ensm5] ? | |
| cse_0[4 * from_el_indices_3[iel_ensm5 + 32 * iface_ensm5] + _pt_data_7[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_3[iel_ensm5 + 32 * iface_ensm5]] | |
| : 0.0) | |
| + (from_el_present[iel_ensm5 + 32 * iface_ensm5] ? | |
| _pt_part_ph_id_0[4 * from_el_indices[iel_ensm5 + 32 * iface_ensm5] + _pt_data_3[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_1_b_face_restr_interior[from_el_indices[iel_ensm5 + 32 * iface_ensm5]] | |
| : 0.0) | |
| + (from_el_present_0[iel_ensm5 + 32 * iface_ensm5] ? | |
| cse[4 * from_el_indices_1[iel_ensm5 + 32 * iface_ensm5] + _pt_data_5[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_1_b_BTAG_PARTITION[from_el_indices_1[iel_ensm5 + 32 * iface_ensm5]] | |
| : 0.0); | |
| _pt_temp[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5] = (from_el_present_6[iel_ensm5 + 32 * iface_ensm5] ? cse_5[4 * from_el_indices_13[iel_ensm5 + 32 * iface_ensm5] + _pt_data_17[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_2_b_BTAG_PARTITION_5[from_el_indices_13[iel_ensm5 + 32 * iface_ensm5]] : 0.0) + (from_el_present_5[iel_ensm5 + 32 * iface_ensm5] ? cse_4[4 * from_el_indices_11[iel_ensm5 + 32 * iface_ensm5] + _pt_data_15[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_2_b_BTAG_PARTITION_4[from_el_indices_11[iel_ensm5 + 32 * iface_ensm5]] : 0.0) + (from_el_present_4[iel_ensm5 + 32 * iface_ensm5] ? cse_3[4 * from_el_indices_9[iel_ensm5 + 32 * iface_ensm5] + _pt_data_13[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_2_b_BTAG_PARTITION_3[from_el_indices_9[iel_ensm5 + 32 * iface_ensm5]] : 0.0) + (from_el_present_3[iel_ensm5 + 32 * iface_ensm5] ? cse_2[4 * from_el_indices_7[iel_ensm5 + 32 * iface_ensm5] + _pt_data_11[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_2_b_BTAG_PARTITION_2[from_el_indices_7[iel_ensm5 + 32 * iface_ensm5]] : 0.0) + (from_el_present_2[iel_ensm5 + 32 * iface_ensm5] ? cse_1[4 * from_el_indices_5[iel_ensm5 + 32 * iface_ensm5] + _pt_data_9[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_2_b_BTAG_PARTITION_1[from_el_indices_5[iel_ensm5 + 32 * iface_ensm5]] : 0.0) + (from_el_present_1[iel_ensm5 + 32 * iface_ensm5] ? cse_0[4 * from_el_indices_3[iel_ensm5 + 32 * iface_ensm5] + _pt_data_7[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_3[iel_ensm5 + 32 * iface_ensm5]] : 0.0) + (from_el_present[iel_ensm5 + 32 * iface_ensm5] ? _pt_part_ph_id_0[4 * from_el_indices[iel_ensm5 + 32 * iface_ensm5] + _pt_data_3[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_2_b_face_restr_interior[from_el_indices[iel_ensm5 + 32 * iface_ensm5]] : 0.0) + (from_el_present_0[iel_ensm5 + 32 * iface_ensm5] ? cse[4 * from_el_indices_1[iel_ensm5 + 32 * iface_ensm5] + _pt_data_5[idof_ensm5 + 4 * iel_ensm5 + 128 * iface_ensm5]] * normal_2_b_BTAG_PARTITION[from_el_indices_1[iel_ensm5 + 32 * iface_ensm5]] : 0.0); | |
| } | |
| for (int idof_ensm6_0 = 0; idof_ensm6_0 <= 9; ++idof_ensm6_0) | |
| for (int iel_ensm6 = 0; iel_ensm6 <= 31; ++iel_ensm6) | |
| { | |
| acc_iface_ensm6_idof_ensm6_0 = 0.0; | |
| acc_iface_ensm6_idof_ensm6 = 0.0; | |
| for (int idof_ensm6 = 0; idof_ensm6 <= 3; ++idof_ensm6) | |
| for (int iface_ensm6 = 0; iface_ensm6 <= 2; ++iface_ensm6) | |
| { | |
| acc_iface_ensm6_idof_ensm6_0 = acc_iface_ensm6_idof_ensm6_0 + _pt_data_1[idof_ensm6 + 4 * iface_ensm6 + 12 * idof_ensm6_0] * _pt_data_2[32 * iface_ensm6 + iel_ensm6] * _pt_temp_1[idof_ensm6 + 128 * iface_ensm6 + 4 * iel_ensm6]; | |
| acc_iface_ensm6_idof_ensm6 = acc_iface_ensm6_idof_ensm6 + _pt_data_1[idof_ensm6 + 4 * iface_ensm6 + 12 * idof_ensm6_0] * _pt_data_2[32 * iface_ensm6 + iel_ensm6] * _pt_temp[idof_ensm6 + 128 * iface_ensm6 + 4 * iel_ensm6]; | |
| } | |
| face_mass_0[idof_ensm6_0 + 10 * iel_ensm6] = acc_iface_ensm6_idof_ensm6_0; | |
| face_mass[idof_ensm6_0 + 10 * iel_ensm6] = acc_iface_ensm6_idof_ensm6; | |
| } | |
| for (int idof_ensm8 = 0; idof_ensm8 <= 9; ++idof_ensm8) | |
| for (int iel_ensm8 = 0; iel_ensm8 <= 31; ++iel_ensm8) | |
| { | |
| mass_inv_inp_0[idof_ensm8 + 10 * iel_ensm8] = _pt_part_ph_id[idof_ensm8 + 10 * iel_ensm8] + -1.0 * face_mass_0[idof_ensm8 + 10 * iel_ensm8]; | |
| mass_inv_inp[idof_ensm8 + 10 * iel_ensm8] = _pt_part_ph_id[320 + idof_ensm8 + 10 * iel_ensm8] + -1.0 * face_mass[idof_ensm8 + 10 * iel_ensm8]; | |
| } | |
| for (int idof_ensm9 = 0; idof_ensm9 <= 9; ++idof_ensm9) | |
| for (int iel_ensm9 = 0; iel_ensm9 <= 31; ++iel_ensm9) | |
| { | |
| acc_idof_ensm9_0 = 0.0; | |
| acc_idof_ensm9_1 = 0.0; | |
| for (int idof_ensm9_0 = 0; idof_ensm9_0 <= 9; ++idof_ensm9_0) | |
| { | |
| acc_idof_ensm9_0 = acc_idof_ensm9_0 + _pt_data[iel_ensm9] * _pt_data_0[idof_ensm9_0 + 10 * idof_ensm9] * mass_inv_inp[idof_ensm9_0 + 10 * iel_ensm9]; | |
| acc_idof_ensm9_1 = acc_idof_ensm9_1 + _pt_data[iel_ensm9] * _pt_data_0[idof_ensm9_0 + 10 * idof_ensm9] * mass_inv_inp_0[idof_ensm9_0 + 10 * iel_ensm9]; | |
| } | |
| _pt_temp_0[idof_ensm9 + 10 * iel_ensm9] = acc_idof_ensm9_0; | |
| _pt_temp_2[idof_ensm9 + 10 * iel_ensm9] = acc_idof_ensm9_1; | |
| } | |
| for (int idof_ensm7 = 0; idof_ensm7 <= 9; ++idof_ensm7) | |
| for (int iel_ensm7 = 0; iel_ensm7 <= 31; ++iel_ensm7) | |
| { | |
| _pt_out_0_0[idof_ensm7 + 10 * iel_ensm7] = -1.0 * _pt_temp_2[idof_ensm7 + 10 * iel_ensm7]; | |
| _pt_out_1_0[idof_ensm7 + 10 * iel_ensm7] = -1.0 * _pt_temp_0[idof_ensm7 + 10 * iel_ensm7]; | |
| } | |
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs(double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_all, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_all, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm5_0; | |
| double acc_idof_ensm5_1; | |
| double acc_iface_ensm2_idof_ensm2; | |
| double acc_iface_ensm2_idof_ensm2_0; | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| for (int iel_ensm1 = 0; iel_ensm1 <= 4; ++iel_ensm1) | |
| { | |
| cse_0[idof_ensm1 + 4 * iel_ensm1] = 0.5 * (_pt_part_ph_id_3[idof_ensm1 + 4 * iel_ensm1] + _pt_dist_id_0[4 * from_el_indices_3[iel_ensm1] + _pt_data_7[idof_ensm1 + 4 * iel_ensm1]]); | |
| cse[idof_ensm1 + 4 * iel_ensm1] = 0.5 * (_pt_part_ph_id_2[idof_ensm1 + 4 * iel_ensm1] + _pt_dist_id[4 * from_el_indices_1[iel_ensm1] + _pt_data_5[idof_ensm1 + 4 * iel_ensm1]]); | |
| } | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| for (int iel_ensm0 = 0; iel_ensm0 <= 31; ++iel_ensm0) | |
| for (int iface_ensm0 = 0; iface_ensm0 <= 2; ++iface_ensm0) | |
| { | |
| _pt_temp_1[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0] = (from_el_present[iel_ensm0 + 32 * iface_ensm0] ? normal_1_b_all[from_el_indices[iel_ensm0 + 32 * iface_ensm0]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]]) : 0.0) + (from_el_present_2[iel_ensm0 + 32 * iface_ensm0] ? cse_0[4 * from_el_indices_4[iel_ensm0 + 32 * iface_ensm0] + _pt_data_8[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_0[iel_ensm0 + 32 * iface_ensm0] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm0 + 32 * iface_ensm0] + _pt_data_4[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_face_restr_interior[from_el_indices_0[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_1[iel_ensm0 + 32 * iface_ensm0] ? cse[4 * from_el_indices_2[iel_ensm0 + 32 * iface_ensm0] + _pt_data_6[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_1_b_BTAG_PARTITION[from_el_indices_2[iel_ensm0 + 32 * iface_ensm0]] : 0.0); | |
| _pt_temp[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0] = (from_el_present[iel_ensm0 + 32 * iface_ensm0] ? normal_2_b_all[from_el_indices[iel_ensm0 + 32 * iface_ensm0]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm0 + 32 * iface_ensm0] + _pt_data_3[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]]) : 0.0) + (from_el_present_2[iel_ensm0 + 32 * iface_ensm0] ? cse_0[4 * from_el_indices_4[iel_ensm0 + 32 * iface_ensm0] + _pt_data_8[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_0[iel_ensm0 + 32 * iface_ensm0] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm0 + 32 * iface_ensm0] + _pt_data_4[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_face_restr_interior[from_el_indices_0[iel_ensm0 + 32 * iface_ensm0]] : 0.0) + (from_el_present_1[iel_ensm0 + 32 * iface_ensm0] ? cse[4 * from_el_indices_2[iel_ensm0 + 32 * iface_ensm0] + _pt_data_6[idof_ensm0 + 4 * iel_ensm0 + 128 * iface_ensm0]] * normal_2_b_BTAG_PARTITION[from_el_indices_2[iel_ensm0 + 32 * iface_ensm0]] : 0.0); | |
| } | |
| for (int idof_ensm2_0 = 0; idof_ensm2_0 <= 9; ++idof_ensm2_0) | |
| for (int iel_ensm2 = 0; iel_ensm2 <= 31; ++iel_ensm2) | |
| { | |
| acc_iface_ensm2_idof_ensm2_0 = 0.0; | |
| acc_iface_ensm2_idof_ensm2 = 0.0; | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| for (int iface_ensm2 = 0; iface_ensm2 <= 2; ++iface_ensm2) | |
| { | |
| acc_iface_ensm2_idof_ensm2_0 = acc_iface_ensm2_idof_ensm2_0 + _pt_data_1[idof_ensm2 + 4 * iface_ensm2 + 12 * idof_ensm2_0] * _pt_data_2[32 * iface_ensm2 + iel_ensm2] * _pt_temp_1[idof_ensm2 + 128 * iface_ensm2 + 4 * iel_ensm2]; | |
| acc_iface_ensm2_idof_ensm2 = acc_iface_ensm2_idof_ensm2 + _pt_data_1[idof_ensm2 + 4 * iface_ensm2 + 12 * idof_ensm2_0] * _pt_data_2[32 * iface_ensm2 + iel_ensm2] * _pt_temp[idof_ensm2 + 128 * iface_ensm2 + 4 * iel_ensm2]; | |
| } | |
| face_mass_0[idof_ensm2_0 + 10 * iel_ensm2] = acc_iface_ensm2_idof_ensm2_0; | |
| face_mass[idof_ensm2_0 + 10 * iel_ensm2] = acc_iface_ensm2_idof_ensm2; | |
| } | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 9; ++idof_ensm3) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 31; ++iel_ensm3) | |
| { | |
| mass_inv_inp_0[idof_ensm3 + 10 * iel_ensm3] = _pt_part_ph_id[idof_ensm3 + 10 * iel_ensm3] + -1.0 * face_mass_0[idof_ensm3 + 10 * iel_ensm3]; | |
| mass_inv_inp[idof_ensm3 + 10 * iel_ensm3] = _pt_part_ph_id[320 + idof_ensm3 + 10 * iel_ensm3] + -1.0 * face_mass[idof_ensm3 + 10 * iel_ensm3]; | |
| } | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 9; ++idof_ensm5) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| { | |
| acc_idof_ensm5_0 = 0.0; | |
| acc_idof_ensm5_1 = 0.0; | |
| for (int idof_ensm5_0 = 0; idof_ensm5_0 <= 9; ++idof_ensm5_0) | |
| { | |
| acc_idof_ensm5_0 = acc_idof_ensm5_0 + _pt_data[iel_ensm5] * _pt_data_0[idof_ensm5_0 + 10 * idof_ensm5] * mass_inv_inp[idof_ensm5_0 + 10 * iel_ensm5]; | |
| acc_idof_ensm5_1 = acc_idof_ensm5_1 + _pt_data[iel_ensm5] * _pt_data_0[idof_ensm5_0 + 10 * idof_ensm5] * mass_inv_inp_0[idof_ensm5_0 + 10 * iel_ensm5]; | |
| } | |
| _pt_temp_0[idof_ensm5 + 10 * iel_ensm5] = acc_idof_ensm5_0; | |
| _pt_temp_2[idof_ensm5 + 10 * iel_ensm5] = acc_idof_ensm5_1; | |
| } | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 9; ++idof_ensm4) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 31; ++iel_ensm4) | |
| { | |
| _pt_out_0_0[idof_ensm4 + 10 * iel_ensm4] = -1.0 * _pt_temp_2[idof_ensm4 + 10 * iel_ensm4]; | |
| _pt_out_1_0[idof_ensm4 + 10 * iel_ensm4] = -1.0 * _pt_temp_0[idof_ensm4 + 10 * iel_ensm4]; | |
| } | |
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| void rhs(double const *__restrict__ _pt_data, double const *__restrict__ _pt_data_0, double const *__restrict__ _pt_part_ph_id, double const *__restrict__ _pt_data_1, double const *__restrict__ _pt_data_2, char const *__restrict__ from_el_present, double const *__restrict__ _pt_part_ph_id_0, double const *__restrict__ normal_2_b_all, int const *__restrict__ from_el_indices, int const *__restrict__ _pt_data_3, char const *__restrict__ from_el_present_0, double const *__restrict__ _pt_part_ph_id_1, double const *__restrict__ normal_2_b_face_restr_interior, int const *__restrict__ from_el_indices_0, int const *__restrict__ _pt_data_4, char const *__restrict__ from_el_present_1, double const *__restrict__ _pt_part_ph_id_2, double const *__restrict__ _pt_dist_id, int const *__restrict__ from_el_indices_1, int const *__restrict__ _pt_data_5, double const *__restrict__ normal_2_b_BTAG_PARTITION, int const *__restrict__ from_el_indices_2, int const *__restrict__ _pt_data_6, char const *__restrict__ from_el_present_2, double const *__restrict__ _pt_part_ph_id_3, double const *__restrict__ _pt_dist_id_0, int const *__restrict__ from_el_indices_3, int const *__restrict__ _pt_data_7, double const *__restrict__ normal_2_b_BTAG_PARTITION_0, int const *__restrict__ from_el_indices_4, int const *__restrict__ _pt_data_8, char const *__restrict__ from_el_present_3, double const *__restrict__ _pt_part_ph_id_4, double const *__restrict__ _pt_dist_id_1, int const *__restrict__ from_el_indices_5, int const *__restrict__ _pt_data_9, double const *__restrict__ normal_2_b_BTAG_PARTITION_1, int const *__restrict__ from_el_indices_6, int const *__restrict__ _pt_data_10, char const *__restrict__ from_el_present_4, double const *__restrict__ _pt_part_ph_id_5, double const *__restrict__ _pt_dist_id_2, int const *__restrict__ from_el_indices_7, int const *__restrict__ _pt_data_11, double const *__restrict__ normal_2_b_BTAG_PARTITION_2, int const *__restrict__ from_el_indices_8, int const *__restrict__ _pt_data_12, double *__restrict__ _pt_out_1_0, double const *__restrict__ normal_1_b_all, double const *__restrict__ normal_1_b_face_restr_interior, double const *__restrict__ normal_1_b_BTAG_PARTITION, double const *__restrict__ normal_1_b_BTAG_PARTITION_0, double const *__restrict__ normal_1_b_BTAG_PARTITION_1, double const *__restrict__ normal_1_b_BTAG_PARTITION_2, double *__restrict__ _pt_out_0_0, double *__restrict__ _pt_temp, double *__restrict__ _pt_temp_0, double *__restrict__ _pt_temp_1, double *__restrict__ _pt_temp_2, double *__restrict__ cse, double *__restrict__ cse_0, double *__restrict__ cse_1, double *__restrict__ cse_2, double *__restrict__ face_mass, double *__restrict__ face_mass_0, double *__restrict__ mass_inv_inp, double *__restrict__ mass_inv_inp_0) | |
| { | |
| double acc_idof_ensm6_0; | |
| double acc_idof_ensm6_1; | |
| double acc_iface_ensm5_idof_ensm5; | |
| double acc_iface_ensm5_idof_ensm5_0; | |
| for (int idof_ensm0 = 0; idof_ensm0 <= 3; ++idof_ensm0) | |
| for (int iel_ensm0 = 0; iel_ensm0 <= 6; ++iel_ensm0) | |
| cse_0[idof_ensm0 + 4 * iel_ensm0] = 0.5 * (_pt_part_ph_id_3[idof_ensm0 + 4 * iel_ensm0] + _pt_dist_id_0[4 * from_el_indices_3[iel_ensm0] + _pt_data_7[idof_ensm0 + 4 * iel_ensm0]]); | |
| for (int idof_ensm1 = 0; idof_ensm1 <= 3; ++idof_ensm1) | |
| { | |
| int const iel_ensm1 = 0; | |
| cse_1[idof_ensm1] = 0.5 * (_pt_part_ph_id_4[idof_ensm1] + _pt_dist_id_1[4 * from_el_indices_5[0] + _pt_data_9[idof_ensm1]]); | |
| } | |
| for (int idof_ensm2 = 0; idof_ensm2 <= 3; ++idof_ensm2) | |
| for (int iel_ensm2 = 0; iel_ensm2 <= 2; ++iel_ensm2) | |
| cse_2[idof_ensm2 + 4 * iel_ensm2] = 0.5 * (_pt_part_ph_id_5[idof_ensm2 + 4 * iel_ensm2] + _pt_dist_id_2[4 * from_el_indices_7[iel_ensm2] + _pt_data_11[idof_ensm2 + 4 * iel_ensm2]]); | |
| for (int idof_ensm3 = 0; idof_ensm3 <= 3; ++idof_ensm3) | |
| for (int iel_ensm3 = 0; iel_ensm3 <= 4; ++iel_ensm3) | |
| cse[idof_ensm3 + 4 * iel_ensm3] = 0.5 * (_pt_part_ph_id_2[idof_ensm3 + 4 * iel_ensm3] + _pt_dist_id[4 * from_el_indices_1[iel_ensm3] + _pt_data_5[idof_ensm3 + 4 * iel_ensm3]]); | |
| for (int idof_ensm4 = 0; idof_ensm4 <= 3; ++idof_ensm4) | |
| for (int iel_ensm4 = 0; iel_ensm4 <= 31; ++iel_ensm4) | |
| for (int iface_ensm4 = 0; iface_ensm4 <= 2; ++iface_ensm4) | |
| { | |
| _pt_temp_1[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4] = (from_el_present[iel_ensm4 + 32 * iface_ensm4] ? normal_1_b_all[from_el_indices[iel_ensm4 + 32 * iface_ensm4]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm4 + 32 * iface_ensm4] + _pt_data_3[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm4 + 32 * iface_ensm4] + _pt_data_3[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4]]) : 0.0) + (from_el_present_4[iel_ensm4 + 32 * iface_ensm4] ? cse_2[4 * from_el_indices_8[iel_ensm4 + 32 * iface_ensm4] + _pt_data_12[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4]] * normal_1_b_BTAG_PARTITION_2[from_el_indices_8[iel_ensm4 + 32 * iface_ensm4]] : 0.0) + (from_el_present_3[iel_ensm4 + 32 * iface_ensm4] ? cse_1[4 * from_el_indices_6[iel_ensm4 + 32 * iface_ensm4] + _pt_data_10[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4]] * normal_1_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm4 + 32 * iface_ensm4]] : 0.0) + (from_el_present_2[iel_ensm4 + 32 * iface_ensm4] ? cse_0[4 * from_el_indices_4[iel_ensm4 + 32 * iface_ensm4] + _pt_data_8[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4]] * normal_1_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm4 + 32 * iface_ensm4]] : 0.0) + (from_el_present_0[iel_ensm4 + 32 * iface_ensm4] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm4 + 32 * iface_ensm4] + _pt_data_4[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4]] * normal_1_b_face_restr_interior[from_el_indices_0[iel_ensm4 + 32 * iface_ensm4]] : 0.0) + (from_el_present_1[iel_ensm4 + 32 * iface_ensm4] ? cse[4 * from_el_indices_2[iel_ensm4 + 32 * iface_ensm4] + _pt_data_6[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4]] * normal_1_b_BTAG_PARTITION[from_el_indices_2[iel_ensm4 + 32 * iface_ensm4]] : 0.0); | |
| _pt_temp[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4] = (from_el_present[iel_ensm4 + 32 * iface_ensm4] ? normal_2_b_all[from_el_indices[iel_ensm4 + 32 * iface_ensm4]] * 0.5 * (_pt_part_ph_id_0[4 * from_el_indices[iel_ensm4 + 32 * iface_ensm4] + _pt_data_3[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4]] + -1.0 * _pt_part_ph_id_0[4 * from_el_indices[iel_ensm4 + 32 * iface_ensm4] + _pt_data_3[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4]]) : 0.0) + (from_el_present_4[iel_ensm4 + 32 * iface_ensm4] ? cse_2[4 * from_el_indices_8[iel_ensm4 + 32 * iface_ensm4] + _pt_data_12[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4]] * normal_2_b_BTAG_PARTITION_2[from_el_indices_8[iel_ensm4 + 32 * iface_ensm4]] : 0.0) + (from_el_present_3[iel_ensm4 + 32 * iface_ensm4] ? cse_1[4 * from_el_indices_6[iel_ensm4 + 32 * iface_ensm4] + _pt_data_10[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4]] * normal_2_b_BTAG_PARTITION_1[from_el_indices_6[iel_ensm4 + 32 * iface_ensm4]] : 0.0) + (from_el_present_2[iel_ensm4 + 32 * iface_ensm4] ? cse_0[4 * from_el_indices_4[iel_ensm4 + 32 * iface_ensm4] + _pt_data_8[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4]] * normal_2_b_BTAG_PARTITION_0[from_el_indices_4[iel_ensm4 + 32 * iface_ensm4]] : 0.0) + (from_el_present_0[iel_ensm4 + 32 * iface_ensm4] ? _pt_part_ph_id_1[4 * from_el_indices_0[iel_ensm4 + 32 * iface_ensm4] + _pt_data_4[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4]] * normal_2_b_face_restr_interior[from_el_indices_0[iel_ensm4 + 32 * iface_ensm4]] : 0.0) + (from_el_present_1[iel_ensm4 + 32 * iface_ensm4] ? cse[4 * from_el_indices_2[iel_ensm4 + 32 * iface_ensm4] + _pt_data_6[idof_ensm4 + 4 * iel_ensm4 + 128 * iface_ensm4]] * normal_2_b_BTAG_PARTITION[from_el_indices_2[iel_ensm4 + 32 * iface_ensm4]] : 0.0); | |
| } | |
| for (int idof_ensm5_0 = 0; idof_ensm5_0 <= 9; ++idof_ensm5_0) | |
| for (int iel_ensm5 = 0; iel_ensm5 <= 31; ++iel_ensm5) | |
| { | |
| acc_iface_ensm5_idof_ensm5_0 = 0.0; | |
| acc_iface_ensm5_idof_ensm5 = 0.0; | |
| for (int idof_ensm5 = 0; idof_ensm5 <= 3; ++idof_ensm5) | |
| for (int iface_ensm5 = 0; iface_ensm5 <= 2; ++iface_ensm5) | |
| { | |
| acc_iface_ensm5_idof_ensm5_0 = acc_iface_ensm5_idof_ensm5_0 + _pt_data_1[idof_ensm5 + 4 * iface_ensm5 + 12 * idof_ensm5_0] * _pt_data_2[32 * iface_ensm5 + iel_ensm5] * _pt_temp_1[idof_ensm5 + 128 * iface_ensm5 + 4 * iel_ensm5]; | |
| acc_iface_ensm5_idof_ensm5 = acc_iface_ensm5_idof_ensm5 + _pt_data_1[idof_ensm5 + 4 * iface_ensm5 + 12 * idof_ensm5_0] * _pt_data_2[32 * iface_ensm5 + iel_ensm5] * _pt_temp[idof_ensm5 + 128 * iface_ensm5 + 4 * iel_ensm5]; | |
| } | |
| face_mass_0[idof_ensm5_0 + 10 * iel_ensm5] = acc_iface_ensm5_idof_ensm5_0; | |
| face_mass[idof_ensm5_0 + 10 * iel_ensm5] = acc_iface_ensm5_idof_ensm5; | |
| } | |
| for (int idof_ensm7 = 0; idof_ensm7 <= 9; ++idof_ensm7) | |
| for (int iel_ensm7 = 0; iel_ensm7 <= 31; ++iel_ensm7) | |
| { | |
| mass_inv_inp_0[idof_ensm7 + 10 * iel_ensm7] = _pt_part_ph_id[idof_ensm7 + 10 * iel_ensm7] + -1.0 * face_mass_0[idof_ensm7 + 10 * iel_ensm7]; | |
| mass_inv_inp[idof_ensm7 + 10 * iel_ensm7] = _pt_part_ph_id[320 + idof_ensm7 + 10 * iel_ensm7] + -1.0 * face_mass[idof_ensm7 + 10 * iel_ensm7]; | |
| } | |
| for (int idof_ensm6 = 0; idof_ensm6 <= 9; ++idof_ensm6) | |
| for (int iel_ensm6 = 0; iel_ensm6 <= 31; ++iel_ensm6) | |
| { | |
| acc_idof_ensm6_0 = 0.0; | |
| acc_idof_ensm6_1 = 0.0; | |
| for (int idof_ensm6_0 = 0; idof_ensm6_0 <= 9; ++idof_ensm6_0) | |
| { | |
| acc_idof_ensm6_0 = acc_idof_ensm6_0 + _pt_data[iel_ensm6] * _pt_data_0[idof_ensm6_0 + 10 * idof_ensm6] * mass_inv_inp[idof_ensm6_0 + 10 * iel_ensm6]; | |
| acc_idof_ensm6_1 = acc_idof_ensm6_1 + _pt_data[iel_ensm6] * _pt_data_0[idof_ensm6_0 + 10 * idof_ensm6] * mass_inv_inp_0[idof_ensm6_0 + 10 * iel_ensm6]; | |
| } | |
| _pt_temp_0[idof_ensm6 + 10 * iel_ensm6] = acc_idof_ensm6_0; | |
| _pt_temp_2[idof_ensm6 + 10 * iel_ensm6] = acc_idof_ensm6_1; | |
| } | |
| for (int idof_ensm8 = 0; idof_ensm8 <= 9; ++idof_ensm8) | |
| for (int iel_ensm8 = 0; iel_ensm8 <= 31; ++iel_ensm8) | |
| { | |
| _pt_out_0_0[idof_ensm8 + 10 * iel_ensm8] = -1.0 * _pt_temp_2[idof_ensm8 + 10 * iel_ensm8]; | |
| _pt_out_1_0[idof_ensm8 + 10 * iel_ensm8] = -1.0 * _pt_temp_0[idof_ensm8 + 10 * iel_ensm8]; | |
| } | |
| } |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment