CoordinateSystemTrafo.h

Go to the documentation of this file.

#ifndef COORDINATESYSTEMTRAFO
#define COORDINATESYSTEMTRAFO
 
#include "Algorithms/Vektor.h"
#include "Algorithms/Quaternion.h"
#include "AppTypes/Tenzor.h"
 
class CoordinateSystemTrafo {
public:
    CoordinateSystemTrafo();
 
    CoordinateSystemTrafo(const CoordinateSystemTrafo &right);
 
    CoordinateSystemTrafo(const Vector_t &origin,
                          const Quaternion &orientation);
 
    Vector_t transformTo(const Vector_t &r) const;
    Vector_t transformFrom(const Vector_t &r) const;
 
    Vector_t rotateTo(const Vector_t &r) const;
    Vector_t rotateFrom(const Vector_t &r) const;
 
    void invert();
    CoordinateSystemTrafo inverted() const;
 
    CoordinateSystemTrafo& operator=(const CoordinateSystemTrafo& right) = default;
    CoordinateSystemTrafo operator*(const CoordinateSystemTrafo &right) const;
    void operator*=(const CoordinateSystemTrafo &right);
 
    Vector_t getOrigin() const;
    Quaternion getRotation() const;
 
    void print(std::ostream&) const;
private:
    Vector_t origin_m;
    Quaternion orientation_m;
    Tenzor<double, 3> rotationMatrix_m;
};
 
inline
std::ostream& operator<<(std::ostream& os, const CoordinateSystemTrafo &trafo) {
    trafo.print(os);
    return os;
}
 
inline
Inform& operator<<(Inform& os, const CoordinateSystemTrafo &trafo) {
    trafo.print(os.getStream());
    return os;
}
 
inline
void CoordinateSystemTrafo::print(std::ostream &os) const {
    os << "Origin: " << origin_m << "\n"
       << "z-axis: " << orientation_m.conjugate().rotate(Vector_t(0,0,1)) << "\n"
       << "x-axis: " << orientation_m.conjugate().rotate(Vector_t(1,0,0));
}
 
inline
Vector_t CoordinateSystemTrafo::getOrigin() const {
    return origin_m;
}
 
inline
Quaternion CoordinateSystemTrafo::getRotation() const {
    return orientation_m;
}
 
inline
CoordinateSystemTrafo CoordinateSystemTrafo::inverted() const {
    CoordinateSystemTrafo result(*this);
    result.invert();
 
    return result;
}
 
inline
Vector_t CoordinateSystemTrafo::transformTo(const Vector_t &r) const {
    return dot(rotationMatrix_m, r - origin_m);
}
 
inline
Vector_t CoordinateSystemTrafo::transformFrom(const Vector_t &r) const {
    return dot(transpose(rotationMatrix_m), r) + origin_m;
}
 
inline
Vector_t CoordinateSystemTrafo::rotateTo(const Vector_t &r) const {
    return dot(rotationMatrix_m, r);
}
 
inline
Vector_t CoordinateSystemTrafo::rotateFrom(const Vector_t &r) const {
    return dot(transpose(rotationMatrix_m), r);
}
 
#endif