Classes
class	BarnesHut
	Класс, определяющий основной алгоритм модификации метода Барнса - Хата More...

class	MortonTree
	Структура, соответствующая мортоновскому дереву More...

class	params
	Класс, содержащий параметры метода Баонса - Хата для CPU. More...

class	PointsCopy
	Класс-обертка для точек и панелей More...

struct	TParticleCode
	Структура, соответствующая частице и ее мортоновскому коду More...

struct	treeCellT
	Структура, соответствующая ячейке мортоновского дерева More...

Functions
void	SizeCheck (std::vector< Point2D > &i00)
	Вспомогательная функция корректировки capacity вектора (при необходимости - удваивает)

Point2D	powz (const Point2D &z, double n)
	Умножение комплексных чисел

template<typename T >
T	sqr (T x)
	Умножение a на комплексно сопряженноe к b.

template<typename T >
int	sign (T val)
	Шаблонная функция знака числа Написана оптимизированная версия, которая работает самым быстрым возможным образом, т.к. не содержит ни одного оператора условного перехода

int	ceilpow2 (unsigned int x, unsigned int p)
	Округление "в потолок" результата деления на степень двойки, эквивалент ceil(x / (2^p))

int	ceilhalf (unsigned int x)
	Округление "в потолок" результата деления пополам, эквивалент ceil(x / 2)

template<typename T >
double	norm (const T &b)
	Шаблонная функция вычисления евклидовой нормы вектора или списка

template<typename T >
std::vector< T >	operator+ (const std::vector< T > &x, const std::vector< T > &y)
	Шаблонная функция сложения двух векторов

template<typename T >
std::vector< T > &	operator+= (std::vector< T > &x, const std::vector< T > &y)
	Шаблонная функция прибавления к одному вектору другого

template<typename T >
std::vector< T >	operator- (const std::vector< T > &x, const std::vector< T > &y)
	Шаблонная функция вычитания векторов

template<typename T >
std::vector< T > &	operator-= (std::vector< T > &x, const std::vector< T > &y)
	Шаблонная функция вычитания из одного вектора другого

template<typename T >
std::vector< T >	operator* (const T lambda, const std::vector< T > &x)
	Шаблонная функция умножения числа на вектор

template<typename T >
T	operator& (const std::vector< T > &x, const std::vector< T > &y)
	Шаблонная функция вычисления скалярного произведения двух векторов

Variables
long long	op
	Глобальная переменная - счетчик количества операций

static const double	PI = 3.1415926535897932384626

static const double	DPI = 2.0 * 3.1415926535897932384626

static const double	IPI = 1.0 / PI

static const double	IDPI = 0.5 / PI

static const int	twoPowCodeLengthVar = (1 << codeLength)
	2 в степени длины мортоновского кода (на каждую координату)

Function Documentation

◆ ceilhalf()

int BH::ceilhalf ( unsigned int x )

inline

Округление "в потолок" результата деления пополам, эквивалент ceil(x / 2)

Definition at line 134 of file defsBH.h.

    {
        return (x >> 1) + (x & 1);
    }

Here is the caller graph for this function:

◆ ceilpow2()

int BH::ceilpow2	(	unsigned int	x,
		unsigned int	p
	)

inline

Округление "в потолок" результата деления на степень двойки, эквивалент ceil(x / (2^p))

Parameters

[in]	x	делимое
[in]	p	показатель степени двойки в делителе

Definition at line 128 of file defsBH.h.

    {
        return (x >> p) + !!(x & ((1 << p) - 1));
    }

Here is the caller graph for this function:

◆ norm()

template<typename T >

double BH::norm ( const T & b )

inline

Шаблонная функция вычисления евклидовой нормы вектора или списка

Parameters

[in] b константная ссылка на вектор или список

Definition at line 143 of file defsBH.h.

    {
        double norm = 0;
#ifndef OLD_OMP
#pragma omp simd reduction(+:norm)
#endif
            for (size_t i = 0; i < b.size(); i++)
                norm += (b[i] * b[i]);
            ADDOP(2*b.size() + 1);
            return sqrt(norm);
    }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ operator&()

template<typename T >

T BH::operator&	(	const std::vector< T > &	x,
		const std::vector< T > &	y
	)

inline

Шаблонная функция вычисления скалярного произведения двух векторов

Definition at line 225 of file defsBH.h.

    {
        T c = 0;
#ifndef OLD_OMP
#pragma omp simd reduction(+:c)
#endif
        for (size_t i = 0; i < x.size(); ++i)
            c += x[i] * y[i];
        ADDOP(x.size());
        return c;
    }

◆ operator*()

template<typename T >

std::vector< T > BH::operator*	(	const T	lambda,
		const std::vector< T > &	x
	)

inline

Шаблонная функция умножения числа на вектор

Definition at line 208 of file defsBH.h.

    {
        std::vector<T> c(x);
        c.resize(x.size());
    
#ifndef OLD_OMP
#pragma omp simd
#endif
        for (size_t i = 0; i < x.size(); ++i)
            c[i] *= lambda;
        ADDOP(x.size());
        return c;
    }

◆ operator+()

template<typename T >

std::vector< T > BH::operator+	(	const std::vector< T > &	x,
		const std::vector< T > &	y
	)

inline

Шаблонная функция сложения двух векторов

Definition at line 157 of file defsBH.h.

    {
        std::vector<T> c(x);
#ifndef OLD_OMP
#pragma omp simd
#endif
        for (size_t i = 0; i < x.size(); ++i)
            c[i] += y[i];
        return c;
    }

◆ operator+=()

template<typename T >

std::vector< T > & BH::operator+=	(	std::vector< T > &	x,
		const std::vector< T > &	y
	)

inline

Шаблонная функция прибавления к одному вектору другого

Definition at line 171 of file defsBH.h.

    {
#ifndef OLD_OMP
#pragma omp simd
#endif
        for (size_t i = 0; i < x.size(); ++i)
            x[i] += y[i];
        return x;
    }

◆ operator-()

template<typename T >

std::vector< T > BH::operator-	(	const std::vector< T > &	x,
		const std::vector< T > &	y
	)

inline

Шаблонная функция вычитания векторов

Definition at line 183 of file defsBH.h.

    {
        std::vector<T> c(x);
#ifndef OLD_OMP
#pragma omp simd
#endif
        for (size_t i = 0; i < x.size(); ++i)
            c[i] -= y[i];
        return c;
    }

◆ operator-=()

template<typename T >

std::vector< T > & BH::operator-=	(	std::vector< T > &	x,
		const std::vector< T > &	y
	)

inline

Шаблонная функция вычитания из одного вектора другого

Definition at line 196 of file defsBH.h.

    {
#ifndef OLD_OMP
#pragma omp simd
#endif
        for (size_t i = 0; i < x.size(); ++i)
            x[i] -= y[i];
        return x;
    }

◆ powz()

Point2D BH::powz	(	const Point2D &	z,
		double	n
	)

inline

Умножение комплексных чисел

Возведение в степень комплексных чисел

Definition at line 83 of file defsBH.h.

    {
        double phi, R;
        ADDOP(10);
        phi = n * atan2(z[1], z[0]);
        R = pow(z.length2(), 0.5*n);
        return Point2D({ R * cos(phi), R * sin(phi) });
    }

Here is the call graph for this function:

◆ sign()

template<typename T >

int BH::sign ( T val )

inline

Шаблонная функция знака числа Написана оптимизированная версия, которая работает самым быстрым возможным образом, т.к. не содержит ни одного оператора условного перехода

Definition at line 111 of file defsBH.h.

    {
        return (T(0) < val) - (val < T(0));
    }

Here is the caller graph for this function:

◆ SizeCheck()

void BH::SizeCheck ( std::vector< Point2D > & i00 )

inline

Вспомогательная функция корректировки capacity вектора (при необходимости - удваивает)

Definition at line 68 of file defsBH.h.

    {
        if (i00.capacity() == i00.size())
            i00.reserve(i00.size() * 2);
    }

◆ sqr()

template<typename T >

T BH::sqr ( T x )

inline

Умножение a на комплексно сопряженноe к b.

Шаблонная функция возведения в квадрат

Definition at line 101 of file defsBH.h.

    {
        ADDOP(1);
        return x * x;
    }

Here is the call graph for this function:

Here is the caller graph for this function:

Variable Documentation

◆ DPI

const double BH::DPI = 2.0 * 3.1415926535897932384626

static

Definition at line 58 of file defsBH.h.

◆ IDPI

const double BH::IDPI = 0.5 / PI

static

Definition at line 60 of file defsBH.h.

◆ IPI

const double BH::IPI = 1.0 / PI

static

Definition at line 59 of file defsBH.h.

◆ op

long long BH::op

extern

Глобальная переменная - счетчик количества операций

Definition at line 46 of file TreeBH.cpp.

◆ PI

const double BH::PI = 3.1415926535897932384626

static

Definition at line 57 of file defsBH.h.

◆ twoPowCodeLengthVar

const int BH::twoPowCodeLengthVar = (1 << codeLength)

static

2 в степени длины мортоновского кода (на каждую координату)

Definition at line 64 of file defsBH.h.

Classes

Functions

Variables

Function Documentation

◆ ceilhalf()

◆ ceilpow2()

◆ norm()

◆ operator&()

◆ operator*()

◆ operator+()

◆ operator+=()

◆ operator-()

◆ operator-=()

◆ powz()

◆ sign()

◆ SizeCheck()

◆ sqr()

Variable Documentation

◆ DPI

◆ IDPI

◆ IPI

◆ op

◆ PI

◆ twoPowCodeLengthVar