bool
HdSt_TriangleIndexBuilderComputation::Resolve()
{
if (!_TryLock()) return false;
HD_TRACE_FUNCTION();
VtVec3iArray trianglesFaceVertexIndices;
VtIntArray primitiveParam;
HdMeshUtil meshUtil(_topology, _id);
meshUtil.ComputeTriangleIndices(
&trianglesFaceVertexIndices,
&primitiveParam);
_SetResult(HdBufferSourceSharedPtr(
new HdVtBufferSource(
HdTokens->indices,
VtValue(trianglesFaceVertexIndices))));
_primitiveParam.reset(new HdVtBufferSource(
HdTokens->primitiveParam,
VtValue(primitiveParam)));
_SetResolved();
return true;
}
bool
HdSt_TriangulateFaceVaryingComputation::Resolve()
{
if (!TF_VERIFY(_source)) return false;
if (!_source->IsResolved()) return false;
if (!_TryLock()) return false;
HD_TRACE_FUNCTION();
HD_PERF_COUNTER_INCR(HdPerfTokens->triangulateFaceVarying);
VtValue result;
HdMeshUtil meshUtil(_topology, _id);
if(meshUtil.ComputeTriangulatedFaceVaryingPrimvar(
_source->GetData(),
_source->GetNumElements(),
_source->GetGLElementDataType(),
&result)) {
_SetResult(HdBufferSourceSharedPtr(
new HdVtBufferSource(
_source->GetName(),
result)));
} else {
_SetResult(_source);
}
_SetResolved();
return true;
}
void
STKSearcher<SpatialDim>::setupSearch()
{
mesh::fem::FEMMetaData& metaData = stk_classic::mesh::fem::FEMMetaData::get(*m_bulk);
mesh::BulkData& bulkData = *m_bulk;
VectorFieldType *coords_field = metaData.get_field<VectorFieldType >("coordinates");
PerceptMesh meshUtil(&metaData, &bulkData);
BBB buildBoundingBoxes(m_boxes, coords_field);
meshUtil.elementOpLoop(buildBoundingBoxes, coords_field);
}
const stk_classic::mesh::Entity *
STKSearcher<SpatialDim>::findElement(MDArray& input_phy_points, MDArray& found_parametric_coordinates,
unsigned& found_it, const mesh::Entity *hint_element )
{
//return 0;
mesh::fem::FEMMetaData& metaData = stk_classic::mesh::fem::FEMMetaData::get(*m_bulk);
mesh::BulkData& bulkData = *m_bulk;
//VectorFieldType *coords_field = metaData.get_field<VectorFieldType >("coordinates");
PerceptMesh meshUtil(&metaData, &bulkData);
double pts[SpatialDim];
for (unsigned iDim = 0; iDim < SpatialDim; iDim++)
{
pts[iDim] = input_phy_points(0, iDim);
}
BPoint pointBoundingBox;
pointBoundingBox.key.ident = 123; // FIXME for multiple points
pointBoundingBox.set_center(pts);
std::vector<BPoint> points(1, pointBoundingBox);
stk_classic::search::FactoryOrder order;
order.m_communicator = bulkData.parallel();
order.m_algorithm = stk_classic::search::FactoryOrder::BIHTREE;
if (0 || EXTRA_PRINT)
{
bool box_p = m_boxes[0].intersect(pointBoundingBox);
bool p_box = pointBoundingBox.intersect(m_boxes[0]);
std::cout << "STKSearcher::findElement: m_boxes[0]= " << m_boxes[0] << std::endl;
std::cout << "STKSearcher::findElement: pointBoundingBox= " << pointBoundingBox << std::endl;
std::cout << "STKSearcher::findElement: box_p= " << box_p << std::endl;
std::cout << "STKSearcher::findElement: p_box= " << p_box << std::endl;
}
if (0 || EXTRA_PRINT) std::cout << "STKSearcher::findElement: nboxes= " << m_boxes.size() << std::endl;
IdentProcRelation relation;
stk_classic::search::coarse_search(relation, m_boxes, points, order);
//stk_classic::search::coarse_search(relation, points, m_boxes, order);
if (0 || EXTRA_PRINT) std::cout << "STKSearcher::findElement: found " << relation.size() << " containing bboxes" << std::endl;
if (relation.size())
{
IsInElement isIn(input_phy_points, found_parametric_coordinates);
for (unsigned i = 0; i < relation.size(); i++)
{
if (0 || EXTRA_PRINT)
std::cout << "relation[ " << i << "]= {" << relation[i].first << "} --> { " << relation[i].second << "}" << std::endl;
mesh::Entity *element = bulkData.get_entity(metaData.element_rank(), relation[i].second.ident);
//bool loop_break = ... intentionally ignoring return value
isIn(*element, bulkData);
if (0 || EXTRA_PRINT) std::cout << "STKSearcher::findElement: found it= " << isIn.m_found_it << std::endl;
if (isIn.m_found_it)
{
found_it = 1;
return isIn.m_foundElement;
}
else
{
found_it = 0;
return 0;
}
}
}
return 0;
}