private static void test01()
//****************************************************************************80
//
// Purpose:
//
// TEST01 tests KEAST_RULE_NUM, KEAST_DEGREE, KEAST_ORDER_NUM.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 14 December 2006
//
// Author:
//
// John Burkardt
//
{
int rule;
Console.WriteLine("");
Console.WriteLine("TEST01");
Console.WriteLine(" KEAST_RULE_NUM returns the number of rules;");
Console.WriteLine(" KEAST_DEGREE returns the degree of a rule;");
Console.WriteLine(" KEAST_ORDER_NUM returns the order of a rule.");
int rule_num = KeastRule.keast_rule_num();
Console.WriteLine("");
Console.WriteLine(" Number of available rules = " + rule_num + "");
Console.WriteLine("");
Console.WriteLine(" Rule Degree Order");
Console.WriteLine("");
for (rule = 1; rule <= rule_num; rule++)
{
int order_num = KeastRule.keast_order_num(rule);
int degree = KeastRule.keast_degree(rule);
Console.WriteLine(" " + rule.ToString().PadLeft(8)
+ " " + degree.ToString().PadLeft(8)
+ " " + order_num.ToString().PadLeft(8) + "");
}
}
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 test03()
//****************************************************************************80
//
// Purpose:
//
// TEST03 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("TEST03");
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, for each");
Console.WriteLine(" quadrature point, the barycentric coordinates");
Console.WriteLine(" add up to 1.");
int rule_num = KeastRule.keast_rule_num();
Console.WriteLine("");
Console.WriteLine(" Rule Suborder Sum of coordinates");
Console.WriteLine("");
for (rule = 1; rule <= rule_num; rule++)
{
int suborder_num = KeastRule.keast_suborder_num(rule);
double[] suborder_xyzz = new double[4 * suborder_num];
double[] suborder_w = new double[suborder_num];
KeastRule.keast_subrule(rule, suborder_num, ref suborder_xyzz, ref suborder_w);
Console.WriteLine("");
Console.WriteLine(" " + rule.ToString().PadLeft(8)
+ " " + suborder_num.ToString().PadLeft(8) + "");
int suborder;
for (suborder = 0; suborder < suborder_num; suborder++)
{
double xyzz_sum = suborder_xyzz[0 + suborder * 4]
+ suborder_xyzz[1 + suborder * 4]
+ suborder_xyzz[2 + suborder * 4]
+ suborder_xyzz[3 + suborder * 4];
Console.WriteLine(" "
+ " " + xyzz_sum.ToString("0.################").PadLeft(25) + "");
}
}
}