@@ -287,45 +287,6 @@ class InsertionAlgorithm : public BaseAlgorithm
287287 }
288288 }
289289 }
290-
291- // ElementIterator::SPtr itTip = l_tipGeom->begin();
292- // auto createTipProximity =
293- // Operations::CreateCenterProximity::Operation::get(itTip->getTypeInfo());
294- // const BaseProximity::SPtr tipProx = createTipProximity(itTip->element());
295- // if (!tipProx) return;
296- //
297- // // 2.1 Check whether coupling point should be added
298- // const type::Vec3 tipToLastCP = m_CPs.back()->getPosition() - tipProx->getPosition();
299- // if (tipToLastCP.norm() > d_tipDistThreshold.getValue())
300- // {
301- // auto findClosestProxOnVol =
302- // Operations::FindClosestProximity::Operation::get(l_volGeom);
303- // auto projectOnVol = Operations::Project::Operation::get(l_volGeom);
304- // const BaseProximity::SPtr volProx =
305- // findClosestProxOnVol(tipProx, l_volGeom.get(), projectOnVol,
306- // getFilterFunc());
307- // // Proximity can be detected before the tip enters the tetra (e.g. near a
308- // boundary
309- // // face) Only accept proximities if the tip is inside the tetra during insertion
310- // if (volProx)
311- // {
312- // TetrahedronProximity::SPtr tetProx =
313- // dynamic_pointer_cast<TetrahedronProximity>(volProx);
314- // if (tetProx)
315- // {
316- // double f0(tetProx->f0()), f1(tetProx->f1()), f2(tetProx->f2()),
317- // f3(tetProx->f3());
318- // bool isInTetra = toolbox::TetrahedronToolBox::isInTetra(
319- // tipProx->getPosition(), tetProx->element()->getTetrahedronInfo(), f0,
320- // f1, f2, f3);
321- // if (isInTetra)
322- // {
323- // volProx->normalize();
324- // m_CPs.push_back(volProx);
325- // }
326- // }
327- // }
328- // }
329290 else // Don't bother with removing the point that was just added
330291 {
331292 // 2.2. Check whether coupling point should be removed
0 commit comments