Commit f4a89168 authored by ghweber's avatar ghweber
Browse files

Testing and debugging of AMR stitch cell operator.

git-svn-id: http://visit.ilight.com/svn/visit/trunk/src@12403 18c085ea-50e0-402c-830e-de6fd14e8384
parent f3e58cfc
......@@ -77,8 +77,8 @@ int tesselationArray3D[5776] = {
-1, 5, 0, 4, 5, 1, 3, 5, 4, 7, 6, 5, 3, 5, 5, 6, 2, 1, 3, -1, 6, 0, 3, 4, 1,
2, 5, 6, 3, 7, 4, 2, 6, 5, -1, 5, 2, 5, 4, 3, 0, 5, 2, 3, 4, 5, 6, -1, 6, 3,
7, 4, 2, 6, 5, 5, 2, 5, 4, 3, 0, 4, 0, 1, 2, 5, -1, 6, 5, 6, 2, 4, 7, 3, 5,
4, 3, 2, 5, 1, 4, 1, 0, 4, 3, -1, 6, 3, 7, 4, 2, 6, 5, 6, 0, 3, 4, 1, 2, 5,
-1, 6, 3, 7, 4, 2, 6, 5, 5, 2, 5, 4, 3, 0, -1, 5, 5, 1, 0, 4, 6, 5, 2, 3, 0,
4, 3, 2, 5, 1, 4, 1, 0, 4, 3, -1, 6, 0, 3, 4, 1, 2, 5, 6, 3, 7, 4, 2, 6, 5,
-1, 6, 0, 3, 4, 1, 2, 5, 5, 2, 3, 4, 5, 7, -1, 5, 5, 1, 0, 4, 6, 5, 2, 3, 0,
1, 6, 5, 4, 0, 3, 7, 6, -1, 4, 5, 1, 0, 6, 4, 5, 0, 4, 6, 4, 1, 2, 0, 6, 4,
0, 2, 3, 6, 4, 6, 7, 3, 0, 4, 4, 7, 6, 0, -1, 5, 0, 3, 7, 4, 1, 4, 5, 4, 6,
1, 4, 4, 7, 6, 1, 4, 6, 3, 2, 1, 4, 6, 7, 3, 1, -1, 5, 1, 0, 4, 5, 6, 5, 1,
......@@ -122,14 +122,14 @@ int tesselationArray3D[5776] = {
6, 5, 3, 6, 5, 0, 1, 4, 1, 2, 3, 6, -1, 5, 5, 4, 7, 6, 3, 5, 2, 3, 0, 1, 5,
4, 5, 0, 4, 3, 4, 2, 6, 3, 5, -1, 6, 6, 7, 3, 5, 4, 0, 5, 5, 0, 3, 6, 2, 4,
2, 1, 5, 0, -1, 5, 2, 1, 5, 6, 3, 5, 1, 0, 4, 5, 3, 5, 5, 4, 7, 6, 3, -1, 6,
0, 4, 5, 3, 7, 6, 6, 1, 0, 5, 2, 3, 6, -1, 6, 0, 4, 5, 3, 7, 6, 5, 3, 6, 5,
0, 1, -1, 4, 6, 4, 7, 3, 7, 2, 1, 5, 6, 3, 0, 4, -1, 5, 0, 3, 7, 4, 6, 5, 5,
1, 0, 5, 2, 3, 6, 6, 0, 4, 5, 3, 7, 6, -1, 6, 1, 0, 5, 2, 3, 6, 5, 3, 0, 5,
6, 4, -1, 4, 6, 4, 7, 3, 7, 2, 1, 5, 6, 3, 0, 4, -1, 5, 0, 3, 7, 4, 6, 5, 5,
6, 2, 1, 0, 4, 0, 2, 3, 6, 4, 5, 4, 6, 0, -1, 4, 4, 0, 3, 1, 4, 4, 3, 7, 6,
4, 1, 6, 2, 3, 4, 5, 6, 1, 4, 4, 6, 3, 1, 4, -1, 5, 5, 1, 0, 4, 6, 5, 1, 2,
3, 0, 6, 5, 0, 3, 7, 4, 6, -1, 4, 6, 5, 4, 2, 5, 5, 1, 0, 4, 2, 5, 6, 7, 3,
2, 4, 4, 2, 3, 0, 4, -1, 6, 4, 7, 6, 0, 3, 2, 5, 0, 2, 6, 4, 5, 4, 5, 1, 0,
2, -1, 5, 5, 6, 2, 1, 4, 5, 6, 7, 3, 2, 4, 5, 2, 3, 0, 1, 4, -1, 6, 4, 7, 6,
0, 3, 2, 6, 5, 4, 6, 1, 0, 2, -1, 6, 4, 7, 6, 0, 3, 2, 5, 0, 2, 6, 4, 5, -1,
2, -1, 5, 5, 6, 2, 1, 4, 5, 6, 7, 3, 2, 4, 5, 2, 3, 0, 1, 4, -1, 6, 5, 4, 6,
1, 0, 2, 6, 4, 7, 6, 0, 3, 2, -1, 6, 5, 4, 6, 1, 0, 2, 5, 0, 4, 6, 2, 7, -1,
4, 7, 6, 4, 3, 4, 6, 5, 4, 3, 5, 6, 2, 1, 5, 3, 5, 5, 1, 0, 4, 3, -1, 6, 5,
1, 0, 6, 2, 3, 5, 6, 3, 0, 5, 4, 4, 4, 7, 6, 3, -1, 4, 4, 0, 3, 1, 4, 4, 3,
7, 6, 4, 1, 6, 2, 3, 4, 5, 6, 1, 4, 4, 6, 3, 1, 4, -1, 5, 5, 1, 0, 4, 3, 4,
......@@ -165,8 +165,8 @@ int tesselationArray3D[5776] = {
1, -1, 5, 0, 3, 7, 4, 1, 5, 3, 2, 6, 7, 1, 5, 7, 6, 5, 4, 1, -1, 4, 3, 7, 4,
2, 5, 7, 6, 5, 4, 2, 5, 3, 0, 1, 2, 4, 4, 2, 1, 5, 4, -1, 5, 3, 7, 4, 0, 2,
5, 7, 6, 5, 4, 2, 5, 4, 5, 1, 0, 2, -1, 6, 2, 6, 7, 1, 5, 4, 5, 1, 4, 7, 2,
3, 4, 3, 0, 1, 4, -1, 6, 4, 5, 1, 7, 6, 2, 6, 0, 4, 1, 3, 7, 2, -1, 6, 4, 5,
1, 7, 6, 2, 5, 7, 2, 1, 4, 0, -1, 5, 1, 2, 3, 0, 5, 4, 4, 0, 7, 5, 4, 0, 3,
3, 4, 3, 0, 1, 4, -1, 6, 0, 4, 1, 3, 7, 2, 6, 4, 5, 1, 7, 6, 2, -1, 6, 0, 4,
1, 3, 7, 2, 5, 7, 4, 1, 2, 5, -1, 5, 1, 2, 3, 0, 5, 4, 4, 0, 7, 5, 4, 0, 3,
7, 5, 4, 7, 2, 6, 5, 4, 7, 3, 2, 5, -1, 4, 0, 1, 2, 5, 4, 0, 2, 3, 7, 4, 5,
7, 6, 2, 4, 4, 7, 5, 0, 4, 7, 2, 5, 0, -1, 5, 5, 1, 0, 4, 7, 5, 5, 6, 2, 1,
7, 4, 7, 0, 3, 1, 4, 1, 2, 3, 7, -1, 4, 4, 5, 0, 7, 4, 5, 1, 0, 7, 5, 5, 6,
......@@ -191,8 +191,8 @@ int tesselationArray3D[5776] = {
1, 0, 4, 5, 7, 4, 0, 2, 3, 7, 4, 1, 2, 0, 7, 4, 1, 5, 2, 7, 4, 5, 6, 2, 7,
-1, 5, 3, 2, 6, 7, 5, 5, 4, 5, 1, 0, 3, 4, 3, 1, 2, 5, 4, 4, 7, 5, 3, -1, 5,
4, 0, 3, 7, 5, 5, 0, 1, 2, 3, 5, 5, 3, 2, 6, 7, 5, -1, 6, 7, 6, 5, 3, 2, 1,
5, 3, 1, 5, 7, 4, 4, 4, 0, 3, 1, -1, 6, 7, 6, 5, 3, 2, 1, 6, 4, 7, 5, 0, 3,
1, -1, 6, 7, 6, 5, 3, 2, 1, 5, 3, 1, 5, 7, 4, -1, 6, 4, 0, 3, 5, 1, 2, 5, 5,
5, 3, 1, 5, 7, 4, 4, 4, 0, 3, 1, -1, 6, 4, 7, 5, 0, 3, 1, 6, 7, 6, 5, 3, 2,
1, -1, 6, 4, 7, 5, 0, 3, 1, 5, 3, 7, 5, 1, 6, -1, 6, 4, 0, 3, 5, 1, 2, 5, 5,
2, 3, 4, 7, 4, 7, 6, 5, 2, -1, 5, 4, 0, 3, 7, 2, 4, 0, 5, 1, 2, 4, 4, 5, 0,
2, 4, 4, 7, 5, 2, 4, 7, 6, 5, 2, -1, 4, 4, 0, 3, 1, 4, 4, 3, 7, 2, 4, 4, 7,
5, 2, 4, 7, 6, 5, 2, 4, 4, 5, 1, 2, 4, 3, 1, 2, 4, -1, 4, 7, 6, 5, 2, 7, 1,
......@@ -209,8 +209,8 @@ int tesselationArray3D[5776] = {
6, 5, 1, 5, 6, 1, 0, 7, 3, 4, 3, 2, 6, 1, -1, 8, 2, 6, 7, 3, 1, 5, 4, 0, -1,
6, 6, 5, 7, 2, 1, 3, -1, 6, 2, 1, 6, 3, 4, 7, -1, 5, 2, 6, 7, 3, 1, -1, 5, 3,
2, 6, 7, 0, 5, 2, 1, 5, 6, 0, 5, 6, 5, 4, 7, 0, -1, 5, 2, 6, 7, 3, 1, 5, 6,
5, 4, 7, 1, 5, 7, 4, 0, 3, 1, -1, 6, 1, 5, 6, 0, 4, 7, 6, 2, 1, 6, 3, 0, 7,
-1, 6, 1, 5, 6, 0, 4, 7, 5, 0, 7, 6, 1, 2, -1, 4, 4, 6, 5, 1, 7, 0, 3, 7, 4,
5, 4, 7, 1, 5, 7, 4, 0, 3, 1, -1, 6, 2, 1, 6, 3, 0, 7, 6, 1, 5, 6, 0, 4, 7,
-1, 6, 2, 1, 6, 3, 0, 7, 5, 0, 1, 6, 7, 5, -1, 4, 4, 6, 5, 1, 7, 0, 3, 7, 4,
1, 2, 6, -1, 4, 4, 7, 6, 0, 5, 7, 3, 2, 6, 0, 5, 4, 5, 1, 0, 6, 4, 0, 1, 2,
6, -1, 4, 5, 4, 6, 1, 4, 4, 7, 6, 1, 5, 4, 0, 3, 7, 1, 5, 7, 3, 2, 6, 1, -1,
6, 0, 3, 7, 1, 2, 6, 5, 1, 6, 7, 0, 4, 4, 4, 5, 1, 6, -1, 5, 4, 0, 3, 7, 2,
......@@ -222,8 +222,8 @@ int tesselationArray3D[5776] = {
0, 1, 2, 3, 6, -1, 4, 6, 2, 1, 3, 4, 6, 1, 5, 4, 4, 3, 4, 0, 1, 4, 7, 4, 3,
6, 4, 4, 1, 3, 6, -1, 5, 0, 3, 7, 4, 6, 4, 3, 1, 2, 6, 4, 0, 1, 3, 6, 4, 0,
4, 1, 6, 4, 4, 5, 1, 6, -1, 5, 7, 3, 2, 6, 4, 5, 3, 0, 1, 2, 4, 5, 2, 1, 5,
6, 4, -1, 6, 0, 1, 2, 4, 5, 6, 6, 3, 0, 2, 7, 4, 6, -1, 6, 0, 1, 2, 4, 5, 6,
5, 4, 6, 2, 0, 3, -1, 5, 2, 6, 7, 3, 4, 4, 6, 1, 5, 4, 4, 2, 1, 6, 4, 4, 2,
6, 4, -1, 6, 3, 0, 2, 7, 4, 6, 6, 0, 1, 2, 4, 5, 6, -1, 6, 3, 0, 2, 7, 4, 6,
5, 4, 0, 2, 6, 1, -1, 5, 2, 6, 7, 3, 4, 4, 6, 1, 5, 4, 4, 2, 1, 6, 4, 4, 2,
3, 1, 4, 4, 3, 0, 1, 4, -1, 4, 4, 5, 1, 6, 7, 2, 3, 0, 1, 6, 7, 4, -1, 4, 7,
5, 4, 0, 7, 3, 2, 6, 7, 0, 1, 5, -1, 4, 4, 7, 5, 0, 4, 7, 6, 5, 0, 5, 7, 3,
2, 6, 0, 5, 6, 2, 1, 5, 0, -1, 5, 7, 3, 2, 6, 1, 5, 5, 1, 0, 4, 7, 4, 7, 0,
......@@ -237,21 +237,21 @@ int tesselationArray3D[5776] = {
1, 3, 7, 5, 6, 4, 6, 2, 1, 3, -1, 4, 1, 2, 3, 6, 4, 1, 3, 0, 7, 4, 1, 0, 5,
7, 4, 0, 4, 5, 7, 4, 1, 5, 6, 7, 4, 3, 6, 7, 1, -1, 5, 6, 7, 3, 2, 5, 5, 7,
4, 0, 3, 5, 5, 3, 0, 1, 2, 5, -1, 5, 3, 2, 6, 7, 5, 4, 2, 0, 1, 5, 4, 3, 0,
2, 5, 4, 3, 7, 0, 5, 4, 7, 4, 0, 5, -1, 6, 5, 4, 7, 1, 0, 3, 6, 6, 5, 7, 2,
1, 3, -1, 6, 5, 4, 7, 1, 0, 3, 5, 1, 3, 7, 5, 6, -1, 4, 0, 4, 5, 7, 7, 6, 2,
2, 5, 4, 3, 7, 0, 5, 4, 7, 4, 0, 5, -1, 6, 6, 5, 7, 2, 1, 3, 6, 5, 4, 7, 1,
0, 3, -1, 6, 6, 5, 7, 2, 1, 3, 5, 1, 5, 7, 3, 4, -1, 4, 0, 4, 5, 7, 7, 6, 2,
1, 5, 7, 3, 0, -1, 8, 5, 4, 7, 6, 1, 0, 3, 2, -1, 6, 5, 1, 4, 6, 3, 7, -1, 6,
4, 0, 7, 5, 1, 6, -1, 5, 5, 4, 7, 6, 1, -1, 5, 4, 7, 6, 5, 0, 5, 7, 3, 2, 6,
0, 5, 6, 2, 1, 5, 0, -1, 5, 5, 4, 7, 6, 1, 5, 4, 0, 3, 7, 1, 5, 7, 3, 2, 6,
1, -1, 6, 6, 2, 1, 7, 3, 0, 6, 5, 6, 1, 4, 7, 0, -1, 6, 6, 2, 1, 7, 3, 0, 5,
7, 0, 1, 6, 5, -1, 5, 6, 5, 4, 7, 2, 5, 5, 1, 0, 4, 2, 5, 4, 0, 3, 7, 2, -1,
1, -1, 6, 5, 6, 1, 4, 7, 0, 6, 6, 2, 1, 7, 3, 0, -1, 6, 5, 6, 1, 4, 7, 0, 5,
7, 6, 1, 0, 2, -1, 5, 6, 5, 4, 7, 2, 5, 5, 1, 0, 4, 2, 5, 4, 0, 3, 7, 2, -1,
5, 4, 7, 6, 5, 2, 4, 7, 0, 3, 2, 4, 4, 0, 7, 2, 4, 4, 5, 0, 2, 4, 5, 1, 0,
2, -1, 6, 7, 3, 2, 4, 0, 1, 6, 6, 7, 2, 5, 4, 1, -1, 6, 7, 3, 2, 4, 0, 1, 5,
4, 1, 2, 7, 6, -1, 4, 2, 3, 0, 7, 7, 4, 5, 1, 0, 7, 6, 2, -1, 5, 7, 6, 5, 4,
3, 5, 6, 2, 1, 5, 3, 5, 5, 1, 0, 4, 3, -1, 6, 5, 1, 0, 6, 2, 3, 6, 4, 5, 0,
7, 6, 3, -1, 6, 5, 1, 0, 6, 2, 3, 5, 6, 3, 0, 5, 4, -1, 5, 5, 4, 7, 6, 3, 4,
2, -1, 6, 6, 7, 2, 5, 4, 1, 6, 7, 3, 2, 4, 0, 1, -1, 6, 6, 7, 2, 5, 4, 1, 5,
4, 7, 2, 1, 3, -1, 4, 2, 3, 0, 7, 7, 4, 5, 1, 0, 7, 6, 2, -1, 5, 7, 6, 5, 4,
3, 5, 6, 2, 1, 5, 3, 5, 5, 1, 0, 4, 3, -1, 6, 4, 5, 0, 7, 6, 3, 6, 5, 1, 0,
6, 2, 3, -1, 6, 4, 5, 0, 7, 6, 3, 5, 6, 5, 0, 3, 1, -1, 5, 5, 4, 7, 6, 3, 4,
4, 1, 0, 3, 4, 5, 1, 4, 3, 4, 5, 6, 1, 3, 4, 6, 2, 1, 3, -1, 4, 6, 2, 1, 3,
7, 0, 4, 5, 1, 3, 7, 6, -1, 6, 4, 0, 3, 5, 1, 2, 6, 7, 4, 3, 6, 5, 2, -1, 6,
4, 0, 3, 5, 1, 2, 5, 5, 2, 3, 4, 7, -1, 4, 5, 1, 0, 2, 7, 3, 7, 4, 0, 2, 6,
7, 0, 4, 5, 1, 3, 7, 6, -1, 6, 7, 4, 3, 6, 5, 2, 6, 4, 0, 3, 5, 1, 2, -1, 6,
7, 4, 3, 6, 5, 2, 5, 5, 4, 3, 2, 0, -1, 4, 5, 1, 0, 2, 7, 3, 7, 4, 0, 2, 6,
5, -1, 4, 1, 0, 4, 3, 7, 7, 6, 5, 4, 3, 2, 1, -1, 8, 0, 1, 2, 3, 4, 5, 6, 7,
-1
};
......
......@@ -1081,11 +1081,19 @@ avtAMRStitchCellFilter::CreateStitchCells(vtkRectilinearGrid *rgrid,
int numVtcs = tesselationArray3D[tessPos++];
vtkIdType vPtId[numVtcs];
int tessPosSave = tessPos;
bool skipCell = false;
for (int vtxNo = 0; vtxNo < numVtcs; ++ vtxNo)
if (skipCellVtx[tesselationArray3D[tessPos++]])
continue; // Skip on to next nexx
tessPos = tessPosSave;
if (skipCellVtx[tesselationArray3D[tessPos+vtxNo]])
{
skipCell = true;
break;
}
if (skipCell)
{
tessPos += numVtcs;
continue; // Skip on to next cell
}
for (int vtxNo =0; vtxNo < numVtcs; ++vtxNo)
{
......@@ -1128,42 +1136,85 @@ avtAMRStitchCellFilter::CreateStitchCells(vtkRectilinearGrid *rgrid,
outputData->InsertNextTuple1(
inputData->GetTuple1(vtkIdType(LOC(vI, vJ, vK))));
#if 0
if ((vI == -1) || (vI == dims[0]-2))
if (((vI == -1) || (vI == dims[0]-2)) &&
(vJ != -1) && (vJ != dims[1]-2) &&
(vK != -1) && (vK != dims[2]-2))
{
if ((vJ != -1) && (vJ != dims[1]-2))
{
if ((vK != -1) && (vK != dims[2]-2))
int baseJ = vJ - (vJ % refinementRatio[1]);
int baseK = vK - (vK % refinementRatio[2]);
for (int refJVtxNo = 0; refJVtxNo < refinementRatio[1];
++refJVtxNo)
for (int refKVtxNo = 0; refKVtxNo < refinementRatio[2];
++refKVtxNo)
{
int baseJ = vJ - (vJ % refinementRatio[1]);
int baseK = vK - (vK % refinementRatio[2]);
for (int refJVtxNo = 0; refJVtxNo < refinementRatio[1];
++refJVtxNo)
for (int refKVtxNo = 0; refKVtxNo < refinementRatio[2];
++refKVtxNo)
{
pointIds[LOC(vI, baseJ+refVtxNo, baseK+refKVtxNo)] =
pointIds[LOC(vI,vJ)];
}
pointIds[LOC(vI, baseJ+refJVtxNo, baseK+refKVtxNo)] =
pointIds[LOC(vI,vJ,vK)];
}
}
else
{
assert ((vJ == -1) || (vJ == dims[1]-2));
}
else if (((vJ == -1) || (vJ == dims[1]-2)) &&
(vI != -1) && (vI != dims[0]-2) &&
(vK != -1) && (vK != dims[2]-2))
{
int baseI = vI - (vI % refinementRatio[0]);
int baseK = vK - (vK % refinementRatio[2]);
if ((vI != -1) && (vI != dims[0]-2))
for (int refIVtxNo = 0; refIVtxNo < refinementRatio[0];
++refIVtxNo)
for (int refKVtxNo = 0; refKVtxNo < refinementRatio[2];
++refKVtxNo)
{
int baseI = vI - (vI % refinementRatio[0]);
for (int refVtxNo = 0; refVtxNo < refinementRatio[0];
++refVtxNo)
{
pointIds[LOC(baseI+refVtxNo,vJ)]=pointIds[LOC(vI,vJ)];
}
pointIds[LOC(baseI+refIVtxNo, vJ, baseK+refKVtxNo)] =
pointIds[LOC(vI,vJ,vK)];
}
}
}
#endif
else if (((vK == -1) || (vK == dims[2]-2)) &&
(vI != -1) && (vI != dims[0]-2) &&
(vJ != -1) && (vJ != dims[1]-2))
{
int baseI = vI - (vI % refinementRatio[0]);
int baseJ = vJ - (vJ % refinementRatio[1]);
for (int refIVtxNo = 0; refIVtxNo < refinementRatio[0];
++refIVtxNo)
for (int refJVtxNo = 0; refJVtxNo < refinementRatio[1];
++refJVtxNo)
{
pointIds[LOC(baseI+refIVtxNo, baseJ+refJVtxNo, vK)] =
pointIds[LOC(vI,vJ,vK)];
}
}
else if (((vI == -1) || (vI == dims[0]-2)) &&
((vJ == -1) || (vJ == dims[1]-2)) &&
(vK != -1) && (vK != dims[2]-2))
{
int baseK = vK - (vK % refinementRatio[2]);
for (int refKVtxNo = 0; refKVtxNo < refinementRatio[2];
++refKVtxNo)
pointIds[LOC(vI, vJ, baseK+refKVtxNo)] =
pointIds[LOC(vI,vJ,vK)];
}
else if (((vI == -1) || (vI == dims[0]-2) ||
(vK == -1) || (vK == dims[2]-2)) &&
(vJ != -1) && (vJ != dims[1]-2))
{
int baseJ = vJ - (vJ % refinementRatio[1]);
for (int refJVtxNo = 0; refJVtxNo < refinementRatio[1];
++refJVtxNo)
pointIds[LOC(vI, baseJ+refJVtxNo, vK)] =
pointIds[LOC(vI,vJ,vK)];
}
else if (((vJ == -1) || (vJ == dims[1]-2) ||
(vK == -1) || (vK == dims[2]-2)) &&
(vI != -1) && (vI != dims[0]-2))
{
int baseI = vI - (vI % refinementRatio[0]);
for (int refIVtxNo = 0; refIVtxNo < refinementRatio[0];
++refIVtxNo)
pointIds[LOC(baseI+refIVtxNo, vJ, vK)] =
pointIds[LOC(vI,vJ,vK)];
}
}
}
......@@ -1171,6 +1222,10 @@ avtAMRStitchCellFilter::CreateStitchCells(vtkRectilinearGrid *rgrid,
vPtId[vtxNo] = pointIds[LOC(vI, vJ, vK)];
}
#if 0
for (int x=0; x<numVtcs; ++x) std::cout << vPtId[x] << " ";
std::cout << std::endl;
#endif
if (numVtcs == 4)
ugrid->InsertNextCell(VTK_TETRA, 4, vPtId);
else if (numVtcs == 5)
......@@ -1213,6 +1268,56 @@ avtAMRStitchCellFilter::CreateStitchCells(vtkRectilinearGrid *rgrid,
}
}
}
#if 0
std::cout << "IJKmin" << std::endl;
for (int i=-1; i<dims[0]-1; ++i)
{
for (int j=-1; j<dims[1]-1; ++j)
std::cout << std::setw(6) << pointIds[LOC(i,j,-1)] << " ";
std::cout << std::endl;
}
std::cout << std::endl;
std::cout << "IJKmax" << std::endl;
for (int i=-1; i<dims[0]-1; ++i)
{
for (int j=-1; j<dims[1]-1; ++j)
std::cout << std::setw(6) << pointIds[LOC(i,j,dims[2]-2)] << " ";
std::cout << std::endl;
}
std::cout << std::endl;
std::cout << "IJminK" << std::endl;
for (int i=-1; i<dims[0]-1; ++i)
{
for (int k=-1; k<dims[2]-1; ++k)
std::cout << std::setw(6) << pointIds[LOC(i,-1,k)] << " ";
std::cout << std::endl;
}
std::cout << std::endl;
std::cout << "IJmaxK" << std::endl;
for (int i=-1; i<dims[0]-1; ++i)
{
for (int k=-1; k<dims[2]-1; ++k)
std::cout << std::setw(6) << pointIds[LOC(i,dims[1]-2,k)] << " ";
std::cout << std::endl;
}
std::cout << std::endl;
std::cout << "IminJK" << std::endl;
for (int j=-1; j<dims[1]-1; ++j)
{
for (int k=-1; k<dims[2]-1; ++k)
std::cout << std::setw(6) << pointIds[LOC(-1,j,k)] << " ";
std::cout << std::endl;
}
std::cout << std::endl;
std::cout << "ImaxJK" << std::endl;
for (int j=-1; j<dims[1]-1; ++j)
{
for (int k=-1; k<dims[2]-1; ++k)
std::cout << std::setw(6) << pointIds[LOC(dims[0]-2,j,k)] << " ";
std::cout << std::endl;
}
std::cout << std::endl;
#endif
}
// Clean-up
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment