Add demo for interpolating expressions #463
Merged
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
jorgensd
commented
Feb 21, 2022
michalhabera
approved these changes
Feb 22, 2022
garth-wells
requested changes
Feb 22, 2022
| @@ -0,0 +1,63 @@ | |||
| # Copyright (C) 2022 Jørgen S. Dokken | |||
There was a problem hiding this comment.
Use the short SPDX license header.
|
|
||
| # Define mesh | ||
| cell = triangle | ||
| v_el = VectorElement("Lagrange", cell, 1) |
There was a problem hiding this comment.
I'd suggest reducing the number of variables to improve readabilty, e.g.
mesh = Mesh(VectorElement("Lagrange", cell, 1))| # Define mixed function space | ||
| el = FiniteElement("CG", cell, 2) | ||
| el_int = VectorElement("Discontinuous Lagrange", cell, 1) | ||
| me = MixedElement([el, el_int]) |
There was a problem hiding this comment.
element = MixedElement([FiniteElement("Lagrange", cell, 2), VectorElement("Discontinuous Lagrange", cell, 1)])| b_cell = basix.cell.string_to_type(cell.cellname()) | ||
|
|
||
| # Find quadrature points for quadrature element | ||
| b_rule = basix.quadrature.string_to_type(q_rule) |
There was a problem hiding this comment.
Can q_rule and q_degree be gotten directly from q_el?
|
|
||
| # Get interpolation points for output space | ||
| family = basix.finite_element.string_to_family("Lagrange", cell.cellname()) | ||
| b_element = basix.create_element(family, b_cell, 4, basix.LagrangeVariant.gll_warped, True) |
| # Get interpolation points for output space | ||
| family = basix.finite_element.string_to_family("Lagrange", cell.cellname()) | ||
| b_element = basix.create_element(family, b_cell, 4, basix.LagrangeVariant.gll_warped, True) | ||
| interpolation_points = b_element.points |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Superseding #453, using
basixto obtain interpolation/quadrature points to generate integration kernels for various expressions that can be used for interpolation in C++.