x 53 old_code/2015linfty/dune_bdf.hxx double x = vec(0), y = vec(1), z = vec(2), t = get_time(); x 55 old_code/2015linfty/dune_bdf.hxx 16.*M_PI*pow(x,3)*cos(2*M_PI*t)/( (pow(y,2) + pow(z,2) + 16.*pow(x,2)/pow(sin(2.*M_PI*t) + 4.,2) )*pow(sin(2.*M_PI*t) + 4.,3) ), x 56 old_code/2015linfty/dune_bdf.hxx 4.*M_PI*pow(x,2)*y*cos(2*M_PI*t)/((pow(y,2) + pow(z,2) + 16.*pow(x,2)/pow(sin(2.*M_PI*t) + 4.,2) )*pow(sin(2.*M_PI*t) + 4.,2)), x 57 old_code/2015linfty/dune_bdf.hxx 4.*M_PI*pow(x,2)*z*cos(2.*M_PI*t)/((pow(y,2) + pow(z,2) + 16.*pow(x,2)/pow(sin(2*M_PI*t) + 4.,2) )*pow(sin(2.*M_PI*t) + 4,2) ) ); x 62 old_code/2015linfty/dune_bdf.hxx double x = vec(0), y = vec(1), z = vec(2), t = get_time(); x 66 old_code/2015linfty/dune_bdf.hxx 48*M_PI*pow(x,2)*cos(2*M_PI*t)/((pow(y,2) + pow(z,2) + 16*pow(x,2)/pow(sin(2*M_PI*t) + 4,2))*pow(sin(2*M_PI*t) + 4,3)) - 512*M_PI*pow(x,4)*cos(2*M_PI*t)/(pow(pow(y,2) + pow(z,2) + 16*pow(x,2)/pow(sin(2*M_PI*t) + 4,2),2)*pow(sin(2*M_PI*t) + 4,5)), x 67 old_code/2015linfty/dune_bdf.hxx -32*M_PI*pow(x,3)*y*cos(2*M_PI*t)/(pow(pow(y,2) + pow(z,2) + 16*pow(x,2)/pow(sin(2*M_PI*t) + 4,2),2)*pow(sin(2*M_PI*t) + 4,3)), x 68 old_code/2015linfty/dune_bdf.hxx -32*M_PI*pow(x,3)*z*cos(2*M_PI*t)/(pow(pow(y,2) + pow(z,2) + 16*pow(x,2)/pow(sin(2*M_PI*t) + 4,2),2)*pow(sin(2*M_PI*t) + 4,3)), x 70 old_code/2015linfty/dune_bdf.hxx 8*M_PI*x*y*cos(2*M_PI*t)/((pow(y,2) + pow(z,2) + 16*pow(x,2)/pow(sin(2*M_PI*t) + 4,2))*pow(sin(2*M_PI*t) + 4,2)) - 128*M_PI*pow(x,3)*y*cos(2*M_PI*t)/(pow(pow(y,2) + pow(z,2) + 16*pow(x,2)/pow(sin(2*M_PI*t) + 4,2),2)*pow(sin(2*M_PI*t) + 4,4)), x 71 old_code/2015linfty/dune_bdf.hxx -8*M_PI*pow(x,2)*pow(y,2)*cos(2*M_PI*t)/(pow(pow(y,2) + pow(z,2) + 16*pow(x,2)/pow(sin(2*M_PI*t) + 4,2),2)*pow(sin(2*M_PI*t) + 4,2)) + 4*M_PI*pow(x,2)*cos(2*M_PI*t)/((pow(y,2) + pow(z,2) + 16*pow(x,2)/pow(sin(2*M_PI*t) + 4,2))*pow(sin(2*M_PI*t) + 4,2)), x 72 old_code/2015linfty/dune_bdf.hxx -8*M_PI*pow(x,2)*y*z*cos(2*M_PI*t)/(pow(pow(y,2) + pow(z,2) + 16*pow(x,2)/pow(sin(2*M_PI*t) + 4,2),2)*pow(sin(2*M_PI*t) + 4,2)), x 74 old_code/2015linfty/dune_bdf.hxx 8*M_PI*x*z*cos(2*M_PI*t)/((pow(y,2) + pow(z,2) + 16*pow(x,2)/pow(sin(2*M_PI*t) + 4,2))*pow(sin(2*M_PI*t) + 4,2)) - 128*M_PI*pow(x,3)*z*cos(2*M_PI*t)/(pow(pow(y,2) + pow(z,2) + 16*pow(x,2)/pow(sin(2*M_PI*t) + 4,2),2)*pow(sin(2*M_PI*t) + 4,4)), x 75 old_code/2015linfty/dune_bdf.hxx -8*M_PI*pow(x,2)*y*z*cos(2*M_PI*t)/(pow(pow(y,2) + pow(z,2) + 16*pow(x,2)/pow(sin(2*M_PI*t) + 4,2),2)*pow(sin(2*M_PI*t) + 4,2)), x 76 old_code/2015linfty/dune_bdf.hxx -8*M_PI*pow(x,2)*pow(z,2)*cos(2*M_PI*t)/(pow(pow(y,2) + pow(z,2) + 16*pow(x,2)/pow(sin(2*M_PI*t) + 4,2),2)*pow(sin(2*M_PI*t) + 4,2)) + 4*M_PI*pow(x,2)*cos(2*M_PI*t)/((pow(y,2) + pow(z,2) + 16*pow(x,2)/pow(sin(2*M_PI*t) + 4,2))*pow(sin(2*M_PI*t) + 4,2)); x 96 old_code/2015linfty/dune_bdf.hxx double x = vec(0), y = vec(1), z = vec(2), t = get_time(); x 98 old_code/2015linfty/dune_bdf.hxx .5 * M_PI * cos(2. * M_PI * t) / ( 2. * sqrt( 1. + .25 * sin(2. * M_PI * t) ) ) *x, 0, 0 ); x 103 old_code/2015linfty/dune_bdf.hxx double x = vec(0), y = vec(1), z = vec(2), t = get_time(); x 32 old_code/tumor_growth/tumor_deformation.h void evaluate (const DomainVector& x, RangeVector& y) const{ x 44 old_code/tumor_growth/tumor_deformation.h std::array<double, 3> domain_node { {x[0], x[1], x[2]} }; x 570 old_code/tumor_growth/tumor_growth.h const typename QuadratureType::CoordinateType& x = quadrature.point(pt); x 571 old_code/tumor_growth/tumor_growth.h const double weight = quadrature.weight(pt) * geometry.integrationElement(x); x 627 old_code/tumor_growth/tumor_growth.h const typename QuadratureType::CoordinateType& x = quadrature.point(pt); x 628 old_code/tumor_growth/tumor_growth.h const double weight = quadrature.weight(pt) * geometry.integrationElement(x); x 680 old_code/tumor_growth/tumor_growth.h const typename QuadratureType::CoordinateType& x = quadrature.point(pt); x 681 old_code/tumor_growth/tumor_growth.h const double weight = quadrature.weight(pt) * geometry.integrationElement(x); x 707 old_code/tumor_growth/tumor_growth.h const typename QuadratureType::CoordinateType& x = bigger_quadrature.point(pt); x 709 old_code/tumor_growth/tumor_growth.h = bigger_quadrature.weight(pt) * geometry.integrationElement(x); x 742 old_code/tumor_growth/tumor_growth.h const typename QuadratureType::CoordinateType& x = quadrature.point(pt); x 743 old_code/tumor_growth/tumor_growth.h const double weight = quadrature.weight(pt) * geometry.integrationElement(x); x 773 old_code/tumor_growth/tumor_growth.h const typename QuadratureType::CoordinateType& x = quadrature.point(pt); x 774 old_code/tumor_growth/tumor_growth.h const double weight = quadrature.weight(pt) * geometry.integrationElement(x); x 37 src/grid_FEfunSet.h Scal_FEfun x; x 44 src/grid_FEfunSet.h :x {name, gt}, y {name, gt}, z {name, gt} {} x 39 src/grid_deformation.cpp static inline void identity(const Domain& x, Range& y) noexcept{ x 40 src/grid_deformation.cpp y[0] = x[0]; x 41 src/grid_deformation.cpp y[1] = x[1]; x 42 src/grid_deformation.cpp y[2] = x[2]; x 52 src/grid_deformation.cpp static inline void logistic_growth(const double t, const Domain& x, Range& y) noexcept{ x 55 src/grid_deformation.cpp y = x; x 60 src/grid_deformation.cpp static inline void dalquist(const double t, const Domain& x, Range& y){ x 62 src/grid_deformation.cpp y = x; x 67 src/grid_deformation.cpp static inline void mcf_sphere(const double t, const Domain& x, Range& y){ x 71 src/grid_deformation.cpp = std::inner_product(&x[0], &x[0] + Domain::dimension, &x[0], 0.); x 73 src/grid_deformation.cpp y = x; x 102 src/grid_deformation.cpp void Esfem::Grid::Deformation::evaluate(const Domain& x, Range& y) const{ x 113 src/grid_deformation.cpp y = d_ptr->hg[x]; x 29 src/grid_evolution.cpp evaluate(const Domain& x, Range& y) const{ x 31 src/grid_evolution.cpp y[0] = x[0]; x 32 src/grid_evolution.cpp y[1] = x[1]; x 33 src/grid_evolution.cpp y[2] = x[2]; x 170 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); x 171 src/secOrd_op_brusselator_impl.cpp const auto weight = q.weight(pt) * g.integrationElement(x); x 187 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); x 188 src/secOrd_op_brusselator_impl.cpp const auto weight = q.weight(pt) * g.integrationElement(x); x 220 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); x 221 src/secOrd_op_brusselator_impl.cpp const auto weight = q.weight(pt) * g.integrationElement(x); x 242 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); x 244 src/secOrd_op_brusselator_impl.cpp = q.weight(pt) * g.integrationElement(x) x 264 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); x 265 src/secOrd_op_brusselator_impl.cpp const auto integral_factor = q.weight(pt) * geometry.integrationElement(x); x 298 src/secOrd_op_brusselator_impl.cpp const auto& x = q.point(pt); x 300 src/secOrd_op_brusselator_impl.cpp = q.weight(pt) * geometry.integrationElement(x) x 92 src/secOrd_op_initData_impl.cpp const double x = d[0]; x 95 src/secOrd_op_initData_impl.cpp r = x * y * std::exp(-6. * t); x 141 src/secOrd_op_initData_impl.cpp void sphere_1EF::evaluate(const domain& x, range& y)const{ x 142 src/secOrd_op_initData_impl.cpp y = x[0]*x[1]*exp(-6*tp.time()); x 157 src/secOrd_op_initData_impl.cpp void sphere_eigenFun::evaluate(const Domain& x, Range& y) const{ x 158 src/secOrd_op_initData_impl.cpp y[0] = x[0]*x[1]; // xy x 159 src/secOrd_op_initData_impl.cpp y[1] = x[1]*x[2]; // yz x 160 src/secOrd_op_initData_impl.cpp y[2] = x[0]*x[2]; // xz x 176 src/secOrd_op_initData_impl.cpp void sphere_mcf_sol::evaluate(const Domain& x, Range& y) const{ x 178 src/secOrd_op_initData_impl.cpp for(int i = 0; i < Domain::dimension; ++i) norm += x[i]*x[i]; x 181 src/secOrd_op_initData_impl.cpp y = x; x 198 src/secOrd_op_initData_impl.cpp void sls_iData::evaluate(const Domain& x, Range& y) const{ x 200 src/secOrd_op_initData_impl.cpp norm = sqrt(inner_product(&x[0], &x[0]+Domain::dimension, &x[0], 0.)), x 203 src/secOrd_op_initData_impl.cpp y = x; x 217 src/secOrd_op_initData_impl.cpp void sls_v_iData::evaluate(const Domain& x, Range& y) const{ x 221 src/secOrd_op_initData_impl.cpp y = x; x 230 src/secOrd_op_initData_impl.cpp void sd_iData::evaluate(const Domain& x, Range& y) const{ x 232 src/secOrd_op_initData_impl.cpp norm = sqrt(inner_product(&x[0], &x[0]+Domain::dimension, &x[0], 0.)), x 234 src/secOrd_op_initData_impl.cpp y = x; x 87 src/secOrd_op_initData_u.cpp const double x = d[0]; x 90 src/secOrd_op_initData_u.cpp r = std::exp(-6.*t)*x*y; x 177 src/secOrd_op_linearHeat.cpp const auto& x = q.point(pt); x 178 src/secOrd_op_linearHeat.cpp const auto weight = q.weight(pt) * g.integrationElement(x); x 196 src/secOrd_op_linearHeat.cpp const auto& x = q.point(pt); x 197 src/secOrd_op_linearHeat.cpp const auto weight = q.weight(pt) * g.integrationElement(x); x 72 src/secOrd_op_rhs_impl.cpp const double x = d[0]; x 76 src/secOrd_op_rhs_impl.cpp r = 2*k*x*y*(rA/(rA*(exp(-k*t) - 1) - rE*exp(-k*t)) + 1)*exp(-6*t) - (k*(pow(x,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1)*(rA/(rA*(exp(-k*t) - 1) - rE*exp(-k*t)) + 1) + k*(pow(y,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1)*(rA/(rA*(exp(-k*t) - 1) - rE*exp(-k*t)) + 1) + k*(pow(z,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1)*(rA/(rA*(exp(-k*t) - 1) - rE*exp(-k*t)) + 1))*x*y*exp(-6*t) + 4*pow(x,3)*y*pow(z,2)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),3) + 4*x*pow(y,3)*pow(z,2)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),3) - 6*x*y*exp(-6*t) + 2*(pow(x,3)*y*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) + (pow(x,3)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - x/(pow(x,2) + pow(y,2) + pow(z,2)))*y*exp(-6*t) - x*y*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)))*(pow(x,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1) + (4*pow(x,2)*pow(y,2)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - pow(x,2)*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)) - (pow(x,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1)*exp(-6*t))*x*y/(pow(x,2) + pow(y,2) + pow(z,2)) + (4*pow(x,2)*pow(y,2)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - pow(y,2)*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)) - (pow(y,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1)*exp(-6*t))*x*y/(pow(x,2) + pow(y,2) + pow(z,2)) + 2*(x*pow(y,3)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) + x*(pow(y,3)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - y/(pow(x,2) + pow(y,2) + pow(z,2)))*exp(-6*t) - x*y*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)))*(pow(y,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1) + 2*(2*pow(x,2)*y*z*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - y*z*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)))*x*z/(pow(x,2) + pow(y,2) + pow(z,2)) + 2*(2*x*pow(y,2)*z*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - x*z*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)))*y*z/(pow(x,2) + pow(y,2) + pow(z,2)) + 2*(2*x*y*pow(z,2)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - x*y*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)))*(pow(z,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1); x 82 src/secOrd_op_rhs_impl.cpp const double x = d[0]; x 86 src/secOrd_op_rhs_impl.cpp r = 2*k*y*z*(rA/(rA*(exp(-k*t) - 1) - rE*exp(-k*t)) + 1)*exp(-6*t) - (k*(pow(x,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1)*(rA/(rA*(exp(-k*t) - 1) - rE*exp(-k*t)) + 1) + k*(pow(y,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1)*(rA/(rA*(exp(-k*t) - 1) - rE*exp(-k*t)) + 1) + k*(pow(z,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1)*(rA/(rA*(exp(-k*t) - 1) - rE*exp(-k*t)) + 1))*y*z*exp(-6*t) - 6*y*z*exp(-6*t) + (4*pow(x,2)*pow(y,3)*z*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),3) + 4*pow(x,2)*y*pow(z,3)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),3) + 2*(2*pow(x,2)*y*z*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - y*z*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)))*(pow(x,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1) + 2*(2*x*pow(y,2)*z*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - x*z*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)))*x*y/(pow(x,2) + pow(y,2) + pow(z,2)) + 2*(pow(y,3)*z*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) + (pow(y,3)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - y/(pow(x,2) + pow(y,2) + pow(z,2)))*z*exp(-6*t) - y*z*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)))*(pow(y,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1) + 2*(2*x*y*pow(z,2)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - x*y*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)))*x*z/(pow(x,2) + pow(y,2) + pow(z,2)) + (4*pow(y,2)*pow(z,2)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - pow(y,2)*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)) - (pow(y,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1)*exp(-6*t))*y*z/(pow(x,2) + pow(y,2) + pow(z,2)) + (4*pow(y,2)*pow(z,2)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - pow(z,2)*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)) - (pow(z,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1)*exp(-6*t))*y*z/(pow(x,2) + pow(y,2) + pow(z,2)) + 2*(y*pow(z,3)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) + y*(pow(z,3)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - z/(pow(x,2) + pow(y,2) + pow(z,2)))*exp(-6*t) - y*z*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)))*(pow(z,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1))*Dc; x 165 src/secOrd_op_rhs_impl.cpp auto u_fun = [&tp = tp](auto x, auto y){ x 166 src/secOrd_op_rhs_impl.cpp return x * y * exp(-6.*tp.time()); x 209 src/secOrd_op_rhs_impl.cpp const auto x = d[0], y = d[1], z = d[2], t = tp.time(); x 210 src/secOrd_op_rhs_impl.cpp return 2*k*x*y*(rA/(rA*(exp(-k*t) - 1) - rE*exp(-k*t)) + 1)*exp(-6*t) - (k*(pow(x,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1)*(rA/(rA*(exp(-k*t) - 1) - rE*exp(-k*t)) + 1) + k*(pow(y,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1)*(rA/(rA*(exp(-k*t) - 1) - rE*exp(-k*t)) + 1) + k*(pow(z,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1)*(rA/(rA*(exp(-k*t) - 1) - rE*exp(-k*t)) + 1))*x*y*exp(-6*t) + 4*pow(x,3)*y*pow(z,2)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),3) + 4*x*pow(y,3)*pow(z,2)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),3) - 6*x*y*exp(-6*t) + 2*(pow(x,3)*y*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) + (pow(x,3)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - x/(pow(x,2) + pow(y,2) + pow(z,2)))*y*exp(-6*t) - x*y*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)))*(pow(x,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1) + (4*pow(x,2)*pow(y,2)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - pow(x,2)*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)) - (pow(x,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1)*exp(-6*t))*x*y/(pow(x,2) + pow(y,2) + pow(z,2)) + (4*pow(x,2)*pow(y,2)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - pow(y,2)*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)) - (pow(y,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1)*exp(-6*t))*x*y/(pow(x,2) + pow(y,2) + pow(z,2)) + 2*(x*pow(y,3)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) + x*(pow(y,3)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - y/(pow(x,2) + pow(y,2) + pow(z,2)))*exp(-6*t) - x*y*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)))*(pow(y,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1) + 2*(2*pow(x,2)*y*z*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - y*z*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)))*x*z/(pow(x,2) + pow(y,2) + pow(z,2)) + 2*(2*x*pow(y,2)*z*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - x*z*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)))*y*z/(pow(x,2) + pow(y,2) + pow(z,2)) + 2*(2*x*y*pow(z,2)*exp(-6*t)/pow(pow(x,2) + pow(y,2) + pow(z,2),2) - x*y*exp(-6*t)/(pow(x,2) + pow(y,2) + pow(z,2)))*(pow(z,2)/(pow(x,2) + pow(y,2) + pow(z,2)) - 1); x 326 src/secOrd_op_rhs_impl.h const auto& x = quad.point(pt); x 327 src/secOrd_op_rhs_impl.h Range fx {rhs(geometry.global(x))}; x 329 src/secOrd_op_rhs_impl.h fx *= quad.weight(pt) * geometry.integrationElement(x); x 114 src/secOrd_op_rhs_u.cpp const double x = d[0]; x 119 src/secOrd_op_rhs_u.cpp const double u = std::exp(-6. * t) * x * y; x 153 src/secOrd_op_rhs_u.cpp const auto& x = q.point(pt); x 154 src/secOrd_op_rhs_u.cpp RHS_data_u::Range fx {rhs_fun(g.global(x))}; x 155 src/secOrd_op_rhs_u.cpp fx *= q.weight(pt) * g.integrationElement(x); x 113 src/secOrd_op_rhs_w.cpp const double x = dom[0]; x 146 src/secOrd_op_rhs_w.cpp const auto& x = q.point(pt); x 147 src/secOrd_op_rhs_w.cpp RHS_data_w::Range fx {rhs_fun(g.global(x))}; x 149 src/secOrd_op_rhs_w.cpp fx *= q.weight(pt) * g.integrationElement(x); x 119 src/secOrd_op_solutionDriven_impl.cpp const auto& x = q.point(pt); x 120 src/secOrd_op_solutionDriven_impl.cpp const auto integral_factor = q.weight(pt) * g.integrationElement(x); x 138 src/secOrd_op_solutionDriven_impl.cpp const auto& x = q.point(pt); x 139 src/secOrd_op_solutionDriven_impl.cpp const auto integral_factor = q.weight(pt) * g.integrationElement(x);