private static void test01()
//****************************************************************************80
//
// Purpose:
//
// TEST01 uses TETRAHEDRON_SAMPLE_01 to compare exact and estimated integrals.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 15 January 2014
//
// Author:
//
// John Burkardt
//
{
int[] e = new int[3];
int i;
const int m = 3;
const int n = 4192;
Console.WriteLine("");
Console.WriteLine("TEST01");
Console.WriteLine(" Estimate monomial integrals using Monte Carlo");
Console.WriteLine(" over the interior of the unit tetrahedron in 3D.");
//
// Get sample points.
//
int seed = 123456789;
double[] x = Integrals.tetrahedron01_sample(n, ref seed);
Console.WriteLine("");
Console.WriteLine(" Number of sample points used is " + n + "");
//
// Run through the exponents.
//
Console.WriteLine("");
Console.WriteLine(" Ex Ey Ez MC-Estimate Exact Error");
Console.WriteLine("");
for (i = 0; i <= 3; i++)
{
e[0] = i;
int j;
for (j = 0; j <= 3; j++)
{
e[1] = j;
int k;
for (k = 0; k <= 3; k++)
{
e[2] = k;
double[] value = Monomial.monomial_value(m, n, e, x);
double result = Integrals.tetrahedron01_volume() * typeMethods.r8vec_sum(n, value)
/ n;
double exact = Integrals.tetrahedron01_monomial_integral(e);
double error = Math.Abs(result - exact);
Console.WriteLine(" " + e[0].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + e[1].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + e[2].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + result.ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + exact.ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + error.ToString(CultureInfo.InvariantCulture).PadLeft(10) + "");
}
}
}
}
