![\"\"](\"https:\/\/123210.net\/wordpress\/tatenaoshi\/wp-content\/uploads\/sites\/3\/2026\/07\/\u5b9f\u884c\u72b6\u6cc110.jpg\")
<\/a><\/p>\n<\/p>\n
\n
C++Builder2007\u30bd\u30fc\u30b9\u30b3\u30fc\u30c9<\/p>\n
\/\/*** Unit1.h ***<\/p>\n
\/\/—————————————————————————<\/p>\n
#ifndef Unit1H
\n#define Unit1H
\n\/\/—————————————————————————
\n#include <Classes.hpp>
\n#include <Controls.hpp>
\n#include <StdCtrls.hpp>
\n#include <Forms.hpp>
\n#include “Chart.hpp”
\n#include “TeEngine.hpp”
\n#include “TeeProcs.hpp”
\n#include <ExtCtrls.hpp><\/p>\n
#include <Series.hpp>
\n#include <ComCtrls.hpp> \/\/ \u2190 \u3053\u308c\u3092\u8ffd\u52a0\uff08\u6700\u91cd\u8981\uff09
\n\/\/—————————————————————————<\/p>\n
\/\/————————————————
\n\/\/ \u5c65\u6b74\u7528\u69cb\u9020\u4f53
\n\/\/————————————————
\nstruct OctoState {
\nint step;
\nint depth;
\ndouble c[8];
\ndouble norm2;
\ndouble N;
\nint lastDeriv;
\n};<\/p>\n
\/\/————————————————
\n\/\/ \u516b\u5143\u6570
\n\/\/————————————————
\nstruct Octonion {
\ndouble e[8];<\/p>\n
Octonion(double a0=0.0, double a1=0.0, double a2=0.0, double a3=0.0,
\ndouble a4=0.0, double a5=0.0, double a6=0.0, double a7=0.0);<\/p>\n
Octonion operator+(const Octonion& o) const;
\nOctonion operator-(const Octonion& o) const; \/\/ \u6e1b\u7b97
\nOctonion operator-() const; \/\/ \u5358\u9805\u30de\u30a4\u30ca\u30b9\uff08\u5ff5\u306e\u305f\u3081\uff09
\nOctonion operator*(const Octonion& o) const; \/\/ \u2190 \u5909\u66f4\uff1a\u5b8c\u5168\u4e57\u7b97
\nOctonion operator*(double scalar) const; \/\/ \u30b9\u30ab\u30e9\u30fc\u4e57\u7b97\u306f\u6b8b\u3059
\ndouble Norm() const;<\/p>\n
Octonion Conjugate() const; \/\/ \u6a19\u6e96\u7684\u306a\u5171\u5f79
\nOctonion PhaseConjugate() const; \/\/ \u4f4d\u76f8\u5171\u5f79\u7528\uff08\u5fc5\u8981\u306b\u5fdc\u3058\u3066\u62e1\u5f35\u53ef\u80fd\uff09
\nOctonion ApplyDerivation(const Octonion& a, const Octonion& b) const; \/\/ D_a,b (this)
\n};<\/p>\n
\/\/————————————————
\n\/\/ \u518d\u5e30\u516b\u5143\u6570
\n\/\/————————————————
\nclass RecursiveOcto {
\npublic:
\nint depth;
\nOctonion value;
\nRecursiveOcto* child; \/\/ \u901a\u5e38\u306e\u5b50
\nRecursiveOcto* childConj; \/\/ \u5171\u5f79\u306e\u5b50\uff08\u8ffd\u52a0\uff09
\nOctonion interference; \/\/ \u5e72\u6e09\u7d50\u679c\u3092\u4fdd\u6301\u3059\u308b\u30e1\u30f3\u30d0\u30fc\u3092\u8ffd\u52a0<\/p>\n
\/\/ \u2605 \u8ffd\u52a0
\ndouble N; \/\/ \u6570\u6f14\u7b97\u5b50\uff08\u52b1\u8d77\u5ea6\uff09
\ndouble Q; \/\/ \u96fb\u8377\u7684\u306a\u91cf = N \/ 3<\/p>\n
RecursiveOcto(int d, Octonion init);
\n~RecursiveOcto();
\nvoid Generate(int maxDepth);
\nvoid Dump(TMemo* memo);<\/p>\n
void ComputeNumberOperator(); \/\/ \u2605 \u8ffd\u52a0<\/p>\n
\/\/ \u89b3\u5bdf\uff08Observe\uff09\u6a5f\u80fd\u3068\u306e\u9023\u643a
\nint preferredDir; \/\/ \u89b3\u5bdf\u3067\u56fa\u5b9a\u3055\u308c\u305f\u65b9\u5411\uff08-1 = \u672a\u89b3\u5bdf\uff09
\nvoid Observe(int dir); \/\/ dir = 0?6\uff08e1?e7\uff09
\n};<\/p>\n
\/\/————————————————
\n\/\/ \u30d5\u30a9\u30fc\u30e0
\n\/\/————————————————
\nclass TForm1 : public TForm
\n{
\n__published: \/\/ IDE-managed Components
\nTButton *ButtonExit;
\nTMemo *Memo1;
\nTButton *Button2;
\nTChart *ChartTime;
\nTChart *ChartNorm;
\nTChart *Chart3D;
\nTPaintBox *PaintBox3D;
\nTTrackBar *TrackBarZoom;
\nTRadioGroup *RadioGroup1;
\nTEdit *Edit1;
\nTUpDown *UpDown1;
\nvoid __fastcall ButtonExitClick(TObject *Sender);
\nvoid __fastcall Button2Click(TObject *Sender);
\nvoid __fastcall FormCreate(TObject *Sender);
\nvoid __fastcall PaintBox3DPaint(TObject *Sender);
\nvoid __fastcall PaintBox3DMouseDown(TObject *Sender, TMouseButton Button,
\nTShiftState Shift, int X, int Y);
\nvoid __fastcall PaintBox3DMouseMove(TObject *Sender, TShiftState Shift, int X,
\nint Y);
\nvoid __fastcall PaintBox3DMouseUp(TObject *Sender, TMouseButton Button,
\nTShiftState Shift, int X, int Y);
\nvoid __fastcall TrackBarZoomChange(TObject *Sender);
\nvoid __fastcall RadioGroup1Click(TObject *Sender);
\nprivate: \/\/ User declarations
\nTList* History;<\/p>\n
void ClearHistory();
\nvoid CollectHistory(RecursiveOcto* node, int& counter);
\nvoid __fastcall InitTimeCharts();
\nvoid Init3DChart();
\nvoid __fastcall UpdateTimeCharts();
\nvoid __fastcall Update3DChart();<\/p>\n
\/\/OpenGL\u98a8
\ndouble rotX, rotY; \/\/ \u56de\u8ee2\u89d2\u5ea6\uff08\u30e9\u30b8\u30a2\u30f3\uff09
\ndouble zoom; \/\/ \u30ba\u30fc\u30e0\uff081.0\u304c\u57fa\u6e96\uff09
\ndouble panX, panY; \/\/ \u753b\u9762\u79fb\u52d5 int oldX, oldY; \/\/ \u30de\u30a6\u30b9\u524d\u56de\u4f4d\u7f6e
\nint oldX, oldY;
\nbool dragging;
\nbool panning; \/\/ \u30d1\u30f3\u4e2d\u304b\u3069\u3046\u304b
\nint projMode; \/\/ 0\uff5e4 \u306a\u3069\u306e\u5c04\u5f71\u756a\u53f7<\/p>\n
void Project3D(double x, double y, double z, int& sx, int& sy);
\nvoid DrawPointCloud();<\/p>\n
public: \/\/ User declarations
\n__fastcall TForm1(TComponent* Owner);
\n};
\n\/\/—————————————————————————
\nextern PACKAGE TForm1 *Form1;
\n\/\/—————————————————————————
\n#endif<\/p>\n
<\/p>\n
\n
\/\/*** Unit1.cpp ***<\/p>\n
\/\/—————————————————————————<\/p>\n
#include <vcl.h>
\n#pragma hdrstop<\/p>\n
#include <math.h>
\n#include “Unit1.h”
\n\/\/—————————————————————————
\n#pragma package(smart_init)
\n#pragma link “Chart”
\n#pragma link “TeEngine”
\n#pragma link “TeeProcs”
\n#pragma resource “*.dfm”
\nTForm1 *Form1;<\/p>\n
\/\/\u30ea\u30f3\u30ab\uff5c\u51fa\u529b\u30aa\u30d7\u30b7\u30e7\u30f3\uff5c\u6700\u5927\u30b9\u30bf\u30c3\u30af\u30b5\u30a4\u30ba\u30000x64000000\uff08\u7d041.56GB\uff09\u5927\u304d\u3059\u304e\u3067\u3042\u3063\u305f\u306e\u3067\u30010x10000000\u306b\u8a2d\u5b9a
\n\/\/—————————————————————————
\n__fastcall TForm1::TForm1(TComponent* Owner)
\n: TForm(Owner)
\n{
\n}
\n\/\/—————————————————————————
\nvoid __fastcall TForm1::ButtonExitClick(TObject *Sender)
\n{
\nClose();
\n}
\n\/\/—————————————————————————
\n\/\/—————————————————————————<\/p>\n
Octonion::Octonion(double a0, double a1, double a2, double a3,
\ndouble a4, double a5, double a6, double a7) {
\ne[0]=a0; e[1]=a1; e[2]=a2; e[3]=a3;
\ne[4]=a4; e[5]=a5; e[6]=a6; e[7]=a7;
\n}<\/p>\n
Octonion Octonion::Conjugate() const {
\nreturn Octonion(
\ne[0],
\n-e[1], -e[2], -e[3], -e[4], -e[5], -e[6], -e[7]
\n);
\n}<\/p>\n
Octonion Octonion::PhaseConjugate() const {
\nreturn Conjugate();
\n}<\/p>\n
Octonion Octonion::operator+(const Octonion& o) const {
\nOctonion res;
\nfor(int i=0; i<8; i++) res.e[i] = e[i] + o.e[i];
\nreturn res;
\n}<\/p>\n
\/\/ \u4e8c\u9805\u6e1b\u7b97
\nOctonion Octonion::operator-(const Octonion& o) const {
\nOctonion res;
\nfor(int i = 0; i < 8; i++) {
\nres.e[i] = e[i] – o.e[i];
\n}
\nreturn res;
\n}<\/p>\n
\/\/ \u5358\u9805\u30de\u30a4\u30ca\u30b9
\nOctonion Octonion::operator-() const {
\nOctonion res;
\nfor(int i = 0; i < 8; i++) {
\nres.e[i] = -e[i];
\n}
\nreturn res;
\n}<\/p>\n
\/\/ \u30b9\u30ab\u30e9\u30fc\u4e57\u7b97\u306f\u6b8b\u3059
\nOctonion Octonion::operator*(double scalar) const {
\nOctonion res;
\nfor(int i=0; i<8; i++) res.e[i] = e[i] * scalar;
\nreturn res;
\n}<\/p>\n
\/\/ \u5b8c\u5168\u306a\u516b\u5143\u6570\u4e57\u7b97
\nOctonion Octonion::operator*(const Octonion& o) const {
\nOctonion res;<\/p>\n
\/\/ e0\uff08\u5b9f\u90e8\uff09
\nres.e[0] = e[0]*o.e[0] – e[1]*o.e[1] – e[2]*o.e[2] – e[3]*o.e[3]
\n– e[4]*o.e[4] – e[5]*o.e[5] – e[6]*o.e[6] – e[7]*o.e[7];<\/p>\n
\/\/ e1
\nres.e[1] = e[0]*o.e[1] + e[1]*o.e[0] + e[2]*o.e[3] – e[3]*o.e[2]
\n+ e[4]*o.e[5] – e[5]*o.e[4] – e[6]*o.e[7] + e[7]*o.e[6];<\/p>\n
\/\/ e2
\nres.e[2] = e[0]*o.e[2] – e[1]*o.e[3] + e[2]*o.e[0] + e[3]*o.e[1]
\n+ e[4]*o.e[6] + e[5]*o.e[7] – e[6]*o.e[4] – e[7]*o.e[5];<\/p>\n
\/\/ e3
\nres.e[3] = e[0]*o.e[3] + e[1]*o.e[2] – e[2]*o.e[1] + e[3]*o.e[0]
\n+ e[4]*o.e[7] – e[5]*o.e[6] + e[6]*o.e[5] – e[7]*o.e[4];<\/p>\n
\/\/ e4
\nres.e[4] = e[0]*o.e[4] – e[1]*o.e[5] – e[2]*o.e[6] – e[3]*o.e[7]
\n+ e[4]*o.e[0] + e[5]*o.e[1] + e[6]*o.e[2] + e[7]*o.e[3];<\/p>\n
\/\/ e5
\nres.e[5] = e[0]*o.e[5] + e[1]*o.e[4] – e[2]*o.e[7] + e[3]*o.e[6]
\n– e[4]*o.e[1] + e[5]*o.e[0] – e[6]*o.e[3] + e[7]*o.e[2];<\/p>\n
\/\/ e6
\nres.e[6] = e[0]*o.e[6] + e[1]*o.e[7] + e[2]*o.e[4] – e[3]*o.e[5]
\n– e[4]*o.e[2] + e[5]*o.e[3] + e[6]*o.e[0] – e[7]*o.e[1];<\/p>\n
\/\/ e7
\nres.e[7] = e[0]*o.e[7] – e[1]*o.e[6] + e[2]*o.e[5] + e[3]*o.e[4]
\n– e[4]*o.e[3] – e[5]*o.e[2] + e[6]*o.e[1] + e[7]*o.e[0];<\/p>\n
return res;
\n}<\/p>\n
double Octonion::Norm() const {
\ndouble sum = 0.0;
\nfor(int i=0; i<8; i++) sum += e[i]*e[i];
\nreturn sqrt(sum);
\n}<\/p>\n
Octonion Octonion::ApplyDerivation(const Octonion& a, const Octonion& b) const
\n{
\n\/\/ 1. \u4ea4\u63db\u5b50 [a,b]
\nOctonion ab = a * b – b * a; \/\/ [a,b]<\/p>\n
\/\/ 2. \u4ea4\u63db\u5b50\u306e\u4f5c\u7528 [[a,b], this]
\nOctonion commutator_term = ab * (*this) – (*this) * ab;<\/p>\n
\/\/ 3. \u7d50\u5408\u5b50 [a,b,this] = (a*b)*this – a*(b*this)
\nOctonion assoc = (a * b) * (*this) – a * (b * (*this));<\/p>\n
\/\/ 4. \u6a19\u6e96\u7684\u306a\u5c0e\u5206\u516c\u5f0f\uff08\u4fc2\u6570\u306f\u8abf\u6574\u53ef\u80fd\uff09
\n\/\/ D = [[a,b],x] – 3[a,b,x]
\nOctonion result = commutator_term – assoc * 3.0;<\/p>\n
\/\/ \u30b9\u30b1\u30fc\u30eb\u8abf\u6574\uff08\u6442\u52d5\u304c\u5927\u304d\u3059\u304e\u306a\u3044\u3088\u3046\u306b\uff09
\n\/\/return result * 0.5;
\n\/\/return result * 0.20; \/\/ \u6442\u52d5\u304c\u5c11\u3057\u5f37\u3081\u306b\u51fa\u3066\u3044\u308b\u3088\u3046\u306a\u306e\u3067\u3001\u30b9\u30b1\u30fc\u30eb\u3092\u5c11\u3057\u5f31\u3081\u30660.5\u21920.20\uff0860%\u6e1b\uff09\u5b89\u5b9a\u6027\u3092\u78ba\u8a8d\u3059\u308b
\n\/\/return result * 0.25; \/\/ \u6442\u52d5\u304c\u5c11\u3057\u5f37\u3081\u306b\u51fa\u3066\u3044\u308b\u3088\u3046\u306a\u306e\u3067\u3001\u30b9\u30b1\u30fc\u30eb\u3092\u5c11\u3057\u5f31\u3081\u30660.5\u21920.25\uff0850%\u6e1b\uff09\u5b89\u5b9a\u6027\u3092\u78ba\u8a8d\u3059\u308b
\n\/\/return result * 0.30; \/\/ \u6442\u52d5\u304c\u5c11\u3057\u5f37\u3081\u306b\u51fa\u3066\u3044\u308b\u3088\u3046\u306a\u306e\u3067\u3001\u30b9\u30b1\u30fc\u30eb\u3092\u5c11\u3057\u5f31\u3081\u30660.5\u21920.30\uff0840%\u6e1b\uff09\u5b89\u5b9a\u6027\u3092\u78ba\u8a8d\u3059\u308b
\nreturn result * 0.35; \/\/ \u6442\u52d5\u304c\u5c11\u3057\u5f37\u3081\u306b\u51fa\u3066\u3044\u308b\u3088\u3046\u306a\u306e\u3067\u3001\u30b9\u30b1\u30fc\u30eb\u3092\u5c11\u3057\u5f31\u3081\u30660.5\u21920.35\uff0830%\u6e1b\uff09\u5b89\u5b9a\u6027\u3092\u78ba\u8a8d\u3059\u308b
\n\/\/return result * 0.40; \/\/ \u6442\u52d5\u304c\u5c11\u3057\u5f37\u3081\u306b\u51fa\u3066\u3044\u308b\u3088\u3046\u306a\u306e\u3067\u3001\u30b9\u30b1\u30fc\u30eb\u3092\u5c11\u3057\u5f31\u3081\u30660.5\u21920.40\uff0820%\u6e1b\uff09\u5b89\u5b9a\u6027\u3092\u78ba\u8a8d\u3059\u308b
\n\/\/return result * 0.45; \/\/ \u6442\u52d5\u304c\u5c11\u3057\u5f37\u3081\u306b\u51fa\u3066\u3044\u308b\u3088\u3046\u306a\u306e\u3067\u3001\u30b9\u30b1\u30fc\u30eb\u3092\u5c11\u3057\u5f31\u3081\u30660.5\u21920.45\uff0810%\u6e1b\uff09\u5b89\u5b9a\u6027\u3092\u78ba\u8a8d\u3059\u308b
\n\/\/*** \u8a55\u4fa1\u7d50\u679c\uff1a0.30\uff1a\u63a7\u3048\u3081\u3060\u304c\u52b9\u679c\u306f\u5341\u5206\u51fa\u308b ***
\n\/\/*** \u8a55\u4fa1\u7d50\u679c\uff1a0.35\uff1a\u5c0e\u5206\u306e\u5f71\u97ff\u304c\u3088\u308a\u660e\u78ba\u306b\u73fe\u308c\u3064\u3064\u3001\u5b89\u5b9a\u6027\u3082\u9ad8\u3044 ***<\/p>\n
}<\/p>\n
\/\/—————————————————————————
\nRecursiveOcto::RecursiveOcto(int d, Octonion v)
\n: depth(d), value(v), child(NULL), childConj(NULL),
\nN(0.0), Q(0.0), preferredDir(-1)
\n{
\n}<\/p>\n
RecursiveOcto::~RecursiveOcto() {
\ndelete child;
\ndelete childConj;
\n}<\/p>\n
void RecursiveOcto::Generate(int maxDepth) {
\nif (depth >= maxDepth) return;
\n\/\/if (depth >= maxDepth || depth > 20) return; \/\/ \u5b89\u5168\u4e0a\u9650\u3092\u8ffd\u52a0<\/p>\n
\/\/ === \u65e2\u5b58\u306e\u56de\u8ee2\u5b50 ===
\nOctonion rotator(cos(depth*0.15), sin(depth*0.15), 0.05, 0,0,0,0,0);<\/p>\n
\/\/ === Fano\u5e73\u9762\u306e\u65b9\u5411\u5b9a\u7fa9 ===
\nOctonion dirs[7][2] = {
\n{ Octonion(0,1,0,0,0,0,0,0), Octonion(0,0,1,0,0,0,0,0) }, \/\/ 0
\n{ Octonion(0,1,0,0,0,0,0,0), Octonion(0,0,0,0,1,0,0,0) }, \/\/ 1
\n{ Octonion(0,1,0,0,0,0,0,0), Octonion(0,0,0,0,0,0,1,0) }, \/\/ 2
\n{ Octonion(0,0,1,0,0,0,0,0), Octonion(0,0,0,0,1,0,0,0) }, \/\/ 3
\n{ Octonion(0,0,1,0,0,0,0,0), Octonion(0,0,0,0,0,1,0,0) }, \/\/ 4
\n{ Octonion(0,0,0,1,0,0,0,0), Octonion(0,0,0,0,1,0,0,0) }, \/\/ 5
\n{ Octonion(0,0,0,1,0,0,0,0), Octonion(0,0,0,0,0,1,0,0) } \/\/ 6
\n};<\/p>\n
int dir = depth % 7;<\/p>\n
\/\/ === \u7c21\u6613\u7684\u306a\u5c40\u6240N\u63a8\u5b9a\uff08Generate\u6642\u70b9\u3067\u4f7f\u3048\u308b\u3088\u3046\u306b\uff09 ===
\ndouble imagStrength = 0.0;
\nfor(int i=1; i<=7; i++) imagStrength += fabs(value.e[i]);<\/p>\n
\/\/ \u7c21\u6613N\uff08\u91cf\u5b50\u5316\u524d\uff09
\ndouble localN = 1.0 + imagStrength * 0.8; \/\/ \u30d9\u30fc\u30b9 + \u865a\u90e8\u306e\u5f37\u3055
\nif (localN > 3.0) localN = 3.0;<\/p>\n
\/\/ N\u306b\u5fdc\u3058\u305f\u30b9\u30b1\u30fc\u30eb\u4fc2\u6570\uff08N\u304c\u5927\u304d\u3044\u307b\u3069\u5f37\u304f\u3059\u308b\uff09
\n\/\/ N=0 \u2192 \u5f31\u3044, N=3 \u2192 \u5f37\u3044
\ndouble nFactor = 0.4 + (localN \/ 3.0) * 0.9; \/\/ 0.4 ? 1.3 \u7a0b\u5ea6\u306e\u7bc4\u56f2<\/p>\n
\/\/*******************************************************************
\n\/\/ Generate \u5185\u306e\u5c0e\u5206\u90e8\u5206\u3067\u3001\u89b3\u5bdf\u6e08\u307f\u306e\u5834\u5408\u306f\u56fa\u5b9a\u65b9\u5411\u3092\u907f\u3051\u308b\uff0f\u5f31\u3081\u308b
\nif (preferredDir >= 0) {
\n\/\/ \u89b3\u5bdf\u6e08\u307f\u306a\u3089\u3001\u56fa\u5b9a\u65b9\u5411\u306b\u95a2\u4fc2\u3059\u308b\u5c0e\u5206\u3092\u5f31\u3081\u308b
\n\/\/ \uff08\u7c21\u6613\u7684\u306b nFactor \u3092\u3055\u3089\u306b\u4e0b\u3052\u308b\uff09
\nnFactor *= 0.5;
\n}
\n\/\/*******************************************************************<\/p>\n
\/\/ === \u8907\u6570\u5c0e\u5206\uff08N\u4f9d\u5b58\u30b9\u30b1\u30fc\u30eb\u3092\u639b\u3051\u308b\uff09 ===
\nOctonion g2_perturb = value.ApplyDerivation(dirs[dir][0], dirs[dir][1]) * (0.07 * nFactor);<\/p>\n
int dir2 = (dir + 1) % 7;
\ng2_perturb = g2_perturb + value.ApplyDerivation(dirs[dir2][0], dirs[dir2][1]) * (0.03 * nFactor);<\/p>\n
int dir3 = (dir + 3) % 7;
\ng2_perturb = g2_perturb + value.ApplyDerivation(dirs[dir3][0], dirs[dir3][1]) * (0.015 * nFactor);<\/p>\n
\/\/ \u6b21\u306e\u72b6\u614b
\nOctonion next = (value + g2_perturb) * rotator;
\nnext = next + Octonion(0.03, 0,0,0,0,0,0,0);<\/p>\n
\/\/ \u5171\u5f79\u7248
\nOctonion nextConj = next.PhaseConjugate();<\/p>\n
\/\/ \u5b50\u751f\u6210 (\uff12\u500b\uff1a\u5171\u5f79\u5bfe\u3092\u751f\u6210)
\nchild = new RecursiveOcto(depth+1, next);
\nchild->Generate(maxDepth);<\/p>\n
childConj = new RecursiveOcto(depth+1, nextConj);
\nchildConj->Generate(maxDepth);<\/p>\n
\/\/ \u5e72\u6e09 + \u30d5\u30a3\u30fc\u30c9\u30d0\u30c3\u30af
\ninterference = (child->value + childConj->value) * 0.5;
\ndouble feedbackStrength = 0.3;
\nvalue = value * (1.0 – feedbackStrength) + interference * feedbackStrength;
\n}<\/p>\n
\/\/ \u6570\u6f14\u7b97\u5b50\u3092\u8a08\u7b97\u3059\u308b\u30e1\u30bd\u30c3\u30c9\uff08\u65b0\u898f\u8ffd\u52a0\uff09
\nvoid RecursiveOcto::ComputeNumberOperator()
\n{
\n\/\/if (depth > 12) return; \/\/ \u8ffd\u52a0<\/p>\n
\/\/ \u57fa\u672c\u52b1\u8d77\uff1a\u81ea\u5206\u81ea\u8eab\u306e\u865a\u6570\u6210\u5206\u306e\u300c\u5f37\u3055\u300d
\n\/\/ \uff08e1?e7 \u306e\u7d76\u5bfe\u5024\u306e\u5408\u8a08\u3092\u9069\u5ea6\u306b\u30b9\u30b1\u30fc\u30eb\uff09
\ndouble imagStrength = 0.0;<\/p>\n
\/\/ \u2605 \u3053\u3053\u3092 e[1] \uff5e e[7] \u306b\u4fee\u6b63
\nimagStrength += fabs(value.e[1]);
\nimagStrength += fabs(value.e[2]);
\nimagStrength += fabs(value.e[3]);
\nimagStrength += fabs(value.e[4]);
\nimagStrength += fabs(value.e[5]);
\nimagStrength += fabs(value.e[6]);
\nimagStrength += fabs(value.e[7]);<\/p>\n
\/\/ \u5b50\u4f9b\u306e\u6709\u7121\u306b\u3088\u308b\u52b1\u8d77
\nint childCount = 0;
\nif (child) childCount++;
\nif (childConj) childCount++;<\/p>\n
\/\/ \u518d\u5e30\u7684\u306b\u5b50\u4f9b\u306e N \u3092\u52a0\u7b97\uff08\u968e\u5c64\u5168\u4f53\u306e\u52b1\u8d77\u3092\u84c4\u7a4d\uff09 \/\/ \u5b50\u4f9b\u306e N \u3092\u518d\u5e30\u7684\u306b\u96c6\u3081\u308b
\ndouble childrenN = 0.0;
\nif (child) {
\nchild->ComputeNumberOperator();
\nchildrenN += child->N;
\n}
\nif (childConj) {
\nchildConj->ComputeNumberOperator();
\nchildrenN += childConj->N;
\n}<\/p>\n
\/\/ \u6700\u7d42\u7684\u306a N \u306e\u5b9a\u7fa9\uff08\u8abf\u6574\u3057\u3084\u3059\u3044\u5f62\uff09
\n\/\/ \u5b50\u4f9b\u306e\u6570 \u00d7 1.0 + \u865a\u6570\u5f37\u5ea6 \u00d7 \u4fc2\u6570 + \u5b50\u4f9b\u305f\u3061\u306e N \u306e\u4e00\u90e8
\n\/\/N = childCount * 1.0 + imagStrength * 0.8 + childrenN * 0.3;<\/p>\n
\/\/ \u2605 \u4fc2\u6570\u3092\u8abf\u6574\u3057\u305f\u30d0\u30fc\u30b8\u30e7\u30f3
\n\/\/ \u5b50\u4f9b\u306e\u5bc4\u4e0e\u3092\u5f37\u3081\u3001\u865a\u6570\u5f37\u5ea6\u306e\u5f71\u97ff\u3092\u6291\u3048\u3081\u306b\u3059\u308b
\n\/*N = childCount * 1.5
\n+ imagStrength * 0.4
\n+ childrenN * 0.25;*\/
\n\/\/ \u2605 \u3053\u3053\u3092\u5927\u304d\u304f\u5909\u66f4
\nN = childCount * 1.0 \/\/ \u5b50\u4f9b\u306e\u6570\u306e\u5f71\u97ff\u3092\u5f31\u3081\u308b
\n+ imagStrength * 0.6
\n+ childrenN * 0.05; \/\/ \u5b50\u5b6b\u304b\u3089\u306e\u84c4\u7a4d\u3092\u304b\u306a\u308a\u5f31\u304f\u3059\u308b<\/p>\n
\/\/ \u3055\u3089\u306b\u3001\u3042\u308b\u7a0b\u5ea6\u6574\u6570\u306b\u8fd1\u3065\u3051\u308b\u305f\u3081\u306e\u4e38\u3081\uff08\u4efb\u610f\uff09
\n\/\/ N = floor(N + 0.5); \/\/ \u6574\u6570\u306b\u3057\u305f\u3044\u5834\u5408\u306f\u3053\u306e\u884c\u3092\u6709\u52b9\u306b\u3059\u308b
\n\/\/ \u2605 \u3053\u3053\u304b\u3089\u91cf\u5b50\u5316\uff08\u96e2\u6563\u5316\uff09<\/p>\n
\/\/ \u91cf\u5b50\u5316
\nif (N < 0.7) N = 0.0;
\nelse if (N < 1.5) N = 1.0;
\nelse if (N < 2.3) N = 2.0;
\nelse N = 3.0;<\/p>\n
\/\/ \u96fb\u8377\u7684\u306a\u91cf
\nQ = N \/ 3.0;
\n}<\/p>\n
void RecursiveOcto::Observe(int dir)
\n{
\nif (depth > 12) return; \/\/ \u8ffd\u52a0<\/p>\n
if (dir < 0 || dir > 6) return;<\/p>\n
preferredDir = dir;<\/p>\n
\/\/ \u6307\u5b9a\u3057\u305f\u865a\u5358\u4f4d\u306e\u6210\u5206\u3092\u5f37\u8abf\u3057\u3001\u4ed6\u3092\u76f8\u5bfe\u7684\u306b\u6291\u5236\u3059\u308b\u7c21\u6613\u6295\u5f71
\n\/\/ \uff08\u5b8c\u5168\u306a\u5c04\u5f71\u3067\u306f\u306a\u304f\u3001\u89b3\u5bdf\u306b\u3088\u308b\u300c\u56fa\u5b9a\u300d\u306e\u30a4\u30e1\u30fc\u30b8\uff09
\ndouble strength = value.e[dir + 1]; \/\/ e1?e7<\/p>\n
\/\/ \u6307\u5b9a\u65b9\u5411\u3092\u6b8b\u3057\u3064\u3064\u3001\u5168\u4f53\u3092\u5c11\u3057\u6b63\u898f\u5316
\nfor(int i = 1; i <= 7; i++) {
\nif (i == dir + 1) {
\n\/\/ \u6307\u5b9a\u65b9\u5411\u306f\u305d\u306e\u307e\u307e\uff08\u307e\u305f\u306f\u5c11\u3057\u5f37\u8abf\uff09
\nvalue.e[i] = strength * 1.1;
\n} else {
\n\/\/ \u4ed6\u306e\u65b9\u5411\u3092\u5f31\u3081\u308b
\nvalue.e[i] *= 0.6;
\n}
\n}<\/p>\n
\/\/ \u5b50\u306b\u3082\u540c\u3058\u89b3\u5bdf\u3092\u4f1d\u64ad\uff08\u518d\u5e30\uff09
\nif (child) child->Observe(dir);
\nif (childConj) childConj->Observe(dir);
\n}<\/p>\n
void RecursiveOcto::Dump(TMemo* memo)
\n{<\/p>\n
if (!memo) return;
\nAnsiString s;
\ns.sprintf(“Depth=%d Norm=%.4f e0=%.4f e1=%.4f N=%.1f Q=%.3f”,
\ndepth, value.Norm(), value.e[0], value.e[1], N, Q);
\nmemo->Lines->Add(s);<\/p>\n
\/\/ \u30aa\u30d7\u30b7\u30e7\u30f3\uff1a\u865a\u90e8\u306e\u5927\u304d\u3055\u3092\u5c11\u3057\u8a73\u3057\u304f
\ndouble imag = 0;
\nfor(int i=1;i<8;i++) imag += fabs(value.e[i]);
\ns.sprintf(” imagStrength=%.4f”, imag);
\nmemo->Lines->Add(s);<\/p>\n
if (child) { memo->Lines->Add(” [\u901a\u5e38\u306e\u5b50]”); child->Dump(memo); }
\nif (childConj) { memo->Lines->Add(” [\u5171\u5f79\u306e\u5b50]”); childConj->Dump(memo); }
\n}<\/p>\n
\/\/—————————————————————————
\n\/\/—————————————————————————<\/p>\n
void __fastcall TForm1::UpdateTimeCharts()
\n{<\/p>\n
if (!ChartTime || !ChartNorm) return;
\nif (ChartTime->SeriesCount() < 8) return;
\nif (ChartNorm->SeriesCount() < 2) return;<\/p>\n
\/\/ \u30af\u30ea\u30a2
\nfor (int i = 0; i < ChartTime->SeriesCount(); i++)
\nChartTime->Series[i]->Clear();
\nChartNorm->Series[0]->Clear();
\nChartNorm->Series[1]->Clear();<\/p>\n
for (int i = 0; i < History->Count; i++) {
\nOctoState* s = (OctoState*)History->Items[i];
\nfor (int k = 0; k < 8; k++)
\nChartTime->Series[k]->AddXY(s->step, s->c[k]);
\nChartNorm->Series[0]->AddXY(s->step, s->norm2);
\nChartNorm->Series[1]->AddXY(s->step, s->N);
\n}<\/p>\n
ChartTime->Refresh();
\nChartNorm->Refresh();
\n}<\/p>\n
void TForm1::ClearHistory()
\n{
\nfor (int i = 0; i < History->Count; i++)
\ndelete (OctoState*)History->Items[i];
\nHistory->Clear();
\n}<\/p>\n
void __fastcall TForm1::Update3DChart()
\n{<\/p>\n
if (!Chart3D) return;
\nif (Chart3D->SeriesCount() == 0) return;<\/p>\n
TPointSeries* ser = dynamic_cast<TPointSeries*>(Chart3D->Series[0]);
\nif (!ser) return;<\/p>\n
ser->Clear();<\/p>\n
for (int i = 0; i < History->Count; i++) {
\nOctoState* s = (OctoState*)History->Items[i];
\ndouble x = s->c[1];
\ndouble y = s->c[2];<\/p>\n
TColor col = clBlack;
\nif (s->N >= 2.5) col = clRed;
\nelse if (s->N >= 1.5) col = clLime;
\nelse if (s->N >= 0.5) col = clBlue;<\/p>\n
ser->AddXY(x, y, “”, col);
\n}<\/p>\n
Chart3D->Axes->Bottom->Automatic = true;
\nChart3D->Axes->Left->Automatic = true;
\nChart3D->Refresh();
\n}<\/p>\n
void TForm1::CollectHistory(RecursiveOcto* node, int& counter)
\n{
\nif (!node) return;<\/p>\n
OctoState* s = new OctoState;
\ns->step = counter++;
\ns->depth = node->depth;
\nfor (int i = 0; i < 8; i++) s->c[i] = node->value.e[i];
\ns->norm2 = node->value.Norm() * node->value.Norm();
\ns->N = node->N;
\ns->lastDeriv = node->depth % 7;<\/p>\n
History->Add(s);<\/p>\n
CollectHistory(node->child, counter);
\nCollectHistory(node->childConj, counter);
\n}<\/p>\n
void __fastcall TForm1::Button2Click(TObject *Sender)
\n{
\nMemo1->Hide();<\/p>\n
Memo1->Clear();
\nClearHistory();<\/p>\n
Octonion rootVal(1.0, 0, 0, 0, 0, 0, 0, 0);
\nRecursiveOcto* root = new RecursiveOcto(0, rootVal);
\n\/\/root->Generate(5); \/\/ \u5f8c\u30677?8\u306b\u4e0a\u3052\u3066\u304f\u3060\u3055\u3044
\nroot->Generate(UpDown1->Position); \/\/ \u5f8c\u30677?8\u306b\u4e0a\u3052\u3066\u304f\u3060\u3055\u3044
\nroot->Observe(2);
\nroot->ComputeNumberOperator();<\/p>\n
int counter = 0;
\nCollectHistory(root, counter);<\/p>\n
root->Dump(Memo1);<\/p>\n
UpdateTimeCharts();
\nUpdate3DChart();<\/p>\n
Memo1->Show();
\ndelete root;<\/p>\n
\/\/ \u65e2\u5b58\u306e UpdateTimeCharts(); Update3DChart(); \u306e\u5f8c\u306b
\nPaintBox3D->Invalidate(); \/\/ \u65b0\u3057\u3044\u70b9\u7fa4\u3092\u63cf\u753b
\n}
\n\/\/—————————————————————————<\/p>\n
void __fastcall TForm1::FormCreate(TObject *Sender)
\n{
\nHistory = new TList();<\/p>\n
\/\/ \u30c1\u30e3\u30fc\u30c8\u304c\u3061\u3083\u3093\u3068\u30d5\u30a9\u30fc\u30e0\u306b\u7d50\u3073\u3064\u3044\u3066\u3044\u308b\u304b\u78ba\u8a8d
\nif (!ChartTime || !ChartNorm || !Chart3D) {
\nShowMessage(“\u30c1\u30e3\u30fc\u30c8\u30b3\u30f3\u30dd\u30fc\u30cd\u30f3\u30c8\u304c\u30d5\u30a9\u30fc\u30e0\u306b\u6b63\u3057\u304f\u7d50\u3073\u3064\u3044\u3066\u3044\u307e\u305b\u3093\u3002\\n”
\n“IDE\u3067\u30b3\u30f3\u30dd\u30fc\u30cd\u30f3\u30c8\u540d\u3092\u78ba\u8a8d\u3057\u3066\u304f\u3060\u3055\u3044\u3002”);
\nreturn;
\n}<\/p>\n
InitTimeCharts();
\nInit3DChart();<\/p>\n
\/\/OpenGL\u98a8
\nrotX = 0.4; \/\/ \u521d\u671f\u306e\u50be\u304d
\nrotY = 0.6;
\nzoom = 1.0;
\npanX = 0.0;
\npanY = 0.0;
\ndragging = false;
\npanning = false;<\/p>\n
\/\/ PaintBox \u306e\u30a4\u30d9\u30f3\u30c8\u3092\u7d50\u3073\u3064\u3051\u308b\uff08IDE\u3067\u8a2d\u5b9a\u3057\u3066\u3082\u53ef\uff09
\nPaintBox3D->OnPaint = PaintBox3DPaint;
\nPaintBox3D->OnMouseDown = PaintBox3DMouseDown;
\nPaintBox3D->OnMouseMove = PaintBox3DMouseMove;
\nPaintBox3D->OnMouseUp = PaintBox3DMouseUp;<\/p>\n
\/\/ \u5c04\u5f71\u8ef8\u306e\u5207\u308a\u66ff\u3048\u7528\u30d1\u30e9\u30e1\u30fc\u30bf
\nprojMode = 0; \/\/ \u521d\u671f\u306f (e1, e2, e4)
\n}
\n\/\/—————————————————————————<\/p>\n
void __fastcall TForm1::InitTimeCharts()
\n{
\nChartTime->View3D = false;
\nChartTime->Legend->Visible = true;
\nChartTime->Title->Text->Text = “Octonion Components”;<\/p>\n
\/\/ 8\u672c\u306e\u7dda\u30b7\u30ea\u30fc\u30ba\u3092\u4f5c\u6210
\nfor (int i = 0; i < 8; i++) {
\nTLineSeries* ser = new TLineSeries(ChartTime);
\nChartTime->AddSeries(ser);
\n\/\/ser->Title = (i == 0) ? “Re” : “e” + IntToStr(i);
\nser->Title = (i == 0) ? AnsiString(“Re”) : AnsiString(“e”) + IntToStr(i);
\nser->LinePen->Width = 1;
\n\/\/ \u8272\u306f\u304a\u597d\u307f\u3067\uff08\u4f8b\uff09
\nstatic TColor cols[8] = {clBlack, clRed, clGreen, clBlue,
\nclFuchsia, clOlive, clNavy, clMaroon};
\nser->SeriesColor = cols[i];
\n}<\/p>\n
\/\/ \u30ce\u30eb\u30e0\u7528\u30c1\u30e3\u30fc\u30c8
\nChartNorm->View3D = false;
\nTLineSeries* sNorm = new TLineSeries(ChartNorm);
\nTLineSeries* sN = new TLineSeries(ChartNorm);
\nChartNorm->AddSeries(sNorm);
\nChartNorm->AddSeries(sN);
\nsNorm->Title = “|O|2”;
\nsN->Title = “N”;
\nsNorm->SeriesColor = clBlack;
\nsN->SeriesColor = clRed;
\n}<\/p>\n
void TForm1::Init3DChart()
\n{<\/p>\n
if (!Chart3D) return;<\/p>\n
\/\/ \u6700\u4f4e\u9650\u306e\u8a2d\u5b9a\u3060\u3051\u306b\u3059\u308b
\nChart3D->View3D = true;
\nChart3D->Legend->Visible = false;<\/p>\n
\/\/ View3DOptions \u3078\u306e\u30a2\u30af\u30bb\u30b9\u3092\u4e00\u65e6\u524a\u9664\uff08\u3053\u308c\u304cAV\u306e\u539f\u56e0\uff09
\n\/\/ Chart3D->View3DOptions->Orthogonal = false;
\n\/\/ Chart3D->View3DOptions->Zoom = 90;<\/p>\n
\/\/ \u65e2\u5b58\u30b7\u30ea\u30fc\u30ba\u3092\u30af\u30ea\u30a2
\nwhile (Chart3D->SeriesCount() > 0)
\nChart3D->Series[0]->Free();<\/p>\n
\/\/ \u30dd\u30a4\u30f3\u30c8\u30b7\u30ea\u30fc\u30ba\u3092\u8ffd\u52a0
\nTPointSeries* ser = new TPointSeries(Chart3D);
\nChart3D->AddSeries(ser);
\nser->Pointer->Style = psCircle;
\nser->Pointer->HorizSize = 3;
\nser->Pointer->VertSize = 3;
\nser->Pointer->Visible = true;
\n}<\/p>\n
<\/p>\n
<\/p>\n
\/\/—————————————————————————
\n\/\/—————————————————————————
\nvoid TForm1::Project3D(double x, double y, double z, int& sx, int& sy)
\n{<\/p>\n
\/\/ \u56de\u8ee2
\ndouble cosy = cos(rotY), siny = sin(rotY);
\ndouble cosx = cos(rotX), sinx = sin(rotX);<\/p>\n
double x1 = x * cosy – z * siny;
\ndouble z1 = x * siny + z * cosy;
\ndouble y1 = y;<\/p>\n
double y2 = y1 * cosx – z1 * sinx;
\ndouble z2 = y1 * sinx + z1 * cosx;
\ndouble x2 = x1;<\/p>\n
\/\/ \u30ba\u30fc\u30e0\u4ed8\u304d\u900f\u8996\u6295\u5f71
\ndouble scale = (250.0 * zoom) \/ (4.0 + z2);
\nsx = (int)(PaintBox3D->Width \/ 2 + panX + x2 * scale);
\nsy = (int)(PaintBox3D->Height \/ 2 + panY – y2 * scale);
\n}<\/p>\n
void TForm1::DrawPointCloud()
\n{
\nTCanvas* c = PaintBox3D->Canvas;<\/p>\n
\/\/ \u80cc\u666f
\nc->Brush->Color = (TColor)0x1A1A2E;
\nc->FillRect(PaintBox3D->ClientRect);<\/p>\n
\/\/ ========== 3D\u30b0\u30ea\u30c3\u30c9\u30fb\u5ea7\u6a19\u8ef8\u3092\u63cf\u753b ==========
\nc->Pen->Width = 1;<\/p>\n
\/\/ \u8ef8\u306e\u9577\u3055\uff08\u30c7\u30fc\u30bf\u306b\u5408\u308f\u305b\u3066\u8abf\u6574\u3057\u3066\u304f\u3060\u3055\u3044\uff09
\ndouble axisLen = 0.8;<\/p>\n
\/\/ \u8272\u3092\u5909\u3048\u30663\u8ef8\u3092\u63cf\u304f
\n\/\/ X\u8ef8\uff08\u8d64\u3063\u307d\u3044\uff09
\nc->Pen->Color = (TColor)0x5555AA;
\nint x1, y1, x2, y2;
\nProject3D(-axisLen, 0, 0, x1, y1);
\nProject3D( axisLen, 0, 0, x2, y2);
\nc->MoveTo(x1, y1); c->LineTo(x2, y2);<\/p>\n
\/\/ Y\u8ef8\uff08\u7dd1\u3063\u307d\u3044\uff09
\nc->Pen->Color = (TColor)0x55AA55;
\nProject3D(0, -axisLen, 0, x1, y1);
\nProject3D(0, axisLen, 0, x2, y2);
\nc->MoveTo(x1, y1); c->LineTo(x2, y2);<\/p>\n
\/\/ Z\u8ef8\uff08\u9752\u3063\u307d\u3044\uff09
\nc->Pen->Color = (TColor)0xAA5555;
\nProject3D(0, 0, -axisLen, x1, y1);
\nProject3D(0, 0, axisLen, x2, y2);
\nc->MoveTo(x1, y1); c->LineTo(x2, y2);<\/p>\n
\/\/ \u7c21\u6613\u30b0\u30ea\u30c3\u30c9\uff08XY\u5e73\u9762\u30fbXZ\u5e73\u9762\u30fbYZ\u5e73\u9762\u306b\u6570\u672c\uff09
\nc->Pen->Color = (TColor)0x333355;
\ndouble g = 0.4; \/\/ \u30b0\u30ea\u30c3\u30c9\u9593\u9694<\/p>\n
\/\/ XY\u5e73\u9762\u306e\u30b0\u30ea\u30c3\u30c9
\nfor (double v = -axisLen; v <= axisLen; v += g) {
\nProject3D(v, -axisLen, 0, x1, y1);
\nProject3D(v, axisLen, 0, x2, y2);
\nc->MoveTo(x1, y1); c->LineTo(x2, y2);<\/p>\n
Project3D(-axisLen, v, 0, x1, y1);
\nProject3D( axisLen, v, 0, x2, y2);
\nc->MoveTo(x1, y1); c->LineTo(x2, y2);
\n}<\/p>\n
\/\/ ========== \u3053\u3053\u304b\u3089\u70b9\u7fa4\u63cf\u753b\uff08\u65e2\u5b58\u306e\u307e\u307e\uff09 ==========
\nif (!History || History->Count == 0) return;<\/p>\n
for (int i = 0; i < History->Count; i++) {
\nOctoState* s = (OctoState*)History->Items[i];<\/p>\n
\/\/double x = s->c[1];
\n\/\/double y = s->c[2];
\n\/\/double z = s->c[4];<\/p>\n
\/\/ ===== \u5c04\u5f71\u8ef8\u306e\u9078\u629e =====
\ndouble x, y, z;<\/p>\n
switch (projMode) {
\ncase 0: \/\/ \u57fa\u672c (e1, e2, e4)
\nx = s->c[1];
\ny = s->c[2];
\nz = s->c[4];
\nbreak;<\/p>\n
case 1: \/\/ Observe\u65b9\u5411\u3092\u542b\u3080 (e1, e3, e7)
\nx = s->c[1];
\ny = s->c[3]; \/\/ e3
\nz = s->c[7];
\nbreak;<\/p>\n
case 2: \/\/ \u5b9f\u90e8\u3068\u4e3b\u306a\u865a\u90e8 (Re, e1, e2)
\nx = s->c[0];
\ny = s->c[1];
\nz = s->c[2];
\nbreak;<\/p>\n
case 3: \/\/ Fano\u4e09\u89d2\u5f62\u98a8
\nx = s->c[1] + s->c[2] + s->c[4];
\ny = s->c[3] + s->c[5] + s->c[6];
\nz = s->c[7];
\nbreak;<\/p>\n
case 4: \/\/ \u52b1\u8d77\u5ea6\u3092Z\u8ef8\u306b (e1, e2, N)
\nx = s->c[1];
\ny = s->c[2];
\nz = s->N * 0.3; \/\/ \u30b9\u30b1\u30fc\u30eb\u8abf\u6574
\nbreak;<\/p>\n
default:
\nx = s->c[1];
\ny = s->c[2];
\nz = s->c[4];
\nbreak;
\n}<\/p>\n
int sx, sy;
\nProject3D(x, y, z, sx, sy);<\/p>\n
if (sx < -20 || sx > PaintBox3D->Width+20 ||
\nsy < -20 || sy > PaintBox3D->Height+20) continue;<\/p>\n
TColor col;
\nif (s->N >= 2.5) col = clRed;
\nelse if (s->N >= 1.5) col = clLime;
\nelse if (s->N >= 0.5) col = clAqua;
\nelse col = clSilver;<\/p>\n
c->Brush->Color = col;
\nc->Pen->Color = col;
\nc->Ellipse(sx-1, sy-1, sx+2, sy+2);
\n}
\n}<\/p>\n
void __fastcall TForm1::PaintBox3DPaint(TObject *Sender)
\n{
\nDrawPointCloud();
\n}
\n\/\/—————————————————————————<\/p>\n
void __fastcall TForm1::PaintBox3DMouseDown(TObject *Sender,
\nTMouseButton Button, TShiftState Shift, int X, int Y)
\n{
\noldX = X;
\noldY = Y;<\/p>\n
if (Button == mbLeft) {
\ndragging = true;
\npanning = false;
\n}
\nelse if (Button == mbRight) {
\npanning = true;
\ndragging = false;
\n}
\n}
\n\/\/—————————————————————————<\/p>\n
void __fastcall TForm1::PaintBox3DMouseMove(TObject *Sender, TShiftState Shift,
\nint X, int Y)
\n{
\nif (dragging) {
\n\/\/ \u56de\u8ee2
\nrotY += (X – oldX) * 0.01;
\nrotX += (Y – oldY) * 0.01;
\noldX = X;
\noldY = Y;
\nPaintBox3D->Invalidate();
\n}
\nelse if (panning) {
\n\/\/ \u79fb\u52d5
\npanX += (X – oldX);
\npanY += (Y – oldY);
\noldX = X;
\noldY = Y;
\nPaintBox3D->Invalidate();
\n}
\n}
\n\/\/—————————————————————————<\/p>\n
void __fastcall TForm1::PaintBox3DMouseUp(TObject *Sender, TMouseButton Button,
\nTShiftState Shift, int X, int Y)
\n{
\n\/\/ dragging = false;<\/p>\n
dragging = false;
\npanning = false;<\/p>\n
}
\n\/\/—————————————————————————<\/p>\n
void __fastcall TForm1::TrackBarZoomChange(TObject *Sender)
\n{
\nint WheelDelta = TrackBarZoom->Position;<\/p>\n
if (WheelDelta > 0)
\nzoom *= 1.15; \/\/ \u62e1\u5927
\nelse
\nzoom \/= 1.15; \/\/ \u7e2e\u5c0f<\/p>\n
\/\/ \u7bc4\u56f2\u5236\u9650
\nif (zoom < 0.1) zoom = 0.1;
\nif (zoom > 30.0) zoom = 30.0;<\/p>\n
PaintBox3D->Invalidate();
\n\/\/Handled = true;
\n}
\n\/\/—————————————————————————<\/p>\n
<\/p>\n
void __fastcall TForm1::RadioGroup1Click(TObject *Sender)
\n{
\nprojMode = RadioGroup1->ItemIndex;
\nPaintBox3D->Invalidate();
\n}
\n\/\/—————————————————————————<\/p>\n
<\/p>\n","protected":false},"excerpt":{"rendered":"
root->Generate(12);\u3068\u3057\u3066 maxDepth\u309212\u306b\u3042\u3052 … \u7d9a\u304d\u3092\u8aad\u3080