Loading...
Searching...
No Matches
EpsilonDivDiv_kernel_gen.py File Reference

Namespaces

namespace  EpsilonDivDiv_kernel_gen
 

Functions

 EpsilonDivDiv_kernel_gen.create_col_assigns (E, col, col_vec)
 
 EpsilonDivDiv_kernel_gen.double_contract (A, B)
 

Variables

 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