#include <EvtComplex.hh>
|
|
|
|
|
|
|
00035 { 00036 00037 double r=_rpart*c._rpart-_ipart*c._ipart; 00038 double i=_rpart*c._ipart+_ipart*c._rpart; 00039 00040 _rpart=r; 00041 _ipart=i; 00042 00043 return *this; 00044 00045 }
|
|
|
|
|
|
00097 { 00098 00099 _rpart+=d; 00100 00101 return *this; 00102 }
|
|
|
|
00104 { 00105 00106 _rpart-=d; 00107 00108 return *this; 00109 }
|
|
|
|
00047 { 00048 00049 double inv=1.0/(c._rpart*c._rpart+c._ipart*c._ipart); 00050 00051 double r=inv*(_rpart*c._rpart+_ipart*c._ipart); 00052 double i=inv*(_rpart*c._ipart-_ipart*c._rpart); 00053 00054 _rpart=r; 00055 _ipart=i; 00056 00057 return *this; 00058 00059 }
|
|
|
|
|
|
|
|
00212 { 00213 00214 double c2=c._rpart*c._rpart+c._ipart*c._ipart; 00215 if (c2<=0.0) return 0.0; 00216 return sqrt(c2); 00217 00218 }
|
|
|
|
00227 { //added by FS 00228 if ((c._rpart==0)&&(c._ipart==0)) {return 0.0;} 00229 if (c._rpart==0){ 00230 if (c._ipart>0){ 00231 return EvtConst::pi/2; 00232 } else { 00233 return -EvtConst::pi/2; 00234 } 00235 } else { 00236 return atan2(c._ipart,c._rpart); 00237 } 00238 }
|
|
00206 { 00207 00208 return EvtComplex(c._rpart,-c._ipart); 00209 00210 }
|
|
|
|
00246 { 00247 00248 return c._ipart; 00249 00250 }
|
|
00168 { 00169 00170 return EvtComplex(c1._rpart*c2._rpart-c1._ipart*c2._ipart, 00171 c1._rpart*c2._ipart+c1._ipart*c2._rpart); 00172 00173 }
|
|
00117 { 00118 00119 return EvtComplex(c._rpart*d,c._ipart*d); 00120 00121 }
|
|
00111 { 00112 00113 return EvtComplex(c._rpart*d,c._ipart*d); 00114 00115 }
|
|
00181 { 00182 00183 return EvtComplex(c1._rpart+c2._rpart,c1._ipart+c2._ipart); 00184 00185 }
|
|
00200 { 00201 00202 return EvtComplex(-c._rpart,-c._ipart); 00203 00204 }
|
|
00175 { 00176 00177 return EvtComplex(c1._rpart-c2._rpart,c1._ipart-c2._ipart); 00178 00179 }
|
|
00159 { 00160 00161 double inv=1.0/(c2._rpart*c2._rpart+c2._ipart*c2._ipart); 00162 00163 return EvtComplex(inv*(c1._rpart*c2._rpart+c1._ipart*c2._ipart), 00164 inv*(c1._ipart*c2._rpart-c1._rpart*c2._ipart)); 00165 00166 }
|
|
00149 { 00150 00151 double Num=d/(c._rpart*c._rpart+c._ipart*c._ipart); 00152 00153 return EvtComplex( Num*c._rpart, -Num*c._ipart ); 00154 00155 }
|
|
00125 { 00126 00127 return EvtComplex(c._rpart/d,c._ipart/d); 00128 00129 }
|
|
|
|
00240 { 00241 00242 return c._rpart; 00243 00244 }
|
|
|
|
|