VM2D::World2D Class Reference

Класс, опеделяющий текущую решаемую задачу More...

#include <World2D.h>

Inheritance diagram for VM2D::World2D:

Collaboration diagram for VM2D::World2D:

Public Member Functions
const Airfoil &	getAirfoil (size_t i) const
	Возврат константной ссылки на объект профиля More...
const Airfoil &	getOldAirfoil (size_t i) const
	Возврат константной ссылки на объект старого профиля More...
Airfoil &	getNonConstAirfoil (size_t i) const
	Возврат неконстантной ссылки на объект профиля More...
size_t	getNumberOfAirfoil () const
	Возврат количества профилей в задаче More...
const Boundary &	getBoundary (size_t i) const
	Возврат константной ссылки на объект граничного условия More...
Boundary &	getNonConstBoundary (size_t i) const
	Возврат неконстантной ссылки на объект граничного условия More...
size_t	getNumberOfBoundary () const
	Возврат количества граничных условий в задаче More...
size_t	getDispBoundaryInSystem (size_t i) const
	Возврат смещения в системе dispBoundaryInSystem. More...
const MeasureVP &	getMeasureVP () const
	Возврат константной ссылки на measureVP. More...
MeasureVP &	getNonConstMeasureVP () const
	Возврат неконстантной ссылки на measureVP. More...
const Mechanics &	getMechanics (size_t i) const
	Возврат константной ссылки на объект механики More...
Mechanics &	getNonConstMechanics (size_t i) const
	Возврат неконстантной ссылки на объект механики More...
const Wake &	getWake () const
	Возврат константной ссылки на вихревой след More...
Wake &	getNonConstWake () const
	Возврат неконстантной ссылки на вихревой след More...
const WakeDataBase &	getSource () const
	Возврат константной ссылки на источники в области течения More...
const Velocity &	getVelocity () const
	Возврат константной ссылки на объект для вычисления скоростей More...
Velocity &	getNonConstVelocity () const
	Возврат неконстантной ссылки на объект для вычисления скоростей More...
const Passport &	getPassport () const
	Возврат константной ссылки на паспорт More...
const Gpu &	getCuda () const
	Возврат константной ссылки на объект, связанный с видеокартой (GPU) More...
Gpu &	getNonConstCuda () const
	Возврат неконстантной ссылки на объект, связанный с видеокартой (GPU) More...
const std::pair< Eigen::MatrixXd, Eigen::MatrixXd > &	getIQ (size_t i, size_t j) const
	Возврат константной ссылки на объект, связанный с матрицей интегралов от (r-xi)/|r-xi|^2. More...
Times &	getTimestat () const
	Возврат ссылки на временную статистику выполнения шага расчета по времени More...
bool	ifDivisible (int val) const
void	SolveLinearSystem ()
	Решение системы линейных алгебраических уравнений More...
void	FillIQ ()
	Заполнение матрицы, состоящей из интегралов от (r-xi) / |r-xi|^2. More...
void	FillMatrixAndRhs ()
	Заполнение матрицы системы для всех профилей More...
void	ReserveMemoryForMatrixAndRhs ()
	Вычисляем размер матрицы и резервируем память под нее и под правую часть More...
void	CalcVortexVelo (bool shiftTime)
	Вычисление скоростей (и конвективных, и диффузионных) вихрей (в пелене и виртуальных), а также в точках вычисления VP. More...
void	CalcPanelsVeloAndAttachedSheets ()
	Вычисление скоростей панелей и интенсивностей присоединенных слоев вихрей и источников More...
void	CalcAndSolveLinearSystem ()
	Набор матрицы, правой части и решение СЛАУ More...
void	MoveVortexes (std::vector< Point2D > &newPos)
	Вычисляем новые положения вихрей (в пелене и виртуальных) More...
void	CheckInside (std::vector< Point2D > &newPos, const std::vector< std::unique_ptr< Airfoil >> &oldAirfoil)
	Проверка проникновения вихрей внутрь профиля More...
void	WakeAndAirfoilsMotion ()
	Перемещение вихрей и профилей на шаге More...
	World2D (const VMlib::PassportGen &passport_)
	Конструктор More...
	~World2D ()
	Деструктор More...
void	GenerateMechanicsHeader (size_t mechanicsNumber)
virtual void	Step () override
	Функция выполнения предварительного шага More...
VMlib::LogStream &	getInfo () const
	Возврат ссылки на объект LogStream Используется в техничеcких целях для организации вывода More...
std::ostream &	getInfo (char x) const
	Возврат ссылки на поток вывода информации Необходимо для вывода телеметрической информации, информации об ошибках и т.п. More...
size_t	getCurrentStep () const
	Возврат константной ссылки на параметры распараллеливания по MPI. More...
const PassportGen &	getPassportGen ()
bool	isFinished () const
	Функция, возвращающая признак завершения счета в решаемой задаче More...

Public Attributes
size_t	currentStep
	Текущий номер шага в решаемой задаче More...

Protected Attributes
LogStream	info
	Поток для вывода логов и сообщений об ошибках More...
const PassportGen &	passportGen
	Константная ссылка на паспорт конкретного расчета More...
std::unique_ptr< TimesGen >	timestat
	Сведения о временах выполнения основных операций More...

Private Attributes
std::vector< std::unique_ptr< Airfoil > >	airfoil
	Список умных указателей на обтекаемые профили More...
std::vector< std::unique_ptr< Airfoil > >	oldAirfoil
	Список умных указателей на обтекаемые профили для сохранения старого положения More...
std::vector< std::unique_ptr< Boundary > >	boundary
	Список умных указателей на формирователи граничных условий на профилях More...
std::vector< size_t >	dispBoundaryInSystem
	Список номеров, с которых начинаются элементы правой части (или матрицы) системы для профилей More...
std::vector< std::unique_ptr< Mechanics > >	mechanics
	Список умных указателей на типы механической системы для каждого профиля More...
std::unique_ptr< Velocity >	velocity
	Умный укзатель на объект, определяющий методику вычисления скоростей More...
std::unique_ptr< Wake >	wake
	Умный указатель на вихревой след More...
std::unique_ptr< WakeDataBase >	source
	Умный указатель на источники More...
std::unique_ptr< MeasureVP >	measureVP
	Умный указатель на алгоритм вычисления полей скоростей и давления (для сохранения в файл) More...
Eigen::MatrixXd	matr
	Матрица системы More...
std::vector< std::vector< std::pair< Eigen::MatrixXd, Eigen::MatrixXd > > >	IQ
	Матрица, состоящая из пар матриц, в которых хранятся касательные и нормальные компоненты интегралов от ядра More...
Eigen::MatrixXd	invMatr
	Обратная матрица More...
bool	useInverseMatrix
	Признак использования обратной матрицы More...
Eigen::VectorXd	rhs
	Правая часть системы More...
Eigen::VectorXd	sol
	Решение системы More...
const Passport &	passport
	Константная ссылка на паспорт конкретного расчета More...
Gpu	cuda
	Объект, управляющий графическим ускорителем More...

Detailed Description

Класс, опеделяющий текущую решаемую задачу

Author

Марчевский Илья Константинович

Сокол Ксения Сергеевна

Рятина Евгения Павловна

Колганова Александра Олеговна 1.12

Date

14 января 2024 г.

Definition at line 68 of file World2D.h.

Constructor & Destructor Documentation

World2D::World2D

(

const VMlib::PassportGen &

passport_

)

Конструктор

Parameters

[in]

passport_

константная ссылка на паспорт расчета

Definition at line 67 of file World2D.cpp.

                                                  :
    WorldGen(passport_),
    passport(dynamic_cast<const Passport&>(passport_)),
    cuda(Gpu(*this))
{
    std::stringstream ss;
    ss << "#" << passport.problemNumber << " (" << passport.problemName << ")";     
    info.assignStream(defaults::defaultWorld2DLogStream, ss.str()); 
        
    passport.physicalProperties.setCurrTime(passport.timeDiscretizationProperties.timeStart);
    currentStep = 0;
        
    timestat.reset(new Times(*this)),
    
    wake.reset(new Wake(*this));
    // загрузка пелены из файла
    if (passport.wakeDiscretizationProperties.fileWake != "")
        wake->ReadFromFile(passport.wakesDir, passport.wakeDiscretizationProperties.fileWake); //Считываем из каталога с пеленой
    
    source.reset(new WakeDataBase(*this));
    // загрузка положений источников из файла
    if (passport.wakeDiscretizationProperties.fileSource != "")
        source->ReadFromFile(passport.dir, passport.wakeDiscretizationProperties.fileSource); //Считываем из текущего каталога

    //считываем массив точек для подсчета и вывода поля скоростей и давлений
    measureVP.reset(new MeasureVP(*this));
    if (passport.timeDiscretizationProperties.saveVPstep != 0)
        measureVP->ReadPointsFromFile(passport.dir);

    switch (passport.numericalSchemes.velocityComputation.second)
    {
    case 0:
        velocity.reset(new VelocityBiotSavart(*this));
        break;
    //case 1:
    //  velocity.reset(new VelocityBarnesHut(*this));
    //  break;
    }

    velocity->virtualVortexesParams.resize(passport.airfoilParams.size());

    for (size_t i = 0; i < passport.airfoilParams.size(); ++i)
    {
        switch (passport.numericalSchemes.panelsType.second)
        {
        case 0:         
            airfoil.emplace_back(new AirfoilRect(*this, i));
            break;
        //case 1:
        //  airfoil.emplace_back(new AirfoilCurv(*this, i));
        //  break;
        }

        airfoil[i]->ReadFromFile(passport.airfoilsDir); //Считываем из каталога с коллекцией профилей
        
        switch (passport.numericalSchemes.boundaryCondition.second)
        {
        case 0:
            boundary.emplace_back(new BoundaryConstLayerAver(*this, i));
            break;

        case 1:
            boundary.emplace_back(new BoundaryLinLayerAver(*this, i));
            break;
        
        case 10:            
            boundary.emplace_back(new BoundaryVortexCollocN(*this, i));
            //info('e') << "BoundaryMDV is not implemented now! " << std::endl;
            //exit(1);
            break;

        default:
            info('e') << "Unknown scheme!" << std::endl;
            exit(1);
        }


        switch (passport.airfoilParams[i].mechanicalSystemType)
        {
        case 0:
            mechanics.emplace_back(new MechanicsRigidImmovable(*this, i));
            break;
        case 1:
            mechanics.emplace_back(new MechanicsRigidGivenLaw(*this, i));
            break;
        case 2:
            mechanics.emplace_back(new MechanicsRigidOscillPart(*this, i));
            break;

        case 3:
            mechanics.emplace_back(new MechanicsRigidRotatePart(*this, i));
            break;
        /*
        case 4:
            mechanics.emplace_back(new MechanicsRigidOscillMon(*this, i));
            break;
        */
        }

    }

    IQ.resize(passport.airfoilParams.size());
    for (size_t i = 0; i < passport.airfoilParams.size(); ++i)
        IQ[i].resize(passport.airfoilParams.size());


    info.endl();
    info('i') << "Start solving problem " << passport.problemName << std::endl;
    info.endl();

    VMlib::PrintLogoToStream(info('_') << std::endl);
}//World2D(...)

Here is the call graph for this function:

Here is the caller graph for this function:

VM2D::World2D::~World2D ( )

inline

Деструктор

Definition at line 306 of file World2D.h.

{ };

Here is the call graph for this function:

Member Function Documentation

void World2D::CalcAndSolveLinearSystem

(

)

Набор матрицы, правой части и решение СЛАУ

Вызывается в Step()

Definition at line 774 of file World2D.cpp.

{
    if (airfoil.size() > 0)
    {
        ReserveMemoryForMatrixAndRhs();

#if defined(__CUDACC__) || defined(USE_CUDA)
        cuda.setAccelCoeff(passport.physicalProperties.accelCft());
        cuda.setMaxGamma(passport.wakeDiscretizationProperties.maxGamma);

        int sch = passport.numericalSchemes.boundaryCondition.second;

        if ((sch == 0) || (sch == 1) || (sch == 10))
            cuda.setSchemeSwitcher(sch + 1);
        else
        {
            info('e') << "schemeSwitcher is not 0, or 1, or 10! " << std::endl;
            exit(1);
        }

        cuda.RefreshWake(1);
        cuda.RefreshAfls(1);
        cuda.RefreshVirtualWakes(1);
#endif

        if (currentStep == 0)
        {
            useInverseMatrix = (
                (mechanics.size() == 1)
                ||
                (mechanics.size() > 1 && !std::any_of(mechanics.begin(), mechanics.end(), [](const std::unique_ptr<Mechanics>& m) { return m->isMoves; }))
                );
        }

        FillIQ();
        FillMatrixAndRhs();

        //{
        //  std::stringstream ss;
        //  ss << "IQ1-" << currentStep;
        //  std::ofstream of(passport.dir + "dbg/" + ss.str());
        //  for (size_t i = 0; i < (IQ[0][0]).first.rows(); ++i)
        //  {
        //      for (size_t j = 0; j < (IQ[0][0]).first.cols(); ++j)
        //          of << (IQ[0][0]).first(i, j) << " ";
        //      of << std::endl;
        //  }
        //  of.close();
        //}

        /*
        {
            std::stringstream ss;
            ss << "matr-" << currentStep;
            std::ofstream of(passport.dir + "dbg/" + ss.str());
            for (size_t i = 0; i < matr.rows(); ++i)
            {
                for (size_t j = 0; j < matr.cols(); ++j)
                    of << matr(i, j) << " ";
                of << std::endl;
            }
            of.close();
        }
        */

        SolveLinearSystem();

        getTimestat().timeOther.first += omp_get_wtime();

        size_t currentRow = 0;
        for (size_t bou = 0; bou < boundary.size(); ++bou)
        {
            size_t nVars = boundary[bou]->GetUnknownsSize();
            Eigen::VectorXd locSol;
            locSol.resize(nVars);
            for (size_t i = 0; i < nVars; ++i)
                locSol(i) = sol(currentRow + i);

            boundary[bou]->SolutionToFreeVortexSheetAndVirtualVortex(locSol);
            currentRow += nVars + 1;
        }

        if (currentStep == 0)
        {
            Point2D addMass = { 0.0, 0.0 };
            for (size_t q = 0; q < airfoil[0]->getNumberOfPanels(); ++q)
            {
                addMass += (sol(q) + boundary[0]->sheets.attachedVortexSheet(q,0)) * 0.5 * (airfoil[0]->getR(q) + airfoil[0]->getR(q + 1)).kcross() * airfoil[0]->len[q];           
            }
            addMass *= passport.physicalProperties.rho;

            std::cout << "AddMass = " << addMass << std::endl;
            //exit(-42);
        }

        getTimestat().timeOther.second += omp_get_wtime();
    }
}//CalcAndSolveLinearSystem()

Here is the call graph for this function:

Here is the caller graph for this function:

void World2D::CalcPanelsVeloAndAttachedSheets

(

)

Вычисление скоростей панелей и интенсивностей присоединенных слоев вихрей и источников

Вызывается в Step()

Definition at line 664 of file World2D.cpp.

{
    getTimestat().timeOther.first += omp_get_wtime();

    //вычисляем скорости панелей
    for (size_t i = 0; i < airfoil.size(); ++i)
        mechanics[i]->VeloOfAirfoilPanels(currentStep * passport.timeDiscretizationProperties.dt);

    //вычисляем интенсивности присоединенных слоев
    for (size_t i = 0; i < airfoil.size(); ++i)
        boundary[i]->ComputeAttachedSheetsIntensity();

    getTimestat().timeOther.second += omp_get_wtime();
}//CalcPanelsVeloAndAttachedSheets(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void World2D::CalcVortexVelo

(

bool

shiftTime

)

Вычисление скоростей (и конвективных, и диффузионных) вихрей (в пелене и виртуальных), а также в точках вычисления VP.

Вызывается в Step()

Definition at line 563 of file World2D.cpp.

{
    getTimestat().timeCalcVortexConvVelo.first += omp_get_wtime();
    
    velocity->ResizeAndZero();
    
    //Конвективные скорости всех вихрей (и в пелене, и виртуальных), индуцируемые вихрями в пелене
    //Подготовка CUDA
#if defined(__CUDACC__) || defined(USE_CUDA)
    cuda.RefreshWake(2);
    cuda.RefreshAfls(2);    
    cuda.RefreshVirtualWakes(2);
#endif
    getTimestat().timeCalcVortexConvVelo.second += omp_get_wtime();

    //Вычисление конвективных скоростей вихрей (и в пелене, и виртуальных), а также в точках wakeVP
    if (shiftTime)
        --currentStep;
    
    velocity->CalcConvVelo();
    
    if (shiftTime)
        ++currentStep;

    //Расчет средних значений eps для каждой панели и их передача на видеокарту
    getTimestat().timeCalcVortexDiffVelo.first += omp_get_wtime();
    for (size_t bou = 0; bou < getNumberOfBoundary(); ++bou)
        getNonConstAirfoil(bou).calcMeanEpsOverPanel();

#if defined(__CUDACC__) || defined(USE_CUDA)
    for (size_t i = 0; i < airfoil.size(); ++i)
        cuda.CopyMemToDev<double, 1>(airfoil[i]->getNumberOfPanels(), airfoil[i]->meanEpsOverPanel.data(), airfoil[i]->devMeanEpsOverPanelPtr);
#endif
    getTimestat().timeCalcVortexDiffVelo.second += omp_get_wtime();
    
    //Вычисление диффузионных скоростей вихрей (и в пелене, и виртуальных)
    velocity->CalcDiffVelo();       
    
    getTimestat().timeSaveKadr.first += omp_get_wtime();
    /*//Сохранение всех параметров для вихрей в пелене

    {
        VMlib::CreateDirectory(passport.dir, "dbg");
        std::ostringstream sss;
        sss << "prmWake";
        sss << currentStep;
        std::ofstream prmtFile(passport.dir + "dbg/" + sss.str());
        prmtFile << "i x y g epsast convVeloX convVeloY diffVeloX diffVeloY I0 I1 I2X I2Y I3X I3Y" << std::endl;
        for (size_t i = 0; i < wake->vtx.size(); ++i)
            prmtFile << i << " " \
            << wake->vtx[i].r()[0] << " " << wake->vtx[i].r()[1] << " " \
            << wake->vtx[i].g() << " " \
            << velocity->wakeVortexesParams.epsastWake[i] << " " \
            << velocity->wakeVortexesParams.convVelo[i][0] << " " << velocity->wakeVortexesParams.convVelo[i][1] << " "\
            << velocity->wakeVortexesParams.diffVelo[i][0] << " " << velocity->wakeVortexesParams.diffVelo[i][1] << " "\
            //<< std::endl;
            << velocity->wakeVortexesParams.I0[i] << " " \
            << velocity->wakeVortexesParams.I1[i] << " " \
            << velocity->wakeVortexesParams.I2[i][0] << " " << velocity->wakeVortexesParams.I2[i][1] << " " \
            << velocity->wakeVortexesParams.I3[i][0] << " " << velocity->wakeVortexesParams.I3[i][1] << " " \
            << std::endl;

        prmtFile.close();
    }
//*/

/*
    //Сохранение всех параметров для виртуальных вихрей
    {
        for (size_t b = 0; b < boundary.size(); ++b)
        {
            std::ostringstream sss;
            sss << "prmVirtual_";
            sss << b << "-";
            sss << currentStep;
            std::ofstream prmFileVirt(passport.dir + "dbg/" + sss.str());
            prmFileVirt << "i x y g epsast convVeloX convVeloY diffVeloX diffVeloY I0 I1 I2X I2Y I3X I3Y" << std::endl;
            for (size_t i = 0; i < boundary[b]->virtualWake.vtx.size(); ++i)
                prmFileVirt << i << " " \
                << boundary[b]->virtualWake.vtx[i].r()[0] << " " << boundary[b]->virtualWake.vtx[i].r()[1] << " " \
                << boundary[b]->virtualWake.vtx[i].g() << " " \
                << velocity->virtualVortexesParams[b].epsastWake[i] << " " \
                << velocity->virtualVortexesParams[b].convVelo[i][0] << " " << velocity->virtualVortexesParams[b].convVelo[i][1] << " "\
                << velocity->virtualVortexesParams[b].diffVelo[i][0] << " " << velocity->virtualVortexesParams[b].diffVelo[i][1] << " "\
                //<< std::endl;             
                << velocity->virtualVortexesParams[b].I0[i] << " " \
                << velocity->virtualVortexesParams[b].I1[i] << " " \
                << velocity->virtualVortexesParams[b].I2[i][0] << " " << velocity->virtualVortexesParams[b].I2[i][1] << " " \
                << velocity->virtualVortexesParams[b].I3[i][0] << " " << velocity->virtualVortexesParams[b].I3[i][1] << " " \
                << std::endl;
            prmFileVirt.close();
        }
        //if (currentStep==2) exit(-123);
    }
*/
    getTimestat().timeSaveKadr.second += omp_get_wtime();
}//CalcVortexVelo()

Here is the call graph for this function:

Here is the caller graph for this function:

void World2D::CheckInside	(	std::vector< Point2D > &	newPos,
		const std::vector< std::unique_ptr< Airfoil >> &	oldAirfoil
	)

Проверка проникновения вихрей внутрь профиля

Вызывается в Step()

Parameters

[in]	newPos	новые позиции вихрей
[in]	oldAirfoil	константная ссылка на вектор из умных указателей на старые положения профилей

Definition at line 294 of file World2D.cpp.

{
    getTimestat().timeCheckInside.first += omp_get_wtime();
        

#if (!defined(USE_CUDA))    
    for (size_t afl = 0; afl < airfoil.size(); ++afl)
        wake->Inside(newPos, *airfoil[afl], mechanics[afl]->isMoves, *oldAirfoil[afl]);
#else   
    if (newPos.size() > 0)
    {
        double* devNewpos_ptr;
        cuReserveDevMem((void*&)devNewpos_ptr, newPos.size() * sizeof(double) * 2, 0);
        cuCopyFixedArray(devNewpos_ptr, newPos.data(), newPos.size() * sizeof(double) * 2, 0);

        for (size_t afl = 0; afl < airfoil.size(); ++afl)
        {
            std::vector<double> gamma(airfoil[afl]->getNumberOfPanels(), 0.0);
            std::vector<unsigned int> hit = lbvh_check_inside((int)currentStep, (int)newPos.size(), devNewpos_ptr, (int)airfoil[afl]->getNumberOfPanels(), airfoil[afl]->devRPtr);
            for (int i = 0; i < newPos.size(); ++i)
            {
                if (hit[i] != (unsigned int)(-1))
                {
                    gamma[hit[i]] += wake->vtx[i].g();
                    wake->vtx[i].g() = 0.0;
                }
            }
            airfoil[afl]->gammaThrough = gamma;
        }//for afl

        cuDeleteFromDev(devNewpos_ptr);
    }
#endif

    getTimestat().timeCheckInside.second += omp_get_wtime();
}

Here is the call graph for this function:

Here is the caller graph for this function:

void World2D::FillIQ

(

)

Заполнение матрицы, состоящей из интегралов от (r-xi) / |r-xi|^2.

Вызывается в Step()

Definition at line 732 of file World2D.cpp.

{
    getTimestat().timeFillMatrixAndRhs.first += omp_get_wtime();
    
    for (size_t bou = 0; bou < boundary.size(); ++bou)
    {
        if (currentStep == 0 || mechanics[bou]->isDeform)
        {
            boundary[bou]->FillIQSelf(IQ[bou][bou]);
        }

        for (size_t oth = 0; oth < boundary.size(); ++oth)
        {           

#ifdef INITIAL
            if (currentStep == 0 || !useInverseMatrix)
            {
                //size_t nVarsOther = boundary[oth]->GetUnknownsSize();
                if (bou != oth)
                {
                    //std::cout << "matr!" << std::endl;
                    boundary[bou]->FillIQFromOther(*boundary[oth], IQ[bou][oth]);
                }
            }// if (currentStep == 0 || !useInverseMatrix)
#endif

#ifdef BRIDGE
            if (currentStep == 0)
            {
                size_t nVarsOther = boundary[oth]->GetUnknownsSize();
                if (bou != oth)
                    boundary[bou]->FillIQFromOther(*boundary[oth], IQ[bou][oth]);
            }// if (currentStep == 0)
#endif

        }// for oth

    }// for bou
       
    getTimestat().timeFillMatrixAndRhs.second += omp_get_wtime();
}//FillIQ()

Here is the call graph for this function:

Here is the caller graph for this function:

void World2D::FillMatrixAndRhs

(

)

Заполнение матрицы системы для всех профилей

Вызывается в Step()

Definition at line 438 of file World2D.cpp.

{
    getTimestat().timeFillMatrixAndRhs.first += omp_get_wtime();

    Eigen::MatrixXd locMatr;
    Eigen::MatrixXd otherMatr;
    Eigen::VectorXd locLastLine, locLastCol;
    
    std::vector<std::vector<Point2D>> locIQ;
    std::vector<std::vector<Point2D>> locOtherIQ;

    //обнуляем матрицу на первом шаге расчета
    if (currentStep == 0)
    {
        for (int i = 0; i < matr.rows(); ++i)
        for (int j = 0; j < matr.cols(); ++j)
            matr(i, j) = 0.0;
    }

    size_t currentRow = 0;

    for (size_t bou = 0; bou < boundary.size(); ++bou)
    {
        size_t nVars = boundary[bou]->GetUnknownsSize();
        if (currentStep == 0)
        {
            locMatr.resize(nVars, nVars);
            locLastLine.resize(nVars);
            locLastCol.resize(nVars);
        }

        if (currentStep == 0 || mechanics[bou]->isDeform)
        {
            boundary[bou]->FillMatrixSelf(locMatr, locLastLine, locLastCol);
        }

        //размазываем матрицу       

        for (size_t i = 0; i < nVars; ++i)
        {
            if (currentStep == 0 || mechanics[bou]->isDeform)
            {
                for (size_t j = 0; j < nVars; ++j)
                    matr(i + currentRow, j + currentRow) = locMatr(i, j);
                matr(currentRow + nVars, i + currentRow) = locLastLine(i);
                matr(i + currentRow, currentRow + nVars) = locLastCol(i);
            }
        }
                
        if ( (currentStep == 0) || (!useInverseMatrix) )
        {
            size_t currentCol = 0;
            for (size_t oth = 0; oth < boundary.size(); ++oth)
            {
                size_t nVarsOther = boundary[oth]->GetUnknownsSize();

                if (bou != oth)
                {
                    otherMatr.resize(nVars, nVarsOther);

                    boundary[bou]->FillMatrixFromOther(*boundary[oth], otherMatr);

                    //размазываем матрицу
                    for (size_t i = 0; i < nVars; ++i)
                    {
                        for (size_t j = 0; j < nVarsOther; ++j)
                            matr(i + currentRow, j + currentCol) = otherMatr(i, j);
                    }
                }// if (bou != oth)
                currentCol += nVarsOther + 1;
            }// for oth
        }// if (currentStep == 0 || mechanics[oth]->isMoves)

        currentRow += nVars + 1;
    }// for bou

    velocity->FillRhs(rhs);
    
    getTimestat().timeFillMatrixAndRhs.second += omp_get_wtime();
}//FillMatrixAndRhs()

Here is the call graph for this function:

Here is the caller graph for this function:

void World2D::GenerateMechanicsHeader

(

size_t

mechanicsNumber

)

Метод-обертка для вызова метода генерации заголовка файла нагрузок и заголовка файла положения (последнее — если профиль движется)

Parameters

[in]

mechanicsNumber

номер профиля, для которого генерируется заголовок файла

Definition at line 725 of file World2D.cpp.

{
    mechanics[mechanicsNumber]->GenerateForcesHeader();
    mechanics[mechanicsNumber]->GeneratePositionHeader();
}//GenerateMechanicsHeader(...)

Here is the caller graph for this function:

const Airfoil& VM2D::World2D::getAirfoil ( size_t  i ) const

inline

Возврат константной ссылки на объект профиля

Parameters

[in]

i

номер профиля, константная ссылка на который возвращается

Returns

константную ссылку на i-й профиль

Definition at line 130 of file World2D.h.

{ return *airfoil[i]; };

Here is the caller graph for this function:

const Boundary& VM2D::World2D::getBoundary ( size_t  i ) const

inline

Возврат константной ссылки на объект граничного условия

Parameters

[in]

i

номер граничного условия, константная ссылка на которое возвращается

Returns

константную ссылку на i-е граничное условие

Definition at line 153 of file World2D.h.

{ return *boundary[i]; };

Here is the caller graph for this function:

const Gpu& VM2D::World2D::getCuda ( ) const

inline

Возврат константной ссылки на объект, связанный с видеокартой (GPU)

Returns

константную ссылку на объект, связанный с видеокартой (GPU)

Definition at line 227 of file World2D.h.

{ return cuda; };

Here is the caller graph for this function:

size_t VMlib::WorldGen::getCurrentStep ( ) const

inlineinherited

Возврат константной ссылки на параметры распараллеливания по MPI.

Returns

константную ссылку на параметры распараллеливания по MPI const Parallel& getParallel() const { return parallel; }; Возврат номера текущего временного шага

номера текущего временного шага

Definition at line 91 of file WorldGen.h.

{ return currentStep; };

Here is the caller graph for this function:

size_t VM2D::World2D::getDispBoundaryInSystem ( size_t  i ) const

inline

Возврат смещения в системе dispBoundaryInSystem.

Parameters

[in]

i

номер граничного условия, константная ссылка на которое возвращается

Returns

константную ссылку на i-е граничное условие

Definition at line 170 of file World2D.h.

{ return dispBoundaryInSystem[i]; };

VMlib::LogStream& VMlib::WorldGen::getInfo ( ) const

inlineinherited

Возврат ссылки на объект LogStream Используется в техничеcких целях для организации вывода

Returns

ссылку на объект LogStream

Definition at line 74 of file WorldGen.h.

{ return info; };

Here is the caller graph for this function:

std::ostream& VMlib::WorldGen::getInfo ( char  x ) const

inlineinherited

Возврат ссылки на поток вывода информации Необходимо для вывода телеметрической информации, информации об ошибках и т.п.

Parameters

[in]

x

символ, определяющий стиль вывода сообщения

Returns

ссылку на поток вывода информации

Definition at line 81 of file WorldGen.h.

{ return info(x); };

const std::pair<Eigen::MatrixXd, Eigen::MatrixXd>& VM2D::World2D::getIQ ( size_t  i, size_t  j  ) const

inline

Возврат константной ссылки на объект, связанный с матрицей интегралов от (r-xi)/|r-xi|^2.

Returns

константную ссылку на объект, связанный с матрицей интегралов от (r-xi)/|r-xi|^2

Definition at line 237 of file World2D.h.

{ return IQ[i][j]; };

Here is the caller graph for this function:

const MeasureVP& VM2D::World2D::getMeasureVP ( ) const

inline

Возврат константной ссылки на measureVP.

Returns

константную ссылку на вихревой след

Definition at line 175 of file World2D.h.

{ return *measureVP; };

Here is the caller graph for this function:

const Mechanics& VM2D::World2D::getMechanics ( size_t  i ) const

inline

Возврат константной ссылки на объект механики

Parameters

[in]

i

номер механики, константная ссылка на который возвращается

Returns

константную ссылку на i-ю механику

Definition at line 186 of file World2D.h.

{ return *mechanics[i]; };

Here is the caller graph for this function:

Airfoil& VM2D::World2D::getNonConstAirfoil ( size_t  i ) const

inline

Возврат неконстантной ссылки на объект профиля

Parameters

[in]

i

номер профиля, неконстантная ссылка на который возвращается

Returns

неконстантную ссылку на i-й профиль

Definition at line 142 of file World2D.h.

{ return *airfoil[i]; };

Here is the caller graph for this function:

Boundary& VM2D::World2D::getNonConstBoundary ( size_t  i ) const

inline

Возврат неконстантной ссылки на объект граничного условия

Parameters

[in]

i

номер граничного условия, неконстантная ссылка на которое возвращается

Returns

неконстантную ссылку на i-е граничное условие

Definition at line 159 of file World2D.h.

{ return *boundary[i]; };

Here is the caller graph for this function:

Gpu& VM2D::World2D::getNonConstCuda ( ) const

inline

Возврат неконстантной ссылки на объект, связанный с видеокартой (GPU)

Returns

неконстантную ссылку на объект, связанный с видеокартой (GPU)

Definition at line 232 of file World2D.h.

{ return cuda; };

Here is the caller graph for this function:

MeasureVP& VM2D::World2D::getNonConstMeasureVP ( ) const

inline

Возврат неконстантной ссылки на measureVP.

Returns

неконстантную ссылку на вихревой след

Definition at line 180 of file World2D.h.

{ return *measureVP; };

Here is the caller graph for this function:

Mechanics& VM2D::World2D::getNonConstMechanics ( size_t  i ) const

inline

Возврат неконстантной ссылки на объект механики

Parameters

[in]

i

номер механики, константная ссылка на который возвращается

Returns

неконстантную ссылку на i-ю механику

Definition at line 192 of file World2D.h.

{ return *mechanics[i]; };

Here is the caller graph for this function:

Velocity& VM2D::World2D::getNonConstVelocity ( ) const

inline

Возврат неконстантной ссылки на объект для вычисления скоростей

Returns

неконстантную ссылку на объект для вычисления скоростей

Definition at line 217 of file World2D.h.

{ return *velocity; };

Here is the caller graph for this function:

Wake& VM2D::World2D::getNonConstWake ( ) const

inline

Возврат неконстантной ссылки на вихревой след

Returns

неконстантную ссылку на вихревой след

Definition at line 202 of file World2D.h.

{ return *wake; };

Here is the caller graph for this function:

size_t VM2D::World2D::getNumberOfAirfoil ( ) const

inline

Возврат количества профилей в задаче

Returns

количество профилей в задаче

Definition at line 147 of file World2D.h.

{ return airfoil.size(); };

Here is the caller graph for this function:

size_t VM2D::World2D::getNumberOfBoundary ( ) const

inline

Возврат количества граничных условий в задаче

Returns

количество граничных условий в задаче

Definition at line 164 of file World2D.h.

{ return boundary.size(); };

Here is the caller graph for this function:

const Airfoil& VM2D::World2D::getOldAirfoil ( size_t  i ) const

inline

Возврат константной ссылки на объект старого профиля

Parameters

[in]

i

номер старого профиля, константная ссылка на который возвращается

Returns

константную ссылку на i-й старый профиль

Definition at line 136 of file World2D.h.

{ return *oldAirfoil[i]; };

Here is the caller graph for this function:

const Passport& VM2D::World2D::getPassport ( ) const

inline

Возврат константной ссылки на паспорт

Returns

константную ссылку на паспорт

Definition at line 222 of file World2D.h.

{ return passport; };

Here is the caller graph for this function:

const PassportGen& VMlib::WorldGen::getPassportGen ( )

inlineinherited

Definition at line 93 of file WorldGen.h.

{ return passportGen; };

Here is the call graph for this function:

const WakeDataBase& VM2D::World2D::getSource ( ) const

inline

Возврат константной ссылки на источники в области течения

Returns

константную ссылку на источники в области течения

Definition at line 207 of file World2D.h.

{ return *source; };

Here is the caller graph for this function:

Times& VM2D::World2D::getTimestat ( ) const

inline

Возврат ссылки на временную статистику выполнения шага расчета по времени

Returns

ссылку на временную статистику выполнения шага расчета по времени

Definition at line 242 of file World2D.h.

{ return dynamic_cast<Times&>(*timestat); };

Here is the caller graph for this function:

const Velocity& VM2D::World2D::getVelocity ( ) const

inline

Возврат константной ссылки на объект для вычисления скоростей

Returns

константную ссылку на объект для вычисления скоростей

Definition at line 212 of file World2D.h.

{ return *velocity; };

Here is the caller graph for this function:

const Wake& VM2D::World2D::getWake ( ) const

inline

Возврат константной ссылки на вихревой след

Returns

константную ссылку на вихревой след

Definition at line 197 of file World2D.h.

{ return *wake; };

Here is the caller graph for this function:

bool VM2D::World2D::ifDivisible ( int  val ) const

inline

Definition at line 244 of file World2D.h.

{ return ((val > 0) && (!(currentStep % val))); };

Here is the call graph for this function:

Here is the caller graph for this function:

bool WorldGen::isFinished ( ) const

inherited

Функция, возвращающая признак завершения счета в решаемой задаче

true, если задача решена и выполнен признак останова; false если требуется еще выполнять шаги по времени

Definition at line 49 of file WorldGen.cpp.

{
    return (passportGen.timeDiscretizationProperties.currTime >= passportGen.timeDiscretizationProperties.timeStop);
}

void World2D::MoveVortexes

(

std::vector< Point2D > &

newPos

)

Вычисляем новые положения вихрей (в пелене и виртуальных)

Вызывается в Step()

Parameters

[out]

newPos

новые позиции вихрей

Definition at line 681 of file World2D.cpp.

{
    getTimestat().timeMoveVortexes.first += omp_get_wtime();

    size_t nvt = wake->vtx.size();
    for (size_t i = 0; i < getNumberOfBoundary(); ++i)
        nvt += boundary[i]->virtualWake.vtx.size();

    //newPos.clear();
    newPos.resize(nvt);

    
#pragma omp parallel for
    for (int i = 0; i < (int)wake->vtx.size(); ++i)
    {
        newPos[i] = wake->vtx[i].r() \
            + (velocity->wakeVortexesParams.convVelo[i] +
                velocity->wakeVortexesParams.diffVelo[i] +
                passport.physicalProperties.V0()) * passport.timeDiscretizationProperties.dt;
    }

    int curCounterNewPos = (int)wake->vtx.size();

    wake->vtx.resize(nvt);

    
    for (size_t bou = 0; bou < boundary.size(); ++bou)
    {
#pragma omp parallel for
        for (int i = 0; i < (int)boundary[bou]->virtualWake.vtx.size(); ++i)
        {
            wake->vtx[curCounterNewPos + i] = (boundary[bou]->virtualWake.vtx[i]);
            newPos[curCounterNewPos + i] = (boundary[bou]->virtualWake.vtx[i].r() \
                + (velocity->virtualVortexesParams[bou].convVelo[i] +
                    velocity->virtualVortexesParams[bou].diffVelo[i] +
                    passport.physicalProperties.V0()) * passport.timeDiscretizationProperties.dt);
        }
        curCounterNewPos += (int)boundary[bou]->virtualWake.vtx.size();
    }

    getTimestat().timeMoveVortexes.second += omp_get_wtime();
}//MoveVortexes(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void World2D::ReserveMemoryForMatrixAndRhs

(

)

Вычисляем размер матрицы и резервируем память под нее и под правую часть

Вызывается в Step()

Definition at line 520 of file World2D.cpp.

{
    getTimestat().timeReserveMemoryForMatrixAndRhs.first += omp_get_wtime();

    
    if (currentStep == 0)
    {
        dispBoundaryInSystem.resize(boundary.size());
        dispBoundaryInSystem[0] = 0;
        
        for (size_t i = 1; i < boundary.size(); ++i)
        {
            dispBoundaryInSystem[i] = dispBoundaryInSystem[i - 1] + boundary[i - 1]->GetUnknownsSize() + 1;
        }

        size_t matrSize = boundary.size();
        for (auto it = boundary.begin(); it != boundary.end(); ++it)
            matrSize += (*it)->GetUnknownsSize();

        matr.resize(matrSize, matrSize);
        matr.setZero();

        for (size_t i = 0; i < getNumberOfAirfoil(); ++i)
        {
            size_t nVari = boundary[i]->GetUnknownsSize();
            for (size_t j = 0; j < getNumberOfAirfoil(); ++j)
            {
                size_t nVarj = boundary[j]->GetUnknownsSize();
                IQ[i][j].first.resize(nVari, nVarj);
                IQ[i][j].second.resize(nVari, nVarj);
            }
        }

        rhs.resize(matrSize);
    }
    rhs.setZero();

    getTimestat().timeReserveMemoryForMatrixAndRhs.second += omp_get_wtime();
}//ReserveMemoryForMatrixAndRhs()

Here is the call graph for this function:

Here is the caller graph for this function:

void World2D::SolveLinearSystem

(

)

Решение системы линейных алгебраических уравнений

Вызывается в Step()

Definition at line 333 of file World2D.cpp.

{
    getTimestat().timeSolveLinearSystem.first += omp_get_wtime();

/*
    if (currentStep == 0)
    {
        invMatr = matr;
        std::ifstream fileMatrix("matrix.txt");
        int nx, ny;
        fileMatrix >> nx;
        fileMatrix >> ny;
        for (int i = 0; i < matr.rows(); ++i)
        {
            for (int j = 0; j < matr.cols(); ++j)
            {
                fileMatrix >> matr(i, j);
            }
        }
        fileMatrix.close();
    }
*/


    if (useInverseMatrix && (currentStep == 0))
    {
        info('t') << "Inverting matrix" << std::endl;

#if (defined(USE_CUDA)) 
        invMatr.resize(matr.rows(), matr.cols());
        for (int i = 0; i < (int)matr.rows(); ++i)
            for (int j = 0; j < (int)matr.cols(); ++j)
                invMatr(i, j) = (i == j) ? 1.0 : 0.0;
        cuInverseMatrix((int)matr.rows(), matr.data(), invMatr.data());
#else
        invMatr = matr.inverse();
#endif
    }



/*
    if (currentStep == 0)
    {
        invMatr = matr;
        std::ifstream fileMatrix("invMatrix.txt");
        int nx, ny;
        fileMatrix >> nx;
        fileMatrix >> ny;
        for (int i = 0; i < invMatr.rows(); ++i)
        {
            for (int j = 0; j < invMatr.cols(); ++j)
            {
                fileMatrix >> invMatr(i, j);
            }
        }
        fileMatrix.close();
    }
*/

/*
    if (currentStep == 0)
    {
        Eigen::MatrixXd mmul = matr * invMatr;
        for (int i = 0; i < invMatr.rows(); ++i)
            mmul(i,i) -= 1.0;
            
        double maxElem = 0.0;
        for (int i = 0; i < mmul.rows(); ++i)
            for (int j = 0; j < mmul.cols(); ++j)
                if (fabs(mmul(i,j)) > maxElem )
                    maxElem = mmul(i,j);

        std::cout << "A*A^-1: MAX ELEMENT = " << maxElem << std::endl;
    }
*/
    
/*
    if (currentStep == 0)
    {
        std::ofstream fileMatrix("invMatrix4x3200.txt");
        fileMatrix << invMatr.rows() << " " << invMatr.cols() << std::endl;
        fileMatrix.precision(18);
        for (int i = 0; i < invMatr.rows(); ++i)
        {
            for (int j = 0; j < invMatr.cols(); ++j)
            {
                fileMatrix << invMatr(i, j) << " ";
            }
            fileMatrix << std::endl;
        }
        fileMatrix.close();
    }
*/


    if (useInverseMatrix)
        sol = invMatr * rhs;
    else
        sol = matr.partialPivLu().solve(rhs);

    getTimestat().timeSolveLinearSystem.second += omp_get_wtime();
}//SolveLinearSystem()

Here is the call graph for this function:

Here is the caller graph for this function:

void World2D::Step ( )

overridevirtual

Функция выполнения предварительного шага

Основная функция выполнения одного шага по времени

Implements VMlib::WorldGen.

Definition at line 196 of file World2D.cpp.

{
    try {
        //Очистка статистики
        getTimestat().ToZero();

        //Засечка времени в начале шага
        getTimestat().timeWholeStep.first += omp_get_wtime();

        CalcPanelsVeloAndAttachedSheets();      
        measureVP->Initialization();
        
        //for (int i = 0; i < wake->vtx.size(); ++i)
        //{
        //  if (std::isnan(wake->vtx[i].r()[0]) || std::isnan(wake->vtx[i].r()[1]) || std::isnan(wake->vtx[i].g()))
        //  {
        //      std::cout << wake->vtx[i].r() << " " << wake->vtx[i].g() << std::endl;
        //      exit(-4);
        //  }
        //}
        
        
        CalcAndSolveLinearSystem();

        //Вычисление скоростей вихрей: для тех, которые в следе, и виртуальных, а также в точках wakeVP
        CalcVortexVelo(false);

        //Расчет и сохранение поля давления
        if (ifDivisible(passport.timeDiscretizationProperties.saveVPstep))
        {
            measureVP->CalcPressure();
            measureVP->SaveVP();
        }

        std::vector<MechanicsRigidOscillPart*> mechOscilPart;
        if (mechanics.size() > 0)
        {
            mechOscilPart.resize(mechanics.size(), nullptr);

            for (size_t s = 0; s < mechanics.size(); ++s)
                mechOscilPart[s] = dynamic_cast<MechanicsRigidOscillPart*>(mechanics[0].get());

            if (getPassport().airfoilParams[0].addedMass.length2() > 0)
            {
                for (size_t s = 0; s < mechanics.size(); ++s)
                    mechanics[s]->hydroDynamForce = { 0.0, 0.0 }; //12-01

                WakeAndAirfoilsMotion();                      //12-01
                CalcPanelsVeloAndAttachedSheets(); //12-01
                CalcAndSolveLinearSystem();        //12-01
            }

        }
        //Вычисление сил, действующих на профиль и сохранение в файл    

        for (auto& mech : mechanics)
        {
            mech->GetHydroDynamForce();         
            mech->GenerateForcesString();
            mech->GeneratePositionString();
        }
        
        for (size_t s = 0; s < mechanics.size(); ++s)
            if (mechOscilPart[s])
            {
                if (getPassport().airfoilParams[0].addedMass.length2() > 0)
                    mechOscilPart[s]->UpdateU();     //12-01
            }

        //Движение вихрей (сброс вихрей + передвижение пелены)
        WakeAndAirfoilsMotion();  //12-01   

        wake->Restruct();

        wake->SaveKadrVtk();

        //Засечка времени в конце шага
        getTimestat().timeWholeStep.second += omp_get_wtime();


        info('i') << "Step = " << currentStep \
            << " PhysTime = " << passport.physicalProperties.getCurrTime() \
            << " StepTime = " << Times::dT(getTimestat().timeWholeStep) << std::endl;
        getTimestat().GenerateStatString();


        passport.physicalProperties.addCurrTime(passport.timeDiscretizationProperties.dt);
        ++currentStep;
    }
    catch (...)
    {
        info('e') << "!!! Exception from unknown source !!!" << std::endl;
        exit(-1);
    }
}//Step()

Here is the call graph for this function:

Here is the caller graph for this function:

void World2D::WakeAndAirfoilsMotion

(

)

Перемещение вихрей и профилей на шаге

Вызывается в Step()

Definition at line 873 of file World2D.cpp.

{
    std::vector<Point2D> newPos;

    MoveVortexes(newPos);

#ifdef BRIDGE
    double totalForce = 0;
    for (size_t afl = 0; afl < airfoil.size(); ++afl)
    {
        totalForce += mechanics[afl]->hydroDynamForce[1];
    }
    totalForce *= 1.2;

    mechanics[0]->hydroDynamForce[1] = totalForce;
#endif

    getTimestat().timeOther.first += omp_get_wtime();

    oldAirfoil.resize(0);
    for (auto& afl : airfoil)
    {
        switch (passport.numericalSchemes.panelsType.second)
        {
        case 0:
            oldAirfoil.emplace_back(new AirfoilRect(*afl));
            break;
        //case 1:
        //  oldAirfoil.emplace_back(new AirfoilCurv(*afl));
        //  break;
        }


#ifdef BRIDGE
        if (afl->numberInPassport == 0)
        {
            mechanics[afl->numberInPassport]->Move();
        }
        else
        {
            Mechanics& mechGen = *mechanics[0];
            MechanicsRigidOscillPart& mech0 = dynamic_cast<MechanicsRigidOscillPart&>(mechGen);
            double dy = mech0.getY() - mech0.getYOld();
            double du = mech0.getU() - mech0.getUOld();
            
            Mechanics& mechGenI = *mechanics[afl->numberInPassport];
            MechanicsRigidOscillPart& mechI = dynamic_cast<MechanicsRigidOscillPart&>(mechGenI);
            
            mechI.getUOld() = mechI.getU();
            mechI.getYOld() = mechI.getY();
            
            //std::cout << "afl = " << afl << ", dy = " << dy << std::endl;
            
            airfoil[afl->numberInPassport]->Move({ 0.0, dy });
            mechI.Vcm[1] += du;
            
            mechI.getY() += dy;
            mechI.getU() += du;
        }
#endif

#ifdef INITIAL
        mechanics[afl->numberInPassport]->Move();       
#endif
    }//for

    for (auto& bou : boundary)
        bou->virtualWake.vtx.clear();

    getTimestat().timeOther.second += omp_get_wtime();

    CheckInside(newPos, oldAirfoil);
    
    getTimestat().timeOther.first += omp_get_wtime();   

    //передача новых положений вихрей в пелену  
    for (size_t i = 0; i < wake->vtx.size(); ++i)
        wake->vtx[i].r() = newPos[i];

//  getWake().SaveKadrVtk();
    getTimestat().timeOther.second += omp_get_wtime();
}//WakeAndAirfoilsMotion()

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

std::vector<std::unique_ptr<Airfoil> > VM2D::World2D::airfoil

private

Список умных указателей на обтекаемые профили

Definition at line 72 of file World2D.h.

std::vector<std::unique_ptr<Boundary> > VM2D::World2D::boundary

private

Список умных указателей на формирователи граничных условий на профилях

Definition at line 78 of file World2D.h.

Gpu VM2D::World2D::cuda

mutableprivate

Объект, управляющий графическим ускорителем

Definition at line 123 of file World2D.h.

size_t VMlib::WorldGen::currentStep

inherited

Текущий номер шага в решаемой задаче

Definition at line 69 of file WorldGen.h.

std::vector<size_t> VM2D::World2D::dispBoundaryInSystem

private

Список номеров, с которых начинаются элементы правой части (или матрицы) системы для профилей

Definition at line 81 of file World2D.h.

LogStream VMlib::WorldGen::info

mutableprotectedinherited

Поток для вывода логов и сообщений об ошибках

Definition at line 59 of file WorldGen.h.

Eigen::MatrixXd VM2D::World2D::invMatr

private

Обратная матрица

Definition at line 108 of file World2D.h.

std::vector<std::vector<std::pair<Eigen::MatrixXd, Eigen::MatrixXd> > > VM2D::World2D::IQ

private

Матрица, состоящая из пар матриц, в которых хранятся касательные и нормальные компоненты интегралов от ядра

Definition at line 105 of file World2D.h.

Eigen::MatrixXd VM2D::World2D::matr

private

Матрица системы

Definition at line 102 of file World2D.h.

std::unique_ptr<MeasureVP> VM2D::World2D::measureVP

private

Умный указатель на алгоритм вычисления полей скоростей и давления (для сохранения в файл)

Definition at line 96 of file World2D.h.

std::vector<std::unique_ptr<Mechanics> > VM2D::World2D::mechanics

private

Список умных указателей на типы механической системы для каждого профиля

Definition at line 84 of file World2D.h.

std::vector<std::unique_ptr<Airfoil> > VM2D::World2D::oldAirfoil

private

Список умных указателей на обтекаемые профили для сохранения старого положения

Definition at line 75 of file World2D.h.

const Passport& VM2D::World2D::passport

private

Константная ссылка на паспорт конкретного расчета

Definition at line 120 of file World2D.h.

const PassportGen& VMlib::WorldGen::passportGen

protectedinherited

Константная ссылка на паспорт конкретного расчета

Definition at line 62 of file WorldGen.h.

Eigen::VectorXd VM2D::World2D::rhs

private

Правая часть системы

Definition at line 114 of file World2D.h.

Eigen::VectorXd VM2D::World2D::sol

private

Решение системы

Definition at line 117 of file World2D.h.

std::unique_ptr<WakeDataBase> VM2D::World2D::source

private

Умный указатель на источники

Definition at line 93 of file World2D.h.

std::unique_ptr<TimesGen> VMlib::WorldGen::timestat

protectedinherited

Сведения о временах выполнения основных операций

Definition at line 65 of file WorldGen.h.

bool VM2D::World2D::useInverseMatrix

private

Признак использования обратной матрицы

Definition at line 111 of file World2D.h.

std::unique_ptr<Velocity> VM2D::World2D::velocity

private

Умный укзатель на объект, определяющий методику вычисления скоростей

Definition at line 87 of file World2D.h.

std::unique_ptr<Wake> VM2D::World2D::wake

private

Умный указатель на вихревой след

Definition at line 90 of file World2D.h.

---

The documentation for this class was generated from the following files:

• VM2D/World2D/World2D.h
• VM2D/World2D/World2D.cpp