<div dir="ltr"><div><span style="font-family:arial,helvetica,sans-serif">Dear Dune,<br><br></span></div><span style="font-family:arial,helvetica,sans-serif">I got to have a question about the quadrature rules implemented in Dune. There seems to be something strange in the Dune tutorial. First, according to the comments, if one aims to do integration with quadrature order p, the rule should be initialized as:<br>
<br></span><div dir="ltr"><span style="font-family:arial,helvetica,sans-serif"><font>const Dune::QuadratureRule <ctype ,dim >& rule = Dune::QuadratureRules <ctype ,dim >::rule(gt,p);<br><br></font></span></div>
<div><span style="font-family:arial,helvetica,sans-serif"><font>However, later for calculation with order two the p value is set to 1. Would you explain a little about this?<br><br></font></span></div><div><span style="font-family:arial,helvetica,sans-serif"><font>My main question is that if I want to consider only one point for the integration (first order quadrature rule), what should I set for the rule order? one or zero?<br>
<br></font></span></div><div><span style="font-family:arial,helvetica,sans-serif"><font>Thank you in advance for your help and best regards,<br></font></span></div><div style="font-family:monospace"><font><span style="font-family:arial,helvetica,sans-serif">Arya</span><br>
</font></div><br></div>