MeasureVP2D.h

Go to the documentation of this file.

/*--------------------------------*- VM2D -*-----------------*---------------*\
| ##  ## ##   ##  ####  #####   |                            | Version 1.14   |
| ##  ## ### ### ##  ## ##  ##  |  VM2D: Vortex Method       | 2026/03/06     |
| ##  ## ## # ##    ##  ##  ##  |  for 2D Flow Simulation    *----------------*
|  ####  ##   ##   ##   ##  ##  |  Open Source Code                           |
|   ##   ##   ## ###### #####   |  https://www.github.com/vortexmethods/VM2D  |
|                                                                             |
| Copyright (C) 2017-2026 I. Marchevsky, K. Sokol, E. Ryatina, A. Kolganova   |
*-----------------------------------------------------------------------------*
| File name: MeasureVP2D.h                                                    |
| Info: Source code of VM2D                                                   |
|                                                                             |
| This file is part of VM2D.                                                  |
| VM2D is free software: you can redistribute it and/or modify it             |
| under the terms of the GNU General Public License as published by           |
| the Free Software Foundation, either version 3 of the License, or           |
| (at your option) any later version.                                         |
|                                                                             |
| VM2D 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 General Public License       |
| for more details.                                                           |
|                                                                             |
| You should have received a copy of the GNU General Public License           |
| along with VM2D.  If not, see <http://www.gnu.org/licenses/>.               |
\*---------------------------------------------------------------------------*/
 
 
#ifndef MEASUREVP2D_H
#define MEASUREVP2D_H
 
#include <memory>
 
#include "defs.h"
 
namespace VM2D
{
    class WakeDataBase;
    class World2D;
 
    class MeasureVP
    {
    private:
        std::vector<Point2D> initialPoints;
 
        std::vector<Point2D> historyPoints;
        
        std::vector<std::ostringstream> sstr;
        int sstrCounter;
        int historyFrequency;
 
 
    public:
        std::vector<Point2D> elasticPoints;
 
    //private:
        std::unique_ptr<WakeDataBase> wakeVP;
 
        std::vector<Point2D> velocity;
 
        std::vector<double> domainRadius;
 
        std::vector<double> pressure;
 
        //const double pInf = 1.0;
 
    protected:
        const World2D& W;
 
    public:
        MeasureVP(const World2D& W_);
 
        ~MeasureVP() { };
 
        void ReadPointsFromFile(const std::string& dir);
        void SetPoints(const std::vector<Point2D>& points, const std::vector<Point2D>& history);
 
        void Initialization();
 
        void CalcPressure();
#ifdef USE_CUDA
        void GPUCalcPressure();
#endif
 
        void SaveVP();
 
        std::vector<std::pair<Point2D, double>> GetVPinElasticPoints();
 
        const WakeDataBase& getWakeVP() const { return *wakeVP; };
 
        const std::vector<Point2D>& getVelocity() const { return velocity; };
 
        std::vector<Point2D>& getNonConstVelocity() { return velocity; };
 
        const std::vector<double>& getDomainRadius() const { return domainRadius; };
 
        std::vector<double>& getNonConstDomainRadius() { return domainRadius; };
 
        size_t getTotalNumberOfRealPoints() const { return initialPoints.size() + historyPoints.size(); };
 
#if defined(USE_CUDA)
        mutable double* devPressurePtr;
#endif
    };
 
}//namespace VM2D
 
#endif