Files
Population-optimization-alg…/MQL5/Include/Math/AOs/TestStand3D.mqh
T

610 lines
45 KiB
Plaintext
Raw Normal View History

2026-03-15 16:01:30 +04:00
//+------------------------------------------------------------------+
//| TestStand3D.mqh |
//| 3D-панель визуализации для тест-стенда метаэвристик |
//| |
//| Порядок #include в скрипте: |
//| 1. TestStandFunctions.mqh — C_TestStand (2D-стенд) |
//| 2. TestFunctions.mqh — C_Function, GenerateDataFixedSize|
//| 3. TestStand3D.mqh — этот файл |
//| 4. алгоритм оптимизации |
//+------------------------------------------------------------------+
#include <Canvas\Canvas3D.mqh>
#include <Canvas\DX\DXSurface.mqh>
#include <Math\AOs\TestFunctions.mqh>
//--- Параметры камеры (можно менять)
#define S3D_DIST_DEF 35.0 // дистанция по умолчанию
#define S3D_DIST_MIN 12.0 // минимальный zoom-in
#define S3D_DIST_MAX 80.0 // максимальный zoom-out
#define S3D_PITCH_DEF 0.45f // наклон камеры по умолчанию (~26°)
#define S3D_FOV 0.5236f // поле зрения = π/6 (30°)
//--- Полуразмеры поверхности в мировых координатах
#define S3D_SX 10.0f
#define S3D_SY 5.0f
#define S3D_SZ 10.0f
//--- Фоновый цвет 3D-панели
#define S3D_BG 0xFF0D0D1A
//+------------------------------------------------------------------+
//| S3DTrailPt — точка трека лучшего решения (мировые координаты) |
//+------------------------------------------------------------------+
struct S3DTrailPt
{
float wx, wy, wz;
};
//+------------------------------------------------------------------+
//| S3DAgent — агент популяции (мировые координаты + цвет) |
//+------------------------------------------------------------------+
struct S3DAgent
{
float wx, wy, wz;
uint clr; // ARGB
bool is_best; // лучший агент — особая отрисовка
};
//+------------------------------------------------------------------+
//| C_TestStand3D |
//| |
//| Отображает левую (H×H) панель тест-стенда как интерактивную |
//| 3D-поверхность функции с проекцией агентов и треком лучшего. |
//| |
//| Управление (EA-контекст, события через OnChartEvent): |
//| LMB-drag — вращение камеры |
//| Scroll — zoom |
//| Double-click — сброс камеры к позиции по умолчанию |
//+------------------------------------------------------------------+
class C_TestStand3D
{
//--- DX-объекты
CCanvas3D m_cv;
CDXSurface m_surf;
//--- параметры канваса
string m_name;
int m_x, m_y, m_w, m_h;
bool m_surf_ok;
bool m_vis;
bool m_shutdown_done;
int m_grid;
//--- параметры камеры
float m_pitch; // наклон (вокруг X)
float m_yaw; // азимут (вокруг Y)
double m_dist; // расстояние от центра сцены
int m_mx, m_my; // предыдущая позиция мыши (локальная)
//--- матрица вид×проекция для ручного проецирования агентов на экран
DXMatrix m_vp;
//--- диапазон домена функции (для маппинга агентов в мировые координаты)
double m_fxMin, m_fxMax, m_fyMin, m_fyMax;
//--- агенты текущей эпохи
S3DAgent m_ag [];
int m_nAg;
//--- трек лучшего решения (кольцевой буфер)
S3DTrailPt m_trail [];
int m_trailLen;
int m_trailMax;
public:
C_TestStand3D ();
~C_TestStand3D () { Shutdown (); }
//--- жизненный цикл
bool Init (const string name, int x, int y, int w, int h);
void Shutdown ();
//--- контент
// grid — вершин по одной оси (10..200); вызывать при смене функции
bool BuildSurface (C_Function &func, int grid = 80);
// args[] = [x0,y0, x1,y1, …] (count пар); best_idx = -1 → без выделения
void SetAgents (double &args [], int count,
C_Function &func, int best_idx = -1);
// один кадр рендера
void Redraw ();
// сброс трека лучшего (вызывать в начале каждого нового теста)
void ResetTrail () { m_trailLen = 0; }
//--- события (маршрутизировать из OnChartEvent в EA)
void OnMouseMove (int chart_x, int chart_y, uint flags);
void OnMouseWheel (double delta);
void OnDblClick ();
// dt в секундах; в скрипте передавать фиксированный виртуальный шаг
void OnTimer (double dt);
//--- видимость и размер
void Show (bool v);
bool IsVisible () const { return m_vis; }
string GetName () const { return m_name; }
void Resize (int new_w, int new_h);
private:
//--- пересчёт матриц камеры и VP
void _CamUpdate ();
//--- проекция мировой точки → экранные пиксели; false = за кадром
bool _Project (float wx, float wy, float wz, int &sx, int &sy);
//--- отрисовка трека лучшего, агентов
void _DrawTrail ();
void _DrawAgents ();
//--- матричные хелперы (portable, без внешних DX-функций)
// Внимание: DXMatrix в MQL5 хранит элементы как m[row][col]
// DXVector3 передаётся только по ссылке
static void _LookAtLH (DXMatrix &M, DXVector3 &eye, DXVector3 &tgt, DXVector3 &up);
static void _PerspLH (DXMatrix &M, float fov, float aspect, float zn, float zf);
static void _MulMat (DXMatrix &C, DXMatrix &A, DXMatrix &B);
//--- радужный цвет по индексу агента: h_deg ∈ [0,360] → ARGB (S=1, L=0.5)
static uint _HSL (int h_deg);
};
//+------------------------------------------------------------------+
C_TestStand3D::C_TestStand3D ()
: m_surf_ok (false), m_vis (false), m_shutdown_done (false),
m_nAg (0), m_grid (80),
m_pitch (S3D_PITCH_DEF), m_yaw (0.0f), m_dist (S3D_DIST_DEF),
m_mx (-1), m_my (-1),
m_fxMin (0), m_fxMax (1), m_fyMin (0), m_fyMax (1),
m_trailLen (0), m_trailMax (600)
{
ArrayResize (m_trail, m_trailMax);
}
//+------------------------------------------------------------------+
bool C_TestStand3D::Init (const string name, int x, int y, int w, int h)
{
m_name = name;
m_x = x; m_y = y; m_w = w; m_h = h;
ResetLastError ();
if (!m_cv.CreateBitmapLabel (m_name, m_x, m_y, m_w, m_h,
COLOR_FORMAT_ARGB_NORMALIZE))
{
PrintFormat ("C_TestStand3D::Init ошибка CreateBitmapLabel: %d",
GetLastError ());
return false;
}
m_cv.ProjectionMatrixSet (S3D_FOV, (float)m_w / m_h, 0.1f, 200.0f);
DXVector3 tgt = DXVector3 (0.0f, 0.0f, 0.0f);
DXVector3 up = DXVector3 (0.0f, 1.0f, 0.0f);
m_cv.ViewTargetSet (tgt);
m_cv.ViewUpDirectionSet (up);
m_cv.LightColorSet (DXColor (1.0f, 0.97f, 0.88f, 0.32f));
m_cv.AmbientColorSet (DXColor (0.82f, 0.88f, 1.0f, 0.72f));
_CamUpdate ();
ObjectSetInteger (0, m_name, OBJPROP_HIDDEN, true);
ObjectSetInteger (0, m_name, OBJPROP_SELECTABLE, false);
return true;
}
//+------------------------------------------------------------------+
void C_TestStand3D::Shutdown ()
{
if (m_shutdown_done) return;
m_shutdown_done = true;
m_surf.Shutdown ();
// m_cv.Destroy() не вызываем — деструктор CCanvas3D сделает сам.
// Двойной вызов Destroy() приводил к Abnormal termination.
ObjectDelete (0, m_name);
}
//+------------------------------------------------------------------+
bool C_TestStand3D::BuildSurface (C_Function &func, int grid)
{
m_grid = (int)MathMax (10.0, MathMin (200.0, (double)grid));
m_fxMin = func.GetMinRangeX (); m_fxMax = func.GetMaxRangeX ();
m_fyMin = func.GetMinRangeY (); m_fyMax = func.GetMaxRangeY ();
double buf [];
if (!GenerateDataFixedSize (m_grid, m_grid, func, buf))
{
Print ("C_TestStand3D::BuildSurface GenerateDataFixedSize failed");
return false;
}
const float VAL_RANGE = 1.0f; // C_Function нормализует значения в [0,1]
DXVector3 bMin = DXVector3 (-S3D_SX, -S3D_SY, -S3D_SZ);
DXVector3 bMax = DXVector3 (S3D_SX, S3D_SY, S3D_SZ);
DXVector2 uv = DXVector2 (1.0f, 1.0f);
if (!m_surf_ok)
{
if (!m_surf.Create (m_cv.DXDispatcher (), m_cv.InputScene (),
buf, (uint)m_grid, (uint)m_grid,
VAL_RANGE, bMin, bMax, uv,
CDXSurface::SF_TWO_SIDED | CDXSurface::SF_USE_NORMALS))
{
Print ("C_TestStand3D::BuildSurface CDXSurface::Create failed");
return false;
}
// CDXSurface::Create не принимает цветовую схему → сразу Update,
// иначе первая функция отображается без градиента (серая).
m_surf.Update (buf, (uint)m_grid, (uint)m_grid,
VAL_RANGE, bMin, bMax, uv,
CDXSurface::SF_TWO_SIDED | CDXSurface::SF_USE_NORMALS,
CDXSurface::CS_COLD_TO_HOT);
m_surf.SpecularColorSet (DXColor (1.0f, 1.0f, 1.0f, 0.55f));
m_cv.ObjectAdd (&m_surf);
m_surf_ok = true;
}
else
{
m_surf.Update (buf, (uint)m_grid, (uint)m_grid,
VAL_RANGE, bMin, bMax, uv,
CDXSurface::SF_TWO_SIDED | CDXSurface::SF_USE_NORMALS,
CDXSurface::CS_COLD_TO_HOT);
}
m_nAg = 0;
m_trailLen = 0; // трек сбрасывается при смене функции
return true;
}
//+------------------------------------------------------------------+
void C_TestStand3D::SetAgents (double &args [], int count,
C_Function &func, int best_idx)
{
m_nAg = MathMax (0, count);
if (m_nAg == 0) return;
ArrayResize (m_ag, m_nAg);
double rx = m_fxMax - m_fxMin; if (rx < 1e-9) rx = 1e-9;
double ry = m_fyMax - m_fyMin; if (ry < 1e-9) ry = 1e-9;
for (int i = 0; i < m_nAg; i++)
{
// клампинг в допустимый диапазон
double ax = MathMax (m_fxMin, MathMin (m_fxMax, args [i * 2]));
double ay = MathMax (m_fyMin, MathMin (m_fyMax, args [i * 2 + 1]));
// маппинг: домен функции → мировые координаты поверхности
float wx = (float)(-S3D_SX + (ax - m_fxMin) / rx * 2.0 * S3D_SX);
float wz = (float)(-S3D_SZ + (ay - m_fyMin) / ry * 2.0 * S3D_SZ);
// высота: нормализованное значение функции [0,1] → мировая Y + отступ
double fv = MathMax (0.0, MathMin (1.0, func.Core (ax, ay)));
float wy = (float)(-S3D_SY + fv * 2.0 * S3D_SY) + 0.22f;
m_ag [i].wx = wx;
m_ag [i].wy = wy;
m_ag [i].wz = wz;
m_ag [i].is_best = (i == best_idx);
// радужный цвет: синий (i=0) → красный (i=last)
m_ag [i].clr = _HSL ((int)(270.0 * i / MathMax (1, m_nAg - 1)));
// добавить позицию лучшего в трек
if (i == best_idx && best_idx >= 0)
{
if (m_trailLen < m_trailMax)
{
m_trail [m_trailLen].wx = wx;
m_trail [m_trailLen].wy = wy;
m_trail [m_trailLen].wz = wz;
m_trailLen++;
}
else
{
// буфер полон — сдвиг на 1, старейшая точка теряется
for (int t = 0; t < m_trailMax - 1; t++) m_trail [t] = m_trail [t + 1];
m_trail [m_trailMax - 1].wx = wx;
m_trail [m_trailMax - 1].wy = wy;
m_trail [m_trailMax - 1].wz = wz;
}
}
}
}
//+------------------------------------------------------------------+
void C_TestStand3D::Redraw ()
{
if (!m_vis || !m_surf_ok) return;
m_cv.Render (DX_CLEAR_COLOR | DX_CLEAR_DEPTH, S3D_BG);
_DrawTrail (); // трек под агентами
_DrawAgents ();
m_cv.Update ();
}
//+------------------------------------------------------------------+
void C_TestStand3D::_CamUpdate ()
{
// Вращаем базовую позицию (0, 0, -dist) вокруг X (pitch), затем Y (yaw)
DXVector4 cam = DXVector4 (0.0f, 0.0f, -(float)m_dist, 1.0f);
DXVector4 lgt = DXVector4 (0.25f, -0.45f, 1.0f, 0.0f);
DXMatrix rot;
DXMatrixRotationX (rot, m_pitch);
DXVec4Transform (cam, cam, rot);
DXVec4Transform (lgt, lgt, rot);
DXMatrixRotationY (rot, m_yaw);
DXVec4Transform (cam, cam, rot);
DXVec4Transform (lgt, lgt, rot);
DXVector3 eye = DXVector3 (cam);
DXVector3 ld = DXVector3 (lgt);
DXVector3 tgt = DXVector3 (0.0f, 0.0f, 0.0f);
DXVector3 up = DXVector3 (0.0f, 1.0f, 0.0f);
// Обновляем все параметры камеры, чтобы CCanvas3D пересчитал view-матрицу
m_cv.ViewPositionSet (eye);
m_cv.ViewTargetSet (tgt);
m_cv.ViewUpDirectionSet (up);
m_cv.LightDirectionSet (ld);
// Строим VP-матрицу для ручного проецирования агентов в _Project()
DXMatrix view, proj;
_LookAtLH (view, eye, tgt, up);
_PerspLH (proj, S3D_FOV, (float)m_w / m_h, 0.1f, 200.0f);
_MulMat (m_vp, view, proj);
}
//+------------------------------------------------------------------+
bool C_TestStand3D::_Project (float wx, float wy, float wz,
int &sx, int &sy)
{
DXVector4 p = DXVector4 (wx, wy, wz, 1.0f);
DXVec4Transform (p, p, m_vp);
if (p.w < 0.01f) return false;
float nx = p.x / p.w;
float ny = p.y / p.w;
if (nx < -1.1f || nx > 1.1f || ny < -1.1f || ny > 1.1f) return false;
sx = (int)((nx + 1.0f) * 0.5f * m_w);
sy = (int)((1.0f - ny) * 0.5f * m_h);
return true;
}
//+------------------------------------------------------------------+
//| _DrawTrail |
//| Золотая линия (#FFC828) по истории позиций лучшего агента. |
//| Градиент прозрачности: старые точки alpha≈0x30, новые — 0xE0. |
//| Двойная линия (смещение +1px по Y) — визуальная толщина. |
//+------------------------------------------------------------------+
void C_TestStand3D::_DrawTrail ()
{
if (m_trailLen < 2) return;
int px = 0, py = 0, cx = 0, cy = 0;
bool hasPrev = false;
for (int t = 0; t < m_trailLen; t++)
{
if (!_Project (m_trail [t].wx, m_trail [t].wy, m_trail [t].wz, cx, cy))
{
hasPrev = false; // разрыв при выходе точки за viewport
continue;
}
if (!hasPrev)
{
px = cx; py = cy;
hasPrev = true;
continue;
}
double ratio = (double)t / (double)(m_trailLen - 1);
uint alpha = (uint)(0x30 + ratio * 0xB0); // 0x30 → 0xE0
uint clr = (alpha << 24) | 0x00FFC828;
uint clrT = ((alpha >> 1) << 24) | 0x00FFC828; // тень +1px
m_cv.LineAA (px, py, cx, cy, clr);
m_cv.LineAA (px, py + 1, cx, cy + 1, clrT);
px = cx; py = cy;
}
}
//+------------------------------------------------------------------+
//| _DrawAgents |
//| Обычный агент: цветной кружок r=4 + чёрная обводка. |
//| Лучший агент: белый кружок r=7 + жирная тень. |
//+------------------------------------------------------------------+
void C_TestStand3D::_DrawAgents ()
{
for (int i = 0; i < m_nAg; i++)
{
int sx, sy;
if (!_Project (m_ag [i].wx, m_ag [i].wy, m_ag [i].wz, sx, sy)) continue;
if (m_ag [i].is_best)
{
m_cv.FillCircle (sx, sy, 7, 0xFFFFFFFF);
m_cv.Circle (sx, sy, 7, 0xFF000000);
m_cv.Circle (sx, sy, 8, 0xFF000000);
m_cv.Circle (sx, sy, 9, 0x88000000);
}
else
{
m_cv.FillCircle (sx, sy, 4, m_ag [i].clr);
m_cv.Circle (sx, sy, 4, 0xFF000000);
m_cv.Circle (sx, sy, 5, 0x55000000);
}
}
}
//+------------------------------------------------------------------+
void C_TestStand3D::OnMouseMove (int chart_x, int chart_y, uint flags)
{
int lx = chart_x - m_x;
int ly = chart_y - m_y;
if ((flags & 1) == 1) // LMB нажата
{
if (m_mx >= 0)
{
m_yaw += (lx - m_mx) / 240.0f;
m_pitch += (ly - m_my) / 240.0f;
// ограничение наклона ±86°
float pMax = (float)(DX_PI * 0.48);
if (m_pitch < -pMax) m_pitch = -pMax;
if (m_pitch > pMax) m_pitch = pMax;
_CamUpdate ();
}
m_mx = lx; m_my = ly;
}
else
{m_mx = -1; m_my = -1; }
}
//+------------------------------------------------------------------+
void C_TestStand3D::OnMouseWheel (double delta)
{
m_dist *= (1.0 - delta * 0.0012);
m_dist = MathMax (S3D_DIST_MIN, MathMin (S3D_DIST_MAX, m_dist));
_CamUpdate ();
}
//+------------------------------------------------------------------+
void C_TestStand3D::OnDblClick ()
{
m_pitch = S3D_PITCH_DEF;
m_yaw = 0.0f;
m_dist = S3D_DIST_DEF;
_CamUpdate ();
}
//+------------------------------------------------------------------+
void C_TestStand3D::OnTimer (double dt)
{
// зарезервировано для будущей реализации авто-вращения
}
//+------------------------------------------------------------------+
void C_TestStand3D::Show (bool v)
{
m_vis = v;
ObjectSetInteger (0, m_name, OBJPROP_HIDDEN, !v);
ChartRedraw ();
}
//+------------------------------------------------------------------+
void C_TestStand3D::Resize (int new_w, int new_h)
{
if (new_w == m_w && new_h == m_h) return;
m_w = new_w; m_h = new_h;
m_cv.Resize (m_w, m_h);
DXContextSetSize (m_cv.DXContext (), m_w, m_h);
m_cv.ProjectionMatrixSet (S3D_FOV, (float)m_w / m_h, 0.1f, 200.0f);
_CamUpdate ();
}
//══════════════════════════════════════════════════════════════════
// Матричные хелперы
//
// DXMatrix в MQL5 хранит элементы как m[row][col], НЕ как _11/_12
// DXVector3 передаётся только по ссылке (объекты — только ref в MQL5)
// Все матрицы — row-major, left-handed (DirectX конвенция)
//══════════════════════════════════════════════════════════════════
//--- LookAt LH: эквивалент D3DXMatrixLookAtLH
void C_TestStand3D::_LookAtLH (DXMatrix &M, DXVector3 &eye, DXVector3 &tgt, DXVector3 &up)
{
// z = normalize(tgt - eye) — ось «вперёд»
float zx = tgt.x - eye.x, zy = tgt.y - eye.y, zz = tgt.z - eye.z;
float zl = (float)MathSqrt (zx * zx + zy * zy + zz * zz);
if (zl < 1e-7f) zl = 1e-7f;
zx /= zl; zy /= zl; zz /= zl;
// x = normalize(cross(up, z)) — ось «вправо»
float xx = up.y * zz - up.z * zy;
float xy = up.z * zx - up.x * zz;
float xz = up.x * zy - up.y * zx;
float xl = (float)MathSqrt (xx * xx + xy * xy + xz * xz);
if (xl < 1e-7f) xl = 1e-7f;
xx /= xl; xy /= xl; xz /= xl;
// y = cross(z, x) — ортогонализованная ось «вверх»
float yx = zy * xz - zz * xy;
float yy = zz * xx - zx * xz;
float yz = zx * xy - zy * xx;
// трансляционная часть (dot-products позиции камеры с осями)
float dx = -(xx * eye.x + xy * eye.y + xz * eye.z);
float dy = -(yx * eye.x + yy * eye.y + yz * eye.z);
float dz = -(zx * eye.x + zy * eye.y + zz * eye.z);
M.m [0] [0] = xx; M.m [0] [1] = yx; M.m [0] [2] = zx; M.m [0] [3] = 0.0f;
M.m [1] [0] = xy; M.m [1] [1] = yy; M.m [1] [2] = zy; M.m [1] [3] = 0.0f;
M.m [2] [0] = xz; M.m [2] [1] = yz; M.m [2] [2] = zz; M.m [2] [3] = 0.0f;
M.m [3] [0] = dx; M.m [3] [1] = dy; M.m [3] [2] = dz; M.m [3] [3] = 1.0f;
}
//--- Perspective LH: эквивалент D3DXMatrixPerspectiveFovLH
void C_TestStand3D::_PerspLH (DXMatrix &M,
float fov, float aspect, float zn, float zf)
{
float ys = 1.0f / (float)MathTan (fov * 0.5);
float xs = ys / aspect;
float zr = zf / (zf - zn);
M.m [0] [0] = xs; M.m [0] [1] = 0.0f; M.m [0] [2] = 0.0f; M.m [0] [3] = 0.0f;
M.m [1] [0] = 0.0f; M.m [1] [1] = ys; M.m [1] [2] = 0.0f; M.m [1] [3] = 0.0f;
M.m [2] [0] = 0.0f; M.m [2] [1] = 0.0f; M.m [2] [2] = zr; M.m [2] [3] = 1.0f;
M.m [3] [0] = 0.0f; M.m [3] [1] = 0.0f; M.m [3] [2] = -zn * zr; M.m [3] [3] = 0.0f;
}
//--- C = A × B (row-major)
void C_TestStand3D::_MulMat (DXMatrix &C, DXMatrix &A, DXMatrix &B)
{
C.m [0] [0] = A.m [0] [0] * B.m [0] [0] + A.m [0] [1] * B.m [1] [0] + A.m [0] [2] * B.m [2] [0] + A.m [0] [3] * B.m [3] [0];
C.m [0] [1] = A.m [0] [0] * B.m [0] [1] + A.m [0] [1] * B.m [1] [1] + A.m [0] [2] * B.m [2] [1] + A.m [0] [3] * B.m [3] [1];
C.m [0] [2] = A.m [0] [0] * B.m [0] [2] + A.m [0] [1] * B.m [1] [2] + A.m [0] [2] * B.m [2] [2] + A.m [0] [3] * B.m [3] [2];
C.m [0] [3] = A.m [0] [0] * B.m [0] [3] + A.m [0] [1] * B.m [1] [3] + A.m [0] [2] * B.m [2] [3] + A.m [0] [3] * B.m [3] [3];
C.m [1] [0] = A.m [1] [0] * B.m [0] [0] + A.m [1] [1] * B.m [1] [0] + A.m [1] [2] * B.m [2] [0] + A.m [1] [3] * B.m [3] [0];
C.m [1] [1] = A.m [1] [0] * B.m [0] [1] + A.m [1] [1] * B.m [1] [1] + A.m [1] [2] * B.m [2] [1] + A.m [1] [3] * B.m [3] [1];
C.m [1] [2] = A.m [1] [0] * B.m [0] [2] + A.m [1] [1] * B.m [1] [2] + A.m [1] [2] * B.m [2] [2] + A.m [1] [3] * B.m [3] [2];
C.m [1] [3] = A.m [1] [0] * B.m [0] [3] + A.m [1] [1] * B.m [1] [3] + A.m [1] [2] * B.m [2] [3] + A.m [1] [3] * B.m [3] [3];
C.m [2] [0] = A.m [2] [0] * B.m [0] [0] + A.m [2] [1] * B.m [1] [0] + A.m [2] [2] * B.m [2] [0] + A.m [2] [3] * B.m [3] [0];
C.m [2] [1] = A.m [2] [0] * B.m [0] [1] + A.m [2] [1] * B.m [1] [1] + A.m [2] [2] * B.m [2] [1] + A.m [2] [3] * B.m [3] [1];
C.m [2] [2] = A.m [2] [0] * B.m [0] [2] + A.m [2] [1] * B.m [1] [2] + A.m [2] [2] * B.m [2] [2] + A.m [2] [3] * B.m [3] [2];
C.m [2] [3] = A.m [2] [0] * B.m [0] [3] + A.m [2] [1] * B.m [1] [3] + A.m [2] [2] * B.m [2] [3] + A.m [2] [3] * B.m [3] [3];
C.m [3] [0] = A.m [3] [0] * B.m [0] [0] + A.m [3] [1] * B.m [1] [0] + A.m [3] [2] * B.m [2] [0] + A.m [3] [3] * B.m [3] [0];
C.m [3] [1] = A.m [3] [0] * B.m [0] [1] + A.m [3] [1] * B.m [1] [1] + A.m [3] [2] * B.m [2] [1] + A.m [3] [3] * B.m [3] [1];
C.m [3] [2] = A.m [3] [0] * B.m [0] [2] + A.m [3] [1] * B.m [1] [2] + A.m [3] [2] * B.m [2] [2] + A.m [3] [3] * B.m [3] [2];
C.m [3] [3] = A.m [3] [0] * B.m [0] [3] + A.m [3] [1] * B.m [1] [3] + A.m [3] [2] * B.m [2] [3] + A.m [3] [3] * B.m [3] [3];
}
//--- HSL→ARGB при S=1, L=0.5 (чистый спектр)
uint C_TestStand3D::_HSL (int h_deg)
{
double hue = h_deg / 360.0;
uint ch [3];
double off [3];
off [0] = 1.0 / 3.0;
off [1] = 0.0;
off [2] = -1.0 / 3.0;
for (int i = 0; i < 3; i++)
{
double vh = hue + off [i];
if (vh < 0.0) vh += 1.0;
if (vh > 1.0) vh -= 1.0;
double v;
if (6.0 * vh < 1.0) v = 6.0 * vh;
else
if (2.0 * vh < 1.0) v = 1.0;
else
if (3.0 * vh < 2.0) v = (2.0 / 3.0 - vh) * 6.0;
else v = 0.0;
ch [i] = (uint)(v * 255.0);
}
return 0xFF000000 | (ch [0] << 16) | (ch [1] << 8) | ch [2];
}
//+------------------------------------------------------------------+