private static void test006()
//****************************************************************************80
//
// Purpose:
//
// TEST006 tests TET_MESH_TET_NEIGHBORS.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 18 August 2009
//
// Author:
//
// John Burkardt
//
{
int node_num = 0;
int tet_num = 0;
const int tet_order = 4;
Console.WriteLine("");
Console.WriteLine("TEST006");
Console.WriteLine(" TET_MESH_TET_NEIGHBORS computes the 4 neighboring");
Console.WriteLine(" tetrahedrons of each tetrahedron in a tet mesh.");
Console.WriteLine(" containing a point.");
//
// Set up the example tetrahedron mesh.
//
TetMesh.tet_mesh_order4_example_size(ref node_num, ref tet_num);
Console.WriteLine("");
Console.WriteLine(" This mesh has tetrahedron order " + tet_order + "");
Console.WriteLine(" The number of tetrahedrons is " + tet_num + "");
double[] node_xyz = new double[3 * node_num];
int[] tet_node = new int[tet_order * tet_num];
TetMesh.tet_mesh_order4_example_set(node_num, tet_num, ref node_xyz, ref tet_node);
//
// Print the tets.
//
typeMethods.i4mat_transpose_print_some(tet_order, tet_num, tet_node,
1, 1, tet_order, 10, " First 10 Tets:");
//
// The TET_NEIGHBOR array is needed by TET_MESH_DELAUNAY_SEARCH.
//
int[] tet_neighbor = TetMesh.tet_mesh_neighbor_tets(tet_order, tet_num, tet_node);
typeMethods.i4mat_transpose_print_some(4, tet_num, tet_neighbor,
1, 1, 4, 10, " First 10 Tet Neighbors:");
}
private static void test007()
//****************************************************************************80
//
// Purpose:
//
// TEST007 tests TET_MESH_SEARCH_NAIVE.
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 19 August 2009
//
// Author:
//
// John Burkardt
//
{
int face = 0;
int node_num = 0;
double[] p = new double[3];
int step_num = 0;
int test;
const int test_num = 5;
int tet_num = 0;
const int tet_order = 4;
double[] tet_xyz = new double[3 * 4];
int seed = 123456789;
Console.WriteLine("");
Console.WriteLine("TEST007");
Console.WriteLine(" TET_MESH_SEARCH_NAIVE uses a naive algorithm");
Console.WriteLine(" to search a tetrahedral mesh for the tetrahedron");
Console.WriteLine(" containing a point.");
//
// Set up the example tetrahedron mesh.
//
TetMesh.tet_mesh_order4_example_size(ref node_num, ref tet_num);
Console.WriteLine("");
Console.WriteLine(" This mesh has tetrahedron order " + tet_order + "");
Console.WriteLine(" The number of tetrahedrons is " + tet_num + "");
double[] node_xyz = new double[3 * node_num];
int[] tet_node = new int[tet_order * tet_num];
TetMesh.tet_mesh_order4_example_set(node_num, tet_num, ref node_xyz, ref tet_node);
//
// The TET_NEIGHBOR array is needed for the Delaunay search.
//
int[] tet_neighbor = TetMesh.tet_mesh_neighbor_tets(tet_order, tet_num, tet_node);
for (test = 1; test <= test_num; test++)
{
//
// Choose a tetrahedron at random.
//
int tet1 = UniformRNG.i4_uniform_ab(0, tet_num - 1, ref seed);
Console.WriteLine("");
Console.WriteLine(" Point was chosen from tetrahedron " + tet1.ToString(CultureInfo.InvariantCulture).PadLeft(8) + "");
int j;
for (j = 0; j < 4; j++)
{
int k = tet_node[j + tet1 * 4];
int i;
for (i = 0; i < 3; i++)
{
tet_xyz[i + j * 3] = node_xyz[i + k * 3];
}
}
//
// Choose a point in the tetrahedron at random.
//
Tetrahedron.tetrahedron_sample(tet_xyz, 1, ref seed, ref p);
//
// Naive search.
//
int tet2 = TetMesh.tet_mesh_search_naive(node_num, node_xyz, tet_order, tet_num,
tet_node, p, ref step_num);
Console.WriteLine(" Naive search ended in tetrahedron " + tet2.ToString(CultureInfo.InvariantCulture).PadLeft(8)
+ ", number of steps = " + step_num + "");
//
// Delaunay search.
//
int tet3 = TetMesh.tet_mesh_search_delaunay(node_num, node_xyz, tet_order,
tet_num, tet_node, tet_neighbor, p, ref face, ref step_num);
Console.WriteLine(" Delaunay search ended in tetrahedron " + tet3.ToString(CultureInfo.InvariantCulture).PadLeft(8)
+ ", number of steps = " + step_num + "");
}
}
private static void Main(string[] args)
//****************************************************************************80
//
// Purpose:
//
// MAIN is the main program for FEM3D_PROJECT.
//
// Discussion:
//
// FEM3D_PROJECT reads files defining a sampling of a (scalar or vector)
// function of 3 arguments, and a list of nodes and tetrahedral elements
// to use for a finite element representation of the data.
//
// It computes a set of finite element coefficients to be associated with
// the given finite element mesh, and writes that information to a file
// so that an FEM representation is formed by the node, element and value
// files.
//
// Usage:
//
// fem3d_project sample_prefix fem_prefix
//
// where 'sample_prefix' is the common prefix for the SAMPLE files:
//
// * sample_prefix_nodes.txt, the node coordinates where samples were taken,
// * sample_prefix_elements.txt, the 4 nodes that make up each element;
// * sample_prefix_values.txt, the sample values.
//
// and 'fem_prefix' is the common prefix for the FEM files:
//
// * fem_prefix_nodes.txt, the node coordinates.
// * fem_prefix_elements.txt, the 4 nodes that make up each element;
// * fem_prefix_values.txt, the values defined at each node (output).
//
// Licensing:
//
// This code is distributed under the GNU LGPL license.
//
// Modified:
//
// 25 August 2009
//
// Author:
//
// John Burkardt
//
{
string fem_prefix;
int i;
int j;
string sample_prefix;
Console.WriteLine("");
Console.WriteLine("FEM3D_PROJECT");
Console.WriteLine("");
Console.WriteLine(" Read files defining a sampling of a function of 3 arguments.");
Console.WriteLine(" Read files defining a finite element mesh.");
Console.WriteLine(" Project the sample data onto the mesh, and");
Console.WriteLine(" write a file of FEM coefficient values.");
//
// Get the number of command line arguments.
//
try
{
sample_prefix = args[0];
}
catch
{
Console.WriteLine("");
Console.WriteLine("Enter the sample file prefix:");
sample_prefix = Console.ReadLine();
}
try
{
fem_prefix = args[1];
}
catch
{
Console.WriteLine("");
Console.WriteLine("Enter the FEM file prefix:");
fem_prefix = Console.ReadLine();
}
//
// Create the filenames.
//
string sample_node_filename = sample_prefix + "_nodes.txt";
string sample_element_filename = sample_prefix + "_elements.txt";
string sample_value_filename = sample_prefix + "_values.txt";
string fem_node_filename = fem_prefix + "_nodes.txt";
string fem_element_filename = fem_prefix + "_elements.txt";
string fem_value_filename = fem_prefix + "_values.txt";
//
// Read the SAMPLE NODE, ELEMENT and VALUE data.
//
TableHeader h = typeMethods.r8mat_header_read(sample_node_filename);
int sample_node_dim = h.m;
int sample_node_num = h.n;
Console.WriteLine("");
Console.WriteLine(" Sample node spatial dimension is " + sample_node_dim + "");
Console.WriteLine(" Sample node number is " + sample_node_num + "");
if (sample_node_dim != 3)
{
Console.WriteLine("");
Console.WriteLine("FEM3D_PROJECT - Fatal error!");
Console.WriteLine(" Spatial dimension of the sample nodes is not 3.");
return;
}
double[] sample_node_xyz = typeMethods.r8mat_data_read(sample_node_filename, sample_node_dim,
sample_node_num);
h = typeMethods.i4mat_header_read(sample_element_filename);
int sample_element_order = h.m;
int sample_element_num = h.n;
Console.WriteLine("");
Console.WriteLine(" Sample element order is " + sample_element_order + "");
Console.WriteLine(" Sample element number is " + sample_element_num + "");
if (sample_element_order != 4)
{
Console.WriteLine("");
Console.WriteLine("FEM3D_PROJECT - Fatal error!");
Console.WriteLine(" The sample element order must be 4.");
return;
}
int[] sample_element_node = typeMethods.i4mat_data_read(sample_element_filename,
sample_element_order, sample_element_num);
int element_min = typeMethods.i4mat_min(sample_element_order, sample_element_num,
sample_element_node);
switch (element_min)
{
case 1:
{
Console.WriteLine("");
Console.WriteLine(" Converting 1-based sample element array to 0 base.");
for (j = 0; j < sample_element_num; j++)
{
for (i = 0; i < sample_element_order; i++)
{
sample_element_node[i + j * sample_element_order] -= 1;
}
}
break;
}
}
h = typeMethods.r8mat_header_read(sample_value_filename);
int sample_value_dim = h.m;
int sample_value_num = h.n;
Console.WriteLine("");
Console.WriteLine(" The sample value dimension is " + sample_value_dim + "");
Console.WriteLine(" The sample value number is " + sample_value_num + "");
if (sample_value_num != sample_node_num)
{
Console.WriteLine("");
Console.WriteLine("FEM3D_PROJECT - Fatal error!");
Console.WriteLine(" Number of sample values and nodes differ.");
return;
}
double[] sample_value = typeMethods.r8mat_data_read(sample_value_filename, sample_value_dim,
sample_value_num);
//
// Create the sample element neighbor array.
//
int[] sample_element_neighbor = TetMesh.tet_mesh_neighbor_tets(sample_element_order,
sample_element_num, sample_element_node);
Console.WriteLine("");
Console.WriteLine(" The element neighbor array has been computed.");
//
// Read the FEM NODE and ELEMENT data.
//
h = typeMethods.r8mat_header_read(fem_node_filename);
int fem_node_dim = h.m;
int fem_node_num = h.n;
Console.WriteLine("");
Console.WriteLine(" The FEM node dimension is " + fem_node_dim + "");
Console.WriteLine(" The FEM node number is " + fem_node_num + "");
if (fem_node_dim != 3)
{
Console.WriteLine("");
Console.WriteLine("FEM3D_PROJECT - Fatal error!");
Console.WriteLine(" Spatial dimension of the nodes is not 3.");
return;
}
double[] fem_node_xyz = typeMethods.r8mat_data_read(fem_node_filename, fem_node_dim, fem_node_num);
h = typeMethods.i4mat_header_read(fem_element_filename);
int fem_element_order = h.m;
int fem_element_num = h.n;
Console.WriteLine(" The FEM element order is " + fem_element_order + "");
Console.WriteLine(" The FEM element number is " + fem_element_num + "");
if (fem_element_order != 4)
{
Console.WriteLine("");
Console.WriteLine("FEM3D_PROJECT - Fatal error!");
Console.WriteLine(" The FEM element order is not 4.");
return;
}
int[] fem_element_node = typeMethods.i4mat_data_read(fem_element_filename, fem_element_order,
fem_element_num);
element_min = typeMethods.i4mat_min(fem_element_order, fem_element_num,
fem_element_node);
switch (element_min)
{
case 1:
{
Console.WriteLine("");
Console.WriteLine(" Converting 1-based FEM element array to 0 base.");
for (j = 0; j < fem_element_num; j++)
{
for (i = 0; i < fem_element_order; i++)
{
fem_element_node[i + j * fem_element_order] -= 1;
}
}
break;
}
}
//
// Compute the FEM values.
//
double[] fem_value = FEM_3D_Transfer.fem3d_transfer(sample_node_num, sample_element_order,
sample_element_num, sample_value_dim, sample_value_num,
sample_node_xyz, sample_element_node, sample_element_neighbor, sample_value,
fem_node_num, fem_element_order,
fem_element_num, sample_value_dim, fem_node_num,
fem_node_xyz, fem_element_node);
//
// Write the FEM values.
//
typeMethods.r8mat_write(fem_value_filename, sample_value_dim, fem_node_num,
fem_value);
Console.WriteLine("");
Console.WriteLine(" FEM value data written to \"" + fem_value_filename + "\"");
//
// Terminate.
//
Console.WriteLine("");
Console.WriteLine("FEM3D_PROJECT");
Console.WriteLine(" Normal end of execution.");
Console.WriteLine("");
}
