geometry/Point.cpp

• View this file on GitHub
• Last update: 2024-08-09 13:08:28+00:00

Required by

• verify/geometry/igor-tests/01.cpp
• verify/geometry/igor-tests/02.cpp
• verify/geometry/igor-tests/03.cpp
• verify/geometry/igor-tests/04.cpp
• verify/geometry/igor-tests/05.cpp
• verify/geometry/igor-tests/06.cpp
• verify/geometry/igor-tests/07.cpp
• verify/geometry/igor-tests/08.cpp
• verify/geometry/igor-tests/09.cpp
• verify/geometry/igor-tests/10.cpp
• verify/geometry/igor-tests/11.cpp
• verify/geometry/igor-tests/12.cpp
• verify/geometry/igor-tests/13.cpp
• verify/geometry/igor-tests/14.cpp
• verify/geometry/igor-tests/15.cpp
• verify/geometry/igor-tests/16.cpp
• verify/geometry/igor-tests/17.cpp
• verify/geometry/igor-tests/20.cpp

Verified with

• verify/geometry/aoj-cgl-1-a.test.cpp
• verify/geometry/aoj-cgl-1-b.test.cpp
• verify/geometry/aoj-cgl-1-c.test.cpp
• verify/geometry/aoj-cgl-2-a.test.cpp
• verify/geometry/aoj-cgl-2-b.test.cpp
• verify/geometry/aoj-cgl-2-c.test.cpp
• verify/geometry/aoj-cgl-2-d.test.cpp
• verify/geometry/aoj-cgl-3-a.test.cpp
• verify/geometry/aoj-cgl-3-b.test.cpp
• verify/geometry/aoj-cgl-3-c.test.cpp
• verify/geometry/aoj-cgl-4-a.test.cpp
• verify/geometry/aoj-cgl-4-b.test.cpp
• verify/geometry/aoj-cgl-4-c.test.cpp
• verify/geometry/aoj/aoj-cgl-1-a.test.cpp
• verify/geometry/aoj/aoj-cgl-1-b.test.cpp
• verify/geometry/aoj/aoj-cgl-1-c.test.cpp
• verify/geometry/aoj/aoj-cgl-2-a.test.cpp
• verify/geometry/aoj/aoj-cgl-2-b.test.cpp
• verify/geometry/aoj/aoj-cgl-2-c.test.cpp
• verify/geometry/aoj/aoj-cgl-2-d.test.cpp
• verify/geometry/aoj/aoj-cgl-3-a.test.cpp
• verify/geometry/aoj/aoj-cgl-3-b.test.cpp
• verify/geometry/aoj/aoj-cgl-3-c.test.cpp
• verify/geometry/aoj/aoj-cgl-4-a.test.cpp
• verify/geometry/aoj/aoj-cgl-4-b.test.cpp
• verify/geometry/aoj/aoj-cgl-4-c.test.cpp
• verify/geometry/igor-tests/18.test.cpp
• verify/geometry/igor-tests/19.test.cpp
• verify/geometry/igor-tests/include-all.test.cpp

Code

/**
 * Author: alexxela12345,daubi,talant
 * Date: 2024-08-03
 * Description: struct Point
 */

const ld EPS = 1e-7;

ld sq(ld x) {
    return x * x;
}

int sign(ld x) {
    if (x < -EPS) {
        return -1;
    }
    if (x > EPS) {
        return 1;
    }
    return 0;
}

#define vec point
struct point {//% - cross, * - dot
    ld x, y;
    auto operator<=>(const point&) const = default;
};
ld operator*(const point &a, const point &b) {
    return a.x * b.x + a.y * b.y;
}
ld operator%(const point &a, const point &b) {
    return a.x * b.y - a.y * b.x;
}
point operator-(const point &a, const point &b) {
    return {a.x - b.x, a.y - b.y};
}
point operator+(const point &a, const point &b) {
    return {a.x + b.x, a.y + b.y};
}
point operator*(const point &a, ld b) {
    return {a.x * b, a.y * b};
}
point operator/(const point &a, ld b) {
    return {a.x / b, a.y / b};
}
bool operator<(const point &a, const point &b)  {
    if (sign(a.y - b.y) != 0) {
        return a.y < b.y;
    } else if (sign(a.x - b.x) != 0) {
        return a.x < b.x;
    }
    return 0;
}
ld len2(const point &a) {
    return sq(a.x) + sq(a.y);
}
ld len(const point &a) {
    return sqrt(len2(a));
}
point norm(point a) {
    return a / len(a);
}
int half(point a) {
    return (sign(a.y) == -1 || (sign(a.y) ==0 && a.x < 0));
}
point ort(point a) {
    return {-a.y, a.x};
}
point turn(point a, ld ang) {
    return {a.x * cos(ang) - a.y * sin(ang), a.x * sin(ang) + a.y * cos(ang)};
}
ld getAngle(point &a, point &b) {
    return atan2(a % b, a * b);
}
bool cmpHalf(const point &a, const point &b) {
    if (half(a) != half(b)) {
        return half(b);
    } else {
        int sgn = sign(a % b);
        if (!sgn) {
            return len2(a) < len2(b);
        } else {
            return sgn == 1;
        }
    }
}

#line 1 "geometry/Point.cpp"
/**
 * Author: alexxela12345,daubi,talant
 * Date: 2024-08-03
 * Description: struct Point
 */

const ld EPS = 1e-7;

ld sq(ld x) {
    return x * x;
}

int sign(ld x) {
    if (x < -EPS) {
        return -1;
    }
    if (x > EPS) {
        return 1;
    }
    return 0;
}

#define vec point
struct point {//% - cross, * - dot
    ld x, y;
    auto operator<=>(const point&) const = default;
};
ld operator*(const point &a, const point &b) {
    return a.x * b.x + a.y * b.y;
}
ld operator%(const point &a, const point &b) {
    return a.x * b.y - a.y * b.x;
}
point operator-(const point &a, const point &b) {
    return {a.x - b.x, a.y - b.y};
}
point operator+(const point &a, const point &b) {
    return {a.x + b.x, a.y + b.y};
}
point operator*(const point &a, ld b) {
    return {a.x * b, a.y * b};
}
point operator/(const point &a, ld b) {
    return {a.x / b, a.y / b};
}
bool operator<(const point &a, const point &b)  {
    if (sign(a.y - b.y) != 0) {
        return a.y < b.y;
    } else if (sign(a.x - b.x) != 0) {
        return a.x < b.x;
    }
    return 0;
}
ld len2(const point &a) {
    return sq(a.x) + sq(a.y);
}
ld len(const point &a) {
    return sqrt(len2(a));
}
point norm(point a) {
    return a / len(a);
}
int half(point a) {
    return (sign(a.y) == -1 || (sign(a.y) ==0 && a.x < 0));
}
point ort(point a) {
    return {-a.y, a.x};
}
point turn(point a, ld ang) {
    return {a.x * cos(ang) - a.y * sin(ang), a.x * sin(ang) + a.y * cos(ang)};
}
ld getAngle(point &a, point &b) {
    return atan2(a % b, a * b);
}
bool cmpHalf(const point &a, const point &b) {
    if (half(a) != half(b)) {
        return half(b);
    } else {
        int sgn = sign(a % b);
        if (!sgn) {
            return len2(a) < len2(b);
        } else {
            return sgn == 1;
        }
    }
}

Back to top page