Skip to content

Body force convergence, convergence issues #103

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.

Sign up for GitHub

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

Closed
jeremyfirst22 wants to merge 8 commits into from
Closed

Body force convergence, convergence issues #103

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
add8c60
bugfix: kernel packing of send buffers must complete before initiatin…
jeremyfirst22 Oct 24, 2025
1681cb6
bugfix: restart to initialize from previous timestep
jeremyfirst22 Oct 17, 2025
38a21f7
bugfix: min_steady_timesteps in lbpm_color_simulator
jeremyfirst22 Oct 24, 2025
a4dee77
lbpm_color_unsteady_simulator for sw convergence
jeremyfirst22 Oct 17, 2025
1ee867e
body force adjustment: scale-then-bisect algorithm
jeremyfirst22 Oct 20, 2025
ebe2ae5
cleanup: fixed point representation not useful for Ca and body force,…
jeremyfirst22 Oct 20, 2025
7e060fc
cleanup: consolidate color model forward propagation steps
jeremyfirst22 Oct 24, 2025
4f3be0c
cleanup: build warnings unused variables
jeremyfirst22 Oct 20, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
16 changes: 15 additions & 1 deletion common/ScaLBL.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2544,6 +2544,8 @@ void ScaLBL_Communicator::SendD3Q7AA(double *Aq, int Component) {
//...Packing for x face(2,8,10,12,14)................................
ScaLBL_D3Q19_Pack(2, dvcSendList_x, 0, sendCount_x, sendbuf_x,
&Aq[Component * 7 * N], N);

ScaLBL_DeviceBarrier();
req1[0] =
MPI_COMM_SCALBL.Isend(sendbuf_x, sendCount_x, rank_x, sendtag + 0);
req2[0] =
Expand All @@ -2552,6 +2554,8 @@ void ScaLBL_Communicator::SendD3Q7AA(double *Aq, int Component) {
//...Packing for X face(1,7,9,11,13)................................
ScaLBL_D3Q19_Pack(1, dvcSendList_X, 0, sendCount_X, sendbuf_X,
&Aq[Component * 7 * N], N);

ScaLBL_DeviceBarrier();
req1[1] =
MPI_COMM_SCALBL.Isend(sendbuf_X, sendCount_X, rank_X, sendtag + 1);
req2[1] =
Expand All @@ -2560,6 +2564,8 @@ void ScaLBL_Communicator::SendD3Q7AA(double *Aq, int Component) {
//...Packing for y face(4,8,9,16,18).................................
ScaLBL_D3Q19_Pack(4, dvcSendList_y, 0, sendCount_y, sendbuf_y,
&Aq[Component * 7 * N], N);

ScaLBL_DeviceBarrier();
req1[2] =
MPI_COMM_SCALBL.Isend(sendbuf_y, sendCount_y, rank_y, sendtag + 2);
req2[2] =
Expand All @@ -2568,6 +2574,8 @@ void ScaLBL_Communicator::SendD3Q7AA(double *Aq, int Component) {
//...Packing for Y face(3,7,10,15,17).................................
ScaLBL_D3Q19_Pack(3, dvcSendList_Y, 0, sendCount_Y, sendbuf_Y,
&Aq[Component * 7 * N], N);

ScaLBL_DeviceBarrier();
req1[3] =
MPI_COMM_SCALBL.Isend(sendbuf_Y, sendCount_Y, rank_Y, sendtag + 3);
req2[3] =
Expand All @@ -2576,6 +2584,8 @@ void ScaLBL_Communicator::SendD3Q7AA(double *Aq, int Component) {
//...Packing for z face(6,12,13,16,17)................................
ScaLBL_D3Q19_Pack(6, dvcSendList_z, 0, sendCount_z, sendbuf_z,
&Aq[Component * 7 * N], N);

ScaLBL_DeviceBarrier();
req1[4] =
MPI_COMM_SCALBL.Isend(sendbuf_z, sendCount_z, rank_z, sendtag + 4);
req2[4] =
Expand All @@ -2584,6 +2594,8 @@ void ScaLBL_Communicator::SendD3Q7AA(double *Aq, int Component) {
//...Packing for Z face(5,11,14,15,18)................................
ScaLBL_D3Q19_Pack(5, dvcSendList_Z, 0, sendCount_Z, sendbuf_Z,
&Aq[Component * 7 * N], N);

ScaLBL_DeviceBarrier();
req1[5] =
MPI_COMM_SCALBL.Isend(sendbuf_Z, sendCount_Z, rank_Z, sendtag + 5);
req2[5] =
Expand Down Expand Up @@ -2697,6 +2709,7 @@ void ScaLBL_Communicator::TriSendD3Q7AA(double *Aq, double *Bq, double *Cq) {

//...................................................................................
// Send all the distributions
ScaLBL_DeviceBarrier();
req1[0] =
MPI_COMM_SCALBL.Isend(sendbuf_x, 3 * sendCount_x, rank_x, sendtag + 0);
req2[0] =
Expand Down Expand Up @@ -2831,10 +2844,11 @@ void ScaLBL_Communicator::SendHalo(double *data) {
ScaLBL_Scalar_Pack(dvcSendList_yZ, sendCount_yZ, sendbuf_yZ, data, N);
ScaLBL_Scalar_Pack(dvcSendList_Yz, sendCount_Yz, sendbuf_Yz, data, N);
ScaLBL_Scalar_Pack(dvcSendList_YZ, sendCount_YZ, sendbuf_YZ, data, N);

//...................................................................................
// Send / Recv all the phase indcator field values
//...................................................................................

ScaLBL_DeviceBarrier();
req1[0] =
MPI_COMM_SCALBL.Isend(sendbuf_x, sendCount_x, rank_x, sendtag + 0);
req2[0] =
Expand Down
Loading