Here's a bug fix for line 258:
Code: Select all
if (bad_vertex < 4)
{
y[bad_vertex - 1] = y[4];
p[bad_vertex - 1] = p[4];
q[bad_vertex - 1] = q[4];
}
Code: Select all
if (bad_vertex < 3)
{
y[bad_vertex - 1] = y[3];
p[bad_vertex - 1] = p[3];
q[bad_vertex - 1] = q[3];
}