Skip to content

ac/adapt-changes capVol fails to swap for cylinder region mesh #510

@bobpaw

Description

@bobpaw

On the ac/adapt-changes branch (b15f7e2), MeshAdapt fails for the capVol cylinder region mesh. The log file is included below. I get the error when I adapt with an MDS mesh as well (the -m flag). During the first iteration, face swap is missing and swapping enters an endless loop.

./capVol -g 1 pumi-meshes/cap/cyl_surf_only.cre out.cre
Initalization of MetricInterpolation ...
Threads : 1
Expiration Date = 2025-8-26
Number of Breps = 1
Loading Data Timer : cpu = 9.155e-03s  wall = 9.57e-03s
GenerateNMesh  : Num Threads = 1
Number of sources : 130
Generating isotropic proxy
For iteration:0 marked cells:0 out of:1210
Max source number:92
Average source number:22.6843
Sizing Field Timing : cpu = 2.806e-02s  wall = 4.01e-01s
Meshing edge 1
Meshing edge 2
Meshing edge 3
Meshing face 1
Meshing face 2
Meshing face 3
Faces Mesher : cpu = 5.432e-02s  wall = 5.86e-02s
Meshing region 1
Missing edges:0 over:2904
Iter=0 npoin=970 nelem=2805
Iter=1 npoin=971 nelem=3203
Iter=2 npoin=971 nelem=3203
Iter=0 npoin=971 nelem=3203 nface=1936
Iter=1 npoin=1560 nelem=6060 nface=3192
Iter=2 npoin=1561 nelem=6088 nface=2750
Iter=3 npoin=1579 nelem=6204 nface=1232
Iter=4 npoin=1871 nelem=8014 nface=1688
Iter=5 npoin=1872 nelem=8045 nface=1380
Iter=6 npoin=1889 nelem=8166 nface=642
Iter=7 npoin=2013 nelem=8933 nface=780
Iter=8 npoin=2017 nelem=8977 nface=622
Iter=9 npoin=2029 nelem=9051 nface=294
Iter=10 npoin=2082 nelem=9375 nface=374
Iter=11 npoin=2085 nelem=9400 nface=266
Iter=12 npoin=2091 nelem=9442 nface=144
Iter=13 npoin=2098 nelem=9488 nface=110
Iter=14 npoin=2099 nelem=9498 nface=58
Iter=15 npoin=2099 nelem=9501 nface=28
Worst quality before optim:362.914
Average quality before optim:1.9096
Iter=0
Edges swapped:612 over:12566 tried:9662
Elements swapped:0 over:9501 tried:35508
Points optimized:644 over:1129 tried:1129
Iter=1
Edges swapped:73 over:12566 tried:7009
Elements swapped:0 over:9501 tried:30098
Points optimized:274 over:1129 tried:1034
Iter=2
Edges swapped:10 over:12566 tried:3478
Elements swapped:0 over:9501 tried:18904
Points optimized:103 over:1129 tried:824
Average quality after optim:1.40664
Worst quality after optim:2.65851
Regions Mesher : cpu = 8.254e-01s  wall = 9.36e-01s
Timing(GenerateNMesh) : cpu = 9.085e-01s  wall = 1.40e+00s
mesh verified in 0.093678 seconds

MeshAdapt: ma::configure:  Based on requested sizefield, MeshAdapt requires at least 2 iterations.
           Setting the number of iteration to 3!

MeshAdapt: version 2.0 !

MeshAdapt: coarsened 89 edges in 0.840241 seconds

MeshAdapt: iteration 0

MeshAdapt: refined 5394 edges in 2.808631 seconds

MeshAdapt: Swap failed: face swap not implemented

MeshAdapt: Swap failed: face swap not implemented

MeshAdapt: Swap failed: face swap not implemented

The last message is repeated, so I hit Ctrl-C. The mesh is the same one in #509.

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugcapstonemesh adaptTopics related to mesh adapt (straight-sided or curved)

    Type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions