/// <inheritdoc /> protected override void OnForward(double[] lp, double[] xy, int startIndex, int numPoints) { for (int i = startIndex; i < startIndex + numPoints; i++) { int phi = i * 2 + PHI; int lam = i * 2 + LAMBDA; int x = i * 2 + X; int y = i * 2 + Y; double[] p = new double[2]; lp[phi] = (lp[phi] - Phi0) * RAD_TO_SEC5; p[R] = _tpsi[NTPSI]; for (int j = NTPSI - 1; j >= 0; j--) { p[R] = _tpsi[j] + lp[phi] * p[R]; } p[R] *= lp[phi]; p[I] = lp[lam]; p = Proj.Zpoly1(p, _bf, NBF); xy[x] = p[I]; xy[y] = p[R]; } }