<div dir="ltr"><div dir="ltr">On Mon, Sep 23, 2019 at 8:46 PM Swarnava Ghosh <<a href="mailto:swarnava89@gmail.com">swarnava89@gmail.com</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">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" target="_blank">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></blockquote><div><br></div><div>Did you give ' -dm_plex_hash_location'? It does not seem possible that you got to this code without giving that option.</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>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" target="_blank">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"><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>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_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>