Skip to content

Conversation

colemanjs
Copy link
Collaborator

Updated logic removes integration points outside of cell bounding box. This is a crude Riemann integration using the midpoint rule, that is second-order accurate for smooth functions and axis-aligned hexahedral cells.

…ls outside of bounding box. This is crude Riemann intergration using the midpoint rule, that is second-order accurate for smooth functions and axis-aligned hex cells.
@colemanjs colemanjs requested a review from kincaidkc February 7, 2025 16:35
@colemanjs colemanjs self-assigned this Feb 7, 2025
Copy link
Collaborator

@kincaidkc kincaidkc left a comment

Choose a reason for hiding this comment

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

changes look good. a couple questions:

  1. do these changes result in any speedup vs. the original version of the code?
  2. would it make sense to add logic to turn off the integration in regions of the mesh that are well-refined?

{
vector dx = cmptDivide(cellBb.span(), vector(nPts));
scalar I = 0.0;
scalar count = 0.0;
Copy link
Collaborator

Choose a reason for hiding this comment

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

this should probably be a label

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.

2 participants