private static void test06()
//****************************************************************************80
//
// Purpose:
//
// TEST06 tests KEAST_RULE.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 26 June 2007
//
// Author:
//
// John Burkardt
//
{
int order;
Console.WriteLine("");
Console.WriteLine("TEST06");
Console.WriteLine(" KEAST_RULE returns the points and weights");
Console.WriteLine(" of a Keast rule for the triangle.");
Console.WriteLine("");
Console.WriteLine(" In this test, we simply print a rule.");
const int rule = 10;
int degree = KeastRule.keast_degree(rule);
int order_num = KeastRule.keast_order_num(rule);
Console.WriteLine("");
Console.WriteLine(" Rule = " + rule + "");
Console.WriteLine(" Degree = " + degree + "");
Console.WriteLine(" Order = " + order_num + "");
Console.WriteLine("");
Console.WriteLine(" I W X Y Z");
Console.WriteLine("");
double[] xyz = new double[3 * order_num];
double[] w = new double[order_num];
KeastRule.keast_rule(rule, order_num, ref xyz, ref w);
for (order = 0; order < order_num; order++)
{
Console.WriteLine(" " + order.ToString().PadLeft(8)
+ " " + w[order].ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + xyz[0 + order * 3].ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + xyz[1 + order * 3].ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + xyz[2 + order * 3].ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
}
}
private static void test02()
//****************************************************************************80
//
// Purpose:
//
// TEST02 tests KEAST_RULE.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 14 December 2006
//
// Author:
//
// John Burkardt
//
{
int rule;
Console.WriteLine("");
Console.WriteLine("TEST02");
Console.WriteLine(" KEAST_RULE returns the points and weights");
Console.WriteLine(" of a Keast rule for the triangle.");
Console.WriteLine("");
Console.WriteLine(" In this test, we simply check that the weights");
Console.WriteLine(" sum to 1.");
int rule_num = KeastRule.keast_rule_num();
Console.WriteLine("");
Console.WriteLine(" Number of available rules = " + rule_num + "");
Console.WriteLine("");
Console.WriteLine(" Rule Order Sum of weights");
Console.WriteLine("");
for (rule = 1; rule <= rule_num; rule++)
{
int order_num = KeastRule.keast_order_num(rule);
double[] xyztab = new double[3 * order_num];
double[] wtab = new double[order_num];
KeastRule.keast_rule(rule, order_num, ref xyztab, ref wtab);
double wtab_sum = 0.0;
int order;
for (order = 0; order < order_num; order++)
{
wtab_sum += wtab[order];
}
Console.WriteLine(" " + rule.ToString().PadLeft(8)
+ " " + order_num.ToString().PadLeft(8)
+ " " + wtab_sum.ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
}
}
private static void test05()
//****************************************************************************80
//
// Purpose:
//
// TEST05 tests KEAST_RULE.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 03 July 2008
//
// Author:
//
// John Burkardt
//
{
int degree;
const int degree_max = 3;
const int dim_num = 3;
int[] expon = new int[3];
Console.WriteLine("");
Console.WriteLine("TEST05");
Console.WriteLine(" Demonstrate the KEAST rules on a sequence of");
Console.WriteLine(" monomial integrands X^A Y^B Z^C");
Console.WriteLine(" on the unit tetrahedron.");
int rule_num = KeastRule.keast_rule_num();
Console.WriteLine("");
Console.WriteLine(" Rule Order Quad");
Console.WriteLine("");
for (degree = 0; degree <= degree_max; degree++)
{
bool more = false;
int h = 0;
int t = 0;
for (;;)
{
Comp.comp_next(degree, dim_num, ref expon, ref more, ref h, ref t);
Console.WriteLine("");
Console.WriteLine(" F(X,Y,Z)"
+ " = X^" + expon[0]
+ " * Y^" + expon[1]
+ " * Z^" + expon[2] + "");
Console.WriteLine("");
int rule;
for (rule = 1; rule <= rule_num; rule++)
{
int order_num = KeastRule.keast_order_num(rule);
double[] xyz = new double[3 * order_num];
double[] w = new double[order_num];
KeastRule.keast_rule(rule, order_num, ref xyz, ref w);
double[] mono = Monomial.monomial_value(dim_num, order_num, xyz, expon);
double quad = typeMethods.r8vec_dot(order_num, w, mono);
Console.WriteLine(" " + rule.ToString().PadLeft(8)
+ " " + order_num.ToString().PadLeft(8)
+ " " + quad.ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
}
if (!more)
{
break;
}
}
}
}
private static void test04()
//****************************************************************************80
//
// Purpose:
//
// TEST04 demonstrates TETRAHEDRON_REFERENCE_TO_PHYSICAL.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 14 December 2006
//
// Author:
//
// John Burkardt
//
{
const int NODE_NUM = 4;
int node;
double[] node_xyz =
{
0.0, 0.0, 0.0,
1.0, 0.0, 0.0,
0.0, 1.0, 0.0,
0.0, 0.0, 1.0
};
double[] node_xyz2 =
{
1.0, 2.0, 3.0,
2.0, 2.0, 3.0,
1.0, 3.0, 3.0,
1.0, 2.0, 9.0
};
int order;
Console.WriteLine("");
Console.WriteLine("TEST04");
Console.WriteLine(" TETRAHEDRON_REFERENCE_TO_PHYSICAL transforms a rule");
Console.WriteLine(" on the unit (reference) tetrahedron to a rule on ");
Console.WriteLine(" an arbitrary (physical) tetrahedron.");
const int rule = 2;
int order_num = KeastRule.keast_order_num(rule);
double[] xyz = new double[3 * order_num];
double[] xyz2 = new double[3 * order_num];
double[] w = new double[order_num];
KeastRule.keast_rule(rule, order_num, ref xyz, ref w);
//
// Here is the reference tetrahedron, and its rule.
//
Console.WriteLine("");
Console.WriteLine(" The reference tetrahedron:");
Console.WriteLine("");
for (node = 0; node < NODE_NUM; node++)
{
Console.WriteLine(" " + (node + 1).ToString().PadLeft(8)
+ " " + node_xyz[0 + node * 3].ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + node_xyz[1 + node * 3].ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + node_xyz[2 + node * 3].ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
}
double volume = Tetrahedron.tetrahedron_volume(node_xyz);
Console.WriteLine("");
Console.WriteLine(" Rule " + rule + " for reference tetrahedron");
Console.WriteLine(" with volume = " + volume + "");
Console.WriteLine("");
Console.WriteLine(" X Y Z W");
Console.WriteLine("");
for (order = 0; order < order_num; order++)
{
Console.WriteLine(" " + order.ToString().PadLeft(8)
+ " " + xyz[0 + order * 3].ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + xyz[1 + order * 3].ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + xyz[2 + order * 3].ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + w[order].ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
}
//
// Transform the rule.
//
Tetrahedron.reference_to_physical_t4(node_xyz2, order_num, xyz, ref xyz2);
//
// Here is the physical tetrahedron, and its transformed rule.
//
Console.WriteLine("");
Console.WriteLine(" The physical tetrahedron:");
Console.WriteLine("");
for (node = 0; node < NODE_NUM; node++)
{
Console.WriteLine(" " + (node + 1).ToString().PadLeft(8)
+ " " + node_xyz2[0 + node * 3].ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + node_xyz2[1 + node * 3].ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + node_xyz2[2 + node * 3].ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
}
double volume2 = Tetrahedron.tetrahedron_volume(node_xyz2);
Console.WriteLine("");
Console.WriteLine(" Rule " + rule + " for physical tetrahedron");
Console.WriteLine(" with volume = " + volume2 + "");
Console.WriteLine("");
Console.WriteLine(" X Y Z W");
Console.WriteLine("");
for (order = 0; order < order_num; order++)
{
Console.WriteLine(" " + order.ToString().PadLeft(8)
+ " " + xyz2[0 + order * 3].ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + xyz2[1 + order * 3].ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + xyz2[2 + order * 3].ToString(CultureInfo.InvariantCulture).PadLeft(14)
+ " " + w[order].ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
}
}