private double[,] readTuple3NameSpecificPointDataArray(vtkUnstructuredGrid unstructuredGrid, string dataArrayName, bool readExtremeForces)
{
double[,] dataArray = new double[, ] {
};
var pointData = unstructuredGrid.GetPointData();
ExtremeDisplacement = new double[3];
for (int i = 0; i < pointData.GetNumberOfArrays(); i++)
{
if (pointData.GetArrayName(i) != dataArrayName)
{
continue;
}
dataArray = new double[pointData.GetArray(i).GetNumberOfTuples(), 3];
for (int j = 0; j < pointData.GetArray(i).GetNumberOfTuples(); j++)
{
var tuple = pointData.GetArray(i).GetTuple3(j);
dataArray[j, 0] = tuple[0]; dataArray[j, 1] = tuple[1]; dataArray[j, 2] = tuple[2];
if (!readExtremeForces)
{
continue;
}
for (int k = 0; k < 3; k++)
{
if (tuple[k] > ExtremeDisplacement[k])
{
ExtremeDisplacement[k] = tuple[k];
}
}
}
}
return(dataArray);
}
public vtkDataArray getNameSpecificDataArrayPointData(vtkUnstructuredGrid unstructuredGrid, string arrayName)
{
vtkDataArray dataArray = null;
var pointData = unstructuredGrid.GetPointData();
for (int i = 0; i < pointData.GetNumberOfArrays(); i++)
{
if (pointData.GetArrayName(i) != arrayName)
{
continue;
}
dataArray = pointData.GetArray(i);
}
return(dataArray);
}
public static UnstructuredModel FromGrid(vtkUnstructuredGrid output, int numArrays)
{
var result = new UnstructuredModel
{
_Desc = DimensionDescription.FromPoints(output.GetPoints()),
_Data = new List <PointData>()
};
for (var i = 0; i < numArrays; ++i)
{
var pdata = PointData.FromArray((vtkFloatArray)output.GetPointData().GetArray(i));
pdata.DimensionDescription = result._Desc;
result._Data.Add(pdata);
}
return(result);
}
private void Simulate(CMshTriangleMesherDouble aMesher, vtkUnstructuredGrid anUnstructuredGrid, vtkRenderer aRenderer)
{
var T = new CPdeFieldDouble();
T.setMesh(aMesher.mesh());
T.setDiscretMethod(EDiscretMethod.DM_FEM);
T.setDiscretOrder(EDiscretOrder.DO_LINEAR);
T.setDiscretLocation(EDiscretLocation.DL_NODE);
T.setSimulationType(ESimulationType.ST_THERMAL);
T.setNbDofs(1);
// Set BC and initial conditions
T.setSize(T.mesh().nbNodes());
for (int i = 0; i < T.mesh().nbNodes(); ++i)
{
int aNodeId = T.mesh().nodeId(i);
var aNode = T.mesh().node(aNodeId);
if (Math.Abs(aNode.x() + 0.0) < 1E-6)
{
T.setFixedValue(i, 0.0);
}
if (Math.Abs(aNode.x() - 1.0) < 1E-6)
{
T.setFixedValue(i, 0.0);
}
if (Math.Abs(aNode.y() + 0.0) < 1E-6)
{
T.setFixedValue(i, 0.0);
}
if (Math.Abs(aNode.y() - 1.0) < 1E-6)
{
T.setFixedValue(i, 0.0);
}
T.setValue(i, 1.0);
}
double rho = 1.0; // density
double Cp = 1.0; // specific heat
double k = 1.0; // conductivity
double time = 0.1;
int nIter = 10;
double dt = time / nIter;
// Unsteady conduction in a plate
for (int i = 0; i < nIter; ++i)
{
Debug.WriteLine("Time = " + dt * (i + 1));
var aTransientTerm = new CSleSystemDouble();
aTransientTerm = ENigMA.ddt(T);
aTransientTerm.Multiply(rho * Cp / dt);
var aDiffusionTerm = new CSleSystemDouble();
aDiffusionTerm = ENigMA.laplacian(T);
aDiffusionTerm.Multiply(k);
var aSystem = new CSleSystemDouble();
aSystem = aTransientTerm;
aSystem.Plus(aDiffusionTerm);
var aPdeEquation = new CPdeEquationDouble(aSystem);
aPdeEquation.setElimination(T);
aPdeEquation.solve(T);
}
var sScalars = vtkFloatArray.New();
sScalars.SetNumberOfTuples(T.mesh().nbNodes());
for (int i = 0; i < T.mesh().nbNodes(); ++i)
{
int aControlVolumeId = T.mesh().nodeId(i);
sScalars.SetTuple1(aControlVolumeId, T.value(i));
}
anUnstructuredGrid.GetPointData().SetScalars(sScalars);
var aGeometryFilter = vtkGeometryFilter.New();
aGeometryFilter.SetInput(anUnstructuredGrid);
aGeometryFilter.Update();
var aBoundedFilter = vtkBandedPolyDataContourFilter.New();
aBoundedFilter.SetInput(aGeometryFilter.GetOutput());
aBoundedFilter.GenerateValues(24, sScalars.GetRange()[0], sScalars.GetRange()[1]);
var aLookupTable = vtkLookupTable.New();
aLookupTable.SetNumberOfColors(256);
aLookupTable.SetHueRange(0.667, 0.0);
aLookupTable.Build();
var aBandedMapper = vtkPolyDataMapper.New();
aBandedMapper.SetInputConnection(aBoundedFilter.GetOutputPort());
aBandedMapper.SetScalarModeToUsePointData();
aBandedMapper.SetScalarRange(sScalars.GetRange()[0], sScalars.GetRange()[1]);
aBandedMapper.SetLookupTable(aLookupTable);
// add actor to the renderer
var MeshActor = vtkActor.New();
aRenderer.AddActor(MeshActor);
var MeshMapper = vtkDataSetMapper.New();
MeshMapper.SetInput(anUnstructuredGrid);
MeshActor.SetMapper(aBandedMapper);
}
/// <summary>
/// 左键释放时获取选择的对象
/// </summary>
/// <param name="sender"></param>
/// <param name="e"></param>
void m_interactor_LeftButtonReleaseEvt(vtkObject sender, vtkObjectEventArgs e)
{
if (m_pickTarget == ePickTarget.None)
{
return;
}
//if (m_iren.GetControlKey() != 0)
// m_pickBoolean = ePickBollean.Add;
//else if (m_iren.GetShiftKey() != 0)
// m_pickBoolean = ePickBollean.Sub;
m_iren.GetEventPosition(ref endx, ref endy);
#region 点选对象
if (m_pickMode == ePickMode.DotPickMode) //如果是点选模式
{
vtkCellPicker cellPicker = new vtkCellPicker();
cellPicker.SetTolerance(0.005);
if (0 == cellPicker.Pick(endx, endy, 0, m_renderer))
{
return; //什么都没有选到就返回
}
FEMActor actor = cellPicker.GetActor() as FEMActor;
if (actor == null)
{
return;
}
if (m_pickTarget == ePickTarget.Actor)
{ //如果选择对象为
actor.IsSelected = true;
m_RWCtrl.Invalidate(); //貌似必须要用这个函数,update 不行
OnObjectsSelected(sender, e);
}
else if (m_pickTarget == ePickTarget.Cell)
{
int cellId = cellPicker.GetCellId();
actor.AddCellId(cellId, m_pickBoolean);
OnObjectsSelected(sender, e);
//if (cellId > 0)
// Console.WriteLine("dotselect cellid:{0}", cellId);
}
else if (m_pickTarget == ePickTarget.Point)
{
int pointId = cellPicker.GetPointId();
if (pointId != -1)
{
actor.AddPointId(pointId, m_pickBoolean);
OnObjectsSelected(sender, e);
// Console.WriteLine("dotselect pointid:{0}", pointId);
}
}
}
#endregion
else if (m_pickMode == ePickMode.RectPickMode) //如果是框选
{
vtkAreaPicker areaPicker = new vtkAreaPicker();
if (0 == areaPicker.AreaPick(startx, starty, endx, endy, m_renderer))
{ //如果什么也没选到就返回
//this.impactViewForm.ClearAllSelectCells();
return;
}
vtkProp3DCollection selectedActors = areaPicker.GetProp3Ds();
selectedActors.InitTraversal();
FEMActor selectedActor;
while ((selectedActor = selectedActors.GetNextProp3D() as FEMActor) != null)
{
vtkSelection selection = new vtkSelection();
vtkSelectionNode selectionNode = new vtkSelectionNode();
selection.AddNode(selectionNode);
vtkIdFilter idFilter = new vtkIdFilter();
idFilter.FieldDataOn();
idFilter.SetInput(selectedActor.GetMapper().GetInput());
idFilter.Update();
vtkExtractSelection frustum = new vtkExtractSelection();
frustum.SetInput(0, idFilter.GetOutput());
frustum.SetInput(1, selection);
selectionNode.Initialize();
selectionNode.SetContentType((int)vtkSelectionNode.SelectionContent.FRUSTUM);
if (m_pickTarget == ePickTarget.Cell)//如果选择的目标是cell
{
selectionNode.SetFieldType((int)vtkSelectionNode.SelectionField.CELL);
}
else if (m_pickTarget == ePickTarget.Point) //如果选择的目标是point
{
selectionNode.SetFieldType((int)vtkSelectionNode.SelectionField.POINT);
}
vtkPoints points = areaPicker.GetClipPoints();
vtkDoubleArray frustumcorners = new vtkDoubleArray();
frustumcorners.SetNumberOfComponents(4);
frustumcorners.SetNumberOfTuples(8);
for (int i = 0; i < 8; ++i)
{
double[] point = points.GetPoint(i);
frustumcorners.SetTuple4(i, point[0], point[1], point[2], 0.0);
}
selectionNode.SetSelectionList(frustumcorners);
frustum.Update();
vtkUnstructuredGrid grid = new vtkUnstructuredGrid();
grid = vtkUnstructuredGrid.SafeDownCast(frustum.GetOutput());
vtkIdTypeArray ids = null;
if (m_pickTarget == ePickTarget.Cell)
{
ids = vtkIdTypeArray.SafeDownCast(grid.GetCellData().GetArray("vtkOriginalCellIds"));
if (ids == null)
{
continue;
}
selectedActor.AddCellId(ids, m_pickBoolean);
OnObjectsSelected(sender, e);
}
else if (m_pickTarget == ePickTarget.Point)
{
ids = vtkIdTypeArray.SafeDownCast(grid.GetPointData().GetArray("vtkOriginalPointIds"));
if (ids == null)
{
continue;
}
selectedActor.AddPointId(ids, m_pickBoolean);
OnObjectsSelected(sender, e);
}
}
}
}
static public void Simulate(CMshTriangleMesherDouble aMesher, vtkUnstructuredGrid anUnstructuredGrid, vtkRenderer aRenderer)
{
var T = new CPdeFieldDouble();
T.setMesh(aMesher.mesh());
T.setDiscretMethod(EDiscretMethod.DM_FEM);
T.setDiscretOrder(EDiscretOrder.DO_LINEAR);
T.setDiscretLocation(EDiscretLocation.DL_NODE);
T.setSimulationType(ESimulationType.ST_THERMAL);
T.setNbDofs(1);
// Set BC and initial conditions
T.setSize(T.mesh().nbNodes());
for (int i = 0; i < T.mesh().nbNodes(); ++i)
{
int aNodeId = T.mesh().nodeId(i);
var aNode = T.mesh().node(aNodeId);
if (Math.Abs(aNode.x() + 0.0) < 1E-6)
{
T.setFixedValue(i, 0.0);
}
if (Math.Abs(aNode.x() - 1.0) < 1E-6)
{
T.setFixedValue(i, 0.0);
}
if (Math.Abs(aNode.y() + 0.0) < 1E-6)
{
T.setFixedValue(i, 0.0);
}
if (Math.Abs(aNode.y() - 1.0) < 1E-6)
{
T.setFixedValue(i, 1.0);
}
T.setValue(i, 0.0);
}
// Steady conduction in a rectangle
var aPdeEquation = new CPdeEquationDouble(ENigMA.laplacian(T));
aPdeEquation.setElimination(T);
aPdeEquation.solve(T);
var sScalars = vtkFloatArray.New();
sScalars.SetNumberOfTuples(T.mesh().nbNodes());
for (int i = 0; i < T.mesh().nbNodes(); ++i)
{
int aNodeId = T.mesh().nodeId(i);
sScalars.SetTuple1(aNodeId, T.value(i));
}
anUnstructuredGrid.GetPointData().SetScalars(sScalars);
var aGeometryFilter = vtkGeometryFilter.New();
aGeometryFilter.SetInput(anUnstructuredGrid);
aGeometryFilter.Update();
var aBoundedFilter = vtkBandedPolyDataContourFilter.New();
aBoundedFilter.SetInput(aGeometryFilter.GetOutput());
aBoundedFilter.GenerateValues(24, sScalars.GetRange()[0], sScalars.GetRange()[1]);
var aLookupTable = vtkLookupTable.New();
aLookupTable.SetNumberOfColors(256);
aLookupTable.SetHueRange(0.667, 0.0);
aLookupTable.Build();
var aBandedMapper = vtkPolyDataMapper.New();
aBandedMapper.SetInputConnection(aBoundedFilter.GetOutputPort());
aBandedMapper.SetScalarModeToUsePointData();
aBandedMapper.SetScalarRange(sScalars.GetRange()[0], sScalars.GetRange()[1]);
aBandedMapper.SetLookupTable(aLookupTable);
// add actor to the renderer
var aMeshActor = vtkActor.New();
aRenderer.AddActor(aMeshActor);
aMeshActor.SetMapper(aBandedMapper);
var aScalarBarActor = vtkScalarBarActor.New();
aRenderer.AddActor(aScalarBarActor);
aScalarBarActor.SetLookupTable(aLookupTable);
aScalarBarActor.SetTitle("Temperature");
aScalarBarActor.SetNumberOfLabels(6);
aScalarBarActor.Modified();
}
/// <summary>
/// The main entry method called by the CSharp driver
/// </summary>
/// <param name="argv"></param>
public static void AVTestCellDerivs(String [] argv)
{
//Prefix Content is: ""
// Demonstrates vtkCellDerivatives for all cell types[]
//[]
// get the interactor ui[]
ren1 = vtkRenderer.New();
renWin = vtkRenderWindow.New();
renWin.AddRenderer((vtkRenderer)ren1);
iren = new vtkRenderWindowInteractor();
iren.SetRenderWindow((vtkRenderWindow)renWin);
// create a scene with one of each cell type[]
// Voxel[]
voxelPoints = new vtkPoints();
voxelPoints.SetNumberOfPoints((int)8);
voxelPoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
voxelPoints.InsertPoint((int)1,(double)1,(double)0,(double)0);
voxelPoints.InsertPoint((int)2,(double)0,(double)1,(double)0);
voxelPoints.InsertPoint((int)3,(double)1,(double)1,(double)0);
voxelPoints.InsertPoint((int)4,(double)0,(double)0,(double)1);
voxelPoints.InsertPoint((int)5,(double)1,(double)0,(double)1);
voxelPoints.InsertPoint((int)6,(double)0,(double)1,(double)1);
voxelPoints.InsertPoint((int)7,(double)1,(double)1,(double)1);
aVoxel = new vtkVoxel();
aVoxel.GetPointIds().SetId((int)0,(int)0);
aVoxel.GetPointIds().SetId((int)1,(int)1);
aVoxel.GetPointIds().SetId((int)2,(int)2);
aVoxel.GetPointIds().SetId((int)3,(int)3);
aVoxel.GetPointIds().SetId((int)4,(int)4);
aVoxel.GetPointIds().SetId((int)5,(int)5);
aVoxel.GetPointIds().SetId((int)6,(int)6);
aVoxel.GetPointIds().SetId((int)7,(int)7);
aVoxelGrid = new vtkUnstructuredGrid();
aVoxelGrid.Allocate((int)1,(int)1);
aVoxelGrid.InsertNextCell((int)aVoxel.GetCellType(),(vtkIdList)aVoxel.GetPointIds());
aVoxelGrid.SetPoints((vtkPoints)voxelPoints);
aVoxelMapper = new vtkDataSetMapper();
aVoxelMapper.SetInput((vtkDataSet)aVoxelGrid);
aVoxelActor = new vtkActor();
aVoxelActor.SetMapper((vtkMapper)aVoxelMapper);
aVoxelActor.GetProperty().BackfaceCullingOn();
// Hexahedron[]
hexahedronPoints = new vtkPoints();
hexahedronPoints.SetNumberOfPoints((int)8);
hexahedronPoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
hexahedronPoints.InsertPoint((int)1,(double)1,(double)0,(double)0);
hexahedronPoints.InsertPoint((int)2,(double)1,(double)1,(double)0);
hexahedronPoints.InsertPoint((int)3,(double)0,(double)1,(double)0);
hexahedronPoints.InsertPoint((int)4,(double)0,(double)0,(double)1);
hexahedronPoints.InsertPoint((int)5,(double)1,(double)0,(double)1);
hexahedronPoints.InsertPoint((int)6,(double)1,(double)1,(double)1);
hexahedronPoints.InsertPoint((int)7,(double)0,(double)1,(double)1);
aHexahedron = new vtkHexahedron();
aHexahedron.GetPointIds().SetId((int)0,(int)0);
aHexahedron.GetPointIds().SetId((int)1,(int)1);
aHexahedron.GetPointIds().SetId((int)2,(int)2);
aHexahedron.GetPointIds().SetId((int)3,(int)3);
aHexahedron.GetPointIds().SetId((int)4,(int)4);
aHexahedron.GetPointIds().SetId((int)5,(int)5);
aHexahedron.GetPointIds().SetId((int)6,(int)6);
aHexahedron.GetPointIds().SetId((int)7,(int)7);
aHexahedronGrid = new vtkUnstructuredGrid();
aHexahedronGrid.Allocate((int)1,(int)1);
aHexahedronGrid.InsertNextCell((int)aHexahedron.GetCellType(),(vtkIdList)aHexahedron.GetPointIds());
aHexahedronGrid.SetPoints((vtkPoints)hexahedronPoints);
aHexahedronMapper = new vtkDataSetMapper();
aHexahedronMapper.SetInput((vtkDataSet)aHexahedronGrid);
aHexahedronActor = new vtkActor();
aHexahedronActor.SetMapper((vtkMapper)aHexahedronMapper);
aHexahedronActor.AddPosition((double)2,(double)0,(double)0);
aHexahedronActor.GetProperty().BackfaceCullingOn();
// Tetra[]
tetraPoints = new vtkPoints();
tetraPoints.SetNumberOfPoints((int)4);
tetraPoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
tetraPoints.InsertPoint((int)1,(double)1,(double)0,(double)0);
tetraPoints.InsertPoint((int)2,(double)0,(double)1,(double)0);
tetraPoints.InsertPoint((int)3,(double)1,(double)1,(double)1);
aTetra = new vtkTetra();
aTetra.GetPointIds().SetId((int)0,(int)0);
aTetra.GetPointIds().SetId((int)1,(int)1);
aTetra.GetPointIds().SetId((int)2,(int)2);
aTetra.GetPointIds().SetId((int)3,(int)3);
aTetraGrid = new vtkUnstructuredGrid();
aTetraGrid.Allocate((int)1,(int)1);
aTetraGrid.InsertNextCell((int)aTetra.GetCellType(),(vtkIdList)aTetra.GetPointIds());
aTetraGrid.SetPoints((vtkPoints)tetraPoints);
aTetraMapper = new vtkDataSetMapper();
aTetraMapper.SetInput((vtkDataSet)aTetraGrid);
aTetraActor = new vtkActor();
aTetraActor.SetMapper((vtkMapper)aTetraMapper);
aTetraActor.AddPosition((double)4,(double)0,(double)0);
aTetraActor.GetProperty().BackfaceCullingOn();
// Wedge[]
wedgePoints = new vtkPoints();
wedgePoints.SetNumberOfPoints((int)6);
wedgePoints.InsertPoint((int)0,(double)0,(double)1,(double)0);
wedgePoints.InsertPoint((int)1,(double)0,(double)0,(double)0);
wedgePoints.InsertPoint((int)2,(double)0,(double).5,(double).5);
wedgePoints.InsertPoint((int)3,(double)1,(double)1,(double)0);
wedgePoints.InsertPoint((int)4,(double)1,(double)0,(double)0);
wedgePoints.InsertPoint((int)5,(double)1,(double).5,(double).5);
aWedge = new vtkWedge();
aWedge.GetPointIds().SetId((int)0,(int)0);
aWedge.GetPointIds().SetId((int)1,(int)1);
aWedge.GetPointIds().SetId((int)2,(int)2);
aWedge.GetPointIds().SetId((int)3,(int)3);
aWedge.GetPointIds().SetId((int)4,(int)4);
aWedge.GetPointIds().SetId((int)5,(int)5);
aWedgeGrid = new vtkUnstructuredGrid();
aWedgeGrid.Allocate((int)1,(int)1);
aWedgeGrid.InsertNextCell((int)aWedge.GetCellType(),(vtkIdList)aWedge.GetPointIds());
aWedgeGrid.SetPoints((vtkPoints)wedgePoints);
aWedgeMapper = new vtkDataSetMapper();
aWedgeMapper.SetInput((vtkDataSet)aWedgeGrid);
aWedgeActor = new vtkActor();
aWedgeActor.SetMapper((vtkMapper)aWedgeMapper);
aWedgeActor.AddPosition((double)6,(double)0,(double)0);
aWedgeActor.GetProperty().BackfaceCullingOn();
// Pyramid[]
pyramidPoints = new vtkPoints();
pyramidPoints.SetNumberOfPoints((int)5);
pyramidPoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
pyramidPoints.InsertPoint((int)1,(double)1,(double)0,(double)0);
pyramidPoints.InsertPoint((int)2,(double)1,(double)1,(double)0);
pyramidPoints.InsertPoint((int)3,(double)0,(double)1,(double)0);
pyramidPoints.InsertPoint((int)4,(double).5,(double).5,(double)1);
aPyramid = new vtkPyramid();
aPyramid.GetPointIds().SetId((int)0,(int)0);
aPyramid.GetPointIds().SetId((int)1,(int)1);
aPyramid.GetPointIds().SetId((int)2,(int)2);
aPyramid.GetPointIds().SetId((int)3,(int)3);
aPyramid.GetPointIds().SetId((int)4,(int)4);
aPyramidGrid = new vtkUnstructuredGrid();
aPyramidGrid.Allocate((int)1,(int)1);
aPyramidGrid.InsertNextCell((int)aPyramid.GetCellType(),(vtkIdList)aPyramid.GetPointIds());
aPyramidGrid.SetPoints((vtkPoints)pyramidPoints);
aPyramidMapper = new vtkDataSetMapper();
aPyramidMapper.SetInput((vtkDataSet)aPyramidGrid);
aPyramidActor = new vtkActor();
aPyramidActor.SetMapper((vtkMapper)aPyramidMapper);
aPyramidActor.AddPosition((double)8,(double)0,(double)0);
aPyramidActor.GetProperty().BackfaceCullingOn();
// Pixel[]
pixelPoints = new vtkPoints();
pixelPoints.SetNumberOfPoints((int)4);
pixelPoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
pixelPoints.InsertPoint((int)1,(double)1,(double)0,(double)0);
pixelPoints.InsertPoint((int)2,(double)0,(double)1,(double)0);
pixelPoints.InsertPoint((int)3,(double)1,(double)1,(double)0);
aPixel = new vtkPixel();
aPixel.GetPointIds().SetId((int)0,(int)0);
aPixel.GetPointIds().SetId((int)1,(int)1);
aPixel.GetPointIds().SetId((int)2,(int)2);
aPixel.GetPointIds().SetId((int)3,(int)3);
aPixelGrid = new vtkUnstructuredGrid();
aPixelGrid.Allocate((int)1,(int)1);
aPixelGrid.InsertNextCell((int)aPixel.GetCellType(),(vtkIdList)aPixel.GetPointIds());
aPixelGrid.SetPoints((vtkPoints)pixelPoints);
aPixelMapper = new vtkDataSetMapper();
aPixelMapper.SetInput((vtkDataSet)aPixelGrid);
aPixelActor = new vtkActor();
aPixelActor.SetMapper((vtkMapper)aPixelMapper);
aPixelActor.AddPosition((double)0,(double)0,(double)2);
aPixelActor.GetProperty().BackfaceCullingOn();
// Quad[]
quadPoints = new vtkPoints();
quadPoints.SetNumberOfPoints((int)4);
quadPoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
quadPoints.InsertPoint((int)1,(double)1,(double)0,(double)0);
quadPoints.InsertPoint((int)2,(double)1,(double)1,(double)0);
quadPoints.InsertPoint((int)3,(double)0,(double)1,(double)0);
aQuad = new vtkQuad();
aQuad.GetPointIds().SetId((int)0,(int)0);
aQuad.GetPointIds().SetId((int)1,(int)1);
aQuad.GetPointIds().SetId((int)2,(int)2);
aQuad.GetPointIds().SetId((int)3,(int)3);
aQuadGrid = new vtkUnstructuredGrid();
aQuadGrid.Allocate((int)1,(int)1);
aQuadGrid.InsertNextCell((int)aQuad.GetCellType(),(vtkIdList)aQuad.GetPointIds());
aQuadGrid.SetPoints((vtkPoints)quadPoints);
aQuadMapper = new vtkDataSetMapper();
aQuadMapper.SetInput((vtkDataSet)aQuadGrid);
aQuadActor = new vtkActor();
aQuadActor.SetMapper((vtkMapper)aQuadMapper);
aQuadActor.AddPosition((double)2,(double)0,(double)2);
aQuadActor.GetProperty().BackfaceCullingOn();
// Triangle[]
trianglePoints = new vtkPoints();
trianglePoints.SetNumberOfPoints((int)3);
trianglePoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
trianglePoints.InsertPoint((int)1,(double)1,(double)0,(double)0);
trianglePoints.InsertPoint((int)2,(double).5,(double).5,(double)0);
triangleTCoords = new vtkFloatArray();
triangleTCoords.SetNumberOfComponents((int)2);
triangleTCoords.SetNumberOfTuples((int)3);
triangleTCoords.InsertTuple2((int)0,(double)1,(double)1);
triangleTCoords.InsertTuple2((int)1,(double)2,(double)2);
triangleTCoords.InsertTuple2((int)2,(double)3,(double)3);
aTriangle = new vtkTriangle();
aTriangle.GetPointIds().SetId((int)0,(int)0);
aTriangle.GetPointIds().SetId((int)1,(int)1);
aTriangle.GetPointIds().SetId((int)2,(int)2);
aTriangleGrid = new vtkUnstructuredGrid();
aTriangleGrid.Allocate((int)1,(int)1);
aTriangleGrid.InsertNextCell((int)aTriangle.GetCellType(),(vtkIdList)aTriangle.GetPointIds());
aTriangleGrid.SetPoints((vtkPoints)trianglePoints);
aTriangleGrid.GetPointData().SetTCoords((vtkDataArray)triangleTCoords);
aTriangleMapper = new vtkDataSetMapper();
aTriangleMapper.SetInput((vtkDataSet)aTriangleGrid);
aTriangleActor = new vtkActor();
aTriangleActor.SetMapper((vtkMapper)aTriangleMapper);
aTriangleActor.AddPosition((double)4,(double)0,(double)2);
aTriangleActor.GetProperty().BackfaceCullingOn();
// Polygon[]
polygonPoints = new vtkPoints();
polygonPoints.SetNumberOfPoints((int)4);
polygonPoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
polygonPoints.InsertPoint((int)1,(double)1,(double)0,(double)0);
polygonPoints.InsertPoint((int)2,(double)1,(double)1,(double)0);
polygonPoints.InsertPoint((int)3,(double)0,(double)1,(double)0);
aPolygon = new vtkPolygon();
aPolygon.GetPointIds().SetNumberOfIds((int)4);
aPolygon.GetPointIds().SetId((int)0,(int)0);
aPolygon.GetPointIds().SetId((int)1,(int)1);
aPolygon.GetPointIds().SetId((int)2,(int)2);
aPolygon.GetPointIds().SetId((int)3,(int)3);
aPolygonGrid = new vtkUnstructuredGrid();
aPolygonGrid.Allocate((int)1,(int)1);
aPolygonGrid.InsertNextCell((int)aPolygon.GetCellType(),(vtkIdList)aPolygon.GetPointIds());
aPolygonGrid.SetPoints((vtkPoints)polygonPoints);
aPolygonMapper = new vtkDataSetMapper();
aPolygonMapper.SetInput((vtkDataSet)aPolygonGrid);
aPolygonActor = new vtkActor();
aPolygonActor.SetMapper((vtkMapper)aPolygonMapper);
aPolygonActor.AddPosition((double)6,(double)0,(double)2);
aPolygonActor.GetProperty().BackfaceCullingOn();
// Triangle strip[]
triangleStripPoints = new vtkPoints();
triangleStripPoints.SetNumberOfPoints((int)5);
triangleStripPoints.InsertPoint((int)0,(double)0,(double)1,(double)0);
triangleStripPoints.InsertPoint((int)1,(double)0,(double)0,(double)0);
triangleStripPoints.InsertPoint((int)2,(double)1,(double)1,(double)0);
triangleStripPoints.InsertPoint((int)3,(double)1,(double)0,(double)0);
triangleStripPoints.InsertPoint((int)4,(double)2,(double)1,(double)0);
triangleStripTCoords = new vtkFloatArray();
triangleStripTCoords.SetNumberOfComponents((int)2);
triangleStripTCoords.SetNumberOfTuples((int)3);
triangleStripTCoords.InsertTuple2((int)0,(double)1,(double)1);
triangleStripTCoords.InsertTuple2((int)1,(double)2,(double)2);
triangleStripTCoords.InsertTuple2((int)2,(double)3,(double)3);
triangleStripTCoords.InsertTuple2((int)3,(double)4,(double)4);
triangleStripTCoords.InsertTuple2((int)4,(double)5,(double)5);
aTriangleStrip = new vtkTriangleStrip();
aTriangleStrip.GetPointIds().SetNumberOfIds((int)5);
aTriangleStrip.GetPointIds().SetId((int)0,(int)0);
aTriangleStrip.GetPointIds().SetId((int)1,(int)1);
aTriangleStrip.GetPointIds().SetId((int)2,(int)2);
aTriangleStrip.GetPointIds().SetId((int)3,(int)3);
aTriangleStrip.GetPointIds().SetId((int)4,(int)4);
aTriangleStripGrid = new vtkUnstructuredGrid();
aTriangleStripGrid.Allocate((int)1,(int)1);
aTriangleStripGrid.InsertNextCell((int)aTriangleStrip.GetCellType(),(vtkIdList)aTriangleStrip.GetPointIds());
aTriangleStripGrid.SetPoints((vtkPoints)triangleStripPoints);
aTriangleStripGrid.GetPointData().SetTCoords((vtkDataArray)triangleStripTCoords);
aTriangleStripMapper = new vtkDataSetMapper();
aTriangleStripMapper.SetInput((vtkDataSet)aTriangleStripGrid);
aTriangleStripActor = new vtkActor();
aTriangleStripActor.SetMapper((vtkMapper)aTriangleStripMapper);
aTriangleStripActor.AddPosition((double)8,(double)0,(double)2);
aTriangleStripActor.GetProperty().BackfaceCullingOn();
// Line[]
linePoints = new vtkPoints();
linePoints.SetNumberOfPoints((int)2);
linePoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
linePoints.InsertPoint((int)1,(double)1,(double)1,(double)0);
aLine = new vtkLine();
aLine.GetPointIds().SetId((int)0,(int)0);
aLine.GetPointIds().SetId((int)1,(int)1);
aLineGrid = new vtkUnstructuredGrid();
aLineGrid.Allocate((int)1,(int)1);
aLineGrid.InsertNextCell((int)aLine.GetCellType(),(vtkIdList)aLine.GetPointIds());
aLineGrid.SetPoints((vtkPoints)linePoints);
aLineMapper = new vtkDataSetMapper();
aLineMapper.SetInput((vtkDataSet)aLineGrid);
aLineActor = new vtkActor();
aLineActor.SetMapper((vtkMapper)aLineMapper);
aLineActor.AddPosition((double)0,(double)0,(double)4);
aLineActor.GetProperty().BackfaceCullingOn();
// Polyline[]
polyLinePoints = new vtkPoints();
polyLinePoints.SetNumberOfPoints((int)3);
polyLinePoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
polyLinePoints.InsertPoint((int)1,(double)1,(double)1,(double)0);
polyLinePoints.InsertPoint((int)2,(double)1,(double)0,(double)0);
aPolyLine = new vtkPolyLine();
aPolyLine.GetPointIds().SetNumberOfIds((int)3);
aPolyLine.GetPointIds().SetId((int)0,(int)0);
aPolyLine.GetPointIds().SetId((int)1,(int)1);
aPolyLine.GetPointIds().SetId((int)2,(int)2);
aPolyLineGrid = new vtkUnstructuredGrid();
aPolyLineGrid.Allocate((int)1,(int)1);
aPolyLineGrid.InsertNextCell((int)aPolyLine.GetCellType(),(vtkIdList)aPolyLine.GetPointIds());
aPolyLineGrid.SetPoints((vtkPoints)polyLinePoints);
aPolyLineMapper = new vtkDataSetMapper();
aPolyLineMapper.SetInput((vtkDataSet)aPolyLineGrid);
aPolyLineActor = new vtkActor();
aPolyLineActor.SetMapper((vtkMapper)aPolyLineMapper);
aPolyLineActor.AddPosition((double)2,(double)0,(double)4);
aPolyLineActor.GetProperty().BackfaceCullingOn();
// Vertex[]
vertexPoints = new vtkPoints();
vertexPoints.SetNumberOfPoints((int)1);
vertexPoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
aVertex = new vtkVertex();
aVertex.GetPointIds().SetId((int)0,(int)0);
aVertexGrid = new vtkUnstructuredGrid();
aVertexGrid.Allocate((int)1,(int)1);
aVertexGrid.InsertNextCell((int)aVertex.GetCellType(),(vtkIdList)aVertex.GetPointIds());
aVertexGrid.SetPoints((vtkPoints)vertexPoints);
aVertexMapper = new vtkDataSetMapper();
aVertexMapper.SetInput((vtkDataSet)aVertexGrid);
aVertexActor = new vtkActor();
aVertexActor.SetMapper((vtkMapper)aVertexMapper);
aVertexActor.AddPosition((double)0,(double)0,(double)6);
aVertexActor.GetProperty().BackfaceCullingOn();
// Polyvertex[]
polyVertexPoints = new vtkPoints();
polyVertexPoints.SetNumberOfPoints((int)3);
polyVertexPoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
polyVertexPoints.InsertPoint((int)1,(double)1,(double)0,(double)0);
polyVertexPoints.InsertPoint((int)2,(double)1,(double)1,(double)0);
aPolyVertex = new vtkPolyVertex();
aPolyVertex.GetPointIds().SetNumberOfIds((int)3);
aPolyVertex.GetPointIds().SetId((int)0,(int)0);
aPolyVertex.GetPointIds().SetId((int)1,(int)1);
aPolyVertex.GetPointIds().SetId((int)2,(int)2);
aPolyVertexGrid = new vtkUnstructuredGrid();
aPolyVertexGrid.Allocate((int)1,(int)1);
aPolyVertexGrid.InsertNextCell((int)aPolyVertex.GetCellType(),(vtkIdList)aPolyVertex.GetPointIds());
aPolyVertexGrid.SetPoints((vtkPoints)polyVertexPoints);
aPolyVertexMapper = new vtkDataSetMapper();
aPolyVertexMapper.SetInput((vtkDataSet)aPolyVertexGrid);
aPolyVertexActor = new vtkActor();
aPolyVertexActor.SetMapper((vtkMapper)aPolyVertexMapper);
aPolyVertexActor.AddPosition((double)2,(double)0,(double)6);
aPolyVertexActor.GetProperty().BackfaceCullingOn();
// Pentagonal prism[]
pentaPoints = new vtkPoints();
pentaPoints.SetNumberOfPoints((int)10);
pentaPoints.InsertPoint((int)0,(double)0.25,(double)0.0,(double)0.0);
pentaPoints.InsertPoint((int)1,(double)0.75,(double)0.0,(double)0.0);
pentaPoints.InsertPoint((int)2,(double)1.0,(double)0.5,(double)0.0);
pentaPoints.InsertPoint((int)3,(double)0.5,(double)1.0,(double)0.0);
pentaPoints.InsertPoint((int)4,(double)0.0,(double)0.5,(double)0.0);
pentaPoints.InsertPoint((int)5,(double)0.25,(double)0.0,(double)1.0);
pentaPoints.InsertPoint((int)6,(double)0.75,(double)0.0,(double)1.0);
pentaPoints.InsertPoint((int)7,(double)1.0,(double)0.5,(double)1.0);
pentaPoints.InsertPoint((int)8,(double)0.5,(double)1.0,(double)1.0);
pentaPoints.InsertPoint((int)9,(double)0.0,(double)0.5,(double)1.0);
aPenta = new vtkPentagonalPrism();
aPenta.GetPointIds().SetId((int)0,(int)0);
aPenta.GetPointIds().SetId((int)1,(int)1);
aPenta.GetPointIds().SetId((int)2,(int)2);
aPenta.GetPointIds().SetId((int)3,(int)3);
aPenta.GetPointIds().SetId((int)4,(int)4);
aPenta.GetPointIds().SetId((int)5,(int)5);
aPenta.GetPointIds().SetId((int)6,(int)6);
aPenta.GetPointIds().SetId((int)7,(int)7);
aPenta.GetPointIds().SetId((int)8,(int)8);
aPenta.GetPointIds().SetId((int)9,(int)9);
aPentaGrid = new vtkUnstructuredGrid();
aPentaGrid.Allocate((int)1,(int)1);
aPentaGrid.InsertNextCell((int)aPenta.GetCellType(),(vtkIdList)aPenta.GetPointIds());
aPentaGrid.SetPoints((vtkPoints)pentaPoints);
aPentaMapper = new vtkDataSetMapper();
aPentaMapper.SetInput((vtkDataSet)aPentaGrid);
aPentaActor = new vtkActor();
aPentaActor.SetMapper((vtkMapper)aPentaMapper);
aPentaActor.AddPosition((double)10,(double)0,(double)0);
aPentaActor.GetProperty().BackfaceCullingOn();
// Hexagonal prism[]
hexaPoints = new vtkPoints();
hexaPoints.SetNumberOfPoints((int)12);
hexaPoints.InsertPoint((int)0,(double)0.0,(double)0.0,(double)0.0);
hexaPoints.InsertPoint((int)1,(double)0.5,(double)0.0,(double)0.0);
hexaPoints.InsertPoint((int)2,(double)1.0,(double)0.5,(double)0.0);
hexaPoints.InsertPoint((int)3,(double)1.0,(double)1.0,(double)0.0);
hexaPoints.InsertPoint((int)4,(double)0.5,(double)1.0,(double)0.0);
hexaPoints.InsertPoint((int)5,(double)0.0,(double)0.5,(double)0.0);
hexaPoints.InsertPoint((int)6,(double)0.0,(double)0.0,(double)1.0);
hexaPoints.InsertPoint((int)7,(double)0.5,(double)0.0,(double)1.0);
hexaPoints.InsertPoint((int)8,(double)1.0,(double)0.5,(double)1.0);
hexaPoints.InsertPoint((int)9,(double)1.0,(double)1.0,(double)1.0);
hexaPoints.InsertPoint((int)10,(double)0.5,(double)1.0,(double)1.0);
hexaPoints.InsertPoint((int)11,(double)0.0,(double)0.5,(double)1.0);
aHexa = new vtkHexagonalPrism();
aHexa.GetPointIds().SetId((int)0,(int)0);
aHexa.GetPointIds().SetId((int)1,(int)1);
aHexa.GetPointIds().SetId((int)2,(int)2);
aHexa.GetPointIds().SetId((int)3,(int)3);
aHexa.GetPointIds().SetId((int)4,(int)4);
aHexa.GetPointIds().SetId((int)5,(int)5);
aHexa.GetPointIds().SetId((int)6,(int)6);
aHexa.GetPointIds().SetId((int)7,(int)7);
aHexa.GetPointIds().SetId((int)8,(int)8);
aHexa.GetPointIds().SetId((int)9,(int)9);
aHexa.GetPointIds().SetId((int)10,(int)10);
aHexa.GetPointIds().SetId((int)11,(int)11);
aHexaGrid = new vtkUnstructuredGrid();
aHexaGrid.Allocate((int)1,(int)1);
aHexaGrid.InsertNextCell((int)aHexa.GetCellType(),(vtkIdList)aHexa.GetPointIds());
aHexaGrid.SetPoints((vtkPoints)hexaPoints);
aHexaMapper = new vtkDataSetMapper();
aHexaMapper.SetInput((vtkDataSet)aHexaGrid);
aHexaActor = new vtkActor();
aHexaActor.SetMapper((vtkMapper)aHexaMapper);
aHexaActor.AddPosition((double)12,(double)0,(double)0);
aHexaActor.GetProperty().BackfaceCullingOn();
ren1.SetBackground((double)1,(double)1,(double)1);
ren1.AddActor((vtkProp)aVoxelActor);
aVoxelActor.GetProperty().SetDiffuseColor((double)1,(double)0,(double)0);
ren1.AddActor((vtkProp)aHexahedronActor);
aHexahedronActor.GetProperty().SetDiffuseColor((double)1,(double)1,(double)0);
ren1.AddActor((vtkProp)aTetraActor);
aTetraActor.GetProperty().SetDiffuseColor((double)0,(double)1,(double)0);
ren1.AddActor((vtkProp)aWedgeActor);
aWedgeActor.GetProperty().SetDiffuseColor((double)0,(double)1,(double)1);
ren1.AddActor((vtkProp)aPyramidActor);
aPyramidActor.GetProperty().SetDiffuseColor((double)1,(double)0,(double)1);
ren1.AddActor((vtkProp)aPixelActor);
aPixelActor.GetProperty().SetDiffuseColor((double)0,(double)1,(double)1);
ren1.AddActor((vtkProp)aQuadActor);
aQuadActor.GetProperty().SetDiffuseColor((double)1,(double)0,(double)1);
ren1.AddActor((vtkProp)aTriangleActor);
aTriangleActor.GetProperty().SetDiffuseColor((double).3,(double)1,(double).5);
ren1.AddActor((vtkProp)aPolygonActor);
aPolygonActor.GetProperty().SetDiffuseColor((double)1,(double).4,(double).5);
ren1.AddActor((vtkProp)aTriangleStripActor);
aTriangleStripActor.GetProperty().SetDiffuseColor((double).3,(double).7,(double)1);
ren1.AddActor((vtkProp)aLineActor);
aLineActor.GetProperty().SetDiffuseColor((double).2,(double)1,(double)1);
ren1.AddActor((vtkProp)aPolyLineActor);
aPolyLineActor.GetProperty().SetDiffuseColor((double)1,(double)1,(double)1);
ren1.AddActor((vtkProp)aVertexActor);
aVertexActor.GetProperty().SetDiffuseColor((double)1,(double)1,(double)1);
ren1.AddActor((vtkProp)aPolyVertexActor);
aPolyVertexActor.GetProperty().SetDiffuseColor((double)1,(double)1,(double)1);
ren1.AddActor((vtkProp)aPentaActor);
aPentaActor.GetProperty().SetDiffuseColor((double)1,(double)1,(double)0);
ren1.AddActor((vtkProp)aHexaActor);
aHexaActor.GetProperty().SetDiffuseColor((double)1,(double)1,(double)0);
//[]
// get the cell center of each type and put a glyph there[]
//[]
ball = new vtkSphereSource();
ball.SetRadius((double).2);
bool tryWorked = false;
aVoxelScalars = new vtkFloatArray();
N = aVoxelGrid.GetNumberOfPoints();
aVoxelScalar = new vtkFloatArray();
aVoxelScalar.SetNumberOfTuples((int)N);
aVoxelScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aVoxelScalar.SetValue(i,0);
i = i + 1;
}
aVoxelScalar.SetValue(0,4);
aVoxelGrid.GetPointData().SetScalars(aVoxelScalar);
aHexahedronScalars = new vtkFloatArray();
N = aHexahedronGrid.GetNumberOfPoints();
aHexahedronScalar = new vtkFloatArray();
aHexahedronScalar.SetNumberOfTuples((int)N);
aHexahedronScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aHexahedronScalar.SetValue(i,0);
i = i + 1;
}
aHexahedronScalar.SetValue(0,4);
aHexahedronGrid.GetPointData().SetScalars(aHexahedronScalar);
aWedgeScalars = new vtkFloatArray();
N = aWedgeGrid.GetNumberOfPoints();
aWedgeScalar = new vtkFloatArray();
aWedgeScalar.SetNumberOfTuples((int)N);
aWedgeScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aWedgeScalar.SetValue(i,0);
i = i + 1;
}
aWedgeScalar.SetValue(0,4);
aWedgeGrid.GetPointData().SetScalars(aWedgeScalar);
aPyramidScalars = new vtkFloatArray();
N = aPyramidGrid.GetNumberOfPoints();
aPyramidScalar = new vtkFloatArray();
aPyramidScalar.SetNumberOfTuples((int)N);
aPyramidScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aPyramidScalar.SetValue(i,0);
i = i + 1;
}
aPyramidScalar.SetValue(0,4);
aPyramidGrid.GetPointData().SetScalars(aPyramidScalar);
aTetraScalars = new vtkFloatArray();
N = aTetraGrid.GetNumberOfPoints();
aTetraScalar = new vtkFloatArray();
aTetraScalar.SetNumberOfTuples((int)N);
aTetraScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aTetraScalar.SetValue(i,0);
i = i + 1;
}
aTetraScalar.SetValue(0,4);
aTetraGrid.GetPointData().SetScalars(aTetraScalar);
aQuadScalars = new vtkFloatArray();
N = aQuadGrid.GetNumberOfPoints();
aQuadScalar = new vtkFloatArray();
aQuadScalar.SetNumberOfTuples((int)N);
aQuadScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aQuadScalar.SetValue(i,0);
i = i + 1;
}
aQuadScalar.SetValue(0,4);
aQuadGrid.GetPointData().SetScalars(aQuadScalar);
aTriangleScalars = new vtkFloatArray();
N = aTriangleGrid.GetNumberOfPoints();
aTriangleScalar = new vtkFloatArray();
aTriangleScalar.SetNumberOfTuples((int)N);
aTriangleScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aTriangleScalar.SetValue(i,0);
i = i + 1;
}
aTriangleScalar.SetValue(0,4);
aTriangleGrid.GetPointData().SetScalars(aTriangleScalar);
aTriangleStripScalars = new vtkFloatArray();
N = aTriangleStripGrid.GetNumberOfPoints();
aTriangleStripScalar = new vtkFloatArray();
aTriangleStripScalar.SetNumberOfTuples((int)N);
aTriangleStripScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aTriangleStripScalar.SetValue(i,0);
i = i + 1;
}
aTriangleStripScalar.SetValue(0,4);
aTriangleStripGrid.GetPointData().SetScalars(aTriangleStripScalar);
aLineScalars = new vtkFloatArray();
N = aLineGrid.GetNumberOfPoints();
aLineScalar = new vtkFloatArray();
aLineScalar.SetNumberOfTuples((int)N);
aLineScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aLineScalar.SetValue(i,0);
i = i + 1;
}
aLineScalar.SetValue(0,4);
aLineGrid.GetPointData().SetScalars(aLineScalar);
aPolyLineScalars = new vtkFloatArray();
N = aPolyLineGrid.GetNumberOfPoints();
aPolyLineScalar = new vtkFloatArray();
aPolyLineScalar.SetNumberOfTuples((int)N);
aPolyLineScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aPolyLineScalar.SetValue(i,0);
i = i + 1;
}
aPolyLineScalar.SetValue(0,4);
aPolyLineGrid.GetPointData().SetScalars(aPolyLineScalar);
aVertexScalars = new vtkFloatArray();
N = aVertexGrid.GetNumberOfPoints();
aVertexScalar = new vtkFloatArray();
aVertexScalar.SetNumberOfTuples((int)N);
aVertexScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aVertexScalar.SetValue(i,0);
i = i + 1;
}
aVertexScalar.SetValue(0,4);
aVertexGrid.GetPointData().SetScalars(aVertexScalar);
aPolyVertexScalars = new vtkFloatArray();
N = aPolyVertexGrid.GetNumberOfPoints();
aPolyVertexScalar = new vtkFloatArray();
aPolyVertexScalar.SetNumberOfTuples((int)N);
aPolyVertexScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aPolyVertexScalar.SetValue(i,0);
i = i + 1;
}
aPolyVertexScalar.SetValue(0,4);
aPolyVertexGrid.GetPointData().SetScalars(aPolyVertexScalar);
aPixelScalars = new vtkFloatArray();
N = aPixelGrid.GetNumberOfPoints();
aPixelScalar = new vtkFloatArray();
aPixelScalar.SetNumberOfTuples((int)N);
aPixelScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aPixelScalar.SetValue(i,0);
i = i + 1;
}
aPixelScalar.SetValue(0,4);
aPixelGrid.GetPointData().SetScalars(aPixelScalar);
aPolygonScalars = new vtkFloatArray();
N = aPolygonGrid.GetNumberOfPoints();
aPolygonScalar = new vtkFloatArray();
aPolygonScalar.SetNumberOfTuples((int)N);
aPolygonScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aPolygonScalar.SetValue(i,0);
i = i + 1;
}
aPolygonScalar.SetValue(0,4);
aPolygonGrid.GetPointData().SetScalars(aPolygonScalar);
aPentaScalars = new vtkFloatArray();
N = aPentaGrid.GetNumberOfPoints();
aPentaScalar = new vtkFloatArray();
aPentaScalar.SetNumberOfTuples((int)N);
aPentaScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aPentaScalar.SetValue(i,0);
i = i + 1;
}
aPentaScalar.SetValue(0,4);
aPentaGrid.GetPointData().SetScalars(aPentaScalar);
aHexaScalars = new vtkFloatArray();
N = aHexaGrid.GetNumberOfPoints();
aHexaScalar = new vtkFloatArray();
aHexaScalar.SetNumberOfTuples((int)N);
aHexaScalar.SetNumberOfComponents(1);
i = 0;
while((i) < N)
{
aHexaScalar.SetValue(i,0);
i = i + 1;
}
aHexaScalar.SetValue(0,4);
aHexaGrid.GetPointData().SetScalars(aHexaScalar);
// write to the temp directory if possible, otherwise use .[]
dir = ".";
dir = TclToCsScriptTestDriver.GetTempDirectory();
aVoxelderivs = new vtkCellDerivatives();
aVoxelderivs.SetInput(aVoxelGrid);
aVoxelderivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aVoxel";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aVoxelWriter = new vtkUnstructuredGridWriter();
aVoxelWriter.SetInputConnection(aVoxelderivs.GetOutputPort());
aVoxelWriter.SetFileName(FileName);
aVoxelWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aVoxelCenters = new vtkCellCenters();
aVoxelCenters.SetInputConnection(aVoxelderivs.GetOutputPort());
aVoxelCenters.VertexCellsOn();
aVoxelhog = new vtkHedgeHog();
aVoxelhog.SetInputConnection(aVoxelCenters.GetOutputPort());
aVoxelmapHog = vtkPolyDataMapper.New();
aVoxelmapHog.SetInputConnection(aVoxelhog.GetOutputPort());
aVoxelmapHog.SetScalarModeToUseCellData();
aVoxelmapHog.ScalarVisibilityOff();
aVoxelhogActor = new vtkActor();
aVoxelhogActor.SetMapper(aVoxelmapHog);
aVoxelhogActor.GetProperty().SetColor(0,1,0);
aVoxelGlyph3D = new vtkGlyph3D();
aVoxelGlyph3D.SetInputConnection(aVoxelCenters.GetOutputPort());
aVoxelGlyph3D.SetSource(ball.GetOutput());
aVoxelCentersMapper = vtkPolyDataMapper.New();
aVoxelCentersMapper.SetInputConnection(aVoxelGlyph3D.GetOutputPort());
aVoxelCentersActor = new vtkActor();
aVoxelCentersActor.SetMapper(aVoxelCentersMapper);
aVoxelhogActor.SetPosition(aVoxelActor.GetPosition()[0],aVoxelActor.GetPosition()[1],aVoxelActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aVoxelhogActor);
aVoxelhogActor.GetProperty().SetRepresentationToWireframe();
aHexahedronderivs = new vtkCellDerivatives();
aHexahedronderivs.SetInput(aHexahedronGrid);
aHexahedronderivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aHexahedron";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aHexahedronWriter = new vtkUnstructuredGridWriter();
aHexahedronWriter.SetInputConnection(aHexahedronderivs.GetOutputPort());
aHexahedronWriter.SetFileName(FileName);
aHexahedronWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aHexahedronCenters = new vtkCellCenters();
aHexahedronCenters.SetInputConnection(aHexahedronderivs.GetOutputPort());
aHexahedronCenters.VertexCellsOn();
aHexahedronhog = new vtkHedgeHog();
aHexahedronhog.SetInputConnection(aHexahedronCenters.GetOutputPort());
aHexahedronmapHog = vtkPolyDataMapper.New();
aHexahedronmapHog.SetInputConnection(aHexahedronhog.GetOutputPort());
aHexahedronmapHog.SetScalarModeToUseCellData();
aHexahedronmapHog.ScalarVisibilityOff();
aHexahedronhogActor = new vtkActor();
aHexahedronhogActor.SetMapper(aHexahedronmapHog);
aHexahedronhogActor.GetProperty().SetColor(0,1,0);
aHexahedronGlyph3D = new vtkGlyph3D();
aHexahedronGlyph3D.SetInputConnection(aHexahedronCenters.GetOutputPort());
aHexahedronGlyph3D.SetSource(ball.GetOutput());
aHexahedronCentersMapper = vtkPolyDataMapper.New();
aHexahedronCentersMapper.SetInputConnection(aHexahedronGlyph3D.GetOutputPort());
aHexahedronCentersActor = new vtkActor();
aHexahedronCentersActor.SetMapper(aHexahedronCentersMapper);
aHexahedronhogActor.SetPosition(aHexahedronActor.GetPosition()[0],aHexahedronActor.GetPosition()[1],aHexahedronActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aHexahedronhogActor);
aHexahedronhogActor.GetProperty().SetRepresentationToWireframe();
aWedgederivs = new vtkCellDerivatives();
aWedgederivs.SetInput(aWedgeGrid);
aWedgederivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aWedge";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aWedgeWriter = new vtkUnstructuredGridWriter();
aWedgeWriter.SetInputConnection(aWedgederivs.GetOutputPort());
aWedgeWriter.SetFileName(FileName);
aWedgeWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aWedgeCenters = new vtkCellCenters();
aWedgeCenters.SetInputConnection(aWedgederivs.GetOutputPort());
aWedgeCenters.VertexCellsOn();
aWedgehog = new vtkHedgeHog();
aWedgehog.SetInputConnection(aWedgeCenters.GetOutputPort());
aWedgemapHog = vtkPolyDataMapper.New();
aWedgemapHog.SetInputConnection(aWedgehog.GetOutputPort());
aWedgemapHog.SetScalarModeToUseCellData();
aWedgemapHog.ScalarVisibilityOff();
aWedgehogActor = new vtkActor();
aWedgehogActor.SetMapper(aWedgemapHog);
aWedgehogActor.GetProperty().SetColor(0,1,0);
aWedgeGlyph3D = new vtkGlyph3D();
aWedgeGlyph3D.SetInputConnection(aWedgeCenters.GetOutputPort());
aWedgeGlyph3D.SetSource(ball.GetOutput());
aWedgeCentersMapper = vtkPolyDataMapper.New();
aWedgeCentersMapper.SetInputConnection(aWedgeGlyph3D.GetOutputPort());
aWedgeCentersActor = new vtkActor();
aWedgeCentersActor.SetMapper(aWedgeCentersMapper);
aWedgehogActor.SetPosition(aWedgeActor.GetPosition()[0],aWedgeActor.GetPosition()[1],aWedgeActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aWedgehogActor);
aWedgehogActor.GetProperty().SetRepresentationToWireframe();
aPyramidderivs = new vtkCellDerivatives();
aPyramidderivs.SetInput(aPyramidGrid);
aPyramidderivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aPyramid";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aPyramidWriter = new vtkUnstructuredGridWriter();
aPyramidWriter.SetInputConnection(aPyramidderivs.GetOutputPort());
aPyramidWriter.SetFileName(FileName);
aPyramidWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aPyramidCenters = new vtkCellCenters();
aPyramidCenters.SetInputConnection(aPyramidderivs.GetOutputPort());
aPyramidCenters.VertexCellsOn();
aPyramidhog = new vtkHedgeHog();
aPyramidhog.SetInputConnection(aPyramidCenters.GetOutputPort());
aPyramidmapHog = vtkPolyDataMapper.New();
aPyramidmapHog.SetInputConnection(aPyramidhog.GetOutputPort());
aPyramidmapHog.SetScalarModeToUseCellData();
aPyramidmapHog.ScalarVisibilityOff();
aPyramidhogActor = new vtkActor();
aPyramidhogActor.SetMapper(aPyramidmapHog);
aPyramidhogActor.GetProperty().SetColor(0,1,0);
aPyramidGlyph3D = new vtkGlyph3D();
aPyramidGlyph3D.SetInputConnection(aPyramidCenters.GetOutputPort());
aPyramidGlyph3D.SetSource(ball.GetOutput());
aPyramidCentersMapper = vtkPolyDataMapper.New();
aPyramidCentersMapper.SetInputConnection(aPyramidGlyph3D.GetOutputPort());
aPyramidCentersActor = new vtkActor();
aPyramidCentersActor.SetMapper(aPyramidCentersMapper);
aPyramidhogActor.SetPosition(aPyramidActor.GetPosition()[0],aPyramidActor.GetPosition()[1],aPyramidActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aPyramidhogActor);
aPyramidhogActor.GetProperty().SetRepresentationToWireframe();
aTetraderivs = new vtkCellDerivatives();
aTetraderivs.SetInput(aTetraGrid);
aTetraderivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aTetra";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aTetraWriter = new vtkUnstructuredGridWriter();
aTetraWriter.SetInputConnection(aTetraderivs.GetOutputPort());
aTetraWriter.SetFileName(FileName);
aTetraWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aTetraCenters = new vtkCellCenters();
aTetraCenters.SetInputConnection(aTetraderivs.GetOutputPort());
aTetraCenters.VertexCellsOn();
aTetrahog = new vtkHedgeHog();
aTetrahog.SetInputConnection(aTetraCenters.GetOutputPort());
aTetramapHog = vtkPolyDataMapper.New();
aTetramapHog.SetInputConnection(aTetrahog.GetOutputPort());
aTetramapHog.SetScalarModeToUseCellData();
aTetramapHog.ScalarVisibilityOff();
aTetrahogActor = new vtkActor();
aTetrahogActor.SetMapper(aTetramapHog);
aTetrahogActor.GetProperty().SetColor(0,1,0);
aTetraGlyph3D = new vtkGlyph3D();
aTetraGlyph3D.SetInputConnection(aTetraCenters.GetOutputPort());
aTetraGlyph3D.SetSource(ball.GetOutput());
aTetraCentersMapper = vtkPolyDataMapper.New();
aTetraCentersMapper.SetInputConnection(aTetraGlyph3D.GetOutputPort());
aTetraCentersActor = new vtkActor();
aTetraCentersActor.SetMapper(aTetraCentersMapper);
aTetrahogActor.SetPosition(aTetraActor.GetPosition()[0],aTetraActor.GetPosition()[1],aTetraActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aTetrahogActor);
aTetrahogActor.GetProperty().SetRepresentationToWireframe();
aQuadderivs = new vtkCellDerivatives();
aQuadderivs.SetInput(aQuadGrid);
aQuadderivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aQuad";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aQuadWriter = new vtkUnstructuredGridWriter();
aQuadWriter.SetInputConnection(aQuadderivs.GetOutputPort());
aQuadWriter.SetFileName(FileName);
aQuadWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aQuadCenters = new vtkCellCenters();
aQuadCenters.SetInputConnection(aQuadderivs.GetOutputPort());
aQuadCenters.VertexCellsOn();
aQuadhog = new vtkHedgeHog();
aQuadhog.SetInputConnection(aQuadCenters.GetOutputPort());
aQuadmapHog = vtkPolyDataMapper.New();
aQuadmapHog.SetInputConnection(aQuadhog.GetOutputPort());
aQuadmapHog.SetScalarModeToUseCellData();
aQuadmapHog.ScalarVisibilityOff();
aQuadhogActor = new vtkActor();
aQuadhogActor.SetMapper(aQuadmapHog);
aQuadhogActor.GetProperty().SetColor(0,1,0);
aQuadGlyph3D = new vtkGlyph3D();
aQuadGlyph3D.SetInputConnection(aQuadCenters.GetOutputPort());
aQuadGlyph3D.SetSource(ball.GetOutput());
aQuadCentersMapper = vtkPolyDataMapper.New();
aQuadCentersMapper.SetInputConnection(aQuadGlyph3D.GetOutputPort());
aQuadCentersActor = new vtkActor();
aQuadCentersActor.SetMapper(aQuadCentersMapper);
aQuadhogActor.SetPosition(aQuadActor.GetPosition()[0],aQuadActor.GetPosition()[1],aQuadActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aQuadhogActor);
aQuadhogActor.GetProperty().SetRepresentationToWireframe();
aTrianglederivs = new vtkCellDerivatives();
aTrianglederivs.SetInput(aTriangleGrid);
aTrianglederivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aTriangle";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aTriangleWriter = new vtkUnstructuredGridWriter();
aTriangleWriter.SetInputConnection(aTrianglederivs.GetOutputPort());
aTriangleWriter.SetFileName(FileName);
aTriangleWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aTriangleCenters = new vtkCellCenters();
aTriangleCenters.SetInputConnection(aTrianglederivs.GetOutputPort());
aTriangleCenters.VertexCellsOn();
aTrianglehog = new vtkHedgeHog();
aTrianglehog.SetInputConnection(aTriangleCenters.GetOutputPort());
aTrianglemapHog = vtkPolyDataMapper.New();
aTrianglemapHog.SetInputConnection(aTrianglehog.GetOutputPort());
aTrianglemapHog.SetScalarModeToUseCellData();
aTrianglemapHog.ScalarVisibilityOff();
aTrianglehogActor = new vtkActor();
aTrianglehogActor.SetMapper(aTrianglemapHog);
aTrianglehogActor.GetProperty().SetColor(0,1,0);
aTriangleGlyph3D = new vtkGlyph3D();
aTriangleGlyph3D.SetInputConnection(aTriangleCenters.GetOutputPort());
aTriangleGlyph3D.SetSource(ball.GetOutput());
aTriangleCentersMapper = vtkPolyDataMapper.New();
aTriangleCentersMapper.SetInputConnection(aTriangleGlyph3D.GetOutputPort());
aTriangleCentersActor = new vtkActor();
aTriangleCentersActor.SetMapper(aTriangleCentersMapper);
aTrianglehogActor.SetPosition(aTriangleActor.GetPosition()[0],aTriangleActor.GetPosition()[1],aTriangleActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aTrianglehogActor);
aTrianglehogActor.GetProperty().SetRepresentationToWireframe();
aTriangleStripderivs = new vtkCellDerivatives();
aTriangleStripderivs.SetInput(aTriangleStripGrid);
aTriangleStripderivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aTriangleStrip";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aTriangleStripWriter = new vtkUnstructuredGridWriter();
aTriangleStripWriter.SetInputConnection(aTriangleStripderivs.GetOutputPort());
aTriangleStripWriter.SetFileName(FileName);
aTriangleStripWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aTriangleStripCenters = new vtkCellCenters();
aTriangleStripCenters.SetInputConnection(aTriangleStripderivs.GetOutputPort());
aTriangleStripCenters.VertexCellsOn();
aTriangleStriphog = new vtkHedgeHog();
aTriangleStriphog.SetInputConnection(aTriangleStripCenters.GetOutputPort());
aTriangleStripmapHog = vtkPolyDataMapper.New();
aTriangleStripmapHog.SetInputConnection(aTriangleStriphog.GetOutputPort());
aTriangleStripmapHog.SetScalarModeToUseCellData();
aTriangleStripmapHog.ScalarVisibilityOff();
aTriangleStriphogActor = new vtkActor();
aTriangleStriphogActor.SetMapper(aTriangleStripmapHog);
aTriangleStriphogActor.GetProperty().SetColor(0,1,0);
aTriangleStripGlyph3D = new vtkGlyph3D();
aTriangleStripGlyph3D.SetInputConnection(aTriangleStripCenters.GetOutputPort());
aTriangleStripGlyph3D.SetSource(ball.GetOutput());
aTriangleStripCentersMapper = vtkPolyDataMapper.New();
aTriangleStripCentersMapper.SetInputConnection(aTriangleStripGlyph3D.GetOutputPort());
aTriangleStripCentersActor = new vtkActor();
aTriangleStripCentersActor.SetMapper(aTriangleStripCentersMapper);
aTriangleStriphogActor.SetPosition(aTriangleStripActor.GetPosition()[0],aTriangleStripActor.GetPosition()[1],aTriangleStripActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aTriangleStriphogActor);
aTriangleStriphogActor.GetProperty().SetRepresentationToWireframe();
aLinederivs = new vtkCellDerivatives();
aLinederivs.SetInput(aLineGrid);
aLinederivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aLine";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aLineWriter = new vtkUnstructuredGridWriter();
aLineWriter.SetInputConnection(aLinederivs.GetOutputPort());
aLineWriter.SetFileName(FileName);
aLineWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aLineCenters = new vtkCellCenters();
aLineCenters.SetInputConnection(aLinederivs.GetOutputPort());
aLineCenters.VertexCellsOn();
aLinehog = new vtkHedgeHog();
aLinehog.SetInputConnection(aLineCenters.GetOutputPort());
aLinemapHog = vtkPolyDataMapper.New();
aLinemapHog.SetInputConnection(aLinehog.GetOutputPort());
aLinemapHog.SetScalarModeToUseCellData();
aLinemapHog.ScalarVisibilityOff();
aLinehogActor = new vtkActor();
aLinehogActor.SetMapper(aLinemapHog);
aLinehogActor.GetProperty().SetColor(0,1,0);
aLineGlyph3D = new vtkGlyph3D();
aLineGlyph3D.SetInputConnection(aLineCenters.GetOutputPort());
aLineGlyph3D.SetSource(ball.GetOutput());
aLineCentersMapper = vtkPolyDataMapper.New();
aLineCentersMapper.SetInputConnection(aLineGlyph3D.GetOutputPort());
aLineCentersActor = new vtkActor();
aLineCentersActor.SetMapper(aLineCentersMapper);
aLinehogActor.SetPosition(aLineActor.GetPosition()[0],aLineActor.GetPosition()[1],aLineActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aLinehogActor);
aLinehogActor.GetProperty().SetRepresentationToWireframe();
aPolyLinederivs = new vtkCellDerivatives();
aPolyLinederivs.SetInput(aPolyLineGrid);
aPolyLinederivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aPolyLine";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aPolyLineWriter = new vtkUnstructuredGridWriter();
aPolyLineWriter.SetInputConnection(aPolyLinederivs.GetOutputPort());
aPolyLineWriter.SetFileName(FileName);
aPolyLineWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aPolyLineCenters = new vtkCellCenters();
aPolyLineCenters.SetInputConnection(aPolyLinederivs.GetOutputPort());
aPolyLineCenters.VertexCellsOn();
aPolyLinehog = new vtkHedgeHog();
aPolyLinehog.SetInputConnection(aPolyLineCenters.GetOutputPort());
aPolyLinemapHog = vtkPolyDataMapper.New();
aPolyLinemapHog.SetInputConnection(aPolyLinehog.GetOutputPort());
aPolyLinemapHog.SetScalarModeToUseCellData();
aPolyLinemapHog.ScalarVisibilityOff();
aPolyLinehogActor = new vtkActor();
aPolyLinehogActor.SetMapper(aPolyLinemapHog);
aPolyLinehogActor.GetProperty().SetColor(0,1,0);
aPolyLineGlyph3D = new vtkGlyph3D();
aPolyLineGlyph3D.SetInputConnection(aPolyLineCenters.GetOutputPort());
aPolyLineGlyph3D.SetSource(ball.GetOutput());
aPolyLineCentersMapper = vtkPolyDataMapper.New();
aPolyLineCentersMapper.SetInputConnection(aPolyLineGlyph3D.GetOutputPort());
aPolyLineCentersActor = new vtkActor();
aPolyLineCentersActor.SetMapper(aPolyLineCentersMapper);
aPolyLinehogActor.SetPosition(aPolyLineActor.GetPosition()[0],aPolyLineActor.GetPosition()[1],aPolyLineActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aPolyLinehogActor);
aPolyLinehogActor.GetProperty().SetRepresentationToWireframe();
aVertexderivs = new vtkCellDerivatives();
aVertexderivs.SetInput(aVertexGrid);
aVertexderivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aVertex";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aVertexWriter = new vtkUnstructuredGridWriter();
aVertexWriter.SetInputConnection(aVertexderivs.GetOutputPort());
aVertexWriter.SetFileName(FileName);
aVertexWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aVertexCenters = new vtkCellCenters();
aVertexCenters.SetInputConnection(aVertexderivs.GetOutputPort());
aVertexCenters.VertexCellsOn();
aVertexhog = new vtkHedgeHog();
aVertexhog.SetInputConnection(aVertexCenters.GetOutputPort());
aVertexmapHog = vtkPolyDataMapper.New();
aVertexmapHog.SetInputConnection(aVertexhog.GetOutputPort());
aVertexmapHog.SetScalarModeToUseCellData();
aVertexmapHog.ScalarVisibilityOff();
aVertexhogActor = new vtkActor();
aVertexhogActor.SetMapper(aVertexmapHog);
aVertexhogActor.GetProperty().SetColor(0,1,0);
aVertexGlyph3D = new vtkGlyph3D();
aVertexGlyph3D.SetInputConnection(aVertexCenters.GetOutputPort());
aVertexGlyph3D.SetSource(ball.GetOutput());
aVertexCentersMapper = vtkPolyDataMapper.New();
aVertexCentersMapper.SetInputConnection(aVertexGlyph3D.GetOutputPort());
aVertexCentersActor = new vtkActor();
aVertexCentersActor.SetMapper(aVertexCentersMapper);
aVertexhogActor.SetPosition(aVertexActor.GetPosition()[0],aVertexActor.GetPosition()[1],aVertexActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aVertexhogActor);
aVertexhogActor.GetProperty().SetRepresentationToWireframe();
aPolyVertexderivs = new vtkCellDerivatives();
aPolyVertexderivs.SetInput(aPolyVertexGrid);
aPolyVertexderivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aPolyVertex";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aPolyVertexWriter = new vtkUnstructuredGridWriter();
aPolyVertexWriter.SetInputConnection(aPolyVertexderivs.GetOutputPort());
aPolyVertexWriter.SetFileName(FileName);
aPolyVertexWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aPolyVertexCenters = new vtkCellCenters();
aPolyVertexCenters.SetInputConnection(aPolyVertexderivs.GetOutputPort());
aPolyVertexCenters.VertexCellsOn();
aPolyVertexhog = new vtkHedgeHog();
aPolyVertexhog.SetInputConnection(aPolyVertexCenters.GetOutputPort());
aPolyVertexmapHog = vtkPolyDataMapper.New();
aPolyVertexmapHog.SetInputConnection(aPolyVertexhog.GetOutputPort());
aPolyVertexmapHog.SetScalarModeToUseCellData();
aPolyVertexmapHog.ScalarVisibilityOff();
aPolyVertexhogActor = new vtkActor();
aPolyVertexhogActor.SetMapper(aPolyVertexmapHog);
aPolyVertexhogActor.GetProperty().SetColor(0,1,0);
aPolyVertexGlyph3D = new vtkGlyph3D();
aPolyVertexGlyph3D.SetInputConnection(aPolyVertexCenters.GetOutputPort());
aPolyVertexGlyph3D.SetSource(ball.GetOutput());
aPolyVertexCentersMapper = vtkPolyDataMapper.New();
aPolyVertexCentersMapper.SetInputConnection(aPolyVertexGlyph3D.GetOutputPort());
aPolyVertexCentersActor = new vtkActor();
aPolyVertexCentersActor.SetMapper(aPolyVertexCentersMapper);
aPolyVertexhogActor.SetPosition(aPolyVertexActor.GetPosition()[0],aPolyVertexActor.GetPosition()[1],aPolyVertexActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aPolyVertexhogActor);
aPolyVertexhogActor.GetProperty().SetRepresentationToWireframe();
aPixelderivs = new vtkCellDerivatives();
aPixelderivs.SetInput(aPixelGrid);
aPixelderivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aPixel";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aPixelWriter = new vtkUnstructuredGridWriter();
aPixelWriter.SetInputConnection(aPixelderivs.GetOutputPort());
aPixelWriter.SetFileName(FileName);
aPixelWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aPixelCenters = new vtkCellCenters();
aPixelCenters.SetInputConnection(aPixelderivs.GetOutputPort());
aPixelCenters.VertexCellsOn();
aPixelhog = new vtkHedgeHog();
aPixelhog.SetInputConnection(aPixelCenters.GetOutputPort());
aPixelmapHog = vtkPolyDataMapper.New();
aPixelmapHog.SetInputConnection(aPixelhog.GetOutputPort());
aPixelmapHog.SetScalarModeToUseCellData();
aPixelmapHog.ScalarVisibilityOff();
aPixelhogActor = new vtkActor();
aPixelhogActor.SetMapper(aPixelmapHog);
aPixelhogActor.GetProperty().SetColor(0,1,0);
aPixelGlyph3D = new vtkGlyph3D();
aPixelGlyph3D.SetInputConnection(aPixelCenters.GetOutputPort());
aPixelGlyph3D.SetSource(ball.GetOutput());
aPixelCentersMapper = vtkPolyDataMapper.New();
aPixelCentersMapper.SetInputConnection(aPixelGlyph3D.GetOutputPort());
aPixelCentersActor = new vtkActor();
aPixelCentersActor.SetMapper(aPixelCentersMapper);
aPixelhogActor.SetPosition(aPixelActor.GetPosition()[0],aPixelActor.GetPosition()[1],aPixelActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aPixelhogActor);
aPixelhogActor.GetProperty().SetRepresentationToWireframe();
aPolygonderivs = new vtkCellDerivatives();
aPolygonderivs.SetInput(aPolygonGrid);
aPolygonderivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aPolygon";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aPolygonWriter = new vtkUnstructuredGridWriter();
aPolygonWriter.SetInputConnection(aPolygonderivs.GetOutputPort());
aPolygonWriter.SetFileName(FileName);
aPolygonWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aPolygonCenters = new vtkCellCenters();
aPolygonCenters.SetInputConnection(aPolygonderivs.GetOutputPort());
aPolygonCenters.VertexCellsOn();
aPolygonhog = new vtkHedgeHog();
aPolygonhog.SetInputConnection(aPolygonCenters.GetOutputPort());
aPolygonmapHog = vtkPolyDataMapper.New();
aPolygonmapHog.SetInputConnection(aPolygonhog.GetOutputPort());
aPolygonmapHog.SetScalarModeToUseCellData();
aPolygonmapHog.ScalarVisibilityOff();
aPolygonhogActor = new vtkActor();
aPolygonhogActor.SetMapper(aPolygonmapHog);
aPolygonhogActor.GetProperty().SetColor(0,1,0);
aPolygonGlyph3D = new vtkGlyph3D();
aPolygonGlyph3D.SetInputConnection(aPolygonCenters.GetOutputPort());
aPolygonGlyph3D.SetSource(ball.GetOutput());
aPolygonCentersMapper = vtkPolyDataMapper.New();
aPolygonCentersMapper.SetInputConnection(aPolygonGlyph3D.GetOutputPort());
aPolygonCentersActor = new vtkActor();
aPolygonCentersActor.SetMapper(aPolygonCentersMapper);
aPolygonhogActor.SetPosition(aPolygonActor.GetPosition()[0],aPolygonActor.GetPosition()[1],aPolygonActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aPolygonhogActor);
aPolygonhogActor.GetProperty().SetRepresentationToWireframe();
aPentaderivs = new vtkCellDerivatives();
aPentaderivs.SetInput(aPentaGrid);
aPentaderivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aPenta";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aPentaWriter = new vtkUnstructuredGridWriter();
aPentaWriter.SetInputConnection(aPentaderivs.GetOutputPort());
aPentaWriter.SetFileName(FileName);
aPentaWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aPentaCenters = new vtkCellCenters();
aPentaCenters.SetInputConnection(aPentaderivs.GetOutputPort());
aPentaCenters.VertexCellsOn();
aPentahog = new vtkHedgeHog();
aPentahog.SetInputConnection(aPentaCenters.GetOutputPort());
aPentamapHog = vtkPolyDataMapper.New();
aPentamapHog.SetInputConnection(aPentahog.GetOutputPort());
aPentamapHog.SetScalarModeToUseCellData();
aPentamapHog.ScalarVisibilityOff();
aPentahogActor = new vtkActor();
aPentahogActor.SetMapper(aPentamapHog);
aPentahogActor.GetProperty().SetColor(0,1,0);
aPentaGlyph3D = new vtkGlyph3D();
aPentaGlyph3D.SetInputConnection(aPentaCenters.GetOutputPort());
aPentaGlyph3D.SetSource(ball.GetOutput());
aPentaCentersMapper = vtkPolyDataMapper.New();
aPentaCentersMapper.SetInputConnection(aPentaGlyph3D.GetOutputPort());
aPentaCentersActor = new vtkActor();
aPentaCentersActor.SetMapper(aPentaCentersMapper);
aPentahogActor.SetPosition(aPentaActor.GetPosition()[0],aPentaActor.GetPosition()[1],aPentaActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aPentahogActor);
aPentahogActor.GetProperty().SetRepresentationToWireframe();
aHexaderivs = new vtkCellDerivatives();
aHexaderivs.SetInput(aHexaGrid);
aHexaderivs.SetVectorModeToComputeGradient();
FileName = dir;
FileName += "/aHexa";
FileName += ".vtk";
// make sure the directory is writeable first[]
tryWorked = false;
try
{
channel = new StreamWriter("" + (dir.ToString()) + "/test.tmp");
tryWorked = true;
}
catch(Exception)
{
tryWorked = false;
}
if(tryWorked)
{
channel.Close();
File.Delete("" + (dir.ToString()) + "/test.tmp");
aHexaWriter = new vtkUnstructuredGridWriter();
aHexaWriter.SetInputConnection(aHexaderivs.GetOutputPort());
aHexaWriter.SetFileName(FileName);
aHexaWriter.Write();
// delete the file[]
File.Delete("FileName");
}
aHexaCenters = new vtkCellCenters();
aHexaCenters.SetInputConnection(aHexaderivs.GetOutputPort());
aHexaCenters.VertexCellsOn();
aHexahog = new vtkHedgeHog();
aHexahog.SetInputConnection(aHexaCenters.GetOutputPort());
aHexamapHog = vtkPolyDataMapper.New();
aHexamapHog.SetInputConnection(aHexahog.GetOutputPort());
aHexamapHog.SetScalarModeToUseCellData();
aHexamapHog.ScalarVisibilityOff();
aHexahogActor = new vtkActor();
aHexahogActor.SetMapper(aHexamapHog);
aHexahogActor.GetProperty().SetColor(0,1,0);
aHexaGlyph3D = new vtkGlyph3D();
aHexaGlyph3D.SetInputConnection(aHexaCenters.GetOutputPort());
aHexaGlyph3D.SetSource(ball.GetOutput());
aHexaCentersMapper = vtkPolyDataMapper.New();
aHexaCentersMapper.SetInputConnection(aHexaGlyph3D.GetOutputPort());
aHexaCentersActor = new vtkActor();
aHexaCentersActor.SetMapper(aHexaCentersMapper);
aHexahogActor.SetPosition(aHexaActor.GetPosition()[0],aHexaActor.GetPosition()[1],aHexaActor.GetPosition()[2]);
ren1.AddActor((vtkProp)aHexahogActor);
aHexahogActor.GetProperty().SetRepresentationToWireframe();
ren1.ResetCamera();
ren1.GetActiveCamera().Azimuth((double)30);
ren1.GetActiveCamera().Elevation((double)20);
ren1.GetActiveCamera().Dolly((double)3.0);
ren1.ResetCameraClippingRange();
renWin.SetSize((int)300,(int)150);
renWin.Render();
// render the image[]
//[]
iren.Initialize();
//deleteAllVTKObjects();
}
/// <summary>
/// The main entry method called by the CSharp driver
/// </summary>
/// <param name="argv"></param>
public static void AVLineIntersectQuadraticCells(String [] argv)
{
//Prefix Content is: ""
// Contour every quadratic cell type[]
// Create a scene with one of each cell type.[]
// QuadraticEdge[]
edgePoints = new vtkPoints();
edgePoints.SetNumberOfPoints((int)3);
edgePoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
edgePoints.InsertPoint((int)1,(double)1.0,(double)0,(double)0);
edgePoints.InsertPoint((int)2,(double)0.5,(double)0.25,(double)0);
edgeScalars = new vtkFloatArray();
edgeScalars.SetNumberOfTuples((int)3);
edgeScalars.InsertValue((int)0,(float)0.0);
edgeScalars.InsertValue((int)1,(float)0.0);
edgeScalars.InsertValue((int)2,(float)0.9);
aEdge = new vtkQuadraticEdge();
aEdge.GetPointIds().SetId((int)0,(int)0);
aEdge.GetPointIds().SetId((int)1,(int)1);
aEdge.GetPointIds().SetId((int)2,(int)2);
aEdgeGrid = new vtkUnstructuredGrid();
aEdgeGrid.Allocate((int)1,(int)1);
aEdgeGrid.InsertNextCell((int)aEdge.GetCellType(),(vtkIdList)aEdge.GetPointIds());
aEdgeGrid.SetPoints((vtkPoints)edgePoints);
aEdgeGrid.GetPointData().SetScalars((vtkDataArray)edgeScalars);
aEdgeMapper = new vtkDataSetMapper();
aEdgeMapper.SetInputData((vtkDataSet)aEdgeGrid);
aEdgeMapper.ScalarVisibilityOff();
aEdgeActor = new vtkActor();
aEdgeActor.SetMapper((vtkMapper)aEdgeMapper);
aEdgeActor.GetProperty().SetRepresentationToWireframe();
aEdgeActor.GetProperty().SetAmbient((double)1.0);
// Quadratic triangle[]
triPoints = new vtkPoints();
triPoints.SetNumberOfPoints((int)6);
triPoints.InsertPoint((int)0,(double)0.0,(double)0.0,(double)0.0);
triPoints.InsertPoint((int)1,(double)1.0,(double)0.0,(double)0.0);
triPoints.InsertPoint((int)2,(double)0.5,(double)0.8,(double)0.0);
triPoints.InsertPoint((int)3,(double)0.5,(double)0.0,(double)0.0);
triPoints.InsertPoint((int)4,(double)0.75,(double)0.4,(double)0.0);
triPoints.InsertPoint((int)5,(double)0.25,(double)0.4,(double)0.0);
triScalars = new vtkFloatArray();
triScalars.SetNumberOfTuples((int)6);
triScalars.InsertValue((int)0,(float)0.0);
triScalars.InsertValue((int)1,(float)0.0);
triScalars.InsertValue((int)2,(float)0.0);
triScalars.InsertValue((int)3,(float)1.0);
triScalars.InsertValue((int)4,(float)0.0);
triScalars.InsertValue((int)5,(float)0.0);
aTri = new vtkQuadraticTriangle();
aTri.GetPointIds().SetId((int)0,(int)0);
aTri.GetPointIds().SetId((int)1,(int)1);
aTri.GetPointIds().SetId((int)2,(int)2);
aTri.GetPointIds().SetId((int)3,(int)3);
aTri.GetPointIds().SetId((int)4,(int)4);
aTri.GetPointIds().SetId((int)5,(int)5);
aTriGrid = new vtkUnstructuredGrid();
aTriGrid.Allocate((int)1,(int)1);
aTriGrid.InsertNextCell((int)aTri.GetCellType(),(vtkIdList)aTri.GetPointIds());
aTriGrid.SetPoints((vtkPoints)triPoints);
aTriGrid.GetPointData().SetScalars((vtkDataArray)triScalars);
aTriMapper = new vtkDataSetMapper();
aTriMapper.SetInputData((vtkDataSet)aTriGrid);
aTriMapper.ScalarVisibilityOff();
aTriActor = new vtkActor();
aTriActor.SetMapper((vtkMapper)aTriMapper);
aTriActor.GetProperty().SetRepresentationToWireframe();
aTriActor.GetProperty().SetAmbient((double)1.0);
// Quadratic quadrilateral[]
quadPoints = new vtkPoints();
quadPoints.SetNumberOfPoints((int)8);
quadPoints.InsertPoint((int)0,(double)0.0,(double)0.0,(double)0.0);
quadPoints.InsertPoint((int)1,(double)1.0,(double)0.0,(double)0.0);
quadPoints.InsertPoint((int)2,(double)1.0,(double)1.0,(double)0.0);
quadPoints.InsertPoint((int)3,(double)0.0,(double)1.0,(double)0.0);
quadPoints.InsertPoint((int)4,(double)0.5,(double)0.0,(double)0.0);
quadPoints.InsertPoint((int)5,(double)1.0,(double)0.5,(double)0.0);
quadPoints.InsertPoint((int)6,(double)0.5,(double)1.0,(double)0.0);
quadPoints.InsertPoint((int)7,(double)0.0,(double)0.5,(double)0.0);
quadScalars = new vtkFloatArray();
quadScalars.SetNumberOfTuples((int)8);
quadScalars.InsertValue((int)0,(float)0.0);
quadScalars.InsertValue((int)1,(float)0.0);
quadScalars.InsertValue((int)2,(float)1.0);
quadScalars.InsertValue((int)3,(float)1.0);
quadScalars.InsertValue((int)4,(float)1.0);
quadScalars.InsertValue((int)5,(float)0.0);
quadScalars.InsertValue((int)6,(float)0.0);
quadScalars.InsertValue((int)7,(float)0.0);
aQuad = new vtkQuadraticQuad();
aQuad.GetPointIds().SetId((int)0,(int)0);
aQuad.GetPointIds().SetId((int)1,(int)1);
aQuad.GetPointIds().SetId((int)2,(int)2);
aQuad.GetPointIds().SetId((int)3,(int)3);
aQuad.GetPointIds().SetId((int)4,(int)4);
aQuad.GetPointIds().SetId((int)5,(int)5);
aQuad.GetPointIds().SetId((int)6,(int)6);
aQuad.GetPointIds().SetId((int)7,(int)7);
aQuadGrid = new vtkUnstructuredGrid();
aQuadGrid.Allocate((int)1,(int)1);
aQuadGrid.InsertNextCell((int)aQuad.GetCellType(),(vtkIdList)aQuad.GetPointIds());
aQuadGrid.SetPoints((vtkPoints)quadPoints);
aQuadGrid.GetPointData().SetScalars((vtkDataArray)quadScalars);
aQuadMapper = new vtkDataSetMapper();
aQuadMapper.SetInputData((vtkDataSet)aQuadGrid);
aQuadMapper.ScalarVisibilityOff();
aQuadActor = new vtkActor();
aQuadActor.SetMapper((vtkMapper)aQuadMapper);
aQuadActor.GetProperty().SetRepresentationToWireframe();
aQuadActor.GetProperty().SetAmbient((double)1.0);
// Quadratic tetrahedron[]
tetPoints = new vtkPoints();
tetPoints.SetNumberOfPoints((int)10);
tetPoints.InsertPoint((int)0,(double)0.0,(double)0.0,(double)0.0);
tetPoints.InsertPoint((int)1,(double)1.0,(double)0.0,(double)0.0);
tetPoints.InsertPoint((int)2,(double)0.5,(double)0.8,(double)0.0);
tetPoints.InsertPoint((int)3,(double)0.5,(double)0.4,(double)1.0);
tetPoints.InsertPoint((int)4,(double)0.5,(double)0.0,(double)0.0);
tetPoints.InsertPoint((int)5,(double)0.75,(double)0.4,(double)0.0);
tetPoints.InsertPoint((int)6,(double)0.25,(double)0.4,(double)0.0);
tetPoints.InsertPoint((int)7,(double)0.25,(double)0.2,(double)0.5);
tetPoints.InsertPoint((int)8,(double)0.75,(double)0.2,(double)0.5);
tetPoints.InsertPoint((int)9,(double)0.50,(double)0.6,(double)0.5);
tetScalars = new vtkFloatArray();
tetScalars.SetNumberOfTuples((int)10);
tetScalars.InsertValue((int)0,(float)1.0);
tetScalars.InsertValue((int)1,(float)1.0);
tetScalars.InsertValue((int)2,(float)1.0);
tetScalars.InsertValue((int)3,(float)1.0);
tetScalars.InsertValue((int)4,(float)0.0);
tetScalars.InsertValue((int)5,(float)0.0);
tetScalars.InsertValue((int)6,(float)0.0);
tetScalars.InsertValue((int)7,(float)0.0);
tetScalars.InsertValue((int)8,(float)0.0);
tetScalars.InsertValue((int)9,(float)0.0);
aTet = new vtkQuadraticTetra();
aTet.GetPointIds().SetId((int)0,(int)0);
aTet.GetPointIds().SetId((int)1,(int)1);
aTet.GetPointIds().SetId((int)2,(int)2);
aTet.GetPointIds().SetId((int)3,(int)3);
aTet.GetPointIds().SetId((int)4,(int)4);
aTet.GetPointIds().SetId((int)5,(int)5);
aTet.GetPointIds().SetId((int)6,(int)6);
aTet.GetPointIds().SetId((int)7,(int)7);
aTet.GetPointIds().SetId((int)8,(int)8);
aTet.GetPointIds().SetId((int)9,(int)9);
aTetGrid = new vtkUnstructuredGrid();
aTetGrid.Allocate((int)1,(int)1);
aTetGrid.InsertNextCell((int)aTet.GetCellType(),(vtkIdList)aTet.GetPointIds());
aTetGrid.SetPoints((vtkPoints)tetPoints);
aTetGrid.GetPointData().SetScalars((vtkDataArray)tetScalars);
aTetMapper = new vtkDataSetMapper();
aTetMapper.SetInputData((vtkDataSet)aTetGrid);
aTetMapper.ScalarVisibilityOff();
aTetActor = new vtkActor();
aTetActor.SetMapper((vtkMapper)aTetMapper);
aTetActor.GetProperty().SetRepresentationToWireframe();
aTetActor.GetProperty().SetAmbient((double)1.0);
// Quadratic hexahedron[]
hexPoints = new vtkPoints();
hexPoints.SetNumberOfPoints((int)20);
hexPoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
hexPoints.InsertPoint((int)1,(double)1,(double)0,(double)0);
hexPoints.InsertPoint((int)2,(double)1,(double)1,(double)0);
hexPoints.InsertPoint((int)3,(double)0,(double)1,(double)0);
hexPoints.InsertPoint((int)4,(double)0,(double)0,(double)1);
hexPoints.InsertPoint((int)5,(double)1,(double)0,(double)1);
hexPoints.InsertPoint((int)6,(double)1,(double)1,(double)1);
hexPoints.InsertPoint((int)7,(double)0,(double)1,(double)1);
hexPoints.InsertPoint((int)8,(double)0.5,(double)0,(double)0);
hexPoints.InsertPoint((int)9,(double)1,(double)0.5,(double)0);
hexPoints.InsertPoint((int)10,(double)0.5,(double)1,(double)0);
hexPoints.InsertPoint((int)11,(double)0,(double)0.5,(double)0);
hexPoints.InsertPoint((int)12,(double)0.5,(double)0,(double)1);
hexPoints.InsertPoint((int)13,(double)1,(double)0.5,(double)1);
hexPoints.InsertPoint((int)14,(double)0.5,(double)1,(double)1);
hexPoints.InsertPoint((int)15,(double)0,(double)0.5,(double)1);
hexPoints.InsertPoint((int)16,(double)0,(double)0,(double)0.5);
hexPoints.InsertPoint((int)17,(double)1,(double)0,(double)0.5);
hexPoints.InsertPoint((int)18,(double)1,(double)1,(double)0.5);
hexPoints.InsertPoint((int)19,(double)0,(double)1,(double)0.5);
hexScalars = new vtkFloatArray();
hexScalars.SetNumberOfTuples((int)20);
hexScalars.InsertValue((int)0,(float)1.0);
hexScalars.InsertValue((int)1,(float)1.0);
hexScalars.InsertValue((int)2,(float)1.0);
hexScalars.InsertValue((int)3,(float)1.0);
hexScalars.InsertValue((int)4,(float)1.0);
hexScalars.InsertValue((int)5,(float)1.0);
hexScalars.InsertValue((int)6,(float)1.0);
hexScalars.InsertValue((int)7,(float)1.0);
hexScalars.InsertValue((int)8,(float)0.0);
hexScalars.InsertValue((int)9,(float)0.0);
hexScalars.InsertValue((int)10,(float)0.0);
hexScalars.InsertValue((int)11,(float)0.0);
hexScalars.InsertValue((int)12,(float)0.0);
hexScalars.InsertValue((int)13,(float)0.0);
hexScalars.InsertValue((int)14,(float)0.0);
hexScalars.InsertValue((int)15,(float)0.0);
hexScalars.InsertValue((int)16,(float)0.0);
hexScalars.InsertValue((int)17,(float)0.0);
hexScalars.InsertValue((int)18,(float)0.0);
hexScalars.InsertValue((int)19,(float)0.0);
aHex = new vtkQuadraticHexahedron();
aHex.GetPointIds().SetId((int)0,(int)0);
aHex.GetPointIds().SetId((int)1,(int)1);
aHex.GetPointIds().SetId((int)2,(int)2);
aHex.GetPointIds().SetId((int)3,(int)3);
aHex.GetPointIds().SetId((int)4,(int)4);
aHex.GetPointIds().SetId((int)5,(int)5);
aHex.GetPointIds().SetId((int)6,(int)6);
aHex.GetPointIds().SetId((int)7,(int)7);
aHex.GetPointIds().SetId((int)8,(int)8);
aHex.GetPointIds().SetId((int)9,(int)9);
aHex.GetPointIds().SetId((int)10,(int)10);
aHex.GetPointIds().SetId((int)11,(int)11);
aHex.GetPointIds().SetId((int)12,(int)12);
aHex.GetPointIds().SetId((int)13,(int)13);
aHex.GetPointIds().SetId((int)14,(int)14);
aHex.GetPointIds().SetId((int)15,(int)15);
aHex.GetPointIds().SetId((int)16,(int)16);
aHex.GetPointIds().SetId((int)17,(int)17);
aHex.GetPointIds().SetId((int)18,(int)18);
aHex.GetPointIds().SetId((int)19,(int)19);
aHexGrid = new vtkUnstructuredGrid();
aHexGrid.Allocate((int)1,(int)1);
aHexGrid.InsertNextCell((int)aHex.GetCellType(),(vtkIdList)aHex.GetPointIds());
aHexGrid.SetPoints((vtkPoints)hexPoints);
aHexGrid.GetPointData().SetScalars((vtkDataArray)hexScalars);
aHexMapper = new vtkDataSetMapper();
aHexMapper.SetInputData((vtkDataSet)aHexGrid);
aHexMapper.ScalarVisibilityOff();
aHexActor = new vtkActor();
aHexActor.SetMapper((vtkMapper)aHexMapper);
aHexActor.GetProperty().SetRepresentationToWireframe();
aHexActor.GetProperty().SetAmbient((double)1.0);
// Quadratic wedge[]
wedgePoints = new vtkPoints();
wedgePoints.SetNumberOfPoints((int)15);
wedgePoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
wedgePoints.InsertPoint((int)1,(double)1,(double)0,(double)0);
wedgePoints.InsertPoint((int)2,(double)0,(double)1,(double)0);
wedgePoints.InsertPoint((int)3,(double)0,(double)0,(double)1);
wedgePoints.InsertPoint((int)4,(double)1,(double)0,(double)1);
wedgePoints.InsertPoint((int)5,(double)0,(double)1,(double)1);
wedgePoints.InsertPoint((int)6,(double)0.5,(double)0,(double)0);
wedgePoints.InsertPoint((int)7,(double)0.5,(double)0.5,(double)0);
wedgePoints.InsertPoint((int)8,(double)0,(double)0.5,(double)0);
wedgePoints.InsertPoint((int)9,(double)0.5,(double)0,(double)1);
wedgePoints.InsertPoint((int)10,(double)0.5,(double)0.5,(double)1);
wedgePoints.InsertPoint((int)11,(double)0,(double)0.5,(double)1);
wedgePoints.InsertPoint((int)12,(double)0,(double)0,(double)0.5);
wedgePoints.InsertPoint((int)13,(double)1,(double)0,(double)0.5);
wedgePoints.InsertPoint((int)14,(double)0,(double)1,(double)0.5);
wedgeScalars = new vtkFloatArray();
wedgeScalars.SetNumberOfTuples((int)15);
wedgeScalars.InsertValue((int)0,(float)1.0);
wedgeScalars.InsertValue((int)1,(float)1.0);
wedgeScalars.InsertValue((int)2,(float)1.0);
wedgeScalars.InsertValue((int)3,(float)1.0);
wedgeScalars.InsertValue((int)4,(float)1.0);
wedgeScalars.InsertValue((int)5,(float)1.0);
wedgeScalars.InsertValue((int)6,(float)1.0);
wedgeScalars.InsertValue((int)7,(float)1.0);
wedgeScalars.InsertValue((int)8,(float)0.0);
wedgeScalars.InsertValue((int)9,(float)0.0);
wedgeScalars.InsertValue((int)10,(float)0.0);
wedgeScalars.InsertValue((int)11,(float)0.0);
wedgeScalars.InsertValue((int)12,(float)0.0);
wedgeScalars.InsertValue((int)13,(float)0.0);
wedgeScalars.InsertValue((int)14,(float)0.0);
aWedge = new vtkQuadraticWedge();
aWedge.GetPointIds().SetId((int)0,(int)0);
aWedge.GetPointIds().SetId((int)1,(int)1);
aWedge.GetPointIds().SetId((int)2,(int)2);
aWedge.GetPointIds().SetId((int)3,(int)3);
aWedge.GetPointIds().SetId((int)4,(int)4);
aWedge.GetPointIds().SetId((int)5,(int)5);
aWedge.GetPointIds().SetId((int)6,(int)6);
aWedge.GetPointIds().SetId((int)7,(int)7);
aWedge.GetPointIds().SetId((int)8,(int)8);
aWedge.GetPointIds().SetId((int)9,(int)9);
aWedge.GetPointIds().SetId((int)10,(int)10);
aWedge.GetPointIds().SetId((int)11,(int)11);
aWedge.GetPointIds().SetId((int)12,(int)12);
aWedge.GetPointIds().SetId((int)13,(int)13);
aWedge.GetPointIds().SetId((int)14,(int)14);
aWedgeGrid = new vtkUnstructuredGrid();
aWedgeGrid.Allocate((int)1,(int)1);
aWedgeGrid.InsertNextCell((int)aWedge.GetCellType(),(vtkIdList)aWedge.GetPointIds());
aWedgeGrid.SetPoints((vtkPoints)wedgePoints);
aWedgeGrid.GetPointData().SetScalars((vtkDataArray)wedgeScalars);
wedgeContours = new vtkClipDataSet();
wedgeContours.SetInputData((vtkDataObject)aWedgeGrid);
wedgeContours.SetValue((double)0.5);
aWedgeContourMapper = new vtkDataSetMapper();
aWedgeContourMapper.SetInputConnection((vtkAlgorithmOutput)wedgeContours.GetOutputPort());
aWedgeContourMapper.ScalarVisibilityOff();
aWedgeMapper = new vtkDataSetMapper();
aWedgeMapper.SetInputData((vtkDataSet)aWedgeGrid);
aWedgeMapper.ScalarVisibilityOff();
aWedgeActor = new vtkActor();
aWedgeActor.SetMapper((vtkMapper)aWedgeMapper);
aWedgeActor.GetProperty().SetRepresentationToWireframe();
aWedgeActor.GetProperty().SetAmbient((double)1.0);
aWedgeContourActor = new vtkActor();
aWedgeContourActor.SetMapper((vtkMapper)aWedgeContourMapper);
aWedgeContourActor.GetProperty().SetAmbient((double)1.0);
// Quadratic pyramid[]
pyraPoints = new vtkPoints();
pyraPoints.SetNumberOfPoints((int)13);
pyraPoints.InsertPoint((int)0,(double)0,(double)0,(double)0);
pyraPoints.InsertPoint((int)1,(double)1,(double)0,(double)0);
pyraPoints.InsertPoint((int)2,(double)1,(double)1,(double)0);
pyraPoints.InsertPoint((int)3,(double)0,(double)1,(double)0);
pyraPoints.InsertPoint((int)4,(double)0,(double)0,(double)1);
pyraPoints.InsertPoint((int)5,(double)0.5,(double)0,(double)0);
pyraPoints.InsertPoint((int)6,(double)1,(double)0.5,(double)0);
pyraPoints.InsertPoint((int)7,(double)0.5,(double)1,(double)0);
pyraPoints.InsertPoint((int)8,(double)0,(double)0.5,(double)0);
pyraPoints.InsertPoint((int)9,(double)0,(double)0,(double)0.5);
pyraPoints.InsertPoint((int)10,(double)0.5,(double)0,(double)0.5);
pyraPoints.InsertPoint((int)11,(double)0.5,(double)0.5,(double)0.5);
pyraPoints.InsertPoint((int)12,(double)0,(double)0.5,(double)0.5);
pyraScalars = new vtkFloatArray();
pyraScalars.SetNumberOfTuples((int)13);
pyraScalars.InsertValue((int)0,(float)1.0);
pyraScalars.InsertValue((int)1,(float)1.0);
pyraScalars.InsertValue((int)2,(float)1.0);
pyraScalars.InsertValue((int)3,(float)1.0);
pyraScalars.InsertValue((int)4,(float)1.0);
pyraScalars.InsertValue((int)5,(float)1.0);
pyraScalars.InsertValue((int)6,(float)1.0);
pyraScalars.InsertValue((int)7,(float)1.0);
pyraScalars.InsertValue((int)8,(float)0.0);
pyraScalars.InsertValue((int)9,(float)0.0);
pyraScalars.InsertValue((int)10,(float)0.0);
pyraScalars.InsertValue((int)11,(float)0.0);
pyraScalars.InsertValue((int)12,(float)0.0);
aPyramid = new vtkQuadraticPyramid();
aPyramid.GetPointIds().SetId((int)0,(int)0);
aPyramid.GetPointIds().SetId((int)1,(int)1);
aPyramid.GetPointIds().SetId((int)2,(int)2);
aPyramid.GetPointIds().SetId((int)3,(int)3);
aPyramid.GetPointIds().SetId((int)4,(int)4);
aPyramid.GetPointIds().SetId((int)5,(int)5);
aPyramid.GetPointIds().SetId((int)6,(int)6);
aPyramid.GetPointIds().SetId((int)7,(int)7);
aPyramid.GetPointIds().SetId((int)8,(int)8);
aPyramid.GetPointIds().SetId((int)9,(int)9);
aPyramid.GetPointIds().SetId((int)10,(int)10);
aPyramid.GetPointIds().SetId((int)11,(int)11);
aPyramid.GetPointIds().SetId((int)12,(int)12);
aPyramidGrid = new vtkUnstructuredGrid();
aPyramidGrid.Allocate((int)1,(int)1);
aPyramidGrid.InsertNextCell((int)aPyramid.GetCellType(),(vtkIdList)aPyramid.GetPointIds());
aPyramidGrid.SetPoints((vtkPoints)pyraPoints);
aPyramidGrid.GetPointData().SetScalars((vtkDataArray)pyraScalars);
pyraContours = new vtkClipDataSet();
pyraContours.SetInputData((vtkDataObject)aPyramidGrid);
pyraContours.SetValue((double)0.5);
aPyramidContourMapper = new vtkDataSetMapper();
aPyramidContourMapper.SetInputConnection((vtkAlgorithmOutput)pyraContours.GetOutputPort());
aPyramidContourMapper.ScalarVisibilityOff();
aPyramidMapper = new vtkDataSetMapper();
aPyramidMapper.SetInputData((vtkDataSet)aPyramidGrid);
aPyramidMapper.ScalarVisibilityOff();
aPyramidActor = new vtkActor();
aPyramidActor.SetMapper((vtkMapper)aPyramidMapper);
aPyramidActor.GetProperty().SetRepresentationToWireframe();
aPyramidActor.GetProperty().SetAmbient((double)1.0);
aPyramidContourActor = new vtkActor();
aPyramidContourActor.SetMapper((vtkMapper)aPyramidContourMapper);
aPyramidContourActor.GetProperty().SetAmbient((double)1.0);
// Create the rendering related stuff.[]
// Since some of our actors are a single vertex, we need to remove all[]
// cullers so the single vertex actors will render[]
ren1 = vtkRenderer.New();
ren1.GetCullers().RemoveAllItems();
renWin = vtkRenderWindow.New();
renWin.SetMultiSamples(0);
renWin.AddRenderer((vtkRenderer)ren1);
iren = new vtkRenderWindowInteractor();
iren.SetRenderWindow((vtkRenderWindow)renWin);
ren1.SetBackground((double).1,(double).2,(double).3);
renWin.SetSize((int)400,(int)200);
// specify properties[]
ren1.AddActor((vtkProp)aEdgeActor);
ren1.AddActor((vtkProp)aTriActor);
ren1.AddActor((vtkProp)aQuadActor);
ren1.AddActor((vtkProp)aTetActor);
ren1.AddActor((vtkProp)aHexActor);
ren1.AddActor((vtkProp)aWedgeActor);
ren1.AddActor((vtkProp)aPyramidActor);
// places everyone!![]
aTriActor.AddPosition((double)2,(double)0,(double)0);
aQuadActor.AddPosition((double)4,(double)0,(double)0);
aTetActor.AddPosition((double)6,(double)0,(double)0);
aHexActor.AddPosition((double)8,(double)0,(double)0);
aWedgeActor.AddPosition((double)10,(double)0,(double)0);
aPyramidActor.AddPosition((double)12,(double)0,(double)0);
BuildBackdrop(-1, 15, -1, 4, -1, 2, .1);
ren1.AddActor((vtkProp)base1);
base1.GetProperty().SetDiffuseColor((double).2,(double).2,(double).2);
ren1.AddActor((vtkProp)left);
left.GetProperty().SetDiffuseColor((double).2,(double).2,(double).2);
ren1.AddActor((vtkProp)back);
back.GetProperty().SetDiffuseColor((double).2,(double).2,(double).2);
ren1.ResetCamera();
ren1.GetActiveCamera().Dolly((double)2.5);
ren1.ResetCameraClippingRange();
renWin.Render();
// create a little scorecard above each of the cells. These are displayed[]
// if a ray cast hits the cell, otherwise they are not shown.[]
pm = new vtkPlaneSource();
pm.SetXResolution((int)1);
pm.SetYResolution((int)1);
pmapper = vtkPolyDataMapper.New();
pmapper.SetInputConnection((vtkAlgorithmOutput)pm.GetOutputPort());
// now try intersecting rays with the cell[]
cellPicker = new vtkCellPicker();
edgeCheck = new vtkActor();
edgeCheck.SetMapper((vtkMapper)pmapper);
edgeCheck.AddPosition((double)0.5,(double)2.5,(double)0);
cellPicker.Pick((double)87,(double)71,(double)0,(vtkRenderer)ren1);
if ((cellPicker.GetCellId()) != -1)
{
ren1.AddActor((vtkProp)edgeCheck);
}
triCheck = new vtkActor();
triCheck.SetMapper((vtkMapper)pmapper);
triCheck.AddPosition((double)2.5,(double)2.5,(double)0);
cellPicker.Pick((double)139,(double)72,(double)0,(vtkRenderer)ren1);
if ((cellPicker.GetCellId()) != -1)
{
ren1.AddActor((vtkProp)triCheck);
}
quadCheck = new vtkActor();
quadCheck.SetMapper((vtkMapper)pmapper);
quadCheck.AddPosition((double)4.5,(double)2.5,(double)0);
cellPicker.Pick((double)192,(double)78,(double)0,(vtkRenderer)ren1);
if ((cellPicker.GetCellId()) != -1)
{
ren1.AddActor((vtkProp)quadCheck);
}
tetCheck = new vtkActor();
tetCheck.SetMapper((vtkMapper)pmapper);
tetCheck.AddPosition((double)6.5,(double)2.5,(double)0);
cellPicker.Pick((double)233,(double)70,(double)0,(vtkRenderer)ren1);
if ((cellPicker.GetCellId()) != -1)
{
ren1.AddActor((vtkProp)tetCheck);
}
hexCheck = new vtkActor();
hexCheck.SetMapper((vtkMapper)pmapper);
hexCheck.AddPosition((double)8.5,(double)2.5,(double)0);
cellPicker.Pick((double)287,(double)80,(double)0,(vtkRenderer)ren1);
if ((cellPicker.GetCellId()) != -1)
{
ren1.AddActor((vtkProp)hexCheck);
}
wedgeCheck = new vtkActor();
wedgeCheck.SetMapper((vtkMapper)pmapper);
wedgeCheck.AddPosition((double)10.5,(double)2.5,(double)0);
cellPicker.Pick((double)287,(double)80,(double)0,(vtkRenderer)ren1);
if ((cellPicker.GetCellId()) != -1)
{
ren1.AddActor((vtkProp)wedgeCheck);
}
pyraCheck = new vtkActor();
pyraCheck.SetMapper((vtkMapper)pmapper);
pyraCheck.AddPosition((double)12.5,(double)2.5,(double)0);
cellPicker.Pick((double)287,(double)80,(double)0,(vtkRenderer)ren1);
if ((cellPicker.GetCellId()) != -1)
{
ren1.AddActor((vtkProp)pyraCheck);
}
// render the image[]
//[]
iren.Initialize();
//deleteAllVTKObjects();
}
public void WriteSimpleVTUExample()
{
string filePath = "C:\\DatenE\\02Studium\\02WiSe1718\\06IndividualProjekt\\03Daten\\testData\\simpleTest4Points.vtu";
vtkPoints points = vtkPoints.New();
points.InsertNextPoint(0, 0, 0);
points.InsertNextPoint(1, 0, 0);
points.InsertNextPoint(1, 1, 0);
points.InsertNextPoint(0, 1, 1);
points.InsertNextPoint(2, 0, 0);
//points.InsertNextPoint(2, 2, 0);
vtkTetra tetra = vtkTetra.New();
tetra.GetPointIds().SetId(0, 0);
tetra.GetPointIds().SetId(1, 1);
tetra.GetPointIds().SetId(2, 2);
tetra.GetPointIds().SetId(3, 3);
vtkCellArray cellArray = vtkCellArray.New();
cellArray.InsertNextCell(tetra);
tetra.GetPointIds().SetId(0, 1);
tetra.GetPointIds().SetId(1, 0);
tetra.GetPointIds().SetId(2, 3);
tetra.GetPointIds().SetId(3, 2);
cellArray.InsertNextCell(tetra);
tetra.GetPointIds().SetId(0, 0);
tetra.GetPointIds().SetId(1, 2);
tetra.GetPointIds().SetId(2, 3);
tetra.GetPointIds().SetId(3, 1);
cellArray.InsertNextCell(tetra);
tetra.GetPointIds().SetId(0, 3);
tetra.GetPointIds().SetId(1, 1);
tetra.GetPointIds().SetId(2, 0);
tetra.GetPointIds().SetId(3, 2);
cellArray.InsertNextCell(tetra);
tetra.GetPointIds().SetId(0, 3);
tetra.GetPointIds().SetId(1, 0);
tetra.GetPointIds().SetId(2, 2);
tetra.GetPointIds().SetId(3, 4);
cellArray.InsertNextCell(tetra);
//tetra.GetPointIds().SetId(0, 4);
//tetra.GetPointIds().SetId(1, 3);
//tetra.GetPointIds().SetId(2, 0);
//tetra.GetPointIds().SetId(3, 2);
//cellArray.InsertNextCell(tetra);
vtkUnstructuredGrid unstructuredGrid = vtkUnstructuredGrid.New();
unstructuredGrid.SetPoints(points);
const int VTK_TETRA = 10;
unstructuredGrid.SetCells(VTK_TETRA, cellArray);
// vX
vtkDoubleArray scalarsX = new vtkDoubleArray();
scalarsX.SetNumberOfValues(5);
scalarsX.SetValue(0, 4);
scalarsX.SetValue(1, 1);
scalarsX.SetValue(2, 2);
scalarsX.SetValue(3, 3);
scalarsX.SetValue(4, 1);
//scalarsX.SetValue(5, 2);
scalarsX.SetName("Vx");
unstructuredGrid.GetPointData().AddArray(scalarsX);
// vY
vtkDoubleArray scalarsY = new vtkDoubleArray();
scalarsY.SetNumberOfValues(5);
scalarsY.SetValue(0, 1);
scalarsY.SetValue(1, 2);
scalarsY.SetValue(2, 3);
scalarsY.SetValue(3, 4);
scalarsY.SetValue(4, 1);
//scalarsY.SetValue(5, 2);
scalarsY.SetName("Vy");
unstructuredGrid.GetPointData().AddArray(scalarsY);
// vZ
vtkDoubleArray scalarsZ = new vtkDoubleArray();
scalarsZ.SetNumberOfValues(5);
scalarsZ.SetValue(0, 3);
scalarsZ.SetValue(1, 1);
scalarsZ.SetValue(2, 4);
scalarsZ.SetValue(3, 2);
scalarsZ.SetValue(4, 1);
//scalarsZ.SetValue(5, 2);
scalarsZ.SetName("Vz");
unstructuredGrid.GetPointData().AddArray(scalarsZ);
// Write file
vtkXMLUnstructuredGridWriter writer = vtkXMLUnstructuredGridWriter.New();
writer.SetFileName(filePath);
writer.SetInput(unstructuredGrid);
writer.Write();
}
public void WriteVTUTriangle(string filePath, List <Vector3d> point_data, List <Vector3> trianglepoints, Dictionary <String, vtkDataArray> scalar_dataArray, List <String> arrayNames, int numberOfPoints, int numberOfTetraPoints)
{
vtkPoints points = vtkPoints.New();
for (int i = 0; i < numberOfPoints; i++)
{
points.InsertNextPoint(point_data[i].X, point_data[i].Y, point_data[i].Z);
}
vtkUnstructuredGrid unstructuredGrid = vtkUnstructuredGrid.New();
vtkCellArray cellArrayOne = vtkCellArray.New();
for (int i = 0; i < numberOfTetraPoints; i++)
{
vtkTriangle triangle = vtkTriangle.New();
triangle.GetPointIds().SetId(0, (long)trianglepoints[i][0]);
triangle.GetPointIds().SetId(1, (long)trianglepoints[i][1]);
triangle.GetPointIds().SetId(2, (long)trianglepoints[i][2]);
cellArrayOne.InsertNextCell(triangle);
}
unstructuredGrid.SetPoints(points);
const int VTK_TRIANGLE = 5;
unstructuredGrid.SetCells(VTK_TRIANGLE, cellArrayOne);
int numberOfScalarData = scalar_dataArray.Count();
for (int i = 0; i < numberOfScalarData; i++)
{
scalar_dataArray.TryGetValue(arrayNames[i], out vtkDataArray scalars);
unstructuredGrid.GetPointData().AddArray(scalars);
}
// add file ending if it is not existent
string suffix = (".vtu");
if (!(filePath.EndsWith(suffix)))
{
filePath += suffix;
}
// Write file
vtkXMLUnstructuredGridWriter writer = vtkXMLUnstructuredGridWriter.New();
writer.SetFileName(filePath);
writer.SetInput(unstructuredGrid);
writer.Write();
// Read and display file for verification that it was written correctly
vtkXMLUnstructuredGridReader reader = vtkXMLUnstructuredGridReader.New();
if (reader.CanReadFile(filePath) == 0)
{
//MessageBox.Show("Cannot read file \"" + filePath + "\"", "Error", MessageBoxButtons.OK);
return;
}
Console.WriteLine("VTU file was writen and is saved at {0}", filePath);
}
public void WriteVTM(string filePath, List <Vector3d> point_data, List <Vector4> tetrapoints, Dictionary <String, vtkDataArray> scalar_dataArray, List <String> arrayNames, int numberOfPoints, int numberOfTetraPoints)
{
vtkPoints points = vtkPoints.New();
for (int i = 0; i < numberOfPoints; i++)
{
points.InsertNextPoint(point_data[i].X, point_data[i].Y, point_data[i].Z);
}
vtkMultiBlockDataSet multiBlockDataSet = vtkMultiBlockDataSet.New();
double countBlocks = Math.Ceiling((double)((1.0 * numberOfTetraPoints) / 1000000));
multiBlockDataSet.SetNumberOfBlocks((uint)1);
for (int j = 0; j < 1; j++)
{
vtkUnstructuredGrid unstructuredGrid = vtkUnstructuredGrid.New();
vtkCellArray cellArrayOne = vtkCellArray.New();
int countTetras = 0;
int startTetras = j * 1000000;
if (numberOfTetraPoints < (j * 1000000) + 1000000)
{
countTetras = numberOfTetraPoints;
}
else
{
countTetras = (j * 1000000) + 1000000;
}
for (int i = startTetras; i < tetrapoints.Count(); i++)
{
vtkTetra tetra = vtkTetra.New();
tetra.GetPointIds().SetId(0, (long)tetrapoints[i][0]);
tetra.GetPointIds().SetId(1, (long)tetrapoints[i][1]);
tetra.GetPointIds().SetId(2, (long)tetrapoints[i][2]);
tetra.GetPointIds().SetId(3, (long)tetrapoints[i][3]);
cellArrayOne.InsertNextCell(tetra);
}
//for (int i = 0; i < tetrapoints2.Count(); i++)
//{
// vtkTetra tetra = vtkTetra.New();
// tetra.GetPointIds().SetId(0, (long)tetrapoints2[i][0]);
// tetra.GetPointIds().SetId(1, (long)tetrapoints2[i][1]);
// tetra.GetPointIds().SetId(2, (long)tetrapoints2[i][2]);
// tetra.GetPointIds().SetId(3, (long)tetrapoints2[i][3]);
// cellArrayOne.InsertNextCell(tetra);
//}
unstructuredGrid.SetPoints(points);
const int VTK_TETRA = 10;
unstructuredGrid.SetCells(VTK_TETRA, cellArrayOne);
int numberOfScalarData = scalar_dataArray.Count();
for (int i = 0; i < numberOfScalarData; i++)
{
scalar_dataArray.TryGetValue(arrayNames[i], out vtkDataArray scalars);
unstructuredGrid.GetPointData().AddArray(scalars);
}
multiBlockDataSet.SetBlock((uint)j, unstructuredGrid);
}
// add file ending if it is not existent
string suffix = (".vtm");
if (!(filePath.EndsWith(suffix)))
{
filePath += suffix;
}
// Write file
vtkXMLMultiBlockDataWriter writer = vtkXMLMultiBlockDataWriter.New();
writer.SetFileName(filePath);
writer.SetInput(multiBlockDataSet);
writer.Write();
// Read and display file for verification that it was written correctly
vtkXMLMultiBlockDataReader reader = vtkXMLMultiBlockDataReader.New();
if (reader.CanReadFile(filePath) == 0)
{
//MessageBox.Show("Cannot read file \"" + filePath + "\"", "Error", MessageBoxButtons.OK);
return;
}
Console.WriteLine("VTU file was writen and is saved at {0}", filePath);
}
public void Read_Unstructured_Grid_File(string filename)
{
// Initalize VTK Reader
vtkXMLUnstructuredGridReader reader = new vtkXMLUnstructuredGridReader();
reader.SetFileName(filename);
reader.Update();
vtkUnstructuredGrid grid = reader.GetOutput();
vtkCellArray data = grid.GetCells();
vtkIdList idList = vtkIdList.New();
int numCells = (int)data.GetNumberOfCells();
if (numCells != 0)
{
if (grid.GetCellType(0) == 10)
{
isTetra = true;
Console.WriteLine("Celltype is tetra");
}
else if (grid.GetCellType(0) == 5)
{
isTetra = false;
Console.WriteLine("Celltype is triangle");
}
else
{
Console.WriteLine("No valid celltype");
}
for (int i = 0; i < numCells; i++)
{
long cellTypeID = grid.GetCellType(i);
// alle punkte durchlaufen und in eine Variable für jeden Punkt die anderen drei Punkte speichern
if (isTetra)
{
Vector4 tetraPoint = new Vector4();
grid.GetCellPoints(i, idList);
// ueber alle vier punkte iterieren und diese in extra liste für jeden Punkt speichern
for (int j = 0; j < 4; j++)
{
tetraPoint[j] = idList.GetId(j);
}
tetraPoints.Add(tetraPoint);
}
else if (!isTetra)
{
Vector3 triangle = new Vector3();
grid.GetCellPoints(i, idList);
// ueber alle drei punkte iterieren und diese in extra liste für jeden Punkt speichern
for (int j = 0; j < 3; j++)
{
triangle[j] = idList.GetId(j);
}
trianglePoints.Add(triangle);
}
}
}
// Read Point Coordinates
vtkPoints points = grid.GetPoints();
int numPoints = (int)points.GetNumberOfPoints();
List <Vector3d> point_dat = new List <Vector3d>();
if (numPoints != 0)
{
// Read Point Data
double[] pt;
for (int i = 0; i < numPoints; i++)
{
pt = points.GetPoint(i);
point_dat.Add(new Vector3d((float)pt[0], (float)pt[1], (float)pt[2]));
}
if (this.vertex_data.ContainsKey("vertices"))
{
this.vertex_data["vertices"] = point_dat;
}
else
{
this.vertex_data.Add("vertices", point_dat);
}
Console.WriteLine("All points read in correctly!");
}
else
{
Console.ForegroundColor = ConsoleColor.Red;
Console.WriteLine("---------------No Points existent");
}
vtkPointData pointData = grid.GetPointData();
// Load point attributes
this.Load_Point_Attributes(pointData);
}
// ============================= DATA COMPRESSION ===================================
public vtkUnstructuredGrid ExportGrid(Database DB, List <string> Result, int step)
{
vtkUnstructuredGrid Output = vtkUnstructuredGrid.New();
// Update Node Coordinates
UpdateNode(DB, step);
// Create Deep Copy of Gird
Output.SetPoints(Points);
// Recreate Mesh for Output Grid
foreach (Element E in DB.ElemLib.Values)
{
// ============= HEXA =================
if (E.Type.Contains("HEX") && E.PID == ID)
{
// Create Hexahedron
vtkHexahedron Hex = vtkHexahedron.New();
Hex.GetPointIds().SetNumberOfIds(8);
// Set Hexa nodes
for (int i = 0; i < 8; i++)
{
Hex.GetPointIds().SetId(i, Nodes[E.NList[i]]);
}
// Save Element index
Output.InsertNextCell(Hex.GetCellType(), Hex.GetPointIds());
}
// ============= PENTA =================
if (E.Type.Contains("PENTA") && E.PID == ID)
{
// Create Wedge
vtkWedge Penta = vtkWedge.New();
Penta.GetPointIds().SetNumberOfIds(6);
// Set Penta nodes
for (int i = 0; i < 6; i++)
{
Penta.GetPointIds().SetId(i, Nodes[E.NList[i]]);
}
// Save Element index
Output.InsertNextCell(Penta.GetCellType(), Penta.GetPointIds());
}
// ============= TETRA =================
if (E.Type.Contains("TET") && E.PID == ID)
{
// Create Tetra
vtkTetra Tetra = vtkTetra.New();
Tetra.GetPointIds().SetNumberOfIds(4);
// Set Tetra nodes
for (int i = 0; i < 4; i++)
{
Tetra.GetPointIds().SetId(i, Nodes[E.NList[i]]);
}
// Save Element index
Output.InsertNextCell(Tetra.GetCellType(), Tetra.GetPointIds());
}
}
foreach (string res in Result)
{
foreach (vtkFloatArray array in PointScalar)
{
if (array.GetName() == res + " INC " + step.ToString())
{
vtkFloatArray ArrayCopy = vtkFloatArray.New();
ArrayCopy.DeepCopy(array);
ArrayCopy.SetName(res);
Output.GetPointData().AddArray(ArrayCopy);
}
}
}
Output.Update();
return(Output);
}
/// <summary>
/// Load scalar results to Grid Point and Cell Data
/// </summary>
public void Load_Scalar(Database DB)
{
int NumbScalar = 24; // Total number of scalar results avaliable
int NumbInc = DB.AnalysisLib.GetResultStepNo() + 1; // Total number of Analysis increments (including 0)
// Initialize Scalar data
CellScalar = new vtkFloatArray[NumbInc, NumbScalar, 3];
PointScalar = new vtkFloatArray[NumbInc, NumbScalar];
// Loop through all Analysis increments
for (int inc = 0; inc < NumbInc; inc++)
{
// Create Scalar Arrays for new increment
for (int s = 0; s < NumbScalar; s++)
{
vtkFloatArray CellMax = vtkFloatArray.New();
vtkFloatArray CellAverage = vtkFloatArray.New();
vtkFloatArray CellMin = vtkFloatArray.New();
vtkFloatArray PointArray = vtkFloatArray.New();
CellMax.SetNumberOfComponents(1);
CellAverage.SetNumberOfComponents(1);
CellMin.SetNumberOfComponents(1);
PointArray.SetNumberOfComponents(1);
CellScalar[inc, s, 0] = CellMax;
CellScalar[inc, s, 1] = CellAverage;
CellScalar[inc, s, 2] = CellMin;
PointScalar[inc, s] = PointArray;
}
// ---------- CELL (ELEMENT) SCALARS ------------------------------------------------
// Set Array Name
for (int i = 0; i < 3; i++)
{
string prefix = "";
if (i == 0)
{
prefix = "Max";
}
if (i == 1)
{
prefix = "Average";
}
if (i == 2)
{
prefix = "Min";
}
CellScalar[inc, 0, i].SetName(prefix + " Displacement X INC " + inc.ToString());
CellScalar[inc, 1, i].SetName(prefix + " Displacement Y INC " + inc.ToString());
CellScalar[inc, 2, i].SetName(prefix + " Displacement Z INC " + inc.ToString());
CellScalar[inc, 3, i].SetName(prefix + " Total Displacement INC " + inc.ToString());
CellScalar[inc, 4, i].SetName(prefix + " Stress XX INC " + inc.ToString());
CellScalar[inc, 5, i].SetName(prefix + " Stress YY INC " + inc.ToString());
CellScalar[inc, 6, i].SetName(prefix + " Stress ZZ INC " + inc.ToString());
CellScalar[inc, 7, i].SetName(prefix + " Stress XY INC " + inc.ToString());
CellScalar[inc, 8, i].SetName(prefix + " Stress YZ INC " + inc.ToString());
CellScalar[inc, 9, i].SetName(prefix + " Stress XZ INC " + inc.ToString());
CellScalar[inc, 10, i].SetName(prefix + " Stress P1 INC " + inc.ToString());
CellScalar[inc, 11, i].SetName(prefix + " Stress P2 INC " + inc.ToString());
CellScalar[inc, 12, i].SetName(prefix + " Stress P3 INC " + inc.ToString());
CellScalar[inc, 13, i].SetName(prefix + " von Mises Stress INC " + inc.ToString());
CellScalar[inc, 14, i].SetName(prefix + " Strain XX INC " + inc.ToString());
CellScalar[inc, 15, i].SetName(prefix + " Strain YY INC " + inc.ToString());
CellScalar[inc, 16, i].SetName(prefix + " Strain ZZ INC " + inc.ToString());
CellScalar[inc, 17, i].SetName(prefix + " Strain XY INC " + inc.ToString());
CellScalar[inc, 18, i].SetName(prefix + " Strain YZ INC " + inc.ToString());
CellScalar[inc, 19, i].SetName(prefix + " Strain XZ INC " + inc.ToString());
CellScalar[inc, 20, i].SetName(prefix + " Strain P1 INC " + inc.ToString());
CellScalar[inc, 21, i].SetName(prefix + " Strain P2 INC " + inc.ToString());
CellScalar[inc, 22, i].SetName(prefix + " Strain P3 INC " + inc.ToString());
CellScalar[inc, 23, i].SetName(prefix + " Effective Strain INC " + inc.ToString());
}
// For each Element collect Nodal values, take Min/Average/Max and insert to array
Parallel.ForEach(Elements, ID =>
{
// Take element from database
int index = Elements.IndexOf(ID);
Element E = DB.ElemLib[Elements[index]];
// List of results in Nodes
List <double[]> Nval = new List <double[]>();
for (int i = 0; i < NumbScalar; i++)
{
Nval.Add(new double[E.NList.Count]);
}
// Assign Nodal values
for (int i = 0; i < E.NList.Count; i++)
{
// Displacement
Nval[0][i] = DB.NodeLib[E.NList[i]].GetDisp(inc, 0);
Nval[1][i] = DB.NodeLib[E.NList[i]].GetDisp(inc, 1);
Nval[2][i] = DB.NodeLib[E.NList[i]].GetDisp(inc, 2);
Nval[3][i] = Math.Sqrt(Math.Pow(Nval[0][i], 2) + Math.Pow(Nval[1][i], 2) + Math.Pow(Nval[2][i], 2));
// Stress
Matrix <double> S = Matrix <double> .Build.Dense(3, 3);
S[0, 0] = E.GetStress(inc, i, 0);
S[1, 1] = E.GetStress(inc, i, 1);
S[2, 2] = E.GetStress(inc, i, 2);
S[0, 1] = E.GetStress(inc, i, 3);
S[1, 0] = E.GetStress(inc, i, 3);
S[1, 2] = E.GetStress(inc, i, 4);
S[2, 1] = E.GetStress(inc, i, 4);
S[0, 2] = E.GetStress(inc, i, 5);
S[2, 0] = E.GetStress(inc, i, 5);
Evd <double> eigen = S.Evd();
double P1 = eigen.EigenValues[2].Real;
double P2 = eigen.EigenValues[1].Real;
double P3 = eigen.EigenValues[0].Real;
Nval[4][i] = S[0, 0];
Nval[5][i] = S[1, 1];
Nval[6][i] = S[2, 2];
Nval[7][i] = S[0, 1];
Nval[8][i] = S[1, 2];
Nval[9][i] = S[0, 2];
Nval[10][i] = P1;
Nval[11][i] = P2;
Nval[12][i] = P3;
Nval[13][i] = Math.Sqrt((Math.Pow(P1 - P2, 2) + Math.Pow(P2 - P3, 2) + Math.Pow(P3 - P1, 2)) / 2);
// Strain
S = Matrix <double> .Build.Dense(3, 3);
S[0, 0] = E.GetStrain(inc, i, 0);
S[1, 1] = E.GetStrain(inc, i, 1);
S[2, 2] = E.GetStrain(inc, i, 2);
S[0, 1] = E.GetStrain(inc, i, 3);
S[1, 0] = E.GetStrain(inc, i, 3);
S[1, 2] = E.GetStrain(inc, i, 4);
S[2, 1] = E.GetStrain(inc, i, 4);
S[0, 2] = E.GetStrain(inc, i, 5);
S[2, 0] = E.GetStrain(inc, i, 5);
eigen = S.Evd();
P1 = eigen.EigenValues[2].Real;
P2 = eigen.EigenValues[1].Real;
P3 = eigen.EigenValues[0].Real;
Nval[14][i] = S[0, 0];
Nval[15][i] = S[1, 1];
Nval[16][i] = S[2, 2];
Nval[17][i] = S[0, 1];
Nval[18][i] = S[1, 2];
Nval[19][i] = S[0, 2];
Nval[20][i] = P1;
Nval[21][i] = P2;
Nval[22][i] = P3;
Nval[23][i] = (2.0 / 3.0) * Math.Sqrt((Math.Pow(P1 - P2, 2) + Math.Pow(P2 - P3, 2) + Math.Pow(P3 - P1, 2)) / 2);
}
lock (CellScalar)
{
for (int s = 0; s < NumbScalar; s++)
{
CellScalar[inc, s, 0].InsertTuple1(index, Nval[s].Max());
CellScalar[inc, s, 1].InsertTuple1(index, Nval[s].Average());
CellScalar[inc, s, 2].InsertTuple1(index, Nval[s].Min());
}
}
});
// Add arrays to Grid CellData
foreach (vtkFloatArray array in CellScalar)
{
Grid.GetCellData().AddArray(array);
}
// ---------- POINT (NODE) SCALARS ------------------------------------------------
// Set Array Name
PointScalar[inc, 0].SetName("Displacement X INC " + inc.ToString());
PointScalar[inc, 1].SetName("Displacement Y INC " + inc.ToString());
PointScalar[inc, 2].SetName("Displacement Z INC " + inc.ToString());
PointScalar[inc, 3].SetName("Total Displacement INC " + inc.ToString());
PointScalar[inc, 4].SetName("Stress XX INC " + inc.ToString());
PointScalar[inc, 5].SetName("Stress YY INC " + inc.ToString());
PointScalar[inc, 6].SetName("Stress ZZ INC " + inc.ToString());
PointScalar[inc, 7].SetName("Stress XY INC " + inc.ToString());
PointScalar[inc, 8].SetName("Stress YZ INC " + inc.ToString());
PointScalar[inc, 9].SetName("Stress XZ INC " + inc.ToString());
PointScalar[inc, 10].SetName("Stress P1 INC " + inc.ToString());
PointScalar[inc, 11].SetName("Stress P2 INC " + inc.ToString());
PointScalar[inc, 12].SetName("Stress P3 INC " + inc.ToString());
PointScalar[inc, 13].SetName("von Mises Stress INC " + inc.ToString());
PointScalar[inc, 14].SetName("Strain XX INC " + inc.ToString());
PointScalar[inc, 15].SetName("Strain YY INC " + inc.ToString());
PointScalar[inc, 16].SetName("Strain ZZ INC " + inc.ToString());
PointScalar[inc, 17].SetName("Strain XY INC " + inc.ToString());
PointScalar[inc, 18].SetName("Strain YZ INC " + inc.ToString());
PointScalar[inc, 19].SetName("Strain XZ INC " + inc.ToString());
PointScalar[inc, 20].SetName("Strain P1 INC " + inc.ToString());
PointScalar[inc, 21].SetName("Strain P2 INC " + inc.ToString());
PointScalar[inc, 22].SetName("Strain P3 INC " + inc.ToString());
PointScalar[inc, 23].SetName("Effective Strain INC " + inc.ToString());
// For each Node collect Nodal values, take Min/Average/Max and insert to array
Parallel.ForEach(Nodes, n =>
{
// Take element from database
Node N = DB.NodeLib[n.Key];
// List of results in Nodes
List <double[]> Nval = new List <double[]>();
for (int i = 0; i < NumbScalar; i++)
{
Nval.Add(new double[N.EList.Count()]);
}
// Assign Nodal values
for (int i = 0; i < N.EList.Count; i++)
{
// Index of this node in i-th Element Node List
int NodeIndex = DB.ElemLib[N.EList[i]].NList.IndexOf(N.ID);
// Displacement
Nval[0][i] = N.GetDisp(inc, 0);
Nval[1][i] = N.GetDisp(inc, 1);
Nval[2][i] = N.GetDisp(inc, 2);
Nval[3][i] = Math.Sqrt(Math.Pow(Nval[0][i], 2) + Math.Pow(Nval[1][i], 2) + Math.Pow(Nval[2][i], 2));
// Stress
Matrix <double> S = Matrix <double> .Build.Dense(3, 3);
S[0, 0] = DB.ElemLib[N.EList[i]].GetStress(inc, NodeIndex, 0);
S[1, 1] = DB.ElemLib[N.EList[i]].GetStress(inc, NodeIndex, 1);
S[2, 2] = DB.ElemLib[N.EList[i]].GetStress(inc, NodeIndex, 2);
S[0, 1] = DB.ElemLib[N.EList[i]].GetStress(inc, NodeIndex, 3);
S[1, 0] = DB.ElemLib[N.EList[i]].GetStress(inc, NodeIndex, 3);
S[1, 2] = DB.ElemLib[N.EList[i]].GetStress(inc, NodeIndex, 4);
S[2, 1] = DB.ElemLib[N.EList[i]].GetStress(inc, NodeIndex, 4);
S[0, 2] = DB.ElemLib[N.EList[i]].GetStress(inc, NodeIndex, 5);
S[2, 0] = DB.ElemLib[N.EList[i]].GetStress(inc, NodeIndex, 5);
Evd <double> eigen = S.Evd();
double P1 = eigen.EigenValues[2].Real;
double P2 = eigen.EigenValues[1].Real;
double P3 = eigen.EigenValues[0].Real;
Nval[4][i] = S[0, 0];
Nval[5][i] = S[1, 1];
Nval[6][i] = S[2, 2];
Nval[7][i] = S[0, 1];
Nval[8][i] = S[1, 2];
Nval[9][i] = S[0, 2];
Nval[10][i] = P1;
Nval[11][i] = P2;
Nval[12][i] = P3;
Nval[13][i] = Math.Sqrt((Math.Pow(P1 - P2, 2) + Math.Pow(P2 - P3, 2) + Math.Pow(P3 - P1, 2)) / 2);
// Strain
S = Matrix <double> .Build.Dense(3, 3);
S[0, 0] = DB.ElemLib[N.EList[i]].GetStrain(inc, NodeIndex, 0);
S[1, 1] = DB.ElemLib[N.EList[i]].GetStrain(inc, NodeIndex, 1);
S[2, 2] = DB.ElemLib[N.EList[i]].GetStrain(inc, NodeIndex, 2);
S[0, 1] = DB.ElemLib[N.EList[i]].GetStrain(inc, NodeIndex, 3);
S[1, 0] = DB.ElemLib[N.EList[i]].GetStrain(inc, NodeIndex, 3);
S[1, 2] = DB.ElemLib[N.EList[i]].GetStrain(inc, NodeIndex, 4);
S[2, 1] = DB.ElemLib[N.EList[i]].GetStrain(inc, NodeIndex, 4);
S[0, 2] = DB.ElemLib[N.EList[i]].GetStrain(inc, NodeIndex, 5);
S[2, 0] = DB.ElemLib[N.EList[i]].GetStrain(inc, NodeIndex, 5);
eigen = S.Evd();
P1 = eigen.EigenValues[2].Real;
P2 = eigen.EigenValues[1].Real;
P3 = eigen.EigenValues[0].Real;
Nval[14][i] = S[0, 0];
Nval[15][i] = S[1, 1];
Nval[16][i] = S[2, 2];
Nval[17][i] = S[0, 1];
Nval[18][i] = S[1, 2];
Nval[19][i] = S[0, 2];
Nval[20][i] = P1;
Nval[21][i] = P2;
Nval[22][i] = P3;
Nval[23][i] = (2.0 / 3.0) * Math.Sqrt((Math.Pow(P1 - P2, 2) + Math.Pow(P2 - P3, 2) + Math.Pow(P3 - P1, 2)) / 2);
}
lock (PointScalar)
{
for (int s = 0; s < NumbScalar; s++)
{
PointScalar[inc, s].InsertTuple1(n.Value, Nval[s].Average());
}
}
});
// Add arrays to Grid PointData
foreach (vtkFloatArray array in PointScalar)
{
Grid.GetPointData().AddArray(array);
}
}
}