SegGeom.cpp interpolation bug
new code to be added in SegGeom.cpp
{
m_xmap->BwdTrans(m_coeffs[i], pts);
len += (pts[npts - 1] - pts[0]) * (pts[npts - 1] - pts[0]);
// OLD CODE -----
//xi += (coords[i] - pts[0]) * (coords[i] - pts[0]);
// ---------------
// NEW CODE ------
xi += (coords[i] - pts[0]) * (pts[npts-1] - pts[0]);
// ---------------
}
// OLD CODE ------
//len = sqrt(len);
//xi = sqrt(xi);
// ---------------
Lcoords[0] = 2 * xi / len - 1.0;
}
else