Skip to content

Commit

Permalink
Avoid processing elements which haven't been created
Browse files Browse the repository at this point in the history
  • Loading branch information
mmalinen committed Feb 6, 2025
1 parent 28602d3 commit eb7e1ec
Show file tree
Hide file tree
Showing 2 changed files with 4 additions and 4 deletions.
6 changes: 3 additions & 3 deletions fem/src/modules/Isosurface.F90
Original file line number Diff line number Diff line change
Expand Up @@ -363,14 +363,14 @@ SUBROUTINE IsosurfaceSolver( Model,Solver,dt,Transient )
NoEdges = 0
CALL Info('IsosurfaceSolver','Creating mesh edges',Level=9)
IF (.NOT.ASSOCIATED(Mesh % Edges)) THEN
NoOrigBoundaryElements = Mesh % NumberOfBoundaryElements
Mesh % NumberOfBoundaryElements = 0
IF( dim == 2 ) THEN
NoOrigBoundaryElements = Mesh % NumberOfBoundaryElements
Mesh % NumberOfBoundaryElements = 0
CALL FindMeshEdges2D(Mesh)
Mesh % NumberOfBoundaryElements = NoOrigBoundaryElements
ELSE
CALL FindMeshEdges3D(Mesh)
END IF
Mesh % NumberOfBoundaryElements = NoOrigBoundaryElements
END IF
NoEdges = Mesh % NumberOfEdges

Expand Down
2 changes: 1 addition & 1 deletion fem/tests/IsosurfaceSTL/case3d.sif
Original file line number Diff line number Diff line change
Expand Up @@ -116,6 +116,6 @@ Boundary Condition 3
End

Solver 1 :: Reference Norm = 1.13166182E+00
Solver 2 :: Reference Norm = 2.09430000E+04
Solver 2 :: Reference Norm = 2.88120000E+04

!End Of File

0 comments on commit eb7e1ec

Please sign in to comment.