|
| | EpsilonDivDiv_kernel_gen.local_subdomain_id |
| |
| | EpsilonDivDiv_kernel_gen.x_cell |
| |
| | EpsilonDivDiv_kernel_gen.y_cell |
| |
| | EpsilonDivDiv_kernel_gen.r_cell |
| |
| | EpsilonDivDiv_kernel_gen.integer |
| |
| str | EpsilonDivDiv_kernel_gen.scalar_type = 'double' |
| |
| | EpsilonDivDiv_kernel_gen.rads |
| |
| | EpsilonDivDiv_kernel_gen.rad_assignments |
| |
| | EpsilonDivDiv_kernel_gen.src_symbol |
| |
| | EpsilonDivDiv_kernel_gen.src_array_declaration |
| |
| | EpsilonDivDiv_kernel_gen.src_assignments |
| |
| | EpsilonDivDiv_kernel_gen.k_symbol |
| |
| | EpsilonDivDiv_kernel_gen.k_array_declaration |
| |
| | EpsilonDivDiv_kernel_gen.k_assignments |
| |
| str | EpsilonDivDiv_kernel_gen.kernel_code = "\n" |
| |
| | EpsilonDivDiv_kernel_gen.qp_data = quad_points_1x1 |
| |
| | EpsilonDivDiv_kernel_gen.qw_data = quad_weights_1x1 |
| |
| int | EpsilonDivDiv_kernel_gen.num_qps = 1 |
| |
| int | EpsilonDivDiv_kernel_gen.dim = 3 |
| |
| str | EpsilonDivDiv_kernel_gen.qp_array_name = "qp_array" |
| |
| str | EpsilonDivDiv_kernel_gen.qw_array_name = "qw_array" |
| |
| | EpsilonDivDiv_kernel_gen.qp_array = IndexedBase(qp_array_name, shape=(num_qps, 3), real=True) |
| |
| | EpsilonDivDiv_kernel_gen.qw_array = IndexedBase(qw_array_name, shape=(num_qps), real=True) |
| |
| int | EpsilonDivDiv_kernel_gen.num_wedges_per_hex_cell = 2 |
| |
| int | EpsilonDivDiv_kernel_gen.num_nodes_per_wedge_surface = 3 |
| |
| int | EpsilonDivDiv_kernel_gen.num_nodes_per_wedge = 6 |
| |
| | EpsilonDivDiv_kernel_gen.cmb_shift |
| |
| | EpsilonDivDiv_kernel_gen.surface_shift |
| |
| | EpsilonDivDiv_kernel_gen.at_surface_boundary |
| |
| | EpsilonDivDiv_kernel_gen.at_cmb_boundary |
| |
| | EpsilonDivDiv_kernel_gen.max_rad |
| |
| | EpsilonDivDiv_kernel_gen.treat_boundary |
| |
| | EpsilonDivDiv_kernel_gen.diagonal |
| |
| | EpsilonDivDiv_kernel_gen.postloop |
| |
| str | EpsilonDivDiv_kernel_gen.dst_array_name = "dst_array" |
| |
| | EpsilonDivDiv_kernel_gen.dst_symbol |
| |
| list | EpsilonDivDiv_kernel_gen.quadloop_body = [] |
| |
| list | EpsilonDivDiv_kernel_gen.quadloop_exprs = [] |
| |
| | EpsilonDivDiv_kernel_gen.q_symbol = sp.symbols("q", integer=True) |
| |
| list | EpsilonDivDiv_kernel_gen.qp = [qp_array[q_symbol, d] for d in range(dim)] |
| |
| | EpsilonDivDiv_kernel_gen.qw = qw_array[q_symbol] |
| |
| | EpsilonDivDiv_kernel_gen.w_symbol |
| |
| | EpsilonDivDiv_kernel_gen.dimi_symbol |
| |
| | EpsilonDivDiv_kernel_gen.dimj_symbol |
| |
| | EpsilonDivDiv_kernel_gen.dim_diagBC_symbol |
| |
| | EpsilonDivDiv_kernel_gen.k_eval_symbol = sp.symbols("k_eval", real=True) |
| |
| | EpsilonDivDiv_kernel_gen.k_eval |
| |
| | EpsilonDivDiv_kernel_gen.k_eval_replacements |
| |
| | EpsilonDivDiv_kernel_gen.k_eval_reduced_exprs |
| |
| | EpsilonDivDiv_kernel_gen.exprs |
| |
| | EpsilonDivDiv_kernel_gen.symbols |
| |
| | EpsilonDivDiv_kernel_gen.g_symbol = symbols("node_idx", integer=True) |
| |
| str | EpsilonDivDiv_kernel_gen.scalar_grad_name = "scalar_grad" |
| |
| | EpsilonDivDiv_kernel_gen.scalar_grad = IndexedBase(scalar_grad_name, shape=(num_nodes_per_wedge, dim), real=True) |
| |
| bool | EpsilonDivDiv_kernel_gen.jac_laterally_precomputed = False |
| |
| list | EpsilonDivDiv_kernel_gen.jac_exprs = [] |
| |
| | EpsilonDivDiv_kernel_gen.J |
| |
| | EpsilonDivDiv_kernel_gen.J_cse_assignments |
| |
| | EpsilonDivDiv_kernel_gen.J_cse_exprs |
| |
| | EpsilonDivDiv_kernel_gen.J_cse = sp.Matrix(3, 3, J_cse_exprs) |
| |
| | EpsilonDivDiv_kernel_gen.J_cse_replaced_assignments |
| |
| | EpsilonDivDiv_kernel_gen.J_cse_replaced |
| |
| | EpsilonDivDiv_kernel_gen.J_det = J_cse_replaced.det() |
| |
| | EpsilonDivDiv_kernel_gen.J_det_replacements |
| |
| | EpsilonDivDiv_kernel_gen.J_det_reduced_exprs |
| |
| | EpsilonDivDiv_kernel_gen.J_det_symbol = sp.symbols(f"J_det", real=True) |
| |
| | EpsilonDivDiv_kernel_gen.J_abs_det = abs(J_det_symbol) |
| |
| | EpsilonDivDiv_kernel_gen.J_inv = J_cse_replaced.adjugate() / J_det_symbol |
| |
| | EpsilonDivDiv_kernel_gen.J_invT = J_inv.transpose() |
| |
| | EpsilonDivDiv_kernel_gen.J_invT_replacements |
| |
| | EpsilonDivDiv_kernel_gen.J_invT_reduced_exprs |
| |
| | EpsilonDivDiv_kernel_gen.J_invT_cse = sp.Matrix(3, 3, J_invT_reduced_exprs) |
| |
| | EpsilonDivDiv_kernel_gen.J_invT_cse_assignments |
| |
| | EpsilonDivDiv_kernel_gen.J_invT_cse_replaced |
| |
| list | EpsilonDivDiv_kernel_gen.grad_exprs = [] |
| |
| list | EpsilonDivDiv_kernel_gen.scalar_grad_is = [] |
| |
| | EpsilonDivDiv_kernel_gen.scalar_grad_i_replacements |
| |
| | EpsilonDivDiv_kernel_gen.scalar_grad_i_reduced_exprs |
| |
| | EpsilonDivDiv_kernel_gen.r_inv |
| |
| | EpsilonDivDiv_kernel_gen.g3 |
| |
| | EpsilonDivDiv_kernel_gen.grad_r |
| |
| | EpsilonDivDiv_kernel_gen.grad_r_inv |
| |
| | EpsilonDivDiv_kernel_gen.real |
| |
| | EpsilonDivDiv_kernel_gen.g1 = IndexedBase("grad1", shape=(3)) |
| |
| | EpsilonDivDiv_kernel_gen.g2 = IndexedBase("grad2", shape=(3)) |
| |
| | EpsilonDivDiv_kernel_gen.d_symbol = symbols("d", integer=True) |
| |
| | EpsilonDivDiv_kernel_gen.value |
| |
| int | EpsilonDivDiv_kernel_gen.r = 1 / r_inv |
| |
| list | EpsilonDivDiv_kernel_gen.dimloop_j_body = [] |
| |
| list | EpsilonDivDiv_kernel_gen.dimloop_j_exprs = [] |
| |
| list | EpsilonDivDiv_kernel_gen.dimloop_i_body = [] |
| |
| list | EpsilonDivDiv_kernel_gen.dimloop_i_exprs = [] |
| |
| str | EpsilonDivDiv_kernel_gen.E_trial_name = "E_grad_trial" |
| |
| | EpsilonDivDiv_kernel_gen.E_trial = IndexedBase(E_trial_name, shape=(3, 3), real=True) |
| |
| | EpsilonDivDiv_kernel_gen.div_u = symbols("div_u", real=True) |
| |
| list | EpsilonDivDiv_kernel_gen.u_grad_loop_exprs = [] |
| |
| | EpsilonDivDiv_kernel_gen.E_trial_matrix = Matrix(3, 3, [E_trial[i, j] for i in range(dim) for j in range(dim)]) |
| |
| str | EpsilonDivDiv_kernel_gen.grad_u_name = "grad_u" |
| |
| | EpsilonDivDiv_kernel_gen.grad_u = IndexedBase(grad_u_name, shape=(3, 3), real=True) |
| |
| | EpsilonDivDiv_kernel_gen.grad_u_matrix = Matrix(3, 3, [grad_u[i, j] for i in range(dim) for j in range(dim)]) |
| |
| float | EpsilonDivDiv_kernel_gen.symm_grad_j = 0.5 * (E_trial_matrix + (E_trial_matrix).transpose()) |
| |
| | EpsilonDivDiv_kernel_gen.symm_grad_j_replacements |
| |
| | EpsilonDivDiv_kernel_gen.symm_grad_j_reduced_exprs |
| |
| str | EpsilonDivDiv_kernel_gen.E_test_name = "E_grad_test" |
| |
| | EpsilonDivDiv_kernel_gen.E_test = IndexedBase(E_test_name, shape=(3, 3), real=True) |
| |
| list | EpsilonDivDiv_kernel_gen.pairing_loop_exprs = [] |
| |
| | EpsilonDivDiv_kernel_gen.E_test_matrix = Matrix(3, 3, [E_test[i, j] for i in range(dim) for j in range(dim)]) |
| |
| float | EpsilonDivDiv_kernel_gen.symm_grad_i = 0.5 * (E_test_matrix + (E_test_matrix).transpose()) |
| |
| | EpsilonDivDiv_kernel_gen.symm_grad_i_replacements |
| |
| | EpsilonDivDiv_kernel_gen.symm_grad_i_reduced_exprs |
| |
| | EpsilonDivDiv_kernel_gen.pairing_replacements |
| |
| | EpsilonDivDiv_kernel_gen.pairing_reduced_exprs |
| |
| list | EpsilonDivDiv_kernel_gen.boundary_loop_exprs = [] |
| |
| str | EpsilonDivDiv_kernel_gen.grad_u_diag_name = "grad_u_diag" |
| |
| | EpsilonDivDiv_kernel_gen.grad_u_diag = IndexedBase(grad_u_diag_name, shape=(3, 3), real=True) |
| |
| | EpsilonDivDiv_kernel_gen.grad_u_diag_matrix |
| |
| | EpsilonDivDiv_kernel_gen.symm_grad_i_reduced_matrix = Matrix(3, 3, symm_grad_i_reduced_exprs[0]) |
| |
| | EpsilonDivDiv_kernel_gen.diag_pairing_replacements |
| |
| | EpsilonDivDiv_kernel_gen.diag_pairing_reduced_exprs |
| |
| list | EpsilonDivDiv_kernel_gen.dimloop_diagBC_body |
| |
| | EpsilonDivDiv_kernel_gen.encoding |
| |