Skip to content

Reproducibility for EB covered state#702

Merged
jrood-nrel merged 7 commits intoAMReX-Combustion:developmentfrom
baperry2:set-eb-state
Oct 11, 2023
Merged

Reproducibility for EB covered state#702
jrood-nrel merged 7 commits intoAMReX-Combustion:developmentfrom
baperry2:set-eb-state

Conversation

@baperry2
Copy link
Copy Markdown
Contributor

@baperry2 baperry2 commented Oct 6, 2023

The state in cells covered by the EB is currently is set based on the value of the state in a somewhat arbitrarily selected fluid cell. Right now, the lexicographically first fluid cell in the first FAB containing fluid cells on the highest rank processor with a fluid cell is selected to populate the covered state.

The identity of this cell depends on the domain decomposition and therefore the covered state may change for a particular case when the number of processors or amr.max_grid_size is changed. See below, for EB-C8 in 2D with 1 and 4 processors.

With this PR, the lexicographically first fluid cell in the entire domain is chosen, so it is no longer depends on domain decomposition. As a result, both cases below will now be equivalent and match the 1 proc version. The code to do this is not very pretty, but I wasn't sure of a better way.

This PR should have no effect on the actual solutions, but may lead to some diffs due to covered cells changing.

1 proc, current
Screen Shot 2023-10-06 at 3 51 37 PM

4 proc, current
Screen Shot 2023-10-06 at 3 51 09 PM

@baperry2 baperry2 requested a review from marchdf October 6, 2023 22:00
@marchdf
Copy link
Copy Markdown
Contributor

marchdf commented Oct 9, 2023

Ah this is great. Thank you for looking at this. Long been on my todo list... I think there's a better way though. Let's chat.

@baperry2 baperry2 requested a review from marchdf October 10, 2023 00:35
@baperry2 baperry2 enabled auto-merge (squash) October 10, 2023 18:26
@jrood-nrel jrood-nrel disabled auto-merge October 11, 2023 16:34
@jrood-nrel jrood-nrel merged commit c0d01b4 into AMReX-Combustion:development Oct 11, 2023
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.

3 participants