Commit 9be305ac authored by vidhi's avatar vidhi

specify reason for exact comparison

parent 23d5dd81
......@@ -1601,6 +1601,13 @@ namespace Nektar
NekDouble xdiff = z[i] - coord;
NekDouble pval = physvals[i];
/*
* (in this specific case) you actually
* want to do the comparison exactly
* (believe it or not!) See chapter 7 of
* the paper here:
*https://people.maths.ox.ac.uk/trefethen/barycentric.pdf
*/
if (xdiff == 0.0)
{
return pval;
......
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