[Nek5000-users] The number of adjacent element

[nek5000-users at lists.mcs.anl.gov]

Something like:


      real x(lx1,ly1,lz1,lelt)

      integer e

      real      eg


      nxyz = lx1*ly1*lz1

      do e=1,nelv

           eg = lglel(e)

            call cfill(x(1,1,1,e),eg, nxyz)

            call cfill(y(1,1,1,e),eg, nxyz)

       enddo


       call dssum(x,nx1,ny1,nz1)


       do e=1,nelv

            eg= -lglel(e)

            call cadd(x(1,1,1,e),eg,nxyz)

        enddo


Now the faces tell you which faces you're connected to.


To get edges and corners is a bit trickier, but not hard --

use dsop(x,'min') and max etc. to discriminate...  A few

passes through the process will reveal all the connections.


hth


Paul


________________________________
From: Nek5000-users <nek5000-users-bounces at lists.mcs.anl.gov> on behalf of nek5000-users at lists.mcs.anl.gov <nek5000-users at lists.mcs.anl.gov>
Sent: Wednesday, August 29, 2018 10:34:25 PM
To: nek5000-users
Subject: [Nek5000-users] The number of adjacent element

Hi, Neks!

Could I easy identify the number of adjacent element? For example, if I know, that now I am in the "e=25" element, could I know "e" neighbors? I know, that the mesh is globaly unstructured, but anyway?

Best regards,
Valdemar
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180830/51275312/attachment.html>