Open
Conversation
This commit 1. Adds a dyad_1D_t type,, 2. Adds a vector_1D_t operator(/) with a RHS integer operand, 3. Adds operator(*) for vector_1D_t operands & a dyad_1D_t result, 4. Updates the burgers-1D example to compute u*u/2.
This commit adds support expressions of the form .div. uu, where uu is a 1D dyadic tensor and the result is stored at cell centers. To Do: 1. Embed an interpolation operator into the divergence operator to interpolate the result back out to the cell faces. 2. Use the interpolated values to construct the vector_1D_t result.
This commit prepares for disambiguation of two types of mimetic matrix abstractions: one for differential operators and one for interpolation operators.
This commit adds interpolation operator matrix abstractions to go from cell centers to faces and vice versa based on the definitions in Dumett & Castillo (2022). https://doi.org/10.13140/RG.2.2.26630.14400
This commit adds a centers-to-faces interpolator and uses that interpolator when computing the divergence of a dyad.
This commit adds the interpolate_1D_t type with a type-bound procedure that interpolates using values at cell centers plus boundaries to estimate values at faces using the 2nd- and 4th-order schemes of Dumett & Castillo (2022).
This commit adds the faces_to_centers_t type with a type-bound procedure that interpolates using values at the faces to estimate to estimate values at cell centers plus boundaries using the 2nd- and 4th-order schemes of Dumett & Castillo (2022).
This file contains hidden or 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
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
1 participant
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.
This pull request adds
faces_to_centers_1D_ttype witha. a user-defined structures constructor and
b. a
center_valuestype-bound function;centers_to_faces_1D_ttype witha. a user-defined structures constructor and
b. a
face_valuestype-bound function;faces_to_centers_1D_tandcenters_to_faces_1D_t; and.div.to adyad_1D_tobject .Given values at cell faces,
center_valuesestimates the corresponding values at the cell centers augmented to include the celll boundaries; whereasface_valuesis the inverse ofcenter_values. Both estimators use the 2nd- or 4th-order schemes of Dumett & Castillo (2022) based on the order specified when constructing the interpolator object.The experimental dyad-divergence operator uses the
centers_to_faces_1D_tinterpolator to estimate the dyad at cell faces. computing the