Skip to content

Dump the outer-boundary input next to outer.msh #41

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
krystophny wants to merge 4 commits into
base: main
Choose a base branch
from
Open

Dump the outer-boundary input next to outer.msh #41

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
d503a15
Dump the outer-boundary input next to outer.msh
krystophny Aug 20, 2026
f67f4ff
build: pin libneo Fortio 0.2.0
krystophny Aug 25, 2026
65d984a
build: track Fortio lifecycle bridge
krystophny Aug 25, 2026
4982f72
build: use main-line libneo lifecycle bridge
krystophny Aug 25, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion cmake/SetupCODE.cmake
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,6 +5,6 @@
# -DLIBNEO_PATH=<dir> local libneo source directory (empty = fetch from git)
# -DTRIANGLE_DIR=<dir> prebuilt Triangle directory (see SetupTriangle)
# -DMFEM_DIR=<dir> prebuilt MFEM directory (see SetupMFEM)
set(LIBNEO_REF "a4620f8075698dabafd774859e28e0ae92e8f339" CACHE STRING
set(LIBNEO_REF "7262c32bf0406594a8ef27e84654d249aef2cb86" CACHE STRING
"libneo branch, tag, or SHA to fetch")
set(LIBNEO_PATH "" CACHE PATH "local libneo source directory (empty = fetch from git)")
11 changes: 11 additions & 0 deletions src/mephit_mesh.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3108,6 +3108,17 @@ subroutine write_FreeFem_mesh
theta(:) = linspace(0d0, 2d0 * pi, npt_outer, 0, 1)
bdry_R(npt_inner+1:) = R_mid + R_rad * cos(theta)
bdry_Z(npt_inner+1:) = Z_mid + Z_rad * sin(theta)
! dump the triangulator input so the annulus meshing problem can be
! reproduced without a MEPHIT build, e.g. for mesher comparisons
open(newunit = fid, file = decorate_filename('outer_boundary.dat', '', basename_suffix), &
status = 'replace', form = 'formatted', action = 'write')
write (fid, '(2(1x, i0))') npt_inner, npt_outer
! 17 significant digits so the dump round-trips the double precision input
write (fid, '(2(1x, es24.16e3))') R_mid, Z_mid
do kpoi = 1, npt_inner + npt_outer
write (fid, '(2(1x, es24.16e3))') bdry_R(kpoi), bdry_Z(kpoi)
end do
close(fid)
call FEM_triangulate_external(npt_inner, npt_outer, bdry_R, bdry_Z, R_mid, Z_mid, &
decorate_filename('outer.msh', '', basename_suffix) // c_null_char)
deallocate(bdry_R, bdry_Z, theta)
Expand Down