gammaCirc.h File Reference

Go to the source code of this file.

Macros
#define	radiuses   { 0.5, 1.0, 1.5, 2.0, 2.5 }
#define	scales   { 1.0 }

Functions
	if (currentStep % 1==0)

Macro Definition Documentation

radiuses

#define radiuses   { 0.5, 1.0, 1.5, 2.0, 2.5 }

Definition at line 15 of file gammaCirc.h.

scales

#define scales   { 1.0 }

Definition at line 19 of file gammaCirc.h.

Function Documentation

if()

if

(

currentStep %

1 = = 0

)

Definition at line 22 of file gammaCirc.h.

{
 
    std::vector<double> rad = radiuses;
    std::vector<double> scale = scales;
 
 
    std::string circFileName = getPassport().dir + "circ";
 
 
    if (currentStep == 0)
    {
        std::ofstream circFile(circFileName.c_str());
 
        VMlib::PrintLogoToTextFile(circFile, circFileName, "Circulations along the circles of fixed radius");
 
        std::stringstream ssCirc;
        ssCirc << "step  time  scale";
        for (auto r : rad)
            ssCirc << "  R=" << r;
 
        VMlib::PrintHeaderToTextFile(circFile, ssCirc.str());
 
        circFile << std::endl;
        circFile.close();
        circFile.clear();
    }
 
    for (int c = 0; c < scale.size(); ++c)
    {
        std::vector<double> circRad(rad.size(), 0.0);
 
        int nPts = 50; //50
 
        double delta = 6.0; //6.0
 
        double R = 0.0;
        double L = 0.0;
        double epsAst = 0.0;
        double gam = 0.0;
        double curCirc = 0.0;
 
        for (size_t i = 0; i < rad.size(); ++i)
        {
            R = rad[i];
            curCirc = 0.0;
 
#pragma omp parallel for default(none) private(L,gam,epsAst) shared(R,nPts,delta,c,scale) reduction(+:curCirc) schedule(dynamic,100)
            for (int j = 0; j < (int)wake->vtx.size(); ++j)
            {
                L = wake->vtx[j].r().length();
                gam = wake->vtx[j].g();
                epsAst = scale[c] * velocity->wakeVortexesParams.epsastWake[j];
 
                if (R - L > delta * epsAst)
                    curCirc += gam;
                else if (R - L > -delta * epsAst)
                {                   
                    double h = fabs(R - L);
                    if (h / R > 1e-3)
                    {
                        double phi = atan(h / sqrt(delta * delta * epsAst * epsAst - h * h));
                        double deltaPhi = (PI - 2.0 * phi) / nPts;
 
                        double addCirc = 0.0;
 
                        auto func = [h, epsAst, delta](double phi) -> double
                            {
                                return (IDPI / epsAst) * (exp(h / (epsAst * sin(phi))) * (epsAst - h / sin(phi)) - exp(-delta) * (1 + delta));
                            };
 
                        for (int s = 0; s < nPts; ++s)
                        {
                            double valAngle = -phi - deltaPhi * (s + 0.5);
 
                            addCirc += func(valAngle) * deltaPhi;
                        }
                        addCirc *= gam;
 
                        if (L > R)
                            curCirc += addCirc;
                        else
                            curCirc += (gam - addCirc);
                    }
                    else
                        curCirc += 0.5 * gam;
                }//else
            }//for j
 
            circRad[i] = curCirc;
        }//for i
 
 
        std::ofstream circFile(circFileName.c_str(), std::ios::app);
 
        circFile << currentStep << " " << getCurrentTime() << " " << scale[c];
 
        for (auto gamma : circRad)
            circFile << " " << gamma;
        circFile << std::endl;
 
 
        circFile.close();
        circFile.clear();
 
    }// for c
 
}

Here is the call graph for this function: