Попадание точки в фигуру на плоскости

#define _CRT_SECURE_NO_WARNINGS
 
#include <cstdio>
#include <cstdlib>
 
int main()
{
    double spo; // Переменные специально разграничил, чтобы было удобнее понять...
    double r; // задаем радиус;
    printf("Enter radius of the circle: ");
    scanf("%lf", &r);
    double px, py; // задаем точку;
    // далее идет список вспомогательных переменных, для удобства я их разбил на принадлежность;
    double ox = r, oy = -r; // задаем центр окружности, границы четвертей берем из вершин ромба!!!;
    // определяем координаты вершин ромба и 4 переменные произведений векторов и точки. Первой вершиной выбрал левую;
    double rx1 = 0, ry1 = -r, rx2 = r, ry2 = -2 * r, rx3 = 2 * r, ry3 = -r, rx4 = r, ry4 = 0, tmp1, tmp2, tmp3, tmp4;
    //
    double sx1 = 0, sy1 = 0, sx2 = 0, sy2 = -2 * r, sx3 = 2 * r, sy3 = -2 * r, sx4 = 2 * r, sy4 = 0; // задаем вершины описанного квадрата;
    
    printf("Enter coordinates of the point [X Y]: ");
    scanf("%lf %lf", &px, &py);
    
    // См. 1 пункт алгоритма -идем против часовой стрелки - с точкой сравниваются последовательно вектора RX2RY2 RX1RY1 --> RX3RY3 RX2RY2 и т.д.;
    tmp1 = (rx2 - rx1)*(py - ry1) - (ry2 - ry1)*(px - rx1);
    tmp2 = (rx3 - rx2)*(py - ry2) - (ry3 - ry2)*(px - rx2);
    tmp3 = (rx4 - rx3)*(py - ry3) - (ry4 - ry3)*(px - rx3);
    tmp4 = (rx1 - rx4)*(py - ry4) - (ry1 - ry4)*(px - rx4);
 
    if ((tmp1 >= 0 && tmp2 >= 0 && tmp3 >= 0 && tmp4 >= 0) || (tmp1 <= 0 && tmp2 <= 0 && tmp3 <= 0 && tmp4))
    {
        printf("The point is inside the white zone\n");
        spo = px * py;
        printf("SPO = x * y = %.3lf\n", spo);
    }
    // См. 2 пункт алгоритма;
    else if ((r*r >= (oy - py) * (oy - py) + (ox - px) * (ox - px)) && ((px >= r && py >= -r) || (px <= r && py <= -r)))
    {
        printf("The point is inside the red zone\n");
        spo = (px + 10) / (px + py);
        printf("SPO = (x + 10) / (x + y) = %.3lf\n", spo);
    }
    else if ((r*r >= (oy - py) * (oy - py) + (ox - px) * (ox - px)) && ((px < r && py > -r) || (px > r && py < -r)))
    {
        printf("The point is inside the blue zone\n");
        spo = px + (1 / py);
        printf("SPO = x + (1 / y) = %.3lf\n", spo);
    }
    // См. 3 пункт алгоритма;
    else if ((r*r < (oy - py) * (oy - py) + (ox - px) * (ox - px)) && px < r && py > -r && px >= 0 && py <= 0)
    {
        printf("The point is inside the white zone\n");
        spo = px * py;
        printf("SPO = x * y = %.3lf\n", spo);
    }
    // См. 4 пункт алгоритма;
    else
    {
        printf("The point is inside the yellow zone\n");
        spo = px / (py - px);
        printf("SPO = x  / (y - x) = %.3lf\n", spo);
    }
 
    system("pause");
    return 0;
}

void my_fill(HWND hWnd, POINTS MP, COLORREF fill_c)
{
    HDC hdc = GetDC(hWnd);
    COLORREF MPc = GetPixel(hdc, MP.x, MP.y);
    if (MPc == fill_c) return;
 
    const unsigned int Pbufsize = 1000;
    POINTS P0[Pbufsize], *pf_beg = P0, *pf_end = P0, P1[Pbufsize], *pt_beg = P1;
    *pf_beg = MP;
    while (pf_end >= pf_beg) {
        POINTS *pf, *pt = pt_beg;
        for (pf = pf_beg; pf <= pf_end; pf++) {
            for (int dx = -1; dx <= 1; dx++)
            for (int dy = -1; dy <= 1; dy++) {
                POINTS P = *pf;
                P.x += dx;
                P.y += dy;
                COLORREF Pc = GetPixel(hdc, P.x, P.y);
                if (Pc == MPc) {
                    SetPixel(hdc, P.x, P.y, fill_c);
                    *pt++ = P;
                }
            }
        }
        pf = pf_beg;
        pf_beg = pt_beg;
        pf_end = pt - 1;
        pt_beg = pf;
    }
}

if ((tmp1 >= 0 && tmp2 >= 0 && tmp3 >= 0 && tmp4 >= 0) || (tmp1 <= 0 && tmp2 <= 0 && tmp3 <= 0 && tmp4 <= 0))