Home
last modified time | relevance | path

Searched refs:JxW (Results 1 – 25 of 198) sorted by relevance

12345678

/dports/math/libmesh/libmesh-1.6.2/examples/systems_of_equations/systems_of_equations_ex4/
H A Dsystems_of_equations_ex4.C189 const std::vector<Real> & JxW = fe->get_JxW(); in assemble_elasticity() local
241 …Kuu(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
244 …Kuu(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
247 …Kuu(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
250 …Kuu(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
261 …Kuv(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
264 …Kuv(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
267 …Kuv(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
270 …Kuv(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
281 …Kvu(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
[all …]
/dports/math/libmesh/libmesh-1.6.2/src/error_estimation/
H A Dhp_coarsentest.C113 Fe(i) += (*JxW)[qp] * in add_projection()
116 Fe(i) += (*JxW)[qp] * in add_projection()
119 Fe(i) += (*JxW)[qp] * in add_projection()
124 Ke(i,j) += (*JxW)[qp] * in add_projection()
127 Ke(i,j) += (*JxW)[qp] * in add_projection()
130 Ke(i,j) += (*JxW)[qp] * in add_projection()
209 JxW = &(fe->get_JxW()); in select_refinement()
342 Fe(i) += (*JxW)[qp] * in select_refinement()
345 Fe(i) += (*JxW)[qp] * in select_refinement()
348 Fe(i) += (*JxW)[qp] * in select_refinement()
[all …]
H A Dweighted_patch_recovery_error_estimator.C231 const std::vector<Real> & JxW = fe->get_JxW(); in operator() local
359 Kp(i,j) += JxW[qp]*psi[i]*psi[j]; in operator()
373 F(i) = JxW[qp]*u_h*psi[i]; in operator()
392 Fx(i) += JxW[qp]*grad_u_h(0)*psi[i]; in operator()
394 Fy(i) += JxW[qp]*grad_u_h(1)*psi[i]; in operator()
397 Fz(i) += JxW[qp]*grad_u_h(2)*psi[i]; in operator()
416 Fx(i) += JxW[qp]*grad_u_h(0)*psi[i]; in operator()
435 Fy(i) += JxW[qp]*grad_u_h(1)*psi[i]; in operator()
455 Fz(i) += JxW[qp]*grad_u_h(2)*psi[i]; in operator()
476 Fx(i) += JxW[qp]*hess_u_h(0,0)*psi[i]; in operator()
[all …]
H A Dpatch_recovery_error_estimator.C330 const std::vector<Real> & JxW = fe->get_JxW(); in operator() local
457 Kp(i,j) += JxW[qp]*psi[i]*psi[j]; in operator()
471 F(i) += JxW[qp]*u_h*psi[i]; in operator()
488 Fx(i) += JxW[qp]*grad_u_h(0)*psi[i]; in operator()
490 Fy(i) += JxW[qp]*grad_u_h(1)*psi[i]; in operator()
493 Fz(i) += JxW[qp]*grad_u_h(2)*psi[i]; in operator()
510 Fx(i) += JxW[qp]*grad_u_h(0)*psi[i]; in operator()
527 Fy(i) += JxW[qp]*grad_u_h(1)*psi[i]; in operator()
545 Fz(i) += JxW[qp]*grad_u_h(2)*psi[i]; in operator()
564 Fx(i) += JxW[qp]*hess_u_h(0,0)*psi[i]; in operator()
[all …]
/dports/math/libmesh/libmesh-1.6.2/examples/systems_of_equations/systems_of_equations_ex5/
H A Dsystems_of_equations_ex5.C200 const std::vector<Real> & JxW = fe->get_JxW(); in assemble_elasticity() local
260 …Kuu(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
263 …Kuu(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
266 …Kuu(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
269 …Kuu(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
280 …Kuv(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
283 …Kuv(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
286 …Kuv(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
289 …Kuv(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
300 …Kvu(i,j) += JxW[qp]*(eval_elasticity_tensor(C_i, C_j, C_k, C_l) * dphi[i][qp](C_j)*dphi[j][qp](C_l… in assemble_elasticity()
[all …]
/dports/math/libmesh/libmesh-1.6.2/examples/reduced_basis/reduced_basis_ex3/
H A Dassembly.h50 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
63 c.get_elem_jacobian()(i,j) += JxW[qp] *phi[j][qp]*phi[i][qp]; in interior_assembly()
77 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
92 c.get_elem_jacobian()(i,j) += JxW[qp] * dphi[j][qp]*dphi[i][qp]; in interior_assembly()
107 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
122 c.get_elem_jacobian()(i,j) += JxW[qp] * dphi[j][qp](0)*phi[i][qp]; in interior_assembly()
136 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
151 c.get_elem_jacobian()(i,j) += JxW[qp] * dphi[j][qp](1)*phi[i][qp]; in interior_assembly()
165 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
177 c.get_elem_residual()(i) += JxW[qp] * (1.*phi[i][qp]); in interior_assembly()
[all …]
/dports/math/libmesh/libmesh-1.6.2/examples/reduced_basis/reduced_basis_ex2/
H A Dassembly.h49 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
81 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
102 c.get_elem_jacobian()(i,j) += JxW[qp] * dphi[j][qp]*dphi[i][qp]; in interior_assembly()
116 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
137 c.get_elem_jacobian()(i,j) += JxW[qp] * dphi[j][qp]*dphi[i][qp]; in interior_assembly()
151 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
166 c.get_elem_jacobian()(i,j) -= JxW[qp] * dphi[i][qp](0)*phi[j][qp]; in interior_assembly()
181 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
193 c.get_elem_residual()(i) += JxW[qp] * (1.*phi[i][qp]); in interior_assembly()
215 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
[all …]
/dports/math/libmesh/libmesh-1.6.2/examples/reduced_basis/reduced_basis_ex1/
H A Dassembly.h49 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
64 c.get_elem_jacobian()(i,j) += JxW[qp] * dphi[j][qp]*dphi[i][qp]; in interior_assembly()
79 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
94 c.get_elem_jacobian()(i,j) += JxW[qp] * dphi[j][qp](0)*phi[i][qp]; in interior_assembly()
108 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
123 c.get_elem_jacobian()(i,j) += JxW[qp] * dphi[j][qp](1)*phi[i][qp]; in interior_assembly()
137 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
149 c.get_elem_residual()(i) += JxW[qp] * (1.*phi[i][qp]); in interior_assembly()
172 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
189 c.get_elem_residual()(i) += JxW[qp] * (1.*phi[i][qp]) / output_area; in interior_assembly()
/dports/math/deal.ii/dealii-803d21ff957e349b3799cd3ef2c840bc78734305/examples/step-30/
H A Dstep-30.cc194 const std::vector<double> &JxW = fe_v.get_JxW_values(); in assemble_cell_term() local
210 rhs[point] * fe_v.shape_value(i, point) * JxW[point]; in assemble_cell_term()
221 const std::vector<double> & JxW = fe_v.get_JxW_values(); in assemble_boundary_term() local
276 JxW[point]; in assemble_face_term()
290 fe_v_neighbor.shape_value(k, point) * JxW[point]; in assemble_face_term()
773 const std::vector<double> &JxW = in set_anisotropic_flags() local
795 area[face_no / 2] += JxW[x]; in set_anisotropic_flags()
818 const std::vector<double> &JxW = in set_anisotropic_flags() local
827 area[face_no / 2] += JxW[x]; in set_anisotropic_flags()
862 const std::vector<double> &JxW = in set_anisotropic_flags() local
[all …]
/dports/math/libmesh/libmesh-1.6.2/examples/adjoints/adjoints_ex4/
H A DL-shaped.C73 const std::vector<Real> & JxW = elem_fe->get_JxW(); in element_time_derivative() local
100 F(i) += JxW[qp] * (grad_T * dphi[i][qp]); in element_time_derivative()
105 K(i,j) += JxW[qp] * (dphi[i][qp] * dphi[j][qp]); in element_time_derivative()
126 const std::vector<Real> & JxW = side_fe->get_JxW(); in side_constraint() local
155 F(i) += JxW[qp] * penalty * (T - u_dirichlet) * phi[i][qp]; in side_constraint()
161 K(i,j) += JxW[qp] * penalty * phi[i][qp] * phi[j][qp]; in side_constraint()
/dports/math/libmesh/libmesh-1.6.2/examples/adjoints/adjoints_ex1/
H A DL-shaped.C72 const std::vector<Real> & JxW = elem_fe->get_JxW(); in element_time_derivative() local
99 F(i) += JxW[qp] * (grad_T * dphi[i][qp]); in element_time_derivative()
104 K(i,j) += JxW[qp] * (dphi[i][qp] * dphi[j][qp]); in element_time_derivative()
125 const std::vector<Real> & JxW = side_fe->get_JxW(); in side_constraint() local
154 F(i) += JxW[qp] * penalty * (T - u_dirichlet) * phi[i][qp]; in side_constraint()
159 K(i,j) += JxW[qp] * penalty * phi[i][qp] * phi[j][qp]; in side_constraint()
/dports/math/libmesh/libmesh-1.6.2/examples/reduced_basis/reduced_basis_ex4/
H A Dassembly.h77 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
92 c.get_elem_jacobian()(i,j) += JxW[qp] * dphi[j][qp]*dphi[i][qp]; in interior_assembly()
107 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
118 c.get_elem_jacobian()(i,j) += JxW[qp] * phi[j][qp]*phi[i][qp]; in interior_assembly()
141 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
155 c.get_elem_residual()(i) += JxW[qp] * (eim_values[qp]*phi[i][qp]); in interior_assembly()
/dports/math/libmesh/libmesh-1.6.2/examples/adjoints/adjoints_ex2/
H A DL-shaped.C69 const std::vector<Real> & JxW = elem_fe->get_JxW(); in element_time_derivative() local
96 F(i) += JxW[qp] * (parameters[0] + (2.*parameters[1])) * (grad_T * dphi[i][qp]); in element_time_derivative()
101 K(i,j) += JxW[qp] * (parameters[0] + (2.*parameters[1])) * (dphi[i][qp] * dphi[j][qp]); in element_time_derivative()
122 const std::vector<Real> & JxW = side_fe->get_JxW(); in side_constraint() local
151 F(i) += JxW[qp] * penalty * (T - u_dirichlet) * phi[i][qp]; in side_constraint()
156 K(i,j) += JxW[qp] * penalty * phi[i][qp] * phi[j][qp]; in side_constraint()
H A DL-qoi.C36 const std::vector<Real> & JxW = elem_fe->get_JxW(); in element_qoi() local
58 dQoI_0 += JxW[qp] * T; in element_qoi()
79 const std::vector<Real> & JxW = elem_fe->get_JxW(); in element_qoi_derivative() local
106 Q(i) += JxW[qp]*phi[i][qp]; in element_qoi_derivative()
/dports/math/libmesh/libmesh-1.6.2/examples/adjoints/adjoints_ex3/
H A DH-qoi.C58 const std::vector<Real> & JxW = side_fe->get_JxW(); in side_qoi_derivative() local
99 Qu(i) += JxW[qp] * -phi[i][qp] * C; in side_qoi_derivative()
100 QC(i) += JxW[qp] * phi[i][qp] * -u; in side_qoi_derivative()
122 const std::vector<Real> & JxW = side_fe->get_JxW(); in side_qoi() local
147 dQoI_0 += JxW[qp] * -u * C; in side_qoi()
H A Dcoupled_system.C228 const std::vector<Real> & JxW = u_elem_fe->get_JxW(); in element_time_derivative() local
320 F[d](i) += JxW[qp] * in element_time_derivative()
325 FC(i) += JxW[qp] * in element_time_derivative()
341 Kdiag[d](i,j) += JxW[qp] * (-(dphi[i][qp]*dphi[j][qp])); in element_time_derivative()
344 C[d](i,j) += JxW[qp] * ((phi[j][qp]*grad_C(d))*phi[i][qp]); in element_time_derivative()
348 KCC(i,j) += JxW[qp]* in element_time_derivative()
356 B[d](i,j) += JxW[qp]*psi[j][qp]*dphi[i][qp](d); in element_time_derivative()
379 const std::vector<Real> & JxW = u_elem_fe->get_JxW(); in element_constraint() local
421 Fp(i) += JxW[qp] * psi[i][qp] * grad_uv[d](d); in element_constraint()
429 B[d](i,j) += JxW[qp]*psi[i][qp]*dphi[j][qp](d); in element_constraint()
/dports/math/deal.ii/dealii-803d21ff957e349b3799cd3ef2c840bc78734305/examples/step-12b/
H A Dstep-12b.cc346 const std::vector<double> &JxW = fe_values.get_JxW_values(); in integrate_cell_term() local
363 JxW[point]; in integrate_cell_term()
378 const std::vector<double> & JxW = fe_face_values.get_JxW_values(); in integrate_boundary_term() local
398 JxW[point]; in integrate_boundary_term()
404 JxW[point]; in integrate_boundary_term()
443 const std::vector<double> & JxW = fe_face_values.get_JxW_values(); in integrate_face_term() local
460 JxW[point]; in integrate_face_term()
470 JxW[point]; in integrate_face_term()
481 JxW[point]; in integrate_face_term()
491 JxW[point]; in integrate_face_term()
/dports/math/libmesh/libmesh-1.6.2/examples/systems_of_equations/systems_of_equations_ex3/
H A Dsystems_of_equations_ex3.C387 const std::vector<Real> & JxW = fe_vel->get_JxW(); in assemble_stokes() local
574 Fu(i) += JxW[qp]*(u_old*phi[i][qp] - // mass-matrix term in assemble_stokes()
581 Fv(i) += JxW[qp]*(v_old*phi[i][qp] - // mass-matrix term in assemble_stokes()
599 Kuv(i,j) += JxW[qp]*theta*dt*u_y*phi[i][qp]*phi[j][qp]; // Newton term in assemble_stokes()
606 Kvu(i,j) += JxW[qp]*theta*dt*v_x*phi[i][qp]*phi[j][qp]; // Newton term in assemble_stokes()
612 Kup(i,j) += JxW[qp]*(-theta*dt*psi[j][qp]*dphi[i][qp](0)); in assemble_stokes()
613 Kvp(i,j) += JxW[qp]*(-theta*dt*psi[j][qp]*dphi[i][qp](1)); in assemble_stokes()
623 Kp_alpha(i,0) += JxW[qp]*psi[i][qp]; in assemble_stokes()
624 Kalpha_p(0,i) += JxW[qp]*psi[i][qp]; in assemble_stokes()
627 Kpu(i,j) += JxW[qp]*psi[i][qp]*dphi[j][qp](0); in assemble_stokes()
[all …]
/dports/math/deal.ii/dealii-803d21ff957e349b3799cd3ef2c840bc78734305/include/deal.II/integrators/
H A Ddivergence.h68 const double dx = fe.JxW(k) * factor; in JxW() function
106 const double dx = factor * fetest.JxW(k);
139 const double dx = factor * fetest.JxW(k);
172 const double dx = fe.JxW(k) * factor; in JxW() function
210 const double dx = factor * fetest.JxW(k);
243 const double dx = factor * fetest.JxW(k);
274 const Tensor<1, dim> ndx = factor * fe.JxW(k) * fe.normal_vector(k);
338 factor * fetest.normal_vector(k) * fetest.JxW(k);
386 const double dx = factor * fe1.JxW(k);
440 const double dx = factor * fe1.JxW(k);
[all …]
/dports/math/libmesh/libmesh-1.6.2/src/apps/
H A Damr.C158 const std::vector<Real> & JxW = fe->get_JxW(); in assemble() local
213 Fu(i) += JxW[gp]*f*phi[i][gp]; in assemble()
214 Fv(i) += JxW[gp]*f*phi[i][gp]; in assemble()
219 Kuu(i,j) += JxW[gp]*((phi[i][gp])*(phi[j][gp])); in assemble()
221 Kvv(i,j) += JxW[gp]*((phi[i][gp])*(phi[j][gp]) + in assemble()
225 vol += JxW[gp]; in assemble()
/dports/math/libmesh/libmesh-1.6.2/examples/vector_fe/vector_fe_ex3/
H A Dcurl_curl_system.C117 const std::vector<Real> & JxW = fe->get_JxW(); in element_time_derivative() local
159 Fu(i) += (curl_u*curl_phi[i][qp] + u*phi[i][qp] - f*phi[i][qp])*JxW[qp]; in element_time_derivative()
165 phi[j][qp]*phi[i][qp])*JxW[qp]; in element_time_derivative()
186 const std::vector<Real> & JxW = side_fe->get_JxW(); in side_time_derivative() local
220 Fu(i) += penalty*Ncu*(phi[i][qp].cross(N))*JxW[qp]; in side_time_derivative()
224 Kuu(i,j) += penalty*(phi[j][qp].cross(N))*(phi[i][qp].cross(N))*JxW[qp]; in side_time_derivative()
/dports/math/libmesh/libmesh-1.6.2/src/solvers/
H A Dfirst_order_unsteady_solver.C52 const std::vector<Real> & JxW = elem_fe->get_JxW(); in compute_second_order_eqns() local
71 Fu(i) += JxW[qp]*(udot-v)*phi[i][qp]; in compute_second_order_eqns()
75 … Number rate_factor = JxW[qp]*context.get_elem_solution_rate_derivative()*phi[i][qp]; in compute_second_order_eqns()
76 Number soln_factor = JxW[qp]*context.get_elem_solution_derivative()*phi[i][qp]; in compute_second_order_eqns()
/dports/math/libmesh/libmesh-1.6.2/examples/vector_fe/vector_fe_ex4/
H A Dcurl_curl_system.C118 const std::vector<Real> & JxW = fe->get_JxW(); in element_time_derivative() local
160 Fu(i) += (curl_u*curl_phi[i][qp] + u*phi[i][qp] - f*phi[i][qp])*JxW[qp]; in element_time_derivative()
166 phi[j][qp]*phi[i][qp])*JxW[qp]; in element_time_derivative()
187 const std::vector<Real> & JxW = side_fe->get_JxW(); in side_time_derivative() local
221 Fu(i) += penalty*Ncu*(phi[i][qp].cross(N))*JxW[qp]; in side_time_derivative()
225 Kuu(i,j) += penalty*(phi[j][qp].cross(N))*(phi[i][qp].cross(N))*JxW[qp]; in side_time_derivative()
/dports/math/libmesh/libmesh-1.6.2/examples/systems_of_equations/systems_of_equations_ex1/
H A Dsystems_of_equations_ex1.C197 const std::vector<Real> & JxW = fe_vel->get_JxW(); in assemble_stokes() local
312 Kuu(i,j) += JxW[qp]*(dphi[i][qp]*dphi[j][qp]); in assemble_stokes()
317 Kup(i,j) += -JxW[qp]*psi[j][qp]*dphi[i][qp](0); in assemble_stokes()
324 Kvv(i,j) += JxW[qp]*(dphi[i][qp]*dphi[j][qp]); in assemble_stokes()
329 Kvp(i,j) += -JxW[qp]*psi[j][qp]*dphi[i][qp](1); in assemble_stokes()
336 Kpu(i,j) += -JxW[qp]*psi[i][qp]*dphi[j][qp](0); in assemble_stokes()
341 Kpv(i,j) += -JxW[qp]*psi[i][qp]*dphi[j][qp](1); in assemble_stokes()
/dports/math/libmesh/libmesh-1.6.2/examples/reduced_basis/reduced_basis_ex5/
H A Dassembly.C74 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
99 JxW[qp]*(C_ijkl * dphi[i][qp](C_j)*dphi[j][qp](C_l)); in interior_assembly()
114 JxW[qp]*(C_ijkl * dphi[i][qp](C_j)*dphi[j][qp](C_l)); in interior_assembly()
133 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
157 JxW[qp]*(C_ijkl * dphi[i][qp](C_j)*dphi[j][qp](C_l)); in interior_assembly()
175 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
202 JxW[qp]*(C_ijkl * dphi[i][qp](C_j)*dphi[j][qp](C_l)); in interior_assembly()
347 const std::vector<Real> & JxW = elem_fe->get_JxW(); in interior_assembly() local
371 Kdiag[d](i,j) += JxW[qp]*(dphi[i][qp]*dphi[j][qp]); in interior_assembly()

12345678