Offset::to_string () const
{
String s;
- s = String (" (") + ::to_string (coordinate_a_[X_AXIS]) + ", "
+ s = String (" (") + ::to_string (coordinate_a_[X_AXIS]) + ", "
+ ::to_string (coordinate_a_[Y_AXIS]) + ")";
return s;
}
#endif
-
bool
isinf_b (Real r)
{
/*
free bsd fix by John Galbraith
- */
-
+*/
+
Offset
complex_multiply (Offset z1, Offset z2)
{
Offset z;
if (!isinf_b (z2[Y_AXIS]))
- {
+ {
z[X_AXIS] = z1[X_AXIS] * z2[X_AXIS] - z1[Y_AXIS]*z2[Y_AXIS];
z[Y_AXIS] = z1[X_AXIS] * z2[Y_AXIS] + z1[Y_AXIS] * z2[X_AXIS];
- }
+ }
return z;
}
-
Offset
complex_conjugate (Offset o)
{
- o[Y_AXIS] = - o[Y_AXIS];
+ o[Y_AXIS] = -o[Y_AXIS];
return o;
}
{
z2 = complex_conjugate (z2);
Offset z = complex_multiply (z1, z2);
- z *= 1/z2.length ();
+ z *= 1 / z2.length ();
return z;
}
-
-
Offset
complex_exp (Offset o)
{
Real s = sin (o[Y_AXIS]);
Real c = cos (o[Y_AXIS]);
-
+
Real r = exp (o[X_AXIS]);
return Offset (r*c, r*s);
/**
euclidian vector length / complex modulus
- */
+*/
Real
Offset::length () const
{