-
Notifications
You must be signed in to change notification settings - Fork 102
Add Ferrite.write_facetset #1259
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
base: master
Are you sure you want to change the base?
Conversation
Codecov Report✅ All modified and coverable lines are covered by tests. Additional details and impacted files@@ Coverage Diff @@
## master #1259 +/- ##
=======================================
Coverage 94.19% 94.19%
=======================================
Files 40 40
Lines 6662 6671 +9
=======================================
+ Hits 6275 6284 +9
Misses 387 387 ☔ View full report in Codecov by Sentry. 🚀 New features to boost your workflow:
|
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I haven't checked the output, but it looks like it could fail for higher order geometric interpolations?
Nvm, looks correct on 2nd thought !
Edit3: Or no, it does look wrong...
| end | ||
|
|
||
| """ | ||
| write_facetset(vtk::VTKGridFile, grid::AbstractGrid, facetset::String) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Would be nice to match the interface from cellsets allowing a vector/set of keys?
Also, since this is not as trivial as cellsets, perhaps it could even be nice to have a lower-level interface providing just an AbstractVecOrSet{<:FacetIndex}, not requiring it to be stored in the grid?
| function write_facetset(vtk, grid::AbstractGrid, facetset::String) | ||
| z = zeros(getnnodes(grid)) | ||
| for fi in getfacetset(grid, facetset) | ||
| for i in facets(getcells(grid, fi[1]))[fi[2]] |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
To make work for higher order elements, I think this is required.
| for i in facets(getcells(grid, fi[1]))[fi[2]] | |
| cell = getcells(grid, fi[1] | |
| for cellnode in facetdof_indices(geometric_interpolation(cell))[fi[2]] | |
| i = get_node_ids(cell)[cellnode] |
|
Btw will (or can) this close #192, or is that already resolved? |
No description provided.