// -----------------------------------------------------------------------------
//
// -----------------------------------------------------------------------------
int TriangleGeom::writeXdmf(QTextStream& out, QString dcName, QString hdfFileName)
{
herr_t err = 0;
// Always start the grid
out << " <!-- *************** START OF " << dcName << " *************** -->" << "\n";
out << " <Grid Name=\"" << dcName << "\" GridType=\"Uniform\">" << "\n";
out << " <Topology TopologyType=\"Triangle\" NumberOfElements=\"" << getNumberOfTris() << "\">" << "\n";
out << " <DataItem Format=\"HDF\" NumberType=\"Int\" Dimensions=\"" << getNumberOfTris() << " 3\">" << "\n";
out << " " << hdfFileName << ":/DataContainers/" << dcName << "/" << DREAM3D::Geometry::Geometry << "/" << DREAM3D::Geometry::SharedTriList << "\n";
out << " </DataItem>" << "\n";
out << " </Topology>" << "\n";
if (m_VertexList.get() == NULL)
{
out << "<!-- ********************* GEOMETRY ERROR ****************************************\n";
out << "The Geometry with name '" << getName() << "' in DataContainer '" << dcName << "' \n";
out << "did not have any vertices assigned.\n";
out << "The Geometry types will be missing from the Xdmf which will cause issues when\n";
out << "trying to load the file\n";
out << " ********************************************************************************** -->\n";
}
else
{
out << " <Geometry Type=\"XYZ\">" << "\n";
out << " <DataItem Format=\"HDF\" Dimensions=\"" << getNumberOfVertices() << " 3\" NumberType=\"Float\" Precision=\"4\">" << "\n";
out << " " << hdfFileName << ":/DataContainers/" << dcName << "/" << DREAM3D::Geometry::Geometry << "/" << DREAM3D::Geometry::SharedVertexList << "\n";
out << " </DataItem>" << "\n";
out << " </Geometry>" << "\n";
out << "" << "\n";
}
return err;
}
// -----------------------------------------------------------------------------
//
// -----------------------------------------------------------------------------
void TriangleGeom::findDerivatives(DoubleArrayType::Pointer field, DoubleArrayType::Pointer derivatives)
{
int64_t numTris = getNumberOfTris();
int cDims = field->getNumberOfComponents();
double* fieldPtr = field->getPointer(0);
double* derivsPtr = derivatives->getPointer(0);
double values[3];
double derivs[3];
int64_t verts[3];
for (int64_t i = 0; i < numTris; i++)
{
getVertsAtTri(i, verts);
for (int j = 0; j < cDims; j++)
{
for (size_t k = 0; k < 3; k++)
{
values[k] = fieldPtr[cDims * verts[k] + j];
}
DerivativeHelpers::TriangleDeriv()(this, i, values, derivs);
derivsPtr[i * 3 * cDims + j * 3] = derivs[0];
derivsPtr[i * 3 * cDims + j * 3 + 1] = derivs[1];
derivsPtr[i * 3 * cDims + j * 3 + 2] = derivs[2];
}
}
}
// -----------------------------------------------------------------------------
//
// -----------------------------------------------------------------------------
void TriangleGeom::findDerivatives(DoubleArrayType::Pointer field, DoubleArrayType::Pointer derivatives, Observable* observable)
{
m_ProgressCounter = 0;
int64_t numTris = getNumberOfTris();
if (observable)
{
connect(this, SIGNAL(filterGeneratedMessage(const PipelineMessage&)),
observable, SLOT(broadcastPipelineMessage(const PipelineMessage&)));
}
#ifdef SIMPLib_USE_PARALLEL_ALGORITHMS
tbb::task_scheduler_init init;
bool doParallel = true;
#endif
#ifdef SIMPLib_USE_PARALLEL_ALGORITHMS
if (doParallel == true)
{
tbb::parallel_for(tbb::blocked_range<int64_t>(0, numTris),
FindTriangleDerivativesImpl(this, field, derivatives), tbb::auto_partitioner());
}
else
#endif
{
FindTriangleDerivativesImpl serial(this, field, derivatives);
serial.compute(0, numTris);
}
}
// -----------------------------------------------------------------------------
//
// -----------------------------------------------------------------------------
int TriangleGeom::findElementCentroids()
{
QVector<size_t> cDims(1, 3);
m_TriangleCentroids = FloatArrayType::CreateArray(getNumberOfTris(), cDims, DREAM3D::StringConstants::TriangleCentroids);
GeometryHelpers::Topology::FindElementCentroids<int64_t>(m_TriList, m_VertexList, m_TriangleCentroids);
if (m_TriangleCentroids.get() == NULL)
{
return -1;
}
return 1;
}
// -----------------------------------------------------------------------------
//
// -----------------------------------------------------------------------------
QString TriangleGeom::getInfoString(DREAM3D::InfoStringFormat format)
{
QString info;
QTextStream ss (&info);
if(format == DREAM3D::HtmlFormat)
{
ss << "<tr bgcolor=\"#D3D8E0\"><th colspan=2>Triangle Geometry Info</th></tr>";
ss << "<tr bgcolor=\"#C3C8D0\"><th align=\"right\">Number of Triangles</th><td>" << getNumberOfTris() << "</td></tr>";
ss << "<tr bgcolor=\"#C3C8D0\"><th align=\"right\">Number of Vertices</th><td>" << getNumberOfVertices() << "</td></tr>";
}
else
{
}
return info;
}
// -----------------------------------------------------------------------------
//
// -----------------------------------------------------------------------------
int TriangleGeom::writeGeometryToHDF5(hid_t parentId, bool DREAM3D_NOT_USED(writeXdmf))
{
herr_t err = 0;
if (m_VertexList.get() != NULL)
{
err = GeometryHelpers::GeomIO::WriteListToHDF5(parentId, m_VertexList);
if (err < 0)
{
return err;
}
}
if (m_EdgeList.get() != NULL)
{
err = GeometryHelpers::GeomIO::WriteListToHDF5(parentId, m_EdgeList);
if (err < 0)
{
return err;
}
}
if (m_TriList.get() != NULL)
{
err = GeometryHelpers::GeomIO::WriteListToHDF5(parentId, m_TriList);
if (err < 0)
{
return err;
}
}
if (m_UnsharedEdgeList.get() != NULL)
{
err = GeometryHelpers::GeomIO::WriteListToHDF5(parentId, m_UnsharedEdgeList);
if (err < 0)
{
return err;
}
}
if (m_TriangleCentroids.get() != NULL)
{
err = GeometryHelpers::GeomIO::WriteListToHDF5(parentId, m_TriangleCentroids);
if (err < 0)
{
return err;
}
}
if (m_TriangleNeighbors.get() != NULL)
{
size_t numTris = getNumberOfTris();
err = GeometryHelpers::GeomIO::WriteDynamicListToHDF5<uint16_t, int64_t>(parentId, m_TriangleNeighbors, numTris, DREAM3D::StringConstants::TriangleNeighbors);
if (err < 0)
{
return err;
}
}
if (m_TrianglesContainingVert.get() != NULL)
{
size_t numVerts = getNumberOfVertices();
err = GeometryHelpers::GeomIO::WriteDynamicListToHDF5<uint16_t, int64_t>(parentId, m_TrianglesContainingVert, numVerts, DREAM3D::StringConstants::TrianglesContainingVert);
if (err < 0)
{
return err;
}
}
return err;
}
