- dxv = D[xv];
- dyv = D[yv];
- duv = D[uv];
- /* if (dxy <= dxu && dxy <= dxv && dxy <= dyu && dxy <= dyv && dxy <= duv) k=1; */
-/* else if (duv <= dxy && duv <= dxu && duv <= dxv && duv <= dyu && duv <= dyv) k=1; */
-/* else if (dxu <= dxy && dxu <= dxv && dxu <= dyu && dxu <= dyv && dxu <= duv) k=2; */
-/* else if (dyv <= dxy && dyv <= dxu && dyv <= dxv && dyv <= dyu && dyv <= duv) k=2; */
-/* else if (dxv <= dxy && dxv <= dxu && dxv <= dyu && dxv <= dyv && dxv <= duv) k=3; */
-/* else if (dyu <= dxy && dyu <= dxu && dyu <= dxv && dyu <= dyv && dyu <= duv) k=3; */
- //Rprintf("%d\t%d\t%d\t%d\t%d\t%d\t\n", xy, xu, xv, yu, yv, uv);
- //Rprintf("dxy = %f\tdxu = %f\tdyu = %f\n", dxy, dxu, dyu);
- //Rprintf("D[xv] = %f\tD[yv] = %f\tD[uv] = %f\n", D[xv], D[yv], D[uv]);
- /* switch (k) { */
-/* case 1 : A = dxv + dyu; B = dxu + dyv; C = dxy + duv; break; */
-/* case 2 : A = dxv + dyu; B = dxy + duv; C = dxu + dyv; break; */
-/* case 3 : A = dxu + dyv; B = dxy + duv; C = dxv + dyu; break; */
-/* } */
- //Rprintf("A = %f\tB = %f\tC = %f\n", A, B, C);