/dports/math/libmesh/libmesh-1.6.2/examples/systems_of_equations/systems_of_equations_ex4/ |
H A D | systems_of_equations_ex4.C | 189 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 D | hp_coarsentest.C | 113 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 D | weighted_patch_recovery_error_estimator.C | 231 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 D | patch_recovery_error_estimator.C | 330 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 D | systems_of_equations_ex5.C | 200 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 D | assembly.h | 50 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 D | assembly.h | 49 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 D | assembly.h | 49 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 D | step-30.cc | 194 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 D | L-shaped.C | 73 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 D | L-shaped.C | 72 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 D | assembly.h | 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 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 D | L-shaped.C | 69 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 D | L-qoi.C | 36 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 D | H-qoi.C | 58 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 D | coupled_system.C | 228 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 D | step-12b.cc | 346 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 D | systems_of_equations_ex3.C | 387 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 D | divergence.h | 68 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 D | amr.C | 158 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 D | curl_curl_system.C | 117 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 D | first_order_unsteady_solver.C | 52 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 D | curl_curl_system.C | 118 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 D | systems_of_equations_ex1.C | 197 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 D | assembly.C | 74 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()
|