| 1 | /**************************************************************************** |
| 2 | ** |
| 3 | ** Copyright (C) 2016 The Qt Company Ltd. |
| 4 | ** Contact: https://www.qt.io/licensing/ |
| 5 | ** |
| 6 | ** This file is part of the QtPositioning module of the Qt Toolkit. |
| 7 | ** |
| 8 | ** $QT_BEGIN_LICENSE:LGPL$ |
| 9 | ** Commercial License Usage |
| 10 | ** Licensees holding valid commercial Qt licenses may use this file in |
| 11 | ** accordance with the commercial license agreement provided with the |
| 12 | ** Software or, alternatively, in accordance with the terms contained in |
| 13 | ** a written agreement between you and The Qt Company. For licensing terms |
| 14 | ** and conditions see https://www.qt.io/terms-conditions. For further |
| 15 | ** information use the contact form at https://www.qt.io/contact-us. |
| 16 | ** |
| 17 | ** GNU Lesser General Public License Usage |
| 18 | ** Alternatively, this file may be used under the terms of the GNU Lesser |
| 19 | ** General Public License version 3 as published by the Free Software |
| 20 | ** Foundation and appearing in the file LICENSE.LGPL3 included in the |
| 21 | ** packaging of this file. Please review the following information to |
| 22 | ** ensure the GNU Lesser General Public License version 3 requirements |
| 23 | ** will be met: https://www.gnu.org/licenses/lgpl-3.0.html. |
| 24 | ** |
| 25 | ** GNU General Public License Usage |
| 26 | ** Alternatively, this file may be used under the terms of the GNU |
| 27 | ** General Public License version 2.0 or (at your option) the GNU General |
| 28 | ** Public license version 3 or any later version approved by the KDE Free |
| 29 | ** Qt Foundation. The licenses are as published by the Free Software |
| 30 | ** Foundation and appearing in the file LICENSE.GPL2 and LICENSE.GPL3 |
| 31 | ** included in the packaging of this file. Please review the following |
| 32 | ** information to ensure the GNU General Public License requirements will |
| 33 | ** be met: https://www.gnu.org/licenses/gpl-2.0.html and |
| 34 | ** https://www.gnu.org/licenses/gpl-3.0.html. |
| 35 | ** |
| 36 | ** $QT_END_LICENSE$ |
| 37 | ** |
| 38 | ****************************************************************************/ |
| 39 | |
| 40 | #include "qdoublevector3d_p.h" |
| 41 | #include <QtCore/qdatastream.h> |
| 42 | #include <QtCore/qmath.h> |
| 43 | #include <QtCore/qdebug.h> |
| 44 | |
| 45 | QT_BEGIN_NAMESPACE |
| 46 | |
| 47 | QDoubleVector3D QDoubleVector3D::normalized() const |
| 48 | { |
| 49 | // Need some extra precision if the length is very small. |
| 50 | double len = double(xp) * double(xp) + |
| 51 | double(yp) * double(yp) + |
| 52 | double(zp) * double(zp); |
| 53 | if (qFuzzyIsNull(d: len - 1.0)) |
| 54 | return *this; |
| 55 | else if (!qFuzzyIsNull(d: len)) |
| 56 | return *this / (double)qSqrt(v: len); |
| 57 | else |
| 58 | return QDoubleVector3D(); |
| 59 | } |
| 60 | |
| 61 | void QDoubleVector3D::normalize() |
| 62 | { |
| 63 | // Need some extra precision if the length is very small. |
| 64 | double len = double(xp) * double(xp) + |
| 65 | double(yp) * double(yp) + |
| 66 | double(zp) * double(zp); |
| 67 | if (qFuzzyIsNull(d: len - 1.0) || qFuzzyIsNull(d: len)) |
| 68 | return; |
| 69 | |
| 70 | len = qSqrt(v: len); |
| 71 | |
| 72 | xp /= len; |
| 73 | yp /= len; |
| 74 | zp /= len; |
| 75 | } |
| 76 | |
| 77 | QDoubleVector3D QDoubleVector3D::normal(const QDoubleVector3D &v1, const QDoubleVector3D &v2) |
| 78 | { |
| 79 | return crossProduct(v1, v2).normalized(); |
| 80 | } |
| 81 | |
| 82 | QDoubleVector3D QDoubleVector3D::normal |
| 83 | (const QDoubleVector3D &v1, const QDoubleVector3D &v2, const QDoubleVector3D &v3) |
| 84 | { |
| 85 | return crossProduct(v1: (v2 - v1), v2: (v3 - v1)).normalized(); |
| 86 | } |
| 87 | |
| 88 | double QDoubleVector3D::distanceToPlane |
| 89 | (const QDoubleVector3D &plane1, const QDoubleVector3D &plane2, const QDoubleVector3D &plane3) const |
| 90 | { |
| 91 | QDoubleVector3D n = normal(v1: plane2 - plane1, v2: plane3 - plane1); |
| 92 | return dotProduct(v1: *this - plane1, v2: n); |
| 93 | } |
| 94 | |
| 95 | double QDoubleVector3D::distanceToLine |
| 96 | (const QDoubleVector3D &point, const QDoubleVector3D &direction) const |
| 97 | { |
| 98 | if (direction.isNull()) |
| 99 | return (*this - point).length(); |
| 100 | QDoubleVector3D p = point + dotProduct(v1: *this - point, v2: direction) * direction; |
| 101 | return (*this - p).length(); |
| 102 | } |
| 103 | |
| 104 | double QDoubleVector3D::length() const |
| 105 | { |
| 106 | return qSqrt(v: xp * xp + yp * yp + zp * zp); |
| 107 | } |
| 108 | |
| 109 | #ifndef QT_NO_DEBUG_STREAM |
| 110 | |
| 111 | QDebug operator<<(QDebug dbg, const QDoubleVector3D &vector) |
| 112 | { |
| 113 | QDebugStateSaver saver(dbg); |
| 114 | dbg.nospace() << "QDoubleVector3D(" |
| 115 | << vector.x() << ", " << vector.y() << ", " << vector.z() << ')'; |
| 116 | return dbg; |
| 117 | } |
| 118 | |
| 119 | #endif |
| 120 | |
| 121 | #ifndef QT_NO_DATASTREAM |
| 122 | |
| 123 | QDataStream &operator<<(QDataStream &stream, const QDoubleVector3D &vector) |
| 124 | { |
| 125 | stream << double(vector.x()) << double(vector.y()) |
| 126 | << double(vector.z()); |
| 127 | return stream; |
| 128 | } |
| 129 | |
| 130 | QDataStream &operator>>(QDataStream &stream, QDoubleVector3D &vector) |
| 131 | { |
| 132 | double x, y, z; |
| 133 | stream >> x; |
| 134 | stream >> y; |
| 135 | stream >> z; |
| 136 | vector.setX(double(x)); |
| 137 | vector.setY(double(y)); |
| 138 | vector.setZ(double(z)); |
| 139 | return stream; |
| 140 | } |
| 141 | |
| 142 | #endif // QT_NO_DATASTREAM |
| 143 | |
| 144 | QT_END_NAMESPACE |
| 145 | |