17 June 2024 04:39:05 PM

ppc_twb_quad_minimal():
  C version.
  Test ppc_twb_quad()

demo1():
  Retrieve and print a single quadrature rule.
  strength = 7,  number of points = 15

            x                    y                    z                    w

   0      1.0000000000000      0.0000000000000      0.0000000000000      0.0102558174092
   1      0.0000000000000      0.0000000000000      1.0000000000000      0.0102558174092
   2      0.0000000000000      1.0000000000000      0.0000000000000      0.0102558174092
   3      0.7839656651012      0.0421382841642      0.1738960507346      0.1116047046647
   4      0.1738960507345      0.7839656651012      0.0421382841643      0.1116047046647
   5      0.1738960507345      0.0421382841642      0.7839656651013      0.1116047046647
   6      0.0421382841642      0.1738960507345      0.7839656651013      0.1116047046647
   7      0.7839656651012      0.1738960507345      0.0421382841643      0.1116047046647
   8      0.0421382841642      0.7839656651012      0.1738960507346      0.1116047046647
   9      0.4743880861752      0.4743880861752      0.0512238276496      0.1679775595335
  10      0.4743880861752      0.0512238276497      0.4743880861751      0.1679775595335
  11      0.0512238276497      0.4743880861752      0.4743880861751      0.1679775595335
  12      0.2385615300181      0.5228769399639      0.2385615300180      0.2652238803946
  13      0.5228769399639      0.2385615300181      0.2385615300180      0.2652238803946
  14      0.2385615300181      0.2385615300181      0.5228769399638      0.2652238803946

demo2():
  Estimate integral of x^2+y^2 over reference triangle.
  Integral estimate = 0.166667

ppc_twb_quad_minimal():
  Normal end of execution.
17 June 2024 04:39:05 PM