[petsc-users] DMPlex cell number containing a point in space

[Swarnava Ghosh]

Hi Matt,

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:

0]PETSC ERROR: No support for this operation for this object type
[0]PETSC ERROR: I have only coded this for 2D
[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for
trouble shooting.

Is there a work around to this? would you please let me know?

Also here is my code

 PetscSF cellSF=NULL;
  Vec v;
  PetscErrorCode ierr;

// create v;
  VecCreate(PETSC_COMM_SELF,&v);
  VecSetSizes(v,PETSC_DECIDE,3);
   VecSetBlockSize(v,3);
  VecSetFromOptions(v);

  VecSetValue(v,0,0.0,INSERT_VALUES);
  VecSetValue(v,1,0.1,INSERT_VALUES);
  VecSetValue(v,2,0.12,INSERT_VALUES);

  VecAssemblyBegin(v);
  VecAssemblyEnd(v);

  PetscInt bs;
  VecGetBlockSize(v,&bs);

  printf("Block size of v=%d \n",bs);

    ierr=DMLocatePoints(pCgdft->dmatom,v,DM_POINTLOCATION_NEAREST,&cellSF);

  // print vector
  VecView(v,PETSC_VIEWER_STDOUT_SELF);


Sincerely,
SG

On Mon, Sep 16, 2019 at 6:37 AM Matthew Knepley <knepley at gmail.com> wrote:

> On Fri, Sep 6, 2019 at 6:07 PM Swarnava Ghosh via petsc-users <
> petsc-users at mcs.anl.gov> wrote:
>
>> Dear Petsc developers and users,
>>
>> 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?
>>
>
> Sorry, I lost this mail.
>
> In serial, you can just use DMLocatePoint(). If you have some points and
> you are not
> sure which process they might be located on, then you need a
> DMInterpolation context.
>
>   Thanks,
>
>     Matt
>
>
>> Thank you,
>> SG
>>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>
> https://www.cse.buffalo.edu/~knepley/
> <http://www.cse.buffalo.edu/~knepley/>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190923/8cf85a67/attachment.html>