Gmres2D.cpp File Reference

���� ���� � ��������� ������� ������ GMRES. More...

#include "Gmres2D.h"
#include "Airfoil2D.h"
#include "Boundary2D.h"
#include "MeasureVP2D.h"
#include "Mechanics2D.h"
#include "Preprocessor.h"
#include "StreamParser.h"
#include "Velocity2D.h"
#include "Wake2D.h"
#include "World2D.h"

Include dependency graph for Gmres2D.cpp:

Go to the source code of this file.

Typedefs
typedef Eigen::Map< const Eigen::VectorXd >	MapVec

Functions
void	PreCalculateCoef (const Airfoil &aflP, size_t nPanelsP, std::vector< Point2D > &prea, std::vector< Point2D > &prec, std::vector< Point2D > &prea1, std::vector< Point2D > &prec1, bool linScheme)

Detailed Description

���� ���� � ��������� ������� ������ GMRES.

Author

���������� ���� ��������������

����� ������ ���������

������ ������� ��������

��������� ���������� ��������

������� ����� ������� \Version 1.14

Date

6 ����� 2026 �.

Definition in file Gmres2D.cpp.

Typedef Documentation

MapVec

typedef Eigen::Map<const Eigen::VectorXd> MapVec

Definition at line 57 of file Gmres2D.cpp.

Function Documentation

PreCalculateCoef()

void PreCalculateCoef	(	const Airfoil &	aflP,
		size_t	nPanelsP,
		std::vector< Point2D > &	prea,
		std::vector< Point2D > &	prec,
		std::vector< Point2D > &	prea1,
		std::vector< Point2D > &	prec1,
		bool	linScheme
	)

Definition at line 528 of file Gmres2D.cpp.

{
    Point2D p1, s1, p2, s2, di, dj;
    numvector<double, 3> alpha, lambda;
    numvector<Point2D, 3> v00, v11;
    Point2D i00, i11;
 
    for (size_t i = 0; i < nPanelsP; ++i)
    {
        std::array<int, 2> vecV = { ((int)i - 1 >= 0) ? int(i - 1) : int(nPanelsP - 1), (i + 1 < nPanelsP) ? int(i + 1) : 0 };
        for (int j : vecV)
        {
            p1 = aflP.getR(i + 1) - aflP.getR(j + 1);
            s1 = aflP.getR(i + 1) - aflP.getR(j);
            p2 = aflP.getR(i) - aflP.getR(j + 1);
            s2 = aflP.getR(i) - aflP.getR(j);
            di = aflP.getR(i + 1) - aflP.getR(i);
            dj = aflP.getR(j + 1) - aflP.getR(j);
 
            alpha = { \
                (aflP.isAfter(j, i)) ? 0.0 : VMlib::Alpha(s2, s1), \
                VMlib::Alpha(s2, p1), \
                (aflP.isAfter(i, j)) ? 0.0 : VMlib::Alpha(p1, p2) \
            };
 
            lambda = { \
                (aflP.isAfter(j, i)) ? 0.0 : VMlib::Lambda(s2, s1), \
                VMlib::Lambda(s2, p1), \
                (aflP.isAfter(i, j)) ? 0.0 : VMlib::Lambda(p1, p2) \
            };
 
            v00 = {
                VMlib::Omega(s1, di.unit(), dj.unit()),
                -VMlib::Omega(di, di.unit(), dj.unit()),
                VMlib::Omega(p2, di.unit(), dj.unit())
            };
 
            i00 = IDPI / di.length() * (-(alpha[0] * v00[0] + alpha[1] * v00[1] + alpha[2] * v00[2]).kcross() \
                + (lambda[0] * v00[0] + lambda[1] * v00[1] + lambda[2] * v00[2]));
 
            if (aflP.isAfter(j, i))
                prec[i] = (i00.kcross()) * (1.0 / dj.length());
            if (aflP.isAfter(i, j))
                prea[i] = (i00.kcross()) * (1.0 / dj.length());
 
            if (linScheme) {
 
                v11 = { 1.0 / (12.0 * di.length() * dj.length()) * \
                (2.0 * (s1 & Omega(s1 - 3.0 * p2, di.unit(), dj.unit())) * VMlib::Omega(s1, di.unit(), dj.unit()) - \
                    s1.length2() * (s1 - 3.0 * p2)) - 0.25 * VMlib::Omega(s1, di.unit(), dj.unit()),
                -di.length() / (12.0 * dj.length()) * Omega(di, dj.unit(), dj.unit()),
                -dj.length() / (12.0 * di.length()) * Omega(dj, di.unit(), di.unit()) };
 
                i11 = (IDPI / di.length()) * ((alpha[0] + alpha[2]) * v11[0] + (alpha[1] + alpha[2]) * v11[1] + alpha[2] * v11[2]\
                    + ((lambda[0] + lambda[2]) * v11[0] + (lambda[1] + lambda[2]) * v11[1] + lambda[2] * v11[2] \
                        + 1.0 / 12.0 * (dj.length() * di.unit() + di.length() * dj.unit() - 2.0 * VMlib::Omega(s1, di.unit(), dj.unit()))).kcross());
                if (aflP.isAfter(j, i))
                    prec1[i] = (i11) * (1.0 / dj.length());
                if (aflP.isAfter(i, j))
                    prea1[i] = (i11) * (1.0 / dj.length());
            }
        }
    }
}

Here is the call graph for this function: