Skip to content

Commit

Permalink
comments and formatting
Browse files Browse the repository at this point in the history
  • Loading branch information
jlchan committed Jun 28, 2024
1 parent 42308c9 commit 7ec2f0a
Show file tree
Hide file tree
Showing 2 changed files with 2 additions and 4 deletions.
2 changes: 1 addition & 1 deletion src/MeshData.jl
Original file line number Diff line number Diff line change
Expand Up @@ -332,7 +332,7 @@ function MeshData(VX, VY, VZ, EToV, rd::RefElemData{3}; is_periodic=(false, fals
x, y, z = (x -> V1 * x[transpose(EToV)]).((VX, VY, VZ))

#Compute connectivity maps: uP = exterior value used in DG numerical fluxes
(; r, s, t, Vf ) = rd
(; r, s, t, Vf) = rd
xf, yf, zf = (x -> Vf * x).((x, y, z))
mapM, mapP, mapB = build_node_maps(rd, FToF, (xf, yf, zf))
mapM = reshape(mapM, :, K)
Expand Down
4 changes: 1 addition & 3 deletions src/connectivity_functions.jl
Original file line number Diff line number Diff line change
Expand Up @@ -166,7 +166,7 @@ end

"""
make_periodic(md::MeshData{Dim}, is_periodic...) where {Dim}
make_periodic(md::MeshData{Dim}, is_periodic = ntuple(x->true,Dim)) where {Dim}
make_periodic(md::MeshData{Dim}, is_periodic = ntuple(x -> true, Dim)) where {Dim}
make_periodic(md::MeshData, is_periodic = true)
Returns new MeshData such that the node maps `mapP` and face maps `FToF` are now periodic.
Expand Down Expand Up @@ -230,7 +230,6 @@ function build_periodic_boundary_maps!(xf, yf, is_periodic_x, is_periodic_y,
ymin, ymax = extrema(yc)

LX, LY = map((x -> x[2] - x[1]) extrema, (xf, yf))
#NODETOL = 100 * max(eps.((LX, LY))...)
NODETOL = tol * max(LX, LY)
if abs(abs(xmax - xmin) - LX) > NODETOL && is_periodic_x
error("periodicity requested in x, but LX = $LX while abs(xmax-xmin) = $(abs(xmax-xmin))")
Expand Down Expand Up @@ -288,7 +287,6 @@ function build_periodic_boundary_maps!(xf, yf, zf,
Flist = 1:length(FToF)
Bfaces = findall(vec(FToF) .== Flist)

# xc, yc, zc = compute_boundary_centroids(xf, yf, zf)
xb, yb, zb = xf[mapB], yf[mapB], zf[mapB]
Nfp = length(xf) ÷ NfacesTotal
Nbfaces = length(xb) ÷ Nfp
Expand Down

0 comments on commit 7ec2f0a

Please sign in to comment.