[MOAB-dev] AdaptiveKDTree and point search

[Lukasz Kaczmarczyk]

Hello,

I switch recently to MOAB 5.0.0 and I encounter some problems with point search.  This is for tet mesh, I expecting that I doing something wrong. It was working with older MOAB. I can work with leaf, but will be nice to have this ok.

   myTree = boost::shared_ptr<AdaptiveKDTree>(new AdaptiveKDTree(&mwlsMoab));
   PetscPrintf(mField.get_comm(),"Build tree ... ");
   rval = myTree->build_tree(levelTets,&treeRoot); CHKERRQ_MOAB(rval);

   boost::shared_ptr<ElemEvaluator> elemEval(myTree,new ElemEvaluator(&mwlsMoab));
   myTree->set_eval(elemEval.get());

   PetscPrintf(mField.get_comm(),"done\n”);

and using it like that

 rval = myTree->point_search(
         material_coords,leaf_out,maxEdgeL*1e-5,maxEdgeL*1e-3,NULL,NULL,&params
       ); CHKERRQ_MOAB(rval);

I get following error,
ssertion failed: (nverts == 4 && verts), function initFcn, file LocalDiscretization/LinearTet.cpp, line 21.
[likask:32071] *** Process received signal ***
[likask:32071] Signal: Abort trap: 6 (6)
[likask:32071] Signal code:  (0)
[likask:32071] [ 0] 0   libsystem_platform.dylib            0x00007fff9d14f52a _sigtramp + 26
[likask:32071] [ 1] 0   ???                                 0x00007fff5acbdea0 0x0 + 140734716698272
[likask:32071] [ 2] 0   libsystem_c.dylib                   0x00007fff8aca76df abort + 129
[likask:32071] [ 3] 0   libsystem_c.dylib                   0x00007fff8ac6edd8 basename + 0
[likask:32071] [ 4] 0   libMOAB.0.dylib                     0x000000010a4fc1ec _ZN4moab9LinearTet7initFcnEPKdiRPd + 108
[likask:32071] [ 5] 0   libMOAB.0.dylib                     0x000000010a1f5cad _ZN4moab13ElemEvaluator14set_ent_handleEm + 1309
[likask:32071] [ 6] 0   libMOAB.0.dylib                     0x000000010a4f6e45 _ZN4moab13ElemEvaluator22find_containing_entityERNS_5RangeEPKdddRmPdPj + 197
[likask:32071] [ 7] 0   libMOAB.0.dylib                     0x000000010a19c033 _ZN4moab13ElemEvaluator22find_containing_entityEmPKdddRmPdPj + 371
[likask:32071] [ 8] 0   libMOAB.0.dylib                     0x000000010a18ef65 _ZN4moab14AdaptiveKDTree12point_searchEPKdRmddPbPmPNS_8CartVectE + 2069
[likask:32071] [ 9] 0   mwls                                0x00000001051a8035 _ZN17FractureMechanics10MWLSApprox17getInfluenceNodesEPdb + 1077
[likask:32071] [10] 0   mwls                                0x0000000104f47913 main + 22179
[likask:32071] [11] 0   libdyld.dylib                       0x00007fff904af5ad start + 1
[likask:32071] [12] 0   ???                                 0x000000000000000b 0x0 + 11


Can you pleas advise me how I can fix this?

Kind regards,
Lukasz

[University of Glasgow: The Times Scottish University of the Year 2018]