mirror of
https://github.com/n5ac/mmtty.git
synced 2025-12-06 04:12:03 +01:00
83 lines
2.2 KiB
C++
83 lines
2.2 KiB
C++
|
//Copyright+LGPL
|
|||
|
|
|
|||
|
|
//-----------------------------------------------------------------------------------------------------------------------------------------------
|
|||
|
|
// Copyright 2000-2013 Makoto Mori, Nobuyuki Oba
|
|||
|
|
//-----------------------------------------------------------------------------------------------------------------------------------------------
|
|||
|
|
// This file is part of MMTTY.
|
|||
|
|
|
|||
|
|
// MMTTY 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.
|
|||
|
|
|
|||
|
|
// MMTTY 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 MMTTY. If not, see
|
|||
|
|
// <http://www.gnu.org/licenses/>.
|
|||
|
|
//-----------------------------------------------------------------------------------------------------------------------------------------------
|
|||
|
|
|
|||
|
|
|
|||
|
|
|
|||
|
|
//---------------------------------------------------------------------------
|
|||
|
|
#include <vcl.h>
|
|||
|
|
#pragma hdrstop
|
|||
|
|
|
|||
|
|
#include "math.h"
|
|||
|
|
#include "CLX.h"
|
|||
|
|
|
|||
|
|
///----------------------------------------------------------------
|
|||
|
|
/// <09><><EFBFBD>f<EFBFBD><66><EFBFBD><EFBFBD><EFBFBD>Z<EFBFBD>N<EFBFBD><4E><EFBFBD>X
|
|||
|
|
void CLX::Div(const CLX &ref)
|
|||
|
|
{
|
|||
|
|
double R, J, D;
|
|||
|
|
R = r * ref.r + j * ref.j;
|
|||
|
|
J = j * ref.r - r * ref.j;
|
|||
|
|
D = ref.r * ref.r + ref.j * ref.j;
|
|||
|
|
if( D ){
|
|||
|
|
r = R/D;
|
|||
|
|
j = J/D;
|
|||
|
|
}
|
|||
|
|
else {
|
|||
|
|
r = R >= 0 ? 1e100 : -1e100;
|
|||
|
|
j = J >= 0 ? 1e100 : -1e100;
|
|||
|
|
}
|
|||
|
|
}
|
|||
|
|
|
|||
|
|
CLX CLX::Sqrt(void)
|
|||
|
|
{
|
|||
|
|
CLX a;
|
|||
|
|
|
|||
|
|
double R, ST;
|
|||
|
|
R = sqrt(sqrt(r*r + j*j));
|
|||
|
|
ST = 0;
|
|||
|
|
if( r ) ST = atan2(j, r)/2.0;
|
|||
|
|
a.r = R*cos(ST);
|
|||
|
|
a.j = R*sin(ST);
|
|||
|
|
return a;
|
|||
|
|
}
|
|||
|
|
|
|||
|
|
CLX operator+(const double R, const CLX &ref)
|
|||
|
|
{
|
|||
|
|
CLX a(R);
|
|||
|
|
a+=ref;
|
|||
|
|
return a;
|
|||
|
|
}
|
|||
|
|
CLX operator-(const double R, const CLX &ref)
|
|||
|
|
{
|
|||
|
|
CLX a(R);
|
|||
|
|
a-=ref;
|
|||
|
|
return a;
|
|||
|
|
}
|
|||
|
|
CLX operator*(const double R, const CLX &ref)
|
|||
|
|
{
|
|||
|
|
CLX a(R);
|
|||
|
|
a*=ref;
|
|||
|
|
return a;
|
|||
|
|
}
|
|||
|
|
CLX operator/(const double R, const CLX &ref)
|
|||
|
|
{
|
|||
|
|
CLX a(R);
|
|||
|
|
a/=ref;
|
|||
|
|
return a;
|
|||
|
|
}
|
|||
|
|
#pragma package(smart_init)
|