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.

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

Example of using DelaunayTriangulation.jl for mesh generation #1040

Open
wants to merge 13 commits into
base: master
Choose a base branch
from

Conversation

DanielVandH
Copy link

Closes #975

Wasn't too sure where to put the example, so I just put in the howto section.

I explain some things that are also explained in more detail in DelaunayTriangulation.jl's docs since I wanted to
try to make this example somewhat self-contained, but I give links to my docs for more information. The mesh refinement here is pretty simple, but I note in the explanation that it could be easily adapted for things like graded meshes or adaptive refinement.

Copy link

codecov bot commented Aug 4, 2024

Codecov Report

All modified and coverable lines are covered by tests ✅

Project coverage is 0.00%. Comparing base (6eead25) to head (6b7f531).

❗ There is a different number of reports uploaded between BASE (6eead25) and HEAD (6b7f531). Click for more details.

HEAD has 5 uploads less than BASE
Flag BASE (6eead25) HEAD (6b7f531)
6 1
Additional details and impacted files
@@            Coverage Diff             @@
##           master   #1040       +/-   ##
==========================================
- Coverage   93.74%   0.00%   -93.75%     
==========================================
  Files          39      39               
  Lines        6238    6152       -86     
==========================================
- Hits         5848       0     -5848     
- Misses        390    6152     +5762     

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

# The solution can be visualised as in previous tutorials. Alternatively, we can use `tricontourf` which
# knows how to use DelaunayTriangulation.jl's `Triangulation`.
sol = evaluate_at_grid_nodes(dh, u, :u)
tricontourf(tri, sol)
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

IMHO it would be advantageous to plot with FerriteViz

Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I wanted to do that originally but didn't because it doesn't work on Ferrite master. I could dev in Ferrite-FEM/FerriteViz.jl#103 to make it work?

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

thats fine in my opinion but ping @KnutAM @termi-official @fredrikekre

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Seems like a good idea to too!
The only downside I see, is that it requires that branch to be continuously kept up to date with Ferrite#master changes now in the time before releasing. But not sure how hard / much work that is?

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I agree with Maxi, but let us do this after 1.0 release, so we can make a FerriteViz release instead of relying on me keeping the FerriteViz branch up to date.

Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I guess this is just waiting on FerriteViz to be made compatible with Ferrite 1.0 before reviewing further?

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Right. The release is already in review, so we should be able to update FerriteViz to 1.0 in the not too far future.

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

You might be wondering, but we are still preparing the release. Several new "bugs" appeared in FerriteViz which block the release.

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Should we just proceed with Plots jl for now and add some message to not delay this further @koehlerson ?

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Copy link
Member

@KnutAM KnutAM left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Cool!
Just a couple of drive-by comments when looking through :)

@DanielVandH
Copy link
Author

I've addressed the review comments from @KnutAM. Let me know if there's anything else.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Example integrating with DelaunayTriangulation.jl
4 participants