#include "Inertia.h"

#include <cmath>

#include "model/Point.h"

Inertia::Inertia() {
    this->mass = 0;
    this->sx = 0;
    this->sxx = 0;
    this->sxy = 0;
    this->sy = 0;
    this->syy = 0;
}

Inertia::Inertia(const Inertia& inertia) { *this = inertia; }

Inertia::~Inertia() = default;

auto Inertia::centerX() const -> double { return this->sx / this->mass; }

auto Inertia::centerY() const -> double { return this->sy / this->mass; }

auto Inertia::xx() const -> double {
    if (this->mass <= 0.0) {
        return 0.0;
    }
    return (this->sxx - this->sx * this->sx / this->mass) / this->mass;
}

auto Inertia::xy() const -> double {
    if (this->mass <= 0.0) {
        return 0.0;
    }
    return (this->sxy - this->sx * this->sy / this->mass) / this->mass;
}

auto Inertia::yy() const -> double {
    if (this->mass <= 0.0) {
        return 0.0;
    }
    return (this->syy - this->sy * this->sy / this->mass) / this->mass;
}

auto Inertia::rad() const -> double {
    double ixx = this->xx();
    double iyy = this->yy();
    if (ixx + iyy <= 0.0) {
        return 0.0;
    }
    return sqrt(ixx + iyy);
}

auto Inertia::det() const -> double {
    double ixx = this->xx();
    double iyy = this->yy();
    double ixy = this->xy();
    if (this->mass <= 0.0) {
        return 0.0;
    }

    if (ixx + iyy <= 0.0) {
        return 0.0;
    }

    return 4 * (ixx * iyy - ixy * ixy) / (ixx + iyy) / (ixx + iyy);
}

auto Inertia::getMass() const -> double { return mass; }

void Inertia::increase(Point p1, Point p2, int coef) {
    double dm = coef * hypot(p2.x - p1.x, p2.y - p1.y);
    this->mass += dm;
    this->sx += dm * p1.x;
    this->sy += dm * p1.y;
    this->sxx += dm * p1.x * p1.x;
    this->syy += dm * p1.y * p1.y;
    this->sxy += dm * p1.x * p1.y;
}

void Inertia::calc(const Point* pt, int start, int end) {
    this->mass = this->sx = this->sy = this->sxx = this->sxy = this->syy = 0.;
    for (int i = start; i < end - 1; i++) {
        this->increase(pt[i], pt[i + 1], 1);
    }
}