/*
This file is part of Magnum.
Original authors — credit is appreciated but not required:
2010, 2011, 2012, 2013, 2014, 2015, 2016, 2017, 2018, 2019,
2020, 2021, 2022, 2023, 2024 — Vladimír Vondruš <mosra@centrum.cz>
2020 — Nghia Truong <nghiatruong.vn@gmail.com>
This is free and unencumbered software released into the public domain.
Anyone is free to copy, modify, publish, use, compile, sell, or distribute
this software, either in source code form or as a compiled binary, for any
purpose, commercial or non-commercial, and by any means.
In jurisdictions that recognize copyright laws, the author or authors of
this software dedicate any and all copyright interest in the software to
the public domain. We make this dedication for the benefit of the public
at large and to the detriment of our heirs and successors. We intend this
dedication to be an overt act of relinquishment in perpetuity of all
present and future rights to this software under copyright law.
THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
THE AUTHORS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER
IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN
CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
*/
#include "ArcBall.h"
#include <Corrade/Utility/Assert.h>
#include <Magnum/Math/Matrix3.h>
namespace Magnum { namespace Examples {
namespace {
/* Project a point in NDC onto the arcball sphere */
Quaternion ndcToArcBall(const Vector2& p) {
const Float dist = Math::dot(p, p);
/* Point is on sphere */
if(dist <= 1.0f)
return {{p.x(), p.y(), Math::sqrt(1.0f - dist)}, 0.0f};
/* Point is outside sphere */
else
return {{p.normalized(), 0.0f}, 0.0f};
}
}
ArcBall::ArcBall(const Vector3& eye, const Vector3& viewCenter,
const Vector3& upDir, Deg fov, const Vector2i& windowSize):
_fov{fov}, _windowSize{windowSize}
{
setViewParameters(eye, viewCenter, upDir);
}
void ArcBall::setViewParameters(const Vector3& eye, const Vector3& viewCenter,
const Vector3& upDir)
{
const Vector3 dir = viewCenter - eye;
Vector3 zAxis = dir.normalized();
Vector3 xAxis = (Math::cross(zAxis, upDir.normalized())).normalized();
Vector3 yAxis = (Math::cross(xAxis, zAxis)).normalized();
xAxis = (Math::cross(zAxis, yAxis)).normalized();
_targetPosition = -viewCenter;
_targetZooming = -dir.length();
_targetQRotation = Quaternion::fromMatrix(
Matrix3x3{xAxis, yAxis, -zAxis}.transposed()).normalized();
_positionT0 = _currentPosition = _targetPosition;
_zoomingT0 = _currentZooming = _targetZooming;
_qRotationT0 = _currentQRotation = _targetQRotation;
updateInternalTransformations();
}
void ArcBall::reset() {
_targetPosition = _positionT0;
_targetZooming = _zoomingT0;
_targetQRotation = _qRotationT0;
}
void ArcBall::setLagging(const Float lagging) {
CORRADE_INTERNAL_ASSERT(lagging >= 0.0f && lagging < 1.0f);
_lagging = lagging;
}
void ArcBall::initTransformation(const Vector2& pointerPosition) {
_prevPointerPositionNdc = screenCoordToNdc(pointerPosition);
}
void ArcBall::rotate(const Vector2& pointerPosition) {
const Vector2 pointerPositionNdc = screenCoordToNdc(pointerPosition);
const Quaternion currentQRotation = ndcToArcBall(pointerPositionNdc);
const Quaternion prevQRotation = ndcToArcBall(_prevPointerPositionNdc);
_prevPointerPositionNdc = pointerPositionNdc;
_targetQRotation =
(currentQRotation*prevQRotation*_targetQRotation).normalized();
}
void ArcBall::translate(const Vector2& pointerPosition) {
const Vector2 mousePosNdc = screenCoordToNdc(pointerPosition);
const Vector2 translationNdc = mousePosNdc - _prevPointerPositionNdc;
_prevPointerPositionNdc = mousePosNdc;
translateDelta(translationNdc);
}
void ArcBall::translateDelta(const Vector2& translationNdc) {
/* Half size of the screen viewport at the view center and perpendicular
with the viewDir */
const Float hh = Math::abs(_targetZooming)*Math::tan(_fov*0.5f);
const Float hw = hh*Vector2{_windowSize}.aspectRatio();
_targetPosition += _inverseView.transformVector(
{translationNdc.x()*hw, translationNdc.y()*hh, 0.0f});
}
void ArcBall::zoom(const Float delta) {
_targetZooming += delta;
}
bool ArcBall::updateTransformation() {
const Vector3 diffViewCenter = _targetPosition - _currentPosition;
const Quaternion diffRotation = _targetQRotation - _currentQRotation;
const Float diffZooming = _targetZooming - _currentZooming;
const Float dViewCenter = Math::dot(diffViewCenter, diffViewCenter);
const Float dRotation = Math::dot(diffRotation, diffRotation);
const Float dZooming = diffZooming * diffZooming;
/* Nothing change */
if(dViewCenter < 1.0e-10f &&
dRotation < 1.0e-10f &&
dZooming < 1.0e-10f) {
return false;
}
/* Nearly done: just jump directly to the target */
if(dViewCenter < 1.0e-6f &&
dRotation < 1.0e-6f &&
dZooming < 1.0e-6f) {
_currentPosition = _targetPosition;
_currentQRotation = _targetQRotation;
_currentZooming = _targetZooming;
/* Interpolate between the current transformation and the target
transformation */
} else {
const Float t = 1 - _lagging;
_currentPosition = Math::lerp(_currentPosition, _targetPosition, t);
_currentZooming = Math::lerp(_currentZooming, _targetZooming, t);
_currentQRotation = Math::slerpShortestPath(
_currentQRotation, _targetQRotation, t);
}
updateInternalTransformations();
return true;
}
void ArcBall::updateInternalTransformations() {
_view = DualQuaternion::translation(Vector3::zAxis(_currentZooming))*
DualQuaternion{_currentQRotation}*
DualQuaternion::translation(_currentPosition);
_inverseView = _view.inverted();
}
Vector2 ArcBall::screenCoordToNdc(const Vector2& pointerPosition) const {
return Vector2::yScale(-1.0f)*
(pointerPosition*2.0f/Vector2{_windowSize} - Vector2{1.0f});
}
}}