<div dir="ltr"><div><div><div><div><div><div><div>Dear Christian,<br><br></div>Thank you very much for your response. It is now clear for me.<br><br></div>For the tutorial statement I referred to, look at the following lines:<br>
<br></div>page 34, line 25 in the presented code<br></div>page 35, line 37 in the presented code<br></div>page 36, the following statement: "Running the executable integration on a YaspGrid in two space dimensions with a quadrature rule of order two the following output is obtained:"<br>
<br></div>Thanks and best regards,<br></div>Arya<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Feb 25, 2014 at 3:40 PM, Christian Engwer <span dir="ltr"><<a href="mailto:christian.engwer@uni-muenster.de" target="_blank">christian.engwer@uni-muenster.de</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Arya,<br>
<div class=""><br>
> I got to have a question about the quadrature rules implemented in Dune.<br>
> There seems to be something strange in the Dune tutorial. First, according<br>
> to the comments, if one aims to do integration with quadrature order p, the<br>
> rule should be initialized as:<br>
><br>
> const Dune::QuadratureRule <ctype ,dim >& rule = Dune::QuadratureRules<br>
> <ctype ,dim >::rule(gt,p);<br>
><br>
> However, later for calculation with order two the p value is set to 1.<br>
> Would you explain a little about this?<br>
<br>
</div>Please be a bit more specific. Where in the tutorial is this<br>
statement. We need to have a closer look...<br>
<div class=""><br>
> My main question is that if I want to consider only one point for the<br>
> integration (first order quadrature rule), what should I set for the rule<br>
> order? one or zero?<br>
<br>
</div>In general the order you specify is the order for which you compute<br>
the exact integral (modulo fp errors). In particular this means that<br>
integrating a first order function you set the rule to 1. The obtained<br>
rule will be the same as for zero order, but this a detail you don't<br>
have to think about.<br>
<br>
Ciao<br>
<br>
Christian<br>
<div class=""><br>
><br>
> Thank you in advance for your help and best regards,<br>
> Arya<br>
<br>
</div>> _______________________________________________<br>
> Dune mailing list<br>
> <a href="mailto:Dune@dune-project.org">Dune@dune-project.org</a><br>
> <a href="http://lists.dune-project.org/mailman/listinfo/dune" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a><br>
<span class="HOEnZb"><font color="#888888"><br>
<br>
--<br>
Prof. Dr. Christian Engwer<br>
Institut für Numerische und Angewandte Mathematik<br>
Fachbereich Mathematik und Informatik der Universität Münster<br>
Einsteinstrasse 62<br>
48149 Münster<br>
<br>
E-Mail  <a href="mailto:christian.engwer@uni-muenster.de">christian.engwer@uni-muenster.de</a><br>
Telefon <a href="tel:%2B49%20251%2083-35067" value="+492518335067">+49 251 83-35067</a><br>
FAX             <a href="tel:%2B49%20251%2083-32729" value="+492518332729">+49 251 83-32729</a><br>
</font></span></blockquote></div><br></div>