void tarch::plotter::griddata::unstructured::vtk::VTKTextFileWriter::VertexDataWriter::plotVertex( int index, const tarch::la::Vector<3,double>& value ) {
assertion(_lastWriteCommandVertexNumber>=-1);
assertion(3<=_recordsPerVertex);
assertion1( static_cast<float>(value(0)) != std::numeric_limits<float>::infinity(), value(0) );
assertion1( value(0) == value(0), value(0) ); // test for not a number
assertion1( static_cast<float>(value(1)) != std::numeric_limits<float>::infinity(), value(1) );
assertion1( value(1) == value(1), value(1) ); // test for not a number
assertion1( static_cast<float>(value(2)) != std::numeric_limits<float>::infinity(), value(2) );
assertion1( value(2) == value(2), value(2) ); // test for not a number
while (_lastWriteCommandVertexNumber<index-1) {
plotVertex(_lastWriteCommandVertexNumber+1,0.0);
}
_lastWriteCommandVertexNumber = index;
_out << static_cast<float>(value(0)) << " ";
_out << static_cast<float>(value(1)) << " ";
_out << static_cast<float>(value(2)) << " ";
for (int i=3; i<_recordsPerVertex; i++) {
_out << 0.0 << " ";
}
_out << std::endl;
if (value(0)<_minValue) _minValue = value(0);
if (value(0)>_maxValue) _maxValue = value(0);
if (value(1)<_minValue) _minValue = value(1);
if (value(1)>_maxValue) _maxValue = value(1);
if (value(2)<_minValue) _minValue = value(2);
if (value(2)>_maxValue) _maxValue = value(2);
}
void peano::applications::pic::demo2::mappings::RegularGrid2PlotSolution::touchVertexFirstTime(
peano::applications::pic::demo2::RegularGridVertex& vertex,
const tarch::la::Vector<DIMENSIONS,double>& x
) {
logTraceInWith2Arguments( "touchVertexFirstTime()", x, vertex );
// @todo Insert your code here
plotVertex(vertex,x);
logTraceOutWith1Argument( "touchVertexFirstTime()", vertex );
}
void peano::applications::pic::demo2::mappings::RegularGrid2PlotSolution::createBoundaryVertex(
const tarch::la::Vector<DIMENSIONS,double>& x,
const tarch::la::Vector<DIMENSIONS,double>& h,
peano::applications::pic::demo2::RegularGridVertex& vertex
) {
logTraceInWith3Arguments( "createBoundaryVertex()", x, h, vertex );
// @todo Insert your code here
plotVertex( vertex, x );
logTraceOutWith1Argument( "createBoundaryVertex()", vertex );
}
int tarch::plotter::griddata::unstructured::binaryvtu::BinaryVTUTextFileWriter::VertexWriter::plotVertex(const tarch::la::Vector<2,double>& position) {
assertion1( _currentVertexNumber>=0, _currentVertexNumber );
tarch::la::Vector<3,double> p;
p(0) = position(0);
p(1) = position(1);
p(2) = 0.0;
return plotVertex(p);
}
void particles::pit::adapters::InitGridAndPlot2VTKGridVisualiser_1::createHangingVertex(
particles::pit::Vertex& fineGridVertex,
const tarch::la::Vector<DIMENSIONS,double>& fineGridX,
const tarch::la::Vector<DIMENSIONS,double>& fineGridH,
particles::pit::Vertex * const coarseGridVertices,
const peano::grid::VertexEnumerator& coarseGridVerticesEnumerator,
particles::pit::Cell& coarseGridCell,
const tarch::la::Vector<DIMENSIONS,int>& fineGridPositionOfVertex
) {
plotVertex( fineGridVertex, fineGridX );
}
void particles::pit::adapters::SetupExperimentWithRandomParticleDistributionAndPlot2VTKGridVisualiser_0::createHangingVertex(
particles::pit::Vertex& fineGridVertex,
const tarch::la::Vector<DIMENSIONS,double>& fineGridX,
const tarch::la::Vector<DIMENSIONS,double>& fineGridH,
particles::pit::Vertex * const coarseGridVertices,
const peano::grid::VertexEnumerator& coarseGridVerticesEnumerator,
particles::pit::Cell& coarseGridCell,
const tarch::la::Vector<DIMENSIONS,int>& fineGridPositionOfVertex
) {
plotVertex( fineGridVertex, fineGridX );
}
void dem::adapters::TimeStepAndPlotOnReluctantDynamicGrid2VTKGridVisualiser_4::createHangingVertex(
dem::Vertex& fineGridVertex,
const tarch::la::Vector<DIMENSIONS,double>& fineGridX,
const tarch::la::Vector<DIMENSIONS,double>& fineGridH,
dem::Vertex * const coarseGridVertices,
const peano::grid::VertexEnumerator& coarseGridVerticesEnumerator,
dem::Cell& coarseGridCell,
const tarch::la::Vector<DIMENSIONS,int>& fineGridPositionOfVertex
) {
plotVertex( fineGridVertex, fineGridX );
}
void peano::applications::pic::demo2::mappings::SpacetreeGrid2PlotSolution::createBoundaryVertex(
peano::applications::pic::demo2::SpacetreeGridVertex& fineGridVertex,
const tarch::la::Vector<DIMENSIONS,double>& fineGridX,
const tarch::la::Vector<DIMENSIONS,double>& fineGridH,
peano::applications::pic::demo2::SpacetreeGridVertex const * const coarseGridVertices,
const peano::kernel::gridinterface::VertexEnumerator& coarseGridVerticesEnumerator,
const peano::applications::pic::demo2::SpacetreeGridCell& coarseGridCell,
const tarch::la::Vector<DIMENSIONS,int>& fineGridPositionOfVertex
) {
logTraceInWith6Arguments( "createBoundaryVertex(...)", fineGridVertex, fineGridX, fineGridH, coarseGridVerticesEnumerator.toString(), coarseGridCell, fineGridPositionOfVertex );
// @todo Insert your code here
plotVertex( fineGridVertex, fineGridX );
logTraceOutWith1Argument( "createBoundaryVertex(...)", fineGridVertex );
}
void peano::applications::poisson::multigrid::mappings::SpacetreeGrid2PlotSolution::createHangingVertex(
peano::applications::poisson::multigrid::SpacetreeGridVertex& fineGridVertex,
const tarch::la::Vector<DIMENSIONS,double>& fineGridX,
const tarch::la::Vector<DIMENSIONS,double>& fineGridH,
peano::applications::poisson::multigrid::SpacetreeGridVertex const * const coarseGridVertices,
const peano::kernel::gridinterface::VertexEnumerator& coarseGridVerticesEnumerator,
const peano::applications::poisson::multigrid::SpacetreeGridCell& coarseGridCell,
const tarch::la::Vector<DIMENSIONS,int>& fineGridPositionOfVertex
) {
logTraceInWith6Arguments( "createHangingVertex(...)", fineGridVertex, fineGridX, fineGridH, coarseGridVerticesEnumerator.toString(), coarseGridCell, fineGridPositionOfVertex );
plotVertex(fineGridVertex, fineGridX);
logTraceOutWith1Argument( "createHangingVertex(...)", fineGridVertex );
}
void peano::applications::faxen::adapters::SpacetreeGrid2PlotGrid::createHangingVertex(
peano::applications::faxen::SpacetreeGridVertex& fineGridVertex,
const tarch::la::Vector<DIMENSIONS,double>& fineGridX,
const tarch::la::Vector<DIMENSIONS,double>& fineGridH,
peano::applications::faxen::SpacetreeGridVertex const * const coarseGridVertices,
const peano::kernel::gridinterface::VertexEnumerator& coarseGridVerticesEnumerator,
const peano::applications::faxen::SpacetreeGridCell& coarseGridCell,
const tarch::la::Vector<DIMENSIONS,int>& fineGridPositionOfVertex
) {
plotVertex( fineGridVertex, fineGridX );
}
void peano::applications::heatequation::timestepping::adapters::SpacetreeGrid2PlotGrid::touchVertexFirstTime(
peano::applications::heatequation::timestepping::SpacetreeGridSingleStepVertex& fineGridVertex,
const tarch::la::Vector<DIMENSIONS,double>& fineGridX,
const tarch::la::Vector<DIMENSIONS,double>& fineGridH,
peano::applications::heatequation::timestepping::SpacetreeGridSingleStepVertex const * const coarseGridVertices,
const peano::kernel::gridinterface::VertexEnumerator& coarseGridVerticesEnumerator,
const peano::applications::heatequation::timestepping::SpacetreeGridCell& coarseGridCell,
const tarch::la::Vector<DIMENSIONS,int>& fineGridPositionOfVertex
) {
logTraceInWith6Arguments( "touchVertexFirstTime(...)", fineGridVertex, fineGridX, fineGridH, coarseGridVerticesEnumerator.toString(), coarseGridCell, fineGridPositionOfVertex );
plotVertex( fineGridVertex, fineGridX );
logTraceOutWith1Argument( "touchVertexFirstTime(...)", _vertex2IndexMap.size() );
}
void dem::adapters::TimeStepAndPlotOnReluctantDynamicGrid2VTKGridVisualiser_4::touchVertexFirstTime(
dem::Vertex& fineGridVertex,
const tarch::la::Vector<DIMENSIONS,double>& fineGridX,
const tarch::la::Vector<DIMENSIONS,double>& fineGridH,
dem::Vertex * const coarseGridVertices,
const peano::grid::VertexEnumerator& coarseGridVerticesEnumerator,
dem::Cell& coarseGridCell,
const tarch::la::Vector<DIMENSIONS,int>& fineGridPositionOfVertex
) {
if (
fineGridVertex.getRefinementControl()==dem::Vertex::Records::Unrefined ||
fineGridVertex.getRefinementControl()==dem::Vertex::Records::RefinementTriggered ||
fineGridVertex.getRefinementControl()==dem::Vertex::Records::Erasing
) {
plotVertex( fineGridVertex, fineGridX );
}
}
void particles::pit::adapters::InitGridAndPlot2VTKGridVisualiser_1::touchVertexFirstTime(
particles::pit::Vertex& fineGridVertex,
const tarch::la::Vector<DIMENSIONS,double>& fineGridX,
const tarch::la::Vector<DIMENSIONS,double>& fineGridH,
particles::pit::Vertex * const coarseGridVertices,
const peano::grid::VertexEnumerator& coarseGridVerticesEnumerator,
particles::pit::Cell& coarseGridCell,
const tarch::la::Vector<DIMENSIONS,int>& fineGridPositionOfVertex
) {
if (
fineGridVertex.getRefinementControl()==particles::pit::Vertex::Records::Unrefined ||
fineGridVertex.getRefinementControl()==particles::pit::Vertex::Records::RefinementTriggered ||
fineGridVertex.getRefinementControl()==particles::pit::Vertex::Records::Erasing ||
fineGridVertex.getRefinementControl()==particles::pit::Vertex::Records::RefineDueToJoinThoughWorkerIsAlreadyErasing
) {
plotVertex( fineGridVertex, fineGridX );
}
}
void peano::applications::heatequation::timestepping::adapters::SpacetreeGrid2PlotGrid::createHangingVertex(
peano::applications::heatequation::timestepping::SpacetreeGridSingleStepVertex& fineGridVertex,
const tarch::la::Vector<DIMENSIONS,double>& fineGridX,
const tarch::la::Vector<DIMENSIONS,double>& fineGridH,
peano::applications::heatequation::timestepping::SpacetreeGridSingleStepVertex const * const coarseGridVertices,
const peano::kernel::gridinterface::VertexEnumerator& coarseGridVerticesEnumerator,
const peano::applications::heatequation::timestepping::SpacetreeGridCell& coarseGridCell,
const tarch::la::Vector<DIMENSIONS,int>& fineGridPositionOfVertex
) {
#ifdef SharedTBB
Vertex2IndexMapSemaphore::scoped_lock localLock(_vertex2IndexMapSemaphore);
#elif SharedOMP
#pragma omp critical
#endif
plotVertex( fineGridVertex, fineGridX );
}
void tarch::plotter::griddata::unstructured::vtk::VTKTextFileWriter::VertexDataWriter::plotVertex( int index, double value ) {
assertion(_lastWriteCommandVertexNumber>=-1);
assertion(1<=_recordsPerVertex);
float valueAsFloat = static_cast<float>(value);
assertion3( valueAsFloat != std::numeric_limits<float>::infinity(), index, value, valueAsFloat );
assertion3( valueAsFloat == valueAsFloat, index, value, valueAsFloat ); // test for not a number
while (_lastWriteCommandVertexNumber<index-1) {
plotVertex(_lastWriteCommandVertexNumber+1,0.0);
}
_lastWriteCommandVertexNumber = index;
_out << valueAsFloat << " ";
for (int i=1; i<_recordsPerVertex; i++) {
_out << 0.0 << " ";
}
_out << std::endl;
if (value<_minValue) _minValue = value;
if (value>_maxValue) _maxValue = value;
}
