[petsc-users] Inquiry about metric-based adaptation using MMG in petsc

[neil liu]

Dear Petsc developers and users,

I am  trying to explore adaptive mesh refinement (tetradedra) with Petsc.
It seems case 12 ($PETSC DIR/src/snes/tutorials/ex12.c, from paper,
https://urldefense.us/v3/__https://arxiv.org/pdf/2201.02806__;!!G_uCfscf7eWS!ZEwaHkvYc_IrjycFr8LiNi-nZcFqu0ZpsydAAomxptQ6C0xkBs3qhn5ba31Z4vipKf4mTrqDm8A5S65DDD_xFg$ ) is a good example.
However when I tried,
./ex12 -run_type full -dm_plex_dim 3 -dm_distribute -dm_plex_box_faces
10,10,10 -bc_type dirichlet -petscspace_degree 1 -variable_coefficient ball
-snes_converged
_reason ::ascii_info_detail -ksp_type cg -pc_type sor -snes_adapt_sequence
3 -adaptor_target_num 1000 -dm_plex_metric_h_max  0.5  -dm_adaptor mmg
      L_2 Error: 1.55486
      Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 2

it shows the following error,
[0]PETSC ERROR: --------------------- Error Message
--------------------------------------------------------------
[0]PETSC ERROR: Argument out of range
[0]PETSC ERROR: Global metric normalization factor must be in (0, inf). Is
the input metric positive-definite?

Do  you have any suggestions here?

Thanks,

Xiaodong
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20240911/a83cb414/attachment.html>