pt 551 old_code/tumor_growth/tumor_growth.h for(size_t pt = 0; pt < quadrature.nop(); ++pt){ pt 554 old_code/tumor_growth/tumor_growth.h uLocal.evaluate(quadrature[pt], u_x); // initialize pt 557 old_code/tumor_growth/tumor_growth.h uLocal.jacobian(quadrature[pt], du_x); // initialize pt 561 old_code/tumor_growth/tumor_growth.h model.ran_to_ran(entity, quadrature[pt], u_x, u_x_psi); pt 567 old_code/tumor_growth/tumor_growth.h model.jac_ran_to_jac_ran_lhs(entity, quadrature[pt], du_x, du_x_dpsi); pt 570 old_code/tumor_growth/tumor_growth.h const typename QuadratureType::CoordinateType& x = quadrature.point(pt); pt 571 old_code/tumor_growth/tumor_growth.h const double weight = quadrature.weight(pt) * geometry.integrationElement(x); pt 576 old_code/tumor_growth/tumor_growth.h wLocal.axpy(quadrature[pt], u_x_psi, du_x_dpsi); pt 603 old_code/tumor_growth/tumor_growth.h for(size_t pt = 0; pt < quadrature.nop(); ++pt){ pt 605 old_code/tumor_growth/tumor_growth.h uLocal.evaluate(quadrature[pt], u_x); // initialize pt 608 old_code/tumor_growth/tumor_growth.h u_approx_local.evaluate(quadrature[pt], u_approx_x); pt 611 old_code/tumor_growth/tumor_growth.h uLocal.jacobian(quadrature[pt], du_x); // initialize pt 615 old_code/tumor_growth/tumor_growth.h model.ran_to_ran(entity, quadrature[pt], u_x, u_x_psi); // initialize pt 620 old_code/tumor_growth/tumor_growth.h model.jac_ran_to_jac_ran_rhs(entity, quadrature[pt], du_x, du_x_dpsi); pt 625 old_code/tumor_growth/tumor_growth.h model.ran_u_normal_rhs(entity, quadrature[pt], u_approx_x, un_x_psi); pt 627 old_code/tumor_growth/tumor_growth.h const typename QuadratureType::CoordinateType& x = quadrature.point(pt); pt 628 old_code/tumor_growth/tumor_growth.h const double weight = quadrature.weight(pt) * geometry.integrationElement(x); pt 635 old_code/tumor_growth/tumor_growth.h wLocal.axpy(quadrature[pt], u_x_psi, du_x_dpsi); pt 636 old_code/tumor_growth/tumor_growth.h wLocal.axpy(quadrature[pt], un_x_psi, 0); pt 661 old_code/tumor_growth/tumor_growth.h for(size_t pt = 0; pt < quadrature.nop(); ++pt){ pt 663 old_code/tumor_growth/tumor_growth.h uLocal.evaluate(quadrature[pt], u_x); // initialize pt 666 old_code/tumor_growth/tumor_growth.h uLocal.jacobian(quadrature[pt], du_x); // initialize pt 670 old_code/tumor_growth/tumor_growth.h model.ran_to_ran_lhs(entity, quadrature[pt], u_x, u_x_chi); pt 677 old_code/tumor_growth/tumor_growth.h model.jac_ran_to_jac_ran_lhs(entity, quadrature[pt], du_x, du_x_dchi); pt 680 old_code/tumor_growth/tumor_growth.h const typename QuadratureType::CoordinateType& x = quadrature.point(pt); pt 681 old_code/tumor_growth/tumor_growth.h const double weight = quadrature.weight(pt) * geometry.integrationElement(x); pt 686 old_code/tumor_growth/tumor_growth.h wLocal.axpy(quadrature[pt], u_x_chi, du_x_dchi); pt 693 old_code/tumor_growth/tumor_growth.h for(size_t pt = 0; pt < bigger_quadrature.nop(); ++pt){ pt 695 old_code/tumor_growth/tumor_growth.h uLocal.evaluate(bigger_quadrature[pt], u_x); // initialize pt 697 old_code/tumor_growth/tumor_growth.h u_approx_local.evaluate(bigger_quadrature[pt], u_app); // initialize pt 699 old_code/tumor_growth/tumor_growth.h u_or_w_approx_local.evaluate(bigger_quadrature[pt], u_or_w_app); pt 704 old_code/tumor_growth/tumor_growth.h model.tri_ran_to_ran_lhs(entity, bigger_quadrature[pt], pt 707 old_code/tumor_growth/tumor_growth.h const typename QuadratureType::CoordinateType& x = bigger_quadrature.point(pt); pt 709 old_code/tumor_growth/tumor_growth.h = bigger_quadrature.weight(pt) * geometry.integrationElement(x); pt 712 old_code/tumor_growth/tumor_growth.h wLocal.axpy(bigger_quadrature[pt], uwu_x_chi, 0); pt 732 old_code/tumor_growth/tumor_growth.h for(size_t pt = 0; pt < quadrature.nop(); ++pt){ pt 735 old_code/tumor_growth/tumor_growth.h uLocal.evaluate(quadrature[pt], u_x); // initialize pt 739 old_code/tumor_growth/tumor_growth.h model.ran_to_ran(entity, quadrature[pt], u_x, u_x_chi); pt 742 old_code/tumor_growth/tumor_growth.h const typename QuadratureType::CoordinateType& x = quadrature.point(pt); pt 743 old_code/tumor_growth/tumor_growth.h const double weight = quadrature.weight(pt) * geometry.integrationElement(x); pt 747 old_code/tumor_growth/tumor_growth.h wLocal.axpy(quadrature[pt], u_x_chi, 0); pt 766 old_code/tumor_growth/tumor_growth.h for(size_t pt = 0; pt < quadrature.nop(); ++pt){ pt 770 old_code/tumor_growth/tumor_growth.h model.ran_to_ran_ones(entity, quadrature[pt], ones_chi); pt 773 old_code/tumor_growth/tumor_growth.h const typename QuadratureType::CoordinateType& x = quadrature.point(pt); pt 774 old_code/tumor_growth/tumor_growth.h const double weight = quadrature.weight(pt) * geometry.integrationElement(x); pt 778 old_code/tumor_growth/tumor_growth.h wLocal.axpy(quadrature[pt], ones_chi, 0); pt 48 src/secOrd_op_brusselator_impl.cpp inline Range massMatrix_loc(const std::size_t pt, const Quadrature& q, pt 51 src/secOrd_op_brusselator_impl.cpp cf.evaluate(q[pt], u); pt 54 src/secOrd_op_brusselator_impl.cpp inline Jacobian_range stiffnessMatrix_loc(const std::size_t pt, const Quadrature& q, pt 57 src/secOrd_op_brusselator_impl.cpp cf.jacobian(q[pt], nabla_u); pt 60 src/secOrd_op_brusselator_impl.cpp inline Range quad_massMatrix_loc(const std::size_t pt, const Quadrature& q, pt 66 src/secOrd_op_brusselator_impl.cpp rhs_loc.evaluate(q[pt], u); pt 67 src/secOrd_op_brusselator_impl.cpp arg1_loc.evaluate(q[pt], tmp); pt 69 src/secOrd_op_brusselator_impl.cpp arg2_loc.evaluate(q[pt], tmp); pt 73 src/secOrd_op_brusselator_impl.cpp inline Range jacobian_quadMass_loc(const std::size_t pt, const Quadrature& q, pt 78 src/secOrd_op_brusselator_impl.cpp arg1_loc.evaluate(q[pt], u); pt 79 src/secOrd_op_brusselator_impl.cpp arg2_loc.evaluate(q[pt], tmp); pt 169 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ pt 170 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); pt 171 src/secOrd_op_brusselator_impl.cpp const auto weight = q.weight(pt) * g.integrationElement(x); pt 172 src/secOrd_op_brusselator_impl.cpp auto Mu = massMatrix_loc(pt, q, rhs_loc); pt 175 src/secOrd_op_brusselator_impl.cpp lhs_loc.axpy(q[pt], Mu); pt 186 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ pt 187 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); pt 188 src/secOrd_op_brusselator_impl.cpp const auto weight = q.weight(pt) * g.integrationElement(x); pt 191 src/secOrd_op_brusselator_impl.cpp lhs_loc.axpy(q[pt], Mu); pt 218 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ pt 220 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); pt 221 src/secOrd_op_brusselator_impl.cpp const auto weight = q.weight(pt) * g.integrationElement(x); pt 223 src/secOrd_op_brusselator_impl.cpp auto Mu = massMatrix_loc(pt, q, rhs_loc); pt 224 src/secOrd_op_brusselator_impl.cpp auto Au = stiffnessMatrix_loc(pt, q, rhs_loc); pt 229 src/secOrd_op_brusselator_impl.cpp lhs_loc.axpy(q[pt], Mu, Au); pt 241 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ pt 242 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); pt 244 src/secOrd_op_brusselator_impl.cpp = q.weight(pt) * g.integrationElement(x) pt 247 src/secOrd_op_brusselator_impl.cpp auto Mu = quad_massMatrix_loc(pt, q, rhs_loc, arg1_loc, arg2_loc); pt 250 src/secOrd_op_brusselator_impl.cpp lhs_loc.axpy(q[pt], Mu); pt 263 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ pt 264 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); pt 265 src/secOrd_op_brusselator_impl.cpp const auto integral_factor = q.weight(pt) * geometry.integrationElement(x); pt 267 src/secOrd_op_brusselator_impl.cpp base_set.evaluateAll(q[pt], vec_range); pt 268 src/secOrd_op_brusselator_impl.cpp base_set.jacobianAll(q[pt], vec_jacRange); pt 297 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ pt 298 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); pt 300 src/secOrd_op_brusselator_impl.cpp = q.weight(pt) * geometry.integrationElement(x) pt 303 src/secOrd_op_brusselator_impl.cpp base_set.evaluateAll(q[pt], vec_range); pt 304 src/secOrd_op_brusselator_impl.cpp const auto u0 = jacobian_quadMass_loc(pt, q, arg1_loc, arg2_loc); pt 60 src/secOrd_op_linearHeat.cpp Range mass_matrix(const std::size_t pt, const Quadrature& q, pt 63 src/secOrd_op_linearHeat.cpp cf.evaluate(q[pt], u); pt 67 src/secOrd_op_linearHeat.cpp Jacobian_range stiffness_matrix(const std::size_t pt, const Quadrature& q, pt 70 src/secOrd_op_linearHeat.cpp cf.jacobian(q[pt], nabla_u); pt 172 src/secOrd_op_linearHeat.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ pt 174 src/secOrd_op_linearHeat.cpp auto Mu = mass_matrix(pt, q, cf); pt 175 src/secOrd_op_linearHeat.cpp auto Au = stiffness_matrix(pt, q, cf); pt 177 src/secOrd_op_linearHeat.cpp const auto& x = q.point(pt); pt 178 src/secOrd_op_linearHeat.cpp const auto weight = q.weight(pt) * g.integrationElement(x); pt 182 src/secOrd_op_linearHeat.cpp f.axpy(q[pt], Mu, Au); pt 192 src/secOrd_op_linearHeat.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ pt 194 src/secOrd_op_linearHeat.cpp auto Mu = mass_matrix(pt, q, cf); pt 196 src/secOrd_op_linearHeat.cpp const auto& x = q.point(pt); pt 197 src/secOrd_op_linearHeat.cpp const auto weight = q.weight(pt) * g.integrationElement(x); pt 200 src/secOrd_op_linearHeat.cpp f.axpy(q[pt], Mu); pt 325 src/secOrd_op_rhs_impl.h for(std::size_t pt = 0; pt < quad.nop(); ++pt){ pt 326 src/secOrd_op_rhs_impl.h const auto& x = quad.point(pt); pt 329 src/secOrd_op_rhs_impl.h fx *= quad.weight(pt) * geometry.integrationElement(x); pt 330 src/secOrd_op_rhs_impl.h fef_local.axpy(quad[pt], fx); pt 152 src/secOrd_op_rhs_u.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ pt 153 src/secOrd_op_rhs_u.cpp const auto& x = q.point(pt); pt 155 src/secOrd_op_rhs_u.cpp fx *= q.weight(pt) * g.integrationElement(x); pt 156 src/secOrd_op_rhs_u.cpp f_loc.axpy(q[pt], fx); pt 145 src/secOrd_op_rhs_w.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ pt 146 src/secOrd_op_rhs_w.cpp const auto& x = q.point(pt); pt 149 src/secOrd_op_rhs_w.cpp fx *= q.weight(pt) * g.integrationElement(x); pt 150 src/secOrd_op_rhs_w.cpp f_loc.axpy(q[pt], fx); pt 117 src/secOrd_op_solutionDriven_impl.cpp for(size_t pt = 0; pt < q.nop(); ++pt){ pt 119 src/secOrd_op_solutionDriven_impl.cpp const auto& x = q.point(pt); pt 120 src/secOrd_op_solutionDriven_impl.cpp const auto integral_factor = q.weight(pt) * g.integrationElement(x); pt 122 src/secOrd_op_solutionDriven_impl.cpp auto X_p = evaluate(pt, q, cf); pt 125 src/secOrd_op_solutionDriven_impl.cpp auto dX_p = jacobian(pt, q, cf); pt 128 src/secOrd_op_solutionDriven_impl.cpp f.axpy(q[pt], X_p, dX_p); pt 136 src/secOrd_op_solutionDriven_impl.cpp for(size_t pt = 0; pt < q.nop(); ++pt){ pt 138 src/secOrd_op_solutionDriven_impl.cpp const auto& x = q.point(pt); pt 139 src/secOrd_op_solutionDriven_impl.cpp const auto integral_factor = q.weight(pt) * g.integrationElement(x); pt 143 src/secOrd_op_solutionDriven_impl.cpp const auto u_p = evaluate(pt, q, u_loc); pt 146 src/secOrd_op_solutionDriven_impl.cpp auto X_p = evaluate(pt, q, cf); pt 150 src/secOrd_op_solutionDriven_impl.cpp auto dX_p = jacobian(pt, q, cf); pt 153 src/secOrd_op_solutionDriven_impl.cpp f.axpy(q[pt], X_p, dX_p); pt 200 src/secOrd_op_solutionDriven_impl.h evaluate(const std::size_t pt, pt 204 src/secOrd_op_solutionDriven_impl.h cf.evaluate(q[pt], X); pt 209 src/secOrd_op_solutionDriven_impl.h evaluate(const std::size_t pt, pt 213 src/secOrd_op_solutionDriven_impl.h cf.evaluate(q[pt], u); pt 218 src/secOrd_op_solutionDriven_impl.h jacobian(const std::size_t pt, pt 222 src/secOrd_op_solutionDriven_impl.h cf.jacobian(q[pt], nabla_X);