Commit 796ba082 authored by Michael Turner's avatar Michael Turner

Merge branch 'fix/BLPoints' into 'master'

Fix boundary layer points

See merge request !892
parents 45ea4fe7 91da1d19
......@@ -76,7 +76,7 @@ v5.0.0
**IncNavierStokesSolver**
- Replace steady-state check based on difference of norms by check based on
norm of the difference, to be consistent with the compressible solver (!832)
- Updated SVV to allow for the DGKernel extension (!851)
- Updated SVV to allow for the DGKernel extension (!851)
**CompressibleFlowSolver**
- Add 3D regression tests (!567)
......@@ -101,6 +101,7 @@ v4.4.2
- Fix a few memory leaks in polylib (!863)
- Fix a crash when Interpolator is called on an empty field (!869)
- Fix petsc compile without MPI (!873)
- Fix calculation of BLPoints (!892)
- Fix deadlock in DiffusionLDG (!885)
- Fix uninitialised coefficients in DirectFull solver (!898)
......
......@@ -63,12 +63,12 @@ namespace Nektar
}
else
{
NekDouble a = 2.0 * (1.0-r) / (1.0 - pow(r,(double)npts));
NekDouble a = 2.0 * (1.0-r) / (1.0 - pow(r,(double)(npts-1)));
m_points[0][0] = -1.0;
for (unsigned int i = 1; i < npts; ++i)
{
m_points[0][i] = m_points[0][i-1] + a*pow(r,(double)i);
m_points[0][i] = m_points[0][i-1] + a*pow(r,(double)(i-1));
}
m_points[0][npts-1] = 1.0;
......
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