Position.h

Go to the documentation of this file.

/****************************************************************************/
// Eclipse SUMO, Simulation of Urban MObility; see https://eclipse.dev/sumo
// Copyright (C) 2001-2025 German Aerospace Center (DLR) and others.
// This program and the accompanying materials are made available under the
// terms of the Eclipse Public License 2.0 which is available at
// https://www.eclipse.org/legal/epl-2.0/
// This Source Code may also be made available under the following Secondary
// Licenses when the conditions for such availability set forth in the Eclipse
// Public License 2.0 are satisfied: GNU General Public License, version 2
// or later which is available at
// https://www.gnu.org/licenses/old-licenses/gpl-2.0-standalone.html
// SPDX-License-Identifier: EPL-2.0 OR GPL-2.0-or-later
/****************************************************************************/
// A position in the 2D- or 3D-world
/****************************************************************************/
#pragma once
#include <config.h>
#include <iostream>
#include <cmath>
 
#include <config.h>
 
// ===========================================================================
// class definitions
// ===========================================================================
class Position {
public:
    Position() :
        myX(0.0), myY(0.0), myZ(0.0) { }
 
    Position(double x, double y) :
        myX(x), myY(y), myZ(0) { }
 
    Position(double x, double y, double z) :
        myX(x), myY(y), myZ(z) { }
 
    inline double x() const {
        return myX;
    }
 
    inline double y() const {
        return myY;
    }
 
    inline double z() const {
        return myZ;
    }
 
    void setx(double x) {
        myX = x;
    }
 
    void sety(double y) {
        myY = y;
    }
 
    void setz(double z) {
        myZ = z;
    }
 
    void set(double x, double y) {
        myX = x;
        myY = y;
    }
 
    void set(double x, double y, double z) {
        myX = x;
        myY = y;
        myZ = z;
    }
 
    void set(const Position& pos) {
        myX = pos.myX;
        myY = pos.myY;
        myZ = pos.myZ;
    }
 
    void mul(double val) {
        myX *= val;
        myY *= val;
        myZ *= val;
    }
 
    void div(double val) {
        myX /= val;
        myY /= val;
        myZ /= val;
    }
 
    void mul(double mx, double my) {
        myX *= mx;
        myY *= my;
    }
 
    void mul(double mx, double my, double mz) {
        myX *= mx;
        myY *= my;
        myZ *= mz;
    }
 
    void add(const Position& pos) {
        myX += pos.myX;
        myY += pos.myY;
        myZ += pos.myZ;
    }
 
    void add(double dx, double dy) {
        myX += dx;
        myY += dy;
    }
 
    void add(double dx, double dy, double dz) {
        myX += dx;
        myY += dy;
        myZ += dz;
    }
 
    void sub(double dx, double dy) {
        myX -= dx;
        myY -= dy;
    }
 
    void sub(double dx, double dy, double dz) {
        myX -= dx;
        myY -= dy;
        myZ -= dz;
    }
 
    void sub(const Position& pos) {
        myX -= pos.myX;
        myY -= pos.myY;
        myZ -= pos.myZ;
    }
 
    inline double length() const {
        return sqrt(myX * myX + myY * myY + myZ * myZ);
    }
 
    inline double length2D() const {
        return sqrt(myX * myX + myY * myY);
    }
 
    inline void norm2D() {
        const double val = length2D();
        if (val != 0.) {
            myX /= val;
            myY /= val;
        }
    }
 
    friend std::ostream& operator<<(std::ostream& os, const Position& p) {
        os << p.x() << "," << p.y();
        if (p.z() != double(0.0)) {
            os << "," << p.z();
        }
        return os;
    }
 
    Position operator+(const Position& p2) const {
        return Position(myX + p2.myX,  myY + p2.myY, myZ + p2.myZ);
    }
 
    Position operator-(const Position& p2) const {
        return Position(myX - p2.myX,  myY - p2.myY, myZ - p2.myZ);
    }
 
    Position operator*(double scalar) const {
        return Position(myX * scalar, myY * scalar, myZ * scalar);
    }
 
    Position operator/(double scalar) const {
        return Position(myX / scalar, myY / scalar, myZ / scalar);
    }
 
    Position operator+(double offset) const {
        const double length = distanceTo(Position(0, 0, 0));
        if (length == 0) {
            return *this;
        }
        const double scalar = (length + offset) / length;
        return Position(myX * scalar, myY * scalar, myZ * scalar);
    }
 
    Position operator-(double offset) const {
        const double length = distanceTo(Position(0, 0, 0));
        if (length == 0) {
            return *this;
        }
        const double scalar = (length - offset) / length;
        return Position(myX * scalar, myY * scalar, myZ * scalar);
    }
 
    bool operator==(const Position& p2) const {
        return myX == p2.myX && myY == p2.myY && myZ == p2.myZ;
    }
 
    bool operator!=(const Position& p2) const {
        return myX != p2.myX || myY != p2.myY || myZ != p2.myZ;
    }
 
    bool operator<(const Position& p2) const {
        if (myX != p2.myX) {
            return myX < p2.myX;
        }
        if (myY != p2.myY) {
            return myY < p2.myY;
        }
        return myZ < p2.myZ;
    }
 
    bool almostSame(const Position& p2, double maxDiv = POSITION_EPS) const {
        return distanceTo(p2) < maxDiv;
    }
 
    inline double distanceTo(const Position& p2) const {
        return sqrt(distanceSquaredTo(p2));
    }
 
    inline double distanceSquaredTo(const Position& p2) const {
        return (myX - p2.myX) * (myX - p2.myX) + (myY - p2.myY) * (myY - p2.myY) + (myZ - p2.myZ) * (myZ - p2.myZ);
    }
 
    inline double distanceTo2D(const Position& p2) const {
        return sqrt(distanceSquaredTo2D(p2));
    }
 
    inline double distanceSquaredTo2D(const Position& p2) const {
        return (myX - p2.myX) * (myX - p2.myX) + (myY - p2.myY) * (myY - p2.myY);
    }
 
    inline double angleTo2D(const Position& other) const {
        return atan2(other.myY - myY, other.myX - myX);
    }
 
    inline double slopeTo2D(const Position& other) const {
        return atan2(other.myZ - myZ, distanceTo2D(other));
    }
 
    Position crossProduct(const Position& pos) {
        return Position(
                   myY * pos.myZ - myZ * pos.myY,
                   myZ * pos.myX - myX * pos.myZ,
                   myX * pos.myY - myY * pos.myX);
    }
 
    inline double dotProduct(const Position& pos) const {
        return myX * pos.myX + myY * pos.myY + myZ * pos.myZ;
    }
 
    Position rotateAround2D(double rad, const Position& origin);
 
    void swapXY() {
        std::swap(myX, myY);
    }
 
    bool isNAN() const {
        return (std::isnan(myX) || std::isnan(myY) || std::isnan(myZ));
    }
 
    void round(int precision);
 
 
    static const Position INVALID;
 
private:
    double myX;
 
    double myY;
 
    double myZ;
};