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 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) + ""); } }