private static void test01()
//****************************************************************************80
//
// Purpose:
//
// TEST01 tests FEM_BASIS_1D
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 05 October 2010
//
// Author:
//
// John Burkardt
//
{
int i2;
Console.WriteLine("");
Console.WriteLine("TEST01");
Console.WriteLine(" FEM_BASIS_1D evaluates an arbitrary");
Console.WriteLine(" basis function over an interval.");
const int i1 = 2;
const int j1 = 1;
int d = i1 + j1;
double x1 = typeMethods.r8_fraction(i1, d);
Console.WriteLine("");
Console.WriteLine(" I J X L(I,J)(X)");
Console.WriteLine("");
Console.WriteLine(" " + i1.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + j1.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + x1.ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + 1.0.ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
Console.WriteLine("");
for (i2 = 0; i2 <= d; i2++)
{
int j2 = d - i2;
double x2 = typeMethods.r8_fraction(i2, d);
double lij = FEM_basis.fem_basis_1d(i1, j1, x2);
Console.WriteLine(" " + i2.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + j2.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + x2.ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + lij.ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
}
}
private static void test02()
//****************************************************************************80
//
// Purpose:
//
// TEST02 tests FEM_BASIS_2D
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 05 October 2010
//
// Author:
//
// John Burkardt
//
{
int j2;
Console.WriteLine("");
Console.WriteLine("TEST02");
Console.WriteLine(" FEM_BASIS_2D evaluates an arbitrary triangular");
Console.WriteLine(" basis function.");
const int i1 = 1;
const int j1 = 0;
const int k1 = 2;
int d = i1 + j1 + k1;
double x1 = typeMethods.r8_fraction(i1, d);
double y1 = typeMethods.r8_fraction(j1, d);
Console.WriteLine("");
Console.WriteLine(" I J K X Y L(I,J,K)(X,Y)");
Console.WriteLine("");
Console.WriteLine(" " + i1.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + j1.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + k1.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + x1.ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + y1.ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + 1.0.ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
Console.WriteLine("");
for (j2 = 0; j2 <= d; j2++)
{
int i2;
for (i2 = 0; i2 <= d - j2; i2++)
{
int k2 = d - i2 - j2;
double x2 = typeMethods.r8_fraction(i2, d);
double y2 = typeMethods.r8_fraction(j2, d);
double lijk = FEM_basis.fem_basis_2d(i1, j1, k1, x2, y2);
Console.WriteLine(" " + i2.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + j2.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + k2.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + x2.ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + y2.ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + lijk.ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
}
}
}
private static void test08()
//****************************************************************************80
//
// Purpose:
//
// TEST08 tests FEM_BASIS_PRISM_TRIANGLE.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 20 October 2010
//
// Author:
//
// John Burkardt
//
{
int[] i1 =
{
2, 0, 1
}
;
int[] i2 = new int[3];
int i_0;
int[] j1 =
{
1, 0
}
;
int[] j2 = new int[2];
double[] xyz1 = new double[3];
double[] xyz2 = new double[3];
Console.WriteLine("");
Console.WriteLine("TEST08");
Console.WriteLine(" FEM_BASIS_PRISM_TRIANGLE evaluates an arbitrary");
Console.WriteLine(" basis function over a right triangular prism.");
Console.WriteLine("");
Console.WriteLine(" Here, we generate basis functions which can be");
Console.WriteLine(" up to degree 3 in X and Y, and up to degree 1 in Z.");
Console.WriteLine("");
Console.WriteLine(" Choose a node N1, define the basis function associated");
Console.WriteLine(" with that node, and then evaluate it at all other nodes.");
int di = typeMethods.i4vec_sum(3, i1);
xyz1[0] = typeMethods.r8_fraction(i1[0], di);
xyz1[1] = typeMethods.r8_fraction(i1[1], di);
int dj = typeMethods.i4vec_sum(2, j1);
xyz1[2] = typeMethods.r8_fraction(j1[0], dj);
Console.WriteLine("");
Console.WriteLine(" I1 I2 I3 J1 J2 X Y Z B(X,Y,Z)");
Console.WriteLine("");
Console.WriteLine(" " + i1[0].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + i1[1].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + i1[2].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + j1[0].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + j1[1].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + xyz1[0].ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + xyz1[1].ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + xyz1[2].ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + 1.0.ToString(CultureInfo.InvariantCulture).PadLeft(10) + "");
Console.WriteLine("");
for (i_0 = 0; i_0 <= di; i_0++)
{
i2[0] = i_0;
xyz2[0] = typeMethods.r8_fraction(i2[0], di);
int i_1;
for (i_1 = 0; i_1 <= di - i2[0]; i_1++)
{
i2[1] = i_1;
xyz2[1] = typeMethods.r8_fraction(i2[1], di);
i2[2] = di - i2[0] - i2[1];
int j_0;
for (j_0 = 0; j_0 <= dj; j_0++)
{
j2[0] = j_0;
j2[1] = dj - j2[0];
xyz2[2] = typeMethods.r8_fraction(j2[0], dj);
double b = FEM_basis.fem_basis_prism_triangle(i1, j1, xyz2);
Console.WriteLine(" " + i2[0].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + i2[1].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + i2[2].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + j2[0].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + j2[1].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + xyz2[0].ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + xyz2[1].ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + xyz2[2].ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + b.ToString(CultureInfo.InvariantCulture).PadLeft(10) + "");
}
}
}
}
private static void test06()
//****************************************************************************80
//
// Purpose:
//
// TEST06 tests FEM_BASIS_MD, repeating TEST03.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 10 October 2010
//
// Author:
//
// John Burkardt
//
{
int i;
int[] i1 = new int[4];
int[] i2 = new int[4];
const int m = 3;
int p3;
double[] x1 = new double[3];
double[] x2 = new double[3];
Console.WriteLine("");
Console.WriteLine("TEST06");
Console.WriteLine(" FEM_BASIS_MD evaluates an arbitrary");
Console.WriteLine(" basis function over an M-dimensional simplex.");
i1[0] = 1;
i1[1] = 0;
i1[2] = 2;
i1[3] = 1;
int d = typeMethods.i4vec_sum(m + 1, i1);
for (i = 0; i < m; i++)
{
x1[i] = typeMethods.r8_fraction(i1[i], d);
}
Console.WriteLine("");
Console.WriteLine(" I J K L X Y Z L(I,J,K,L)(X,Y,Z)");
Console.WriteLine("");
Console.WriteLine(" " + i1[0].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + i1[1].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + i1[2].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + i1[3].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + x1[0].ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + x1[1].ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + x1[2].ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + 1.0.ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
Console.WriteLine("");
for (p3 = 0; p3 <= d; p3++)
{
i2[2] = p3;
int p2;
for (p2 = 0; p2 <= d - p3; p2++)
{
i2[1] = p2;
int p1;
for (p1 = 0; p1 <= d - p3 - p2; p1++)
{
i2[0] = p1;
i2[3] = d - i2[0] - i2[1] - i2[2];
for (i = 0; i < m; i++)
{
x2[i] = typeMethods.r8_fraction(i2[i], d);
}
double l = FEM_basis.fem_basis_md(m, i1, x2);
Console.WriteLine(" " + i2[0].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + i2[1].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + i2[2].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + i2[3].ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + x2[0].ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + x2[1].ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + x2[2].ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + l.ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
}
}
}
}
private static void test03()
//****************************************************************************80
//
// Purpose:
//
// TEST03 tests FEM_BASIS_3D
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 05 October 2010
//
// Author:
//
// John Burkardt
//
{
int k2;
Console.WriteLine("");
Console.WriteLine("TEST03");
Console.WriteLine(" FEM_BASIS_3D evaluates an arbitrary tetrahedral");
Console.WriteLine(" basis function.");
const int i1 = 1;
const int j1 = 0;
const int k1 = 2;
const int l1 = 1;
int d = i1 + j1 + k1 + l1;
double x1 = typeMethods.r8_fraction(i1, d);
double y1 = typeMethods.r8_fraction(j1, d);
double z1 = typeMethods.r8_fraction(k1, d);
Console.WriteLine("");
Console.WriteLine(" I J K L X Y Z L(I,J,K,L)(X,Y,Z)");
Console.WriteLine("");
Console.WriteLine(" " + i1.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + j1.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + k1.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + l1.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + x1.ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + y1.ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + z1.ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + 1.0.ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
Console.WriteLine("");
for (k2 = 0; k2 <= d; k2++)
{
int j2;
for (j2 = 0; j2 <= d - k2; j2++)
{
int i2;
for (i2 = 0; i2 <= d - j2 - k2; i2++)
{
int l2 = d - i2 - j2 - k2;
double x2 = typeMethods.r8_fraction(i2, d);
double y2 = typeMethods.r8_fraction(j2, d);
double z2 = typeMethods.r8_fraction(k2, d);
double lijkl = FEM_basis.fem_basis_3d(i1, j1, k1, l1, x2, y2, z2);
Console.WriteLine(" " + i2.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + j2.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + k2.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + l2.ToString(CultureInfo.InvariantCulture).PadLeft(2)
+ " " + x2.ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + y2.ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + z2.ToString(CultureInfo.InvariantCulture).PadLeft(10)
+ " " + lijkl.ToString(CultureInfo.InvariantCulture).PadLeft(14) + "");
}
}
}
}