Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
#ifndef VECTOR_H_
#define VECTOR_H_
#include "BaseVector.h"
#include "QuantityVector.h"
#include <phys/units/quantity.hpp>
template <typename dim>
class Vector : public BaseVector<dim>
{
using Quantity = phys::units::quantity<dim, double>;
Vector(CoordinateSystem const& pCS, QuantityVector<dim> pQVector) :
BaseVector<Quantity>(pCS, pQVector)
{
}
public:
Vector(CoordinateSystem const& cs, Quantity x, Quantity y, Quantity z) :
BaseVector<dim>(cs, QuantityVector<dim>(x, y, z))
{
}
auto getComponents() const
{
return BaseVector<dim>::qVector;
}
auto getComponents(CoordinateSystem const& pCS) const
{
if (&pCS == BaseVector<dim>::cs)
{
return BaseVector<dim>::qVector;
}
else
{
return QuantityVector<dim>(CoordinateSystem::getTransformation(*BaseVector<dim>::cs, pCS).linear() * BaseVector<dim>::qVector.eVector);
}
}
void rebase(CoordinateSystem const& pCS)
{
BaseVector<dim>::qVector = getComponents(pCS);
BaseVector<dim>::cs = &pCS;
}
auto norm() const
{
return Quantity(BaseVector<dim>::qVector.eVector.norm());
}
//~ template <typename dim2>
//~ auto operator*(Vector<dim2> const& pVec)
//~ {
//~ auto constexpr resulting_dim = dimension
//~ return Vector<
//~ }
};
#endif