VM2D::WakeDataBase Class Reference

Класс, опеделяющий набор вихрей More...

#include <WakeDataBase2D.h>

Inheritance diagram for VM2D::WakeDataBase:

Collaboration diagram for VM2D::WakeDataBase:

Public Member Functions
	WakeDataBase (const World2D &W_)
	Конструктор инициализации More...
virtual	~WakeDataBase ()
void	ReadFromFile (const std::string &dir, const std::string &fileName)
	Считывание вихревого следа из файла More...
void	SaveKadrVtk (const std::string &filePrefix="Kadr") const
	Сохранение вихревого следа в файл .vtk. More...

Public Attributes
const World2D &	W
	Константная ссылка на решаемую задачу More...
std::vector< Vortex2D >	vtx
	Список вихревых элементов More...

Detailed Description

Класс, опеделяющий набор вихрей

Author

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

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

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

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

Date

14 января 2024 г.

Definition at line 64 of file WakeDataBase2D.h.

Constructor & Destructor Documentation

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

inline

Конструктор инициализации

Parameters

[in]

W_

константная ссылка на решаемую задачу

Definition at line 74 of file WakeDataBase2D.h.

            : W(W_) 
        {};

virtual VM2D::WakeDataBase::~WakeDataBase ( )

inlinevirtual

Definition at line 105 of file WakeDataBase2D.h.

{};

Here is the call graph for this function:

Member Function Documentation

void WakeDataBase::ReadFromFile	(	const std::string &	dir,
		const std::string &	fileName
	)

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

Parameters

[in]	dir	константная ссылка на строку, задающую каталог, где лежит файл с вихревым следом
[in]	fileName	константная ссылка на строку, задающую имя файла с вихревым следом

Definition at line 56 of file WakeDataBase2D.cpp.

{
    std::string filename = dir + fileName;
    std::ifstream wakeFile;
    
    if (fileExistTest(filename, W.getInfo(), { "txt", "TXT" }))
    {
        std::stringstream wakeFile(VMlib::Preprocessor(filename).resultString);
        
        VMlib::LogStream XXX;
        VMlib::StreamParser wakeParser(XXX, "vortex wake file parser", wakeFile);
                
        wakeParser.get("vtx", vtx);
    }
}//ReadFromFile(...)

Here is the call graph for this function:

Here is the caller graph for this function:

void WakeDataBase::SaveKadrVtk

(

const std::string &

filePrefix = "Kadr"

)

const

Сохранение вихревого следа в файл .vtk.

Definition at line 73 of file WakeDataBase2D.cpp.

{
    W.getTimestat().timeSaveKadr.first += omp_get_wtime();

    if (W.ifDivisible(W.getPassport().timeDiscretizationProperties.saveVtxStep))
    {       
        std::ofstream outfile;
        size_t numberNonZero = 0;

        if (vtx.size() > 0)
            numberNonZero += vtx.size();
        else
            for (size_t q = 0; q < W.getNumberOfAirfoil(); ++q)
                numberNonZero += W.getAirfoil(q).getNumberOfPanels();
        VMlib::CreateDirectory(W.getPassport().dir, "snapshots");

        if (W.getPassport().timeDiscretizationProperties.fileTypeVtx.second == 0) //text format vtk
        {
            std::string fname = VMlib::fileNameStep(filePrefix, W.getPassport().timeDiscretizationProperties.nameLength, W.getCurrentStep(), "vtk");
            outfile.open(W.getPassport().dir + "snapshots/" + fname);

            outfile << "# vtk DataFile Version 2.0" << std::endl;
            outfile << "VM2D VTK result: " << (W.getPassport().dir + "snapshots/" + fname).c_str() << " saved " << VMlib::CurrentDataTime() << std::endl;
            outfile << "ASCII" << std::endl;
            outfile << "DATASET UNSTRUCTURED_GRID" << std::endl;
            outfile << "POINTS " << numberNonZero << " float" << std::endl;

            if (vtx.size() > 0)
                for (auto& v : vtx)
                {
                    const Point2D& r = v.r();
                    outfile << r[0] << " " << r[1] << " " << "0.0" << std::endl;
                }//for v        
            else
                for (size_t q = 0; q < W.getNumberOfAirfoil(); ++q)
                    for (size_t s = 0; s < W.getAirfoil(q).getNumberOfPanels(); ++s)
                    {
                        const Point2D& r = W.getAirfoil(q).getR(s);
                        outfile << r[0] << " " << r[1] << " " << "0.0" << std::endl;
                    }

            outfile << "CELLS " << numberNonZero << " " << 2 * numberNonZero << std::endl;
            for (size_t i = 0; i < numberNonZero; ++i)
                outfile << "1 " << i << std::endl;

            outfile << "CELL_TYPES " << numberNonZero << std::endl;
            for (size_t i = 0; i < numberNonZero; ++i)
                outfile << "1" << std::endl;

            outfile << std::endl;
            outfile << "POINT_DATA " << numberNonZero << std::endl;
            outfile << "SCALARS Gamma float 1" << std::endl;
            outfile << "LOOKUP_TABLE default" << std::endl;

            if (vtx.size() > 0)
                for (auto& v : vtx)
                    outfile << v.g() << std::endl;
            else
                for (size_t q = 0; q < W.getNumberOfAirfoil(); ++q)
                    for (size_t s = 0; s < W.getAirfoil(q).getNumberOfPanels(); ++s)
                    {
                        //const Point2D& r = W.getAirfoil(q).getR(s);
                        outfile << "0.0" << std::endl;
                    }

            outfile.close();
        }//if fileType = text
        else if (W.getPassport().timeDiscretizationProperties.fileTypeVtx.second == 1) //binary format vtk
        {
            //Тест способа хранения чисел
            uint16_t x = 0x0001;
            bool littleEndian = (*((uint8_t*)&x));
            const char eolnBIN[] = "\n";

            std::string fname = VMlib::fileNameStep(filePrefix, W.getPassport().timeDiscretizationProperties.nameLength, W.getCurrentStep(), "vtk");
            outfile.open(W.getPassport().dir + "snapshots/" + fname, std::ios::out | std::ios::binary);

            outfile << "# vtk DataFile Version 3.0" << "\r\n" << "VM2D VTK result: " << (W.getPassport().dir + "snapshots/" + fname).c_str() << " saved " << VMlib::CurrentDataTime() << eolnBIN;
            outfile << "BINARY" << eolnBIN;
            outfile << "DATASET UNSTRUCTURED_GRID" << eolnBIN << "POINTS " << numberNonZero << " " << "float" << eolnBIN;


            if (vtx.size() > 0)
            {
                Eigen::VectorXf rData = Eigen::VectorXf::Zero(vtx.size() * 3);
                for (size_t i = 0; i < vtx.size(); ++i)
                {
                    rData(3 * i) = (float)(vtx[i].r())[0];
                    rData(3 * i + 1) = (float)(vtx[i].r())[1];
                }//for i

                if (littleEndian)
                    for (int i = 0; i < vtx.size() * 3; ++i)
                        VMlib::SwapEnd(rData(i));
                outfile.write(reinterpret_cast<char*>(rData.data()), vtx.size() * 3 * sizeof(float));                   
            }
            else
                for (size_t q = 0; q < W.getNumberOfAirfoil(); ++q)
                {
                    Eigen::VectorXf rData = Eigen::VectorXf::Zero(W.getAirfoil(q).getNumberOfPanels() * 3);
                    for (size_t s = 0; s < W.getAirfoil(q).getNumberOfPanels(); ++s)
                    {
                        rData(3 * s) = (float)(W.getAirfoil(q).getR(s))[0];
                        rData(3 * s + 1) = (float)(W.getAirfoil(q).getR(s))[1];
                    }//for i

                    if (littleEndian)
                        for (int i = 0; i < W.getAirfoil(q).getNumberOfPanels() * 3; ++i)
                            VMlib::SwapEnd(rData(i));
                    outfile.write(reinterpret_cast<char*>(rData.data()), W.getAirfoil(q).getNumberOfPanels() * 3 * sizeof(float));
                }

            // CELLS
            std::vector<int> cells(2 * numberNonZero);
            for (size_t i = 0; i < numberNonZero; ++i)
            {
                cells[2 * i] = 1;
                cells[2 * i + 1] = (int)i;
            }

            std::vector<int> cellsTypes;
            cellsTypes.resize(numberNonZero, 1);

            if (littleEndian)
            {
                for (int i = 0; i < numberNonZero * 2; ++i)
                    VMlib::SwapEnd(cells[i]);

                for (int i = 0; i < numberNonZero; ++i)
                    VMlib::SwapEnd(cellsTypes[i]);
            }

            outfile << eolnBIN << "CELLS " << numberNonZero << " " << numberNonZero * 2 << eolnBIN;
            outfile.write(reinterpret_cast<char*>(cells.data()), numberNonZero * 2 * sizeof(int));
            outfile << eolnBIN << "CELL_TYPES " << numberNonZero << eolnBIN;
            outfile.write(reinterpret_cast<char*>(cellsTypes.data()), numberNonZero * sizeof(int));

            //gammas
            outfile << eolnBIN << "POINT_DATA " << numberNonZero << eolnBIN;
            outfile << eolnBIN << "SCALARS Gamma " << "float" << " 1" << eolnBIN;
            outfile << "LOOKUP_TABLE default" << eolnBIN;           

            if (vtx.size() > 0)
            {
                Eigen::VectorXf pData = Eigen::VectorXf::Zero(numberNonZero);
                for (int s = 0; s < vtx.size(); ++s)
                    pData(s) = (float)vtx[s].g();
                
                if (littleEndian)
                    for (int i = 0; i < vtx.size(); ++i)
                        VMlib::SwapEnd(pData(i));
                outfile.write(reinterpret_cast<char*>(pData.data()), vtx.size() * sizeof(float));
            }
            else
                for (size_t q = 0; q < W.getNumberOfAirfoil(); ++q)
                {
                    Eigen::VectorXf pData = Eigen::VectorXf::Zero(W.getAirfoil(q).getNumberOfPanels());
                    for (int s = 0; s < W.getAirfoil(q).getNumberOfPanels(); ++s)
                        pData(s) = 0;

                    if (littleEndian)
                        for (int i = 0; i < W.getAirfoil(q).getNumberOfPanels(); ++i)
                            VMlib::SwapEnd(pData(i));
                    outfile.write(reinterpret_cast<char*>(pData.data()), W.getAirfoil(q).getNumberOfPanels() * sizeof(float));                  
                }
            
            outfile << eolnBIN;
            outfile.close();
        }//if binary
        if (W.getPassport().timeDiscretizationProperties.fileTypeVtx.second == 2) //csv
        {
            std::string fname = VMlib::fileNameStep(filePrefix, W.getPassport().timeDiscretizationProperties.nameLength, W.getCurrentStep(), "csv");
            outfile.open(W.getPassport().dir + "snapshots/" + fname);

            outfile << "point,x,y,G " << std::endl;

            int counter = 0;
            if (vtx.size() > 0)
                for (auto& v : vtx)                                 
                    outfile << counter++ << "," << v.r()[0] << "," << v.r()[1] << "," << v.g() << std::endl;                
            else
                for (size_t q = 0; q < W.getNumberOfAirfoil(); ++q)
                    for (size_t s = 0; s < W.getAirfoil(q).getNumberOfPanels(); ++s)
                    {
                        const Point2D& r = W.getAirfoil(q).getR(s);
                        outfile << counter++ << "," << r[0] << "," << r[1] << "," << "0.0" << std::endl;
                    }
            outfile.close();
        }//if fileType = text

    }

    W.getTimestat().timeSaveKadr.second += omp_get_wtime();
}

Member Data Documentation

std::vector<Vortex2D> VM2D::WakeDataBase::vtx

Список вихревых элементов

Definition at line 76 of file WakeDataBase2D.h.

const World2D& VM2D::WakeDataBase::W

Константная ссылка на решаемую задачу

Definition at line 69 of file WakeDataBase2D.h.

---

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

• VM2D/WakeDataBase2D/WakeDataBase2D.h
• VM2D/WakeDataBase2D/WakeDataBase2D.cpp