q 90 old_code/tumor_growth/tumor_growth.h void evaluate(const Domain& p, Range& q) const{ q 91 old_code/tumor_growth/tumor_growth.h q = random_fun(); q 110 old_code/tumor_growth/tumor_growth.h void evaluate(const Domain& p, Range& q) const{ q 112 old_code/tumor_growth/tumor_growth.h q = Range {base}; q 48 src/secOrd_op_brusselator_impl.cpp inline Range massMatrix_loc(const std::size_t pt, const Quadrature& q, q 51 src/secOrd_op_brusselator_impl.cpp cf.evaluate(q[pt], u); q 54 src/secOrd_op_brusselator_impl.cpp inline Jacobian_range stiffnessMatrix_loc(const std::size_t pt, const Quadrature& q, q 57 src/secOrd_op_brusselator_impl.cpp cf.jacobian(q[pt], nabla_u); q 60 src/secOrd_op_brusselator_impl.cpp inline Range quad_massMatrix_loc(const std::size_t pt, const Quadrature& q, q 66 src/secOrd_op_brusselator_impl.cpp rhs_loc.evaluate(q[pt], u); q 67 src/secOrd_op_brusselator_impl.cpp arg1_loc.evaluate(q[pt], tmp); q 69 src/secOrd_op_brusselator_impl.cpp arg2_loc.evaluate(q[pt], tmp); q 73 src/secOrd_op_brusselator_impl.cpp inline Range jacobian_quadMass_loc(const std::size_t pt, const Quadrature& q, q 78 src/secOrd_op_brusselator_impl.cpp arg1_loc.evaluate(q[pt], u); q 79 src/secOrd_op_brusselator_impl.cpp arg2_loc.evaluate(q[pt], tmp); q 168 src/secOrd_op_brusselator_impl.cpp Quadrature q {entity, lhs_loc.order() + rhs_loc.order()}; q 169 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ q 170 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); q 171 src/secOrd_op_brusselator_impl.cpp const auto weight = q.weight(pt) * g.integrationElement(x); q 172 src/secOrd_op_brusselator_impl.cpp auto Mu = massMatrix_loc(pt, q, rhs_loc); q 175 src/secOrd_op_brusselator_impl.cpp lhs_loc.axpy(q[pt], Mu); q 184 src/secOrd_op_brusselator_impl.cpp Quadrature q {entity, lhs_loc.order()}; q 186 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ q 187 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); q 188 src/secOrd_op_brusselator_impl.cpp const auto weight = q.weight(pt) * g.integrationElement(x); q 191 src/secOrd_op_brusselator_impl.cpp lhs_loc.axpy(q[pt], Mu); q 216 src/secOrd_op_brusselator_impl.cpp Quadrature q {e, lhs_loc.order() + rhs_loc.order()}; q 218 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ q 220 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); q 221 src/secOrd_op_brusselator_impl.cpp const auto weight = q.weight(pt) * g.integrationElement(x); q 223 src/secOrd_op_brusselator_impl.cpp auto Mu = massMatrix_loc(pt, q, rhs_loc); q 224 src/secOrd_op_brusselator_impl.cpp auto Au = stiffnessMatrix_loc(pt, q, rhs_loc); q 229 src/secOrd_op_brusselator_impl.cpp lhs_loc.axpy(q[pt], Mu, Au); q 238 src/secOrd_op_brusselator_impl.cpp Quadrature q q 241 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ q 242 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); q 244 src/secOrd_op_brusselator_impl.cpp = q.weight(pt) * g.integrationElement(x) q 247 src/secOrd_op_brusselator_impl.cpp auto Mu = quad_massMatrix_loc(pt, q, rhs_loc, arg1_loc, arg2_loc); q 250 src/secOrd_op_brusselator_impl.cpp lhs_loc.axpy(q[pt], Mu); q 262 src/secOrd_op_brusselator_impl.cpp Quadrature q {entity, 2 * fef_loc.order()}; q 263 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ q 264 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); q 265 src/secOrd_op_brusselator_impl.cpp const auto integral_factor = q.weight(pt) * geometry.integrationElement(x); q 267 src/secOrd_op_brusselator_impl.cpp base_set.evaluateAll(q[pt], vec_range); q 268 src/secOrd_op_brusselator_impl.cpp base_set.jacobianAll(q[pt], vec_jacRange); q 296 src/secOrd_op_brusselator_impl.cpp Quadrature q {entity, 2 * fef_loc.order() + arg1_loc.order() + arg2_loc.order()}; q 297 src/secOrd_op_brusselator_impl.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ q 298 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); q 300 src/secOrd_op_brusselator_impl.cpp = q.weight(pt) * geometry.integrationElement(x) q 303 src/secOrd_op_brusselator_impl.cpp base_set.evaluateAll(q[pt], vec_range); q 304 src/secOrd_op_brusselator_impl.cpp const auto u0 = jacobian_quadMass_loc(pt, q, arg1_loc, arg2_loc); q 404 src/secOrd_op_initData_impl.h evaluate(const Domain&, Range& q) const{ q 405 src/secOrd_op_initData_impl.h q = random_fun(); q 60 src/secOrd_op_linearHeat.cpp Range mass_matrix(const std::size_t pt, const Quadrature& q, q 63 src/secOrd_op_linearHeat.cpp cf.evaluate(q[pt], u); q 67 src/secOrd_op_linearHeat.cpp Jacobian_range stiffness_matrix(const std::size_t pt, const Quadrature& q, q 70 src/secOrd_op_linearHeat.cpp cf.jacobian(q[pt], nabla_u); q 166 src/secOrd_op_linearHeat.cpp const Quadrature& q, q 172 src/secOrd_op_linearHeat.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ q 174 src/secOrd_op_linearHeat.cpp auto Mu = mass_matrix(pt, q, cf); q 175 src/secOrd_op_linearHeat.cpp auto Au = stiffness_matrix(pt, q, cf); q 177 src/secOrd_op_linearHeat.cpp const auto& x = q.point(pt); q 178 src/secOrd_op_linearHeat.cpp const auto weight = q.weight(pt) * g.integrationElement(x); q 182 src/secOrd_op_linearHeat.cpp f.axpy(q[pt], Mu, Au); q 186 src/secOrd_op_linearHeat.cpp const Quadrature& q, q 192 src/secOrd_op_linearHeat.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ q 194 src/secOrd_op_linearHeat.cpp auto Mu = mass_matrix(pt, q, cf); q 196 src/secOrd_op_linearHeat.cpp const auto& x = q.point(pt); q 197 src/secOrd_op_linearHeat.cpp const auto weight = q.weight(pt) * g.integrationElement(x); q 200 src/secOrd_op_linearHeat.cpp f.axpy(q[pt], Mu); q 149 src/secOrd_op_rhs_u.cpp void massMatrix_for_entity(const Geometry& g, const Quadrature& q, q 152 src/secOrd_op_rhs_u.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ q 153 src/secOrd_op_rhs_u.cpp const auto& x = q.point(pt); q 155 src/secOrd_op_rhs_u.cpp fx *= q.weight(pt) * g.integrationElement(x); q 156 src/secOrd_op_rhs_u.cpp f_loc.axpy(q[pt], fx); q 142 src/secOrd_op_rhs_w.cpp void massMatrix_for_entity(const Geometry& g, const Quadrature& q, q 145 src/secOrd_op_rhs_w.cpp for(std::size_t pt = 0; pt < q.nop(); ++pt){ q 146 src/secOrd_op_rhs_w.cpp const auto& x = q.point(pt); q 149 src/secOrd_op_rhs_w.cpp fx *= q.weight(pt) * g.integrationElement(x); q 150 src/secOrd_op_rhs_w.cpp f_loc.axpy(q[pt], fx); q 114 src/secOrd_op_solutionDriven_impl.cpp const Quadrature& q, q 117 src/secOrd_op_solutionDriven_impl.cpp for(size_t pt = 0; pt < q.nop(); ++pt){ q 119 src/secOrd_op_solutionDriven_impl.cpp const auto& x = q.point(pt); q 120 src/secOrd_op_solutionDriven_impl.cpp const auto integral_factor = q.weight(pt) * g.integrationElement(x); q 122 src/secOrd_op_solutionDriven_impl.cpp auto X_p = evaluate(pt, q, cf); q 125 src/secOrd_op_solutionDriven_impl.cpp auto dX_p = jacobian(pt, q, cf); q 128 src/secOrd_op_solutionDriven_impl.cpp f.axpy(q[pt], X_p, dX_p); q 132 src/secOrd_op_solutionDriven_impl.cpp const Quadrature& q, q 136 src/secOrd_op_solutionDriven_impl.cpp for(size_t pt = 0; pt < q.nop(); ++pt){ q 138 src/secOrd_op_solutionDriven_impl.cpp const auto& x = q.point(pt); q 139 src/secOrd_op_solutionDriven_impl.cpp const auto integral_factor = q.weight(pt) * g.integrationElement(x); q 143 src/secOrd_op_solutionDriven_impl.cpp const auto u_p = evaluate(pt, q, u_loc); q 146 src/secOrd_op_solutionDriven_impl.cpp auto X_p = evaluate(pt, q, cf); q 150 src/secOrd_op_solutionDriven_impl.cpp auto dX_p = jacobian(pt, q, cf); q 153 src/secOrd_op_solutionDriven_impl.cpp f.axpy(q[pt], X_p, dX_p); q 201 src/secOrd_op_solutionDriven_impl.h const MCF_op::Quadrature& q, q 204 src/secOrd_op_solutionDriven_impl.h cf.evaluate(q[pt], X); q 210 src/secOrd_op_solutionDriven_impl.h const MCF_op::Quadrature& q, q 213 src/secOrd_op_solutionDriven_impl.h cf.evaluate(q[pt], u); q 219 src/secOrd_op_solutionDriven_impl.h const MCF_op::Quadrature& q, q 222 src/secOrd_op_solutionDriven_impl.h cf.jacobian(q[pt], nabla_X);