Skip to content
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

Reduce generated Faces #15

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

Reduce generated Faces #15

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
618eeb6
using just one point for bottom center
Cirromulus Apr 11, 2022
bb29132
renamed height to depth. Using better centerpoint calculation.
Cirromulus Apr 12, 2022
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
141 changes: 56 additions & 85 deletions src/create_mesh.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -555,36 +555,39 @@ $(SIGNATURES)
"""
function solidify(inputMesh, offset=100)
width = inputMesh.width
height = inputMesh.height
totalNodes = width * height * 2
depth = inputMesh.height
totalNodes = width * depth * 2
nodeList = Vector{Point3D}(undef, totalNodes)
nodeArrayTop = Matrix{Point3D}(undef, width, height)
nodeArrayBottom = Matrix{Point3D}(undef, width, height)

# imagine a 4x4 image. 4 * 2 + 2 * 2 = 12
numEdgeNodes = width * 2 + (height - 2) * 2
numEdgeNodes = width * 2 + (depth - 2) * 2

numTrianglesTop = (width - 1) * (height - 1) * 2
numTrianglesBottom = numTrianglesTop
numTrianglesTop = (width - 1) * (depth - 1) * 2
numTrianglesBottom = numEdgeNodes
numTrianglesEdges = numEdgeNodes * 2

totalTriangles = numTrianglesBottom + numTrianglesTop + numTrianglesEdges

println("Specs: $(width) $(height) $(totalNodes) $(numEdgeNodes) $(numTrianglesBottom) $(totalTriangles)")
println("Specs: $(width)x$(depth). $(totalNodes) nodes, $(numEdgeNodes) edges, $(numTrianglesBottom), faces on bottom, $(totalTriangles) faces total")

centerpoint_idx = 0
# Build the bottom surface
count = 1
for y = 1:height
for y = 1:depth
for x = 1:width
newPoint = Point3D(x, y, -offset, x, y)
nodeList[count] = newPoint
nodeArrayBottom[x, y] = newPoint
# most of them are not used, but just the centerpoint for bottom
nodeList[count] = Point3D(x, y, -offset, x, y)
if y == floor(depth/2) && x == floor(width/2)
centerpoint_idx = count
end
count += 1
end
end

println("Found a centerpoint at ($(nodeList[centerpoint_idx].x), $(nodeList[centerpoint_idx].y))")

# Copy in the top surface
for y = 1:height
for y = 1:depth
for x = 1:width
node = inputMesh.nodeArray[x, y]
copiedPoint = Point3D(node.x, node.y, node.z, node.ix, node.iy)
Expand All @@ -596,45 +599,23 @@ function solidify(inputMesh, offset=100)
end

nodeList[count] = copiedPoint
nodeArrayTop[x, y] = copiedPoint
count += 1
end
end

println("We now have $(count - 1) valid nodes")

triangles = Vector{Triangle}(undef, totalTriangles)
# Build the triangles for the bottom surface
count = 1
for y = 1:(height - 1)
for x = 1:(width - 1)
# here x and y establish the column of squares we're in
index_ul = (y - 1) * width + x
index_ur = index_ul + 1

index_ll = y * width + x
index_lr = index_ll + 1

triangles[count] = Triangle(index_ul, index_ll, index_ur)
count += 1
triangles[count] = Triangle(index_lr, index_ur, index_ll)
count += 1
end
end

println("We've filled up $(count - 1) triangles")
if count != numTrianglesBottom + 1
println("Hmm aren't count and triangles bottom equal? $(count) vs $(numTrianglesBottom + 1)")
end

# Build the triangles for the top surface
for y = 1:(height - 1)
for y = 1:(depth - 1)
for x = 1:(width - 1)
# here x and y establish the column of squares we're in
index_ul = (y - 1) * width + x + totalNodes / 2
index_ul = (y - 1) * width + x + (width*depth)
index_ur = index_ul + 1

index_ll = y * width + x + totalNodes / 2
index_ll = y * width + x + (width*depth)
index_lr = index_ll + 1

triangles[count] = Triangle(index_ul, index_ur, index_ll)
Expand All @@ -644,58 +625,48 @@ function solidify(inputMesh, offset=100)
end
end

println("We've filled up $(count - 1) triangles")

# Build the triangles to close the mesh
x = 1
for y = 1:(height - 1)
ll = (y - 1) * width + x
ul = ll + totalNodes / 2
lr = y * width + x
ur = lr + totalNodes / 2
triangles[count] = Triangle(ll, ul, ur)
count += 1
triangles[count] = Triangle(ur, lr, ll)
count += 1
end

x = width
for y = 1:(height - 1)
ll = (y - 1) * width + x
ul = ll + totalNodes / 2
lr = y * width + x
ur = lr + totalNodes / 2
triangles[count] = Triangle(ll, ur, ul)
count += 1
triangles[count] = Triangle(ur, ll, lr)
count += 1
println("We've filled up $(count - 1) top triangles")

# Build the Side triangles to close the mesh
for x = [1, width]
for y = 1:(depth - 1)
ll = (y - 1) * width + x
ul = ll + (width * depth)
lr = y * width + x
ur = lr + (width * depth)
triangles[count] = Triangle(ll, ul, ur)
count += 1
triangles[count] = Triangle(ur, lr, ll)
count += 1
triangles[count] = Triangle(lr, ll, centerpoint_idx)
count += 1
end
end

y = 1
for x = 2:width
ll = (y - 1) * width + x
ul = ll + totalNodes / 2
lr = (y - 1) * width + (x - 1)
ur = lr + totalNodes / 2
triangles[count] = Triangle(ll, ul, ur)
count += 1
triangles[count] = Triangle(ur, lr, ll)
count += 1
for y = [1, depth]
for x = 2:width
ll = (y - 1) * width + x
ul = ll + (width * depth)
lr = (y - 1) * width + (x - 1)
ur = lr + (width * depth)
triangles[count] = Triangle(ll, ul, ur)
count += 1
triangles[count] = Triangle(ur, lr, ll)
count += 1
triangles[count] = Triangle(lr, ll, centerpoint_idx)
count += 1
end
end

println("We've filled up $(count - 1) triangles, including bottom")

y = height
for x = 2:width
ll = (y - 1) * width + x
ul = ll + totalNodes / 2
lr = (y - 1) * width + (x - 1)
ur = lr + totalNodes / 2
triangles[count] = Triangle(ll, ur, ul)
count += 1
triangles[count] = Triangle(ur, ll, lr)
count += 1
if count-1 != totalTriangles
println("Hmm aren't count and triangles bottom equal? $(count-1) vs $(totalTriangles)")
end

Mesh(nodeList, nodeArrayBottom, triangles, width, height)

# Array not actually used by calling functions here
sumArray = Matrix{Point3D}(undef, width, depth)
Mesh(nodeList, sumArray, triangles, width, depth)
end


Expand Down