Page 1 of 1

EdgeIndexes

Posted: 23 Feb 2015, 12:07
by amonnier
Hello,
For my project, I need to identify the element which are in a level set fonction and recreate the boundary with the exterior element edge.
I didn't fine such subroutine an LevelSet module.
So I make my own subroutine, but I have some difficulties this the EdgeIndexes.

Code: Select all

     
   DO i=1,nbrEltSolid
        EltSolidB = ElementSolid(i)
        Element =>Solver % Mesh % Elements(EltSolidB)
        EdgeIndexes => Element % EdgeIndexes
        DO j=1,Element %  Type % NumberOfEdges   ! boule sur j, sur le nombre d'edge
              write(*,*) "element n°", EltSolidB, "edge n°", j
              k = Element % EdgeIndexes(j)
        END DO   ! fin boule sur j, sur le nombre d'edge
     END DO   ! fin boule sur i, sur les élements solides
The problem is at the line "k = Element %EgdeIndexes (j)", with k define as a integer.

Thanks,
Antoine

Re: EdgeIndexes

Posted: 23 Feb 2015, 18:08
by raback
Hi

You should probably call "FindMeshEdges" prior to trying to print these out.

-Peter

Re: EdgeIndexes

Posted: 24 Feb 2015, 12:38
by amonnier
Hello Peter,

Thanks for your help, I will test that.
Regards,