42   quadPts = 0.5*(gqSolver.Points().transpose()+Eigen::VectorXd::Ones(numQuadPts));  in MonotoneExpansion()
152   Eigen::VectorXd quadEvals(quadPts.size());  in EvaluateForward()
157     evalPt(i) = x(i) * quadPts(0);  in EvaluateForward()
162     for(int k=1; k<quadPts.size(); ++k){  in EvaluateForward()
163       evalPt(i) = x(i) * quadPts(k);  in EvaluateForward()
201     for(int k=0; k<quadPts.size(); ++k){  in JacobianWrtX()
202       evalPt(i) = x(i) * quadPts(k);  in JacobianWrtX()
206       jac(i,i) += quadWeights(k) * (polyEval * polyEval + 2.0*x(i)*polyGrad(i)*polyEval*quadPts(k));  in JacobianWrtX()
252       for(int k=0; k<quadPts.size(); ++k){  in JacobianImpl()
253         evalPt(i) = x(i) * quadPts(k);  in JacobianImpl()
317     for(int k=0; k<quadPts.size(); ++k){  in GradLogDeterminant()
318       evalPt(i) = x(i) * quadPts(k);  in GradLogDeterminant()
324       part1 += quadWeights(k) * (polyEval * polyEval + 2.0*x(i)*polyGradX(i)*polyEval*quadPts(k));  in GradLogDeterminant()
325 …part2 += 2.0 * quadWeights(k) * (polyEval * polyGradC + x(i) * quadPts(k) * (polyEval*polyD2.row(i…  in GradLogDeterminant()