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
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
|
/********************************************************************************
* *
* S i n g l e - P r e c i s i o n C o m p l e x N u m b e r *
* *
*********************************************************************************
* Copyright (C) 2006,2022 by Jeroen van der Zijp. All Rights Reserved. *
*********************************************************************************
* This library is free software; you can redistribute it and/or modify *
* it under the terms of the GNU Lesser General Public License as published by *
* the Free Software Foundation; either version 3 of the License, or *
* (at your option) any later version. *
* *
* This library is distributed in the hope that it will be useful, *
* but WITHOUT ANY WARRANTY; without even the implied warranty of *
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
* GNU Lesser General Public License for more details. *
* *
* You should have received a copy of the GNU Lesser General Public License *
* along with this program. If not, see <http://www.gnu.org/licenses/> *
********************************************************************************/
#include "xincs.h"
#include "fxver.h"
#include "fxdefs.h"
#include "fxmath.h"
#include "FXArray.h"
#include "FXHash.h"
#include "FXStream.h"
#include "FXObject.h"
#include "FXComplexf.h"
using namespace FX;
/*******************************************************************************/
namespace FX {
// Complex square root
FXComplexf csqrt(const FXComplexf& c){
FXfloat mag=abs(c);
FXfloat rr=Math::sqrt((mag+c.re)*0.5f);
FXfloat ii=Math::sqrt((mag-c.re)*0.5f);
return FXComplexf(rr,Math::copysign(ii,c.im));
}
// Complex sine
FXComplexf csin(const FXComplexf& c){
return FXComplexf(Math::sin(c.re)*Math::cosh(c.im),Math::cos(c.re)*Math::sinh(c.im));
}
// Complex cosine
FXComplexf ccos(const FXComplexf& c){
return FXComplexf(Math::cos(c.re)*Math::cosh(c.im),-Math::sin(c.re)*Math::sinh(c.im));
}
// Complex tangent
FXComplexf ctan(const FXComplexf& c){
FXComplexf em=exp(FXComplexf(c.im,-c.re));
FXComplexf ep=exp(FXComplexf(-c.im,c.re));
FXComplexf t=(em-ep)/(em+ep);
return FXComplexf(-t.im,t.re);
}
// Complex hyperbolic sine
FXComplexf csinh(const FXComplexf& c){
return FXComplexf(Math::cos(c.im)*Math::sinh(c.re),Math::sin(c.im)*Math::cosh(c.re));
}
// Complex hyperbolic cosine
FXComplexf ccosh(const FXComplexf& c){
return FXComplexf(Math::cos(c.im)*Math::cosh(c.re),Math::sin(c.im)*Math::sinh(c.re));
}
// Complex hyperbolic tangent
FXComplexf ctanh(const FXComplexf& c){
return csinh(c)/ccosh(c);
}
FXStream& operator<<(FXStream& store,const FXComplexf& c){
store << c.re << c.im;
return store;
}
FXStream& operator>>(FXStream& store,FXComplexf& c){
store >> c.re >> c.im;
return store;
}
}
|
