<div dir="ltr">Hi Matt,<div><br></div><div>I am trying to get this working. However, It seems that this does not work for 3D. I have tets as elements, and the dmplex is serial. I get the following error:</div><div><br></div><div>0]PETSC ERROR: No support for this operation for this object type<br>[0]PETSC ERROR: I have only coded this for 2D<br>[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.<br></div><div><br></div><div>Is there a work around to this? would you please let me know?</div><div><br></div><div>Also here is my code<br></div><div><br></div><div> PetscSF cellSF=NULL;<br>  Vec v;<br>  PetscErrorCode ierr;<br></div><div><br></div><div>// create v;<br>  VecCreate(PETSC_COMM_SELF,&v);<br>  VecSetSizes(v,PETSC_DECIDE,3);<br>   VecSetBlockSize(v,3);<br>  VecSetFromOptions(v);<br><br>  VecSetValue(v,0,0.0,INSERT_VALUES);<br>  VecSetValue(v,1,0.1,INSERT_VALUES);<br>  VecSetValue(v,2,0.12,INSERT_VALUES);<br>  <br>  VecAssemblyBegin(v);<br>  VecAssemblyEnd(v);<br><br>  PetscInt bs;<br>  VecGetBlockSize(v,&bs);<br><br>  printf("Block size of v=%d \n",bs);<br>  <br>    ierr=DMLocatePoints(pCgdft->dmatom,v,DM_POINTLOCATION_NEAREST,&cellSF);<br><br>  // print vector<br>  VecView(v,PETSC_VIEWER_STDOUT_SELF);<br></div><div><br></div><div><br></div><div>Sincerely,</div><div>SG</div><div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Sep 16, 2019 at 6:37 AM Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">On Fri, Sep 6, 2019 at 6:07 PM Swarnava Ghosh via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear Petsc developers and users,<div><br></div><div>I have a DMPlex mesh in 3D. Given a point with (x,y,z) coordinates, I am trying the find the cell number in which this point lies, and the vertices of the cell. Is there any DMPlex function that will give me the cell number?</div></div></blockquote><div><br></div><div>Sorry, I lost this mail.</div><div><br></div><div>In serial, you can just use DMLocatePoint(). If you have some points and you are not</div><div>sure which process they might be located on, then you need a DMInterpolation context.</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Thank you,</div><div>SG  </div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail-m_-7519753545807991045m_-2142680694173864663gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>
</blockquote></div>