Commit b5c244f4 authored by Dave Moxey's avatar Dave Moxey
Browse files

Fix CompressibleFlowSolver with FR

parent 78557eb0
......@@ -1905,7 +1905,7 @@ namespace Nektar
for (j = 0; j < nquad1; ++j)
{
cnt = (nquad0 - 1) + j*nquad0 - i;
cnt = j*nquad0 + i;
divCFluxE2[cnt] = fluxJumps[i] * m_dGR_xi2[n][j];
}
}
......@@ -1917,7 +1917,7 @@ namespace Nektar
//fluxJumps[i] = -(m_Q2D_e3[n][i]) * fluxJumps[i];
for (j = 0; j < nquad0; ++j)
{
cnt = (nquad0*nquad1 - nquad0) + j - i*nquad0;
cnt = j + i*nquad0;
divCFluxE3[cnt] = fluxJumps[i] * m_dGL_xi1[n][j];
}
}
......@@ -2106,7 +2106,7 @@ namespace Nektar
for (j = 0; j < nquad1; ++j)
{
cnt = (nquad0 - 1) + j*nquad0 - i;
cnt = j*nquad0 + i;
divCFluxE2[cnt] = fluxJumps[i] * m_dGR_xi2[n][j];
}
}
......@@ -2118,7 +2118,7 @@ namespace Nektar
fluxJumps[i] = -(m_Q2D_e3[n][i]) * fluxJumps[i];
for (j = 0; j < nquad0; ++j)
{
cnt = (nquad0*nquad1 - nquad0) + j - i*nquad0;
cnt = j + i*nquad0;
divCFluxE3[cnt] = fluxJumps[i] * m_dGL_xi1[n][j];
}
......@@ -2408,7 +2408,7 @@ namespace Nektar
{
for (j = 0; j < nquad1; ++j)
{
cnt = (nquad0 - 1) + j*nquad0 - i;
cnt = j*nquad0 + i;
divCFluxE2[cnt] =
fluxJumps[i] * m_dGR_xi2[n][j];
}
......@@ -2474,8 +2474,7 @@ namespace Nektar
{
for (j = 0; j < nquad0; ++j)
{
cnt = (nquad0*nquad1 - nquad0) + j
- i*nquad0;
cnt = j + i*nquad0;
divCFluxE3[cnt] =
-fluxJumps[i] * m_dGL_xi1[n][j];
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment