立直し-地上天国

時間と空間の原理 新八元数と相対性理論 改訂第２版 後藤 博茂さん
https://gotoclinic.org/timespace2.pdf

実行結果
Depth=0 Norm=1.0093 e0=1.0093 e1=0.0000
→ 干渉 Norm=1.0309 e0=1.0309 e1=0.0000
[通常の子]
Depth=1 Norm=1.0308 e0=1.0302 e1=0.0000
→ 干渉 Norm=1.0307 e0=1.0307 e1=0.0000
[通常の子]
Depth=2 Norm=1.0244 e0=1.0162 e1=0.1077
→ 干渉 Norm=0.9468 e0=0.9468 e1=0.0000
[通常の子]
Depth=3 Norm=0.9498 e0=0.8883 e1=0.3196
→ 干渉 Norm=0.6774 e0=0.6774 e1=0.0000
[通常の子]
Depth=4 Norm=0.8027 e0=0.5329 e1=0.5868
→ 干渉 Norm=0.1305 e0=0.1305 e1=0.0000
[通常の子]
Depth=5 Norm=1.1149 e0=0.1305 e1=1.0935
[共役の子]
Depth=5 Norm=1.1149 e0=0.1305 e1=-1.0935
[共役の子]
Depth=4 Norm=1.0178 e0=0.8219 e1=-0.5868
→ 干渉 Norm=1.0939 e0=1.0939 e1=0.0000
[通常の子]
Depth=5 Norm=1.1405 e0=1.0939 e1=-0.2969
[共役の子]
Depth=5 Norm=1.1405 e0=1.0939 e1=0.2969
[共役の子]
Depth=3 Norm=1.0601 e0=1.0054 e1=-0.3196
→ 干渉 Norm=1.0676 e0=1.0676 e1=0.0000
[通常の子]
Depth=4 Norm=1.0686 e0=1.0664 e1=0.0092
→ 干渉 Norm=0.9481 e0=0.9481 e1=0.0000
[通常の子]
Depth=5 Norm=1.1475 e0=0.9481 e1=0.6383
[共役の子]
Depth=5 Norm=1.1475 e0=0.9481 e1=-0.6383
[共役の子]
Depth=4 Norm=1.0710 e0=1.0688 e1=-0.0092
→ 干渉 Norm=0.9559 e0=0.9559 e1=0.0000
[通常の子]
Depth=5 Norm=1.1477 e0=0.9559 e1=0.6233
[共役の子]
Depth=5 Norm=1.1477 e0=0.9559 e1=-0.6233
[共役の子]
Depth=2 Norm=1.0531 e0=1.0452 e1=-0.1077
→ 干渉 Norm=1.0434 e0=1.0434 e1=0.0000
[通常の子]
Depth=3 Norm=1.0309 e0=1.0241 e1=0.1132
→ 干渉 Norm=0.8941 e0=0.8941 e1=0.0000
[通常の子]
Depth=4 Norm=0.9000 e0=0.7897 e1=0.4296
→ 干渉 Norm=0.4529 e0=0.4529 e1=0.0000
[通常の子]
Depth=5 Norm=1.1323 e0=0.4529 e1=1.0313
[共役の子]
Depth=5 Norm=1.1323 e0=0.4529 e1=-1.0313
[共役の子]
Depth=4 Norm=1.0878 e0=0.9985 e1=-0.4296
→ 干渉 Norm=1.1489 e0=1.1489 e1=0.0000
[通常の子]
Depth=5 Norm=1.1506 e0=1.1489 e1=0.0235
[共役の子]
Depth=5 Norm=1.1506 e0=1.1489 e1=-0.0235
[共役の子]
Depth=3 Norm=1.0695 e0=1.0628 e1=-0.1132
→ 干渉 Norm=1.0234 e0=1.0234 e1=0.0000
[通常の子]
Depth=4 Norm=0.9958 e0=0.9671 e1=0.2279
→ 干渉 Norm=0.7257 e0=0.7257 e1=0.0000
[通常の子]
Depth=5 Norm=1.1431 e0=0.7257 e1=0.8758
[共役の子]
Depth=5 Norm=1.1431 e0=0.7257 e1=-0.8758
[共役の子]
Depth=4 Norm=1.1055 e0=1.0797 e1=-0.2279
→ 干渉 Norm=1.1013 e0=1.1013 e1=0.0000
[通常の子]
Depth=5 Norm=1.1529 e0=1.1013 e1=0.3331
[共役の子]
Depth=5 Norm=1.1529 e0=1.1013 e1=-0.3331
[共役の子]
Depth=1 Norm=1.0323 e0=1.0317 e1=0.0000
→ 干渉 Norm=1.0356 e0=1.0356 e1=0.0000
[通常の子]
Depth=2 Norm=1.0282 e0=1.0225 e1=0.1077
→ 干渉 Norm=0.9561 e0=0.9561 e1=0.0000
[通常の子]
Depth=3 Norm=0.9549 e0=0.8988 e1=0.3201
→ 干渉 Norm=0.6896 e0=0.6896 e1=0.0000
[通常の子]
Depth=4 Norm=0.8061 e0=0.5456 e1=0.5887
→ 干渉 Norm=0.1427 e0=0.1427 e1=0.0000
[通常の子]
Depth=5 Norm=1.1160 e0=0.1427 e1=1.0991
[共役の子]
Depth=5 Norm=1.1160 e0=0.1427 e1=-1.0991
[共役の子]
Depth=4 Norm=1.0233 e0=0.8337 e1=-0.5887
→ 干渉 Norm=1.1030 e0=1.1030 e1=0.0000
[通常の子]
Depth=5 Norm=1.1415 e0=1.1030 e1=-0.2879
[共役の子]
Depth=5 Norm=1.1415 e0=1.1030 e1=0.2879
[共役の子]
Depth=3 Norm=1.0636 e0=1.0135 e1=-0.3201
→ 干渉 Norm=1.0720 e0=1.0720 e1=0.0000
[通常の子]
Depth=4 Norm=1.0690 e0=1.0689 e1=0.0132
→ 干渉 Norm=0.9455 e0=0.9455 e1=0.0000
[通常の子]
Depth=5 Norm=1.1480 e0=0.9455 e1=0.6495
[共役の子]
Depth=5 Norm=1.1480 e0=0.9455 e1=-0.6495
[共役の子]
Depth=4 Norm=1.0752 e0=1.0750 e1=-0.0132
→ 干渉 Norm=0.9661 e0=0.9661 e1=0.0000
[通常の子]
Depth=5 Norm=1.1485 e0=0.9661 e1=0.6172
[共役の子]
Depth=5 Norm=1.1485 e0=0.9661 e1=-0.6172
[共役の子]
Depth=2 Norm=1.0542 e0=1.0486 e1=-0.1077
→ 干渉 Norm=1.0433 e0=1.0433 e1=0.0000
[通常の子]
Depth=3 Norm=1.0294 e0=1.0223 e1=0.1147
→ 干渉 Norm=0.8888 e0=0.8888 e1=0.0000
[通常の子]
Depth=4 Norm=0.8967 e0=0.7827 e1=0.4325
→ 干渉 Norm=0.4428 e0=0.4428 e1=0.0000
[通常の子]
Depth=5 Norm=1.1320 e0=0.4428 e1=1.0355
[共役の子]
Depth=5 Norm=1.1320 e0=0.4428 e1=-1.0355
[共役の子]
Depth=4 Norm=1.0868 e0=0.9948 e1=-0.4325
→ 干渉 Norm=1.1497 e0=1.1497 e1=0.0000
[通常の子]
Depth=5 Norm=1.1506 e0=1.1497 e1=0.0130
[共役の子]
Depth=5 Norm=1.1506 e0=1.1497 e1=-0.0130
[共役の子]
Depth=3 Norm=1.0710 e0=1.0642 e1=-0.1147
→ 干渉 Norm=1.0284 e0=1.0284 e1=0.0000
[通常の子]
Depth=4 Norm=0.9996 e0=0.9738 e1=0.2249
→ 干渉 Norm=0.7357 e0=0.7357 e1=0.0000
[通常の子]
Depth=5 Norm=1.1437 e0=0.7357 e1=0.8713
[共役の子]
Depth=5 Norm=1.1437 e0=0.7357 e1=-0.8713
[共役の子]
Depth=4 Norm=1.1064 e0=1.0831 e1=-0.2249
→ 干渉 Norm=1.1002 e0=1.1002 e1=0.0000
[通常の子]
Depth=5 Norm=1.1532 e0=1.1002 e1=0.3436
[共役の子]
Depth=5 Norm=1.1532 e0=1.1002 e1=-0.3436

C++Builder2007 ソースコード
//************************
//****  Unit1.h ***********
//************************

//—————————————————————————

#ifndef Unit1H
#define Unit1H
//—————————————————————————
#include <Classes.hpp>
#include <Controls.hpp>
#include <StdCtrls.hpp>
#include <Forms.hpp>
//—————————————————————————
struct Octonion {
double e[8];

Octonion(double a0=0.0, double a1=0.0, double a2=0.0, double a3=0.0,
double a4=0.0, double a5=0.0, double a6=0.0, double a7=0.0);

Octonion operator+(const Octonion& o) const;
Octonion operator*(const Octonion& o) const; // ← 変更：完全乗算
Octonion operator*(double scalar) const; // スカラー乗算は残す
double Norm() const;

Octonion Conjugate() const; // 標準的な共役
Octonion PhaseConjugate() const; // 位相共役用（必要に応じて拡張可能）
};

class RecursiveOcto {
public:
int depth;
Octonion value;
RecursiveOcto* child; // 通常の子
RecursiveOcto* childConj; // 共役の子（追加）
Octonion interference; // 干渉結果を保持するメンバーを追加

RecursiveOcto(int d, Octonion init);
~RecursiveOcto();
void Generate(int maxDepth);
void Dump(TMemo* memo);
};

class TForm1 : public TForm
{
__published: // IDE-managed Components
TButton *ButtonExit;
TButton *Button1;
TMemo *Memo1;
void __fastcall ButtonExitClick(TObject *Sender);
void __fastcall Button1Click(TObject *Sender);
private: // User declarations
public: // User declarations
__fastcall TForm1(TComponent* Owner);
};
//—————————————————————————
extern PACKAGE TForm1 *Form1;
//—————————————————————————
#endif

//************************
//****  Unit1.cpp *********
//************************

//—————————————————————————

#include <vcl.h>
#pragma hdrstop

#include <math.h>
#include “Unit1.h”
//—————————————————————————
#pragma package(smart_init)
#pragma resource “*.dfm”
TForm1 *Form1;
//—————————————————————————
__fastcall TForm1::TForm1(TComponent* Owner)
: TForm(Owner)
{
}
//—————————————————————————
void __fastcall TForm1::ButtonExitClick(TObject *Sender)
{
Close();
}
//—————————————————————————
//—————————————————————————

Octonion::Octonion(double a0, double a1, double a2, double a3,
double a4, double a5, double a6, double a7) {
e[0]=a0; e[1]=a1; e[2]=a2; e[3]=a3;
e[4]=a4; e[5]=a5; e[6]=a6; e[7]=a7;
}

Octonion Octonion::Conjugate() const {
return Octonion(
e[0],
-e[1], -e[2], -e[3], -e[4], -e[5], -e[6], -e[7]
);
}

Octonion Octonion::PhaseConjugate() const {
return Conjugate();
}

Octonion Octonion::operator+(const Octonion& o) const {
Octonion res;
for(int i=0; i<8; i++) res.e[i] = e[i] + o.e[i];
return res;
}

RecursiveOcto::RecursiveOcto(int d, Octonion init)
: depth(d), value(init), child(NULL), childConj(NULL), interference(Octonion()) {}

RecursiveOcto::~RecursiveOcto() {
delete child;
delete childConj;
}

// スカラー乗算は残す
Octonion Octonion::operator*(double scalar) const {
Octonion res;
for(int i=0; i<8; i++) res.e[i] = e[i] * scalar;
return res;
}

// 完全な八元数乗算
Octonion Octonion::operator*(const Octonion& o) const {
Octonion res;

// e0（実部）
res.e[0] = e[0]*o.e[0] – e[1]*o.e[1] – e[2]*o.e[2] – e[3]*o.e[3]
– e[4]*o.e[4] – e[5]*o.e[5] – e[6]*o.e[6] – e[7]*o.e[7];

// e1
res.e[1] = e[0]*o.e[1] + e[1]*o.e[0] + e[2]*o.e[3] – e[3]*o.e[2]
+ e[4]*o.e[5] – e[5]*o.e[4] – e[6]*o.e[7] + e[7]*o.e[6];

// e2
res.e[2] = e[0]*o.e[2] – e[1]*o.e[3] + e[2]*o.e[0] + e[3]*o.e[1]
+ e[4]*o.e[6] + e[5]*o.e[7] – e[6]*o.e[4] – e[7]*o.e[5];

// e3
res.e[3] = e[0]*o.e[3] + e[1]*o.e[2] – e[2]*o.e[1] + e[3]*o.e[0]
+ e[4]*o.e[7] – e[5]*o.e[6] + e[6]*o.e[5] – e[7]*o.e[4];

// e4
res.e[4] = e[0]*o.e[4] – e[1]*o.e[5] – e[2]*o.e[6] – e[3]*o.e[7]
+ e[4]*o.e[0] + e[5]*o.e[1] + e[6]*o.e[2] + e[7]*o.e[3];

// e5
res.e[5] = e[0]*o.e[5] + e[1]*o.e[4] – e[2]*o.e[7] + e[3]*o.e[6]
– e[4]*o.e[1] + e[5]*o.e[0] – e[6]*o.e[3] + e[7]*o.e[2];

// e6
res.e[6] = e[0]*o.e[6] + e[1]*o.e[7] + e[2]*o.e[4] – e[3]*o.e[5]
– e[4]*o.e[2] + e[5]*o.e[3] + e[6]*o.e[0] – e[7]*o.e[1];

// e7
res.e[7] = e[0]*o.e[7] – e[1]*o.e[6] + e[2]*o.e[5] + e[3]*o.e[4]
– e[4]*o.e[3] – e[5]*o.e[2] + e[6]*o.e[1] + e[7]*o.e[0];

return res;
}

double Octonion::Norm() const {
double sum = 0.0;
for(int i=0; i<8; i++) sum += e[i]*e[i];
return sqrt(sum);
}

void RecursiveOcto::Generate(int maxDepth) {
if (depth >= maxDepth) return;

// 通常の次ステップ
Octonion rotator(cos(depth * 0.15), sin(depth * 0.15), 0.05, 0, 0, 0, 0, 0);
Octonion next = value * rotator;
next = next + Octonion(0.03, 0, 0, 0, 0, 0, 0, 0);

// 共役版
Octonion nextConj = next.PhaseConjugate();

// 子を生成
child = new RecursiveOcto(depth + 1, next);
child->Generate(maxDepth);

childConj = new RecursiveOcto(depth + 1, nextConj);
childConj->Generate(maxDepth);

// ===== 干渉計算 =====
interference = (child->value + childConj->value) * 0.5;

// ===== 親へのフィードバック（ここが追加） =====
// 干渉結果を親のvalueに少し混ぜる（フィードバック強度は調整可能）
double feedbackStrength = 0.3; // 0.1?0.5くらいで試すと良い
value = value * (1.0 – feedbackStrength) + interference * feedbackStrength;
}

void RecursiveOcto::Dump(TMemo* memo)
{
if (!memo) return;

Octonion conj = value.PhaseConjugate();
memo->Lines->Add(Format(“Depth=%d Norm=%.4f e0=%.4f e1=%.4f”,
ARRAYOFCONST((depth, value.Norm(), value.e[0], value.e[1]))));

// 干渉結果を表示（子がある場合）
if (child && childConj) {
memo->Lines->Add(Format(” → 干渉 Norm=%.4f e0=%.4f e1=%.4f”,
ARRAYOFCONST((interference.Norm(), interference.e[0], interference.e[1]))));
}

if (child) {
memo->Lines->Add(” [通常の子]”);
child->Dump(memo);
}
if (childConj) {
memo->Lines->Add(” [共役の子]”);
childConj->Dump(memo);
}
}

//—————————————————————————
//—————————————————————————

void __fastcall TForm1::Button1Click(TObject *Sender)
{
RecursiveOcto* root = new RecursiveOcto(0, Octonion(1.0));
root->Generate(5); // ← 2分岐なので深さは浅めに（5?6推奨）

Memo1->Clear();
root->Dump(Memo1);

delete root;
}
//—————————————————————————

https://youtu.be/-IfEcZd4rjA?si=eeEf7_jIIsMZwTZ8

https://youtu.be/xb_ykQwRJro?si=37SqWPJbT3c4hHyO

https://youtu.be/M_A-xgK2suE?si=r0yCzERn9DiTtasv

https://youtu.be/uj-iw8MFQTk?si=tylY9Al7EVhpQDgn

https://youtu.be/toVI4SdE5iU?si=IHPeWw3NNYaqNNh8

https://youtu.be/CZS66WNQW4s?si=JJMcQyhNPzQP6H2Q

https://youtu.be/mEBRdWlCXYM?si=RCSF_1WZ2uNFRVxf

https://youtu.be/x-PpHlcpwzM?si=3XVEO0xUZQ0uGp6w

https://youtu.be/zbafuZXVQhI?si=5ii26vyle2wG0kz6

https://ikuro-kotaro.sakura.ne.jp/koramu/1869_mz.htm

https://youtu.be/7fJ0hRkVvB4?si=hPIAszrSyJfuPFS6

https://youtu.be/xl5oLyz5VyM?si=gpQMEGMPNSWylA0O

https://youtu.be/KiE6nvZYDe4?si=1EvyzMqqTTOeupvb

https://youtu.be/WHSz53mq2_U?si=IQDnATTD7N4bw2JW

https://youtu.be/IY00QTsO8tM?si=fIVmd17g0_hfgBQ-

https://youtu.be/Q4AD0HjnZPw?si=i2HR4u6W-r-JS636

https://youtu.be/rs1Qlir1-O4?si=M55BtcW6PTYHToE5

https://youtu.be/kgK6ohO5HXI?si=ZrpR3DBjp3AhzU7s

https://youtu.be/Po8m8bwg3TU?si=SP9gGLevB3zxEzt2

https://youtu.be/3aVl-_i-4cM?si=30sLCfXQCLjWWukb

https://youtu.be/JLwvyMvonr8?si=xtz_79a0wKsQ6L57

https://youtu.be/qOKRlzciulM?si=l34_H7su5hNV0xIb

https://youtu.be/r1n7DS7UeCw?si=3p3T4NKX8XYsY8_m

https://youtu.be/KKIXHfw9GwA?si=sHs-UrImHJs9EnW5

https://youtu.be/KyT1bor3tt0?si=TCE0dPKxbtagvJTb

https://youtu.be/2kAMbIouVjM?si=e8sOO1YxEjdzyqN0

https://youtu.be/43323VJXry8?si=JZ61-1VOU2uAVaE5

https://youtu.be/-8qKieFfh08?si=ef0iiniV0HgioZjb

https://youtu.be/MDTyciV1sFo?si=BOH8aL2EXxVveo1-

https://youtu.be/vb8kbHbvpPw?si=WuGqWOB6zDb-vHje

https://youtu.be/2kAMbIouVjM?si=yMjFpDOTbmFBr40C

https://youtu.be/jlGvb6qaT7w?si=ZxcJJFHkGKj-vUgo

https://youtu.be/MXTuEwwTUKM?si=JYYLZOH9CewNqH-l

https://youtu.be/KU7mraaXsYI?si=GUkcF5dz8q8i5-57

https://youtu.be/J4ZjrHZF3mU?si=xQzdWzG1bJmMQBrM

https://youtu.be/8R1XTqtWl5o?si=W3kL87gNBAWmhSK5

https://youtu.be/I31NHVvvUSQ?si=pv6FWomBkgkqRiIQ

https://youtu.be/r7aNDwyfuW4?si=ZxNN1O3OUEE_NO2y

https://youtu.be/fm5snCUwl_c?si=kXjjS9Ozpm_Wh40_

https://youtu.be/3TUZZjSrYtg?si=7KTsXL-Gan2BzECl

https://youtu.be/Hr019gUsjKA?si=0lLMuefE1K3t9nF3

https://youtu.be/Zp4J14T1o-A?si=hx-BHKb8Dj-cPw_T

https://youtu.be/Tyca1B_2zrI?si=0OBFAuxncjSX7ib6

https://youtu.be/XVT5yxpGz9w?si=W7ztkhiS68y9ur8u

https://youtu.be/Rgvnq23JhYI?si=nqSCjRQvLeHxqcER

https://youtu.be/_tcyjQfWMFE?si=0NBgRxnoX466gzJq

実行結果

C++Builder2007ソースコード

//************************
//***** Unit1.h ********
//************************

//—————————————————————————

#ifndef Unit1H
#define Unit1H
//—————————————————————————
#include <Classes.hpp>
#include <Controls.hpp>
#include <StdCtrls.hpp>
#include <Forms.hpp>
#include “Chart.hpp”
#include “TeEngine.hpp”
#include “TeeProcs.hpp”
#include <ExtCtrls.hpp>
#include “Series.hpp”

#include <vector>
//—————————————————————————
const int max_depth = 8; // 最大階層数

class RecursiveWave
{
public:
int current_depth; // クラスメンバーとして
double Value;
RecursiveWave *child1, *child2;
RecursiveWave *z_1, *z_2;
int direction; // +1: 順方向、 -1:逆方向

// RecursiveWaveクラスに追加（public） 逆転版（トップダウン用）
bool isTopDownMode;

RecursiveWave(int n, int dir);
~RecursiveWave();
void LinkParents(RecursiveWave* parent1, RecursiveWave* parent2); // 後からリンク
void SetupParents(RecursiveWave* node, RecursiveWave* p1, RecursiveWave* p2);
void Dump(TMemo *memo);
double GetValue(){ return Value; }
};

std::vector<RecursiveWave*> all_waves1;
std::vector<RecursiveWave*> all_waves2;

//—————————————————————————
class TForm1 : public TForm
{
__published: // IDE-managed Components
TChart *Chart1;
TButton *ButtonExit;
TLineSeries *Series1;
TMemo *Memo1;
TLineSeries *Series2;
TButton *Button1;
TRadioGroup *RadioGroup1;
void __fastcall ButtonExitClick(TObject *Sender);
void __fastcall Button1Click(TObject *Sender);
private: // User declarations
public: // User declarations
__fastcall TForm1(TComponent* Owner);

RecursiveWave *Top;
void ChartSetting();

};
//—————————————————————————
extern PACKAGE TForm1 *Form1;
//—————————————————————————
#endif

//************************
//***** Unit1.cpp ******
//************************

//—————————————————————————

#include <vcl.h>
#pragma hdrstop

#include <math.h>
#include “Unit1.h”
//—————————————————————————
#pragma package(smart_init)
#pragma link “Chart”
#pragma link “TeEngine”
#pragma link “TeeProcs”
#pragma link “Series”
#pragma resource “*.dfm”
TForm1 *Form1;
//—————————————————————————
__fastcall TForm1::TForm1(TComponent* Owner)
: TForm(Owner)
{
}
//—————————————————————————
void __fastcall TForm1::ButtonExitClick(TObject *Sender)
{
Close();
}
//—————————————————————————
// ***** ChartSetting *****
void TForm1::ChartSetting()
{
Chart1->SeriesList->ClearValues();

// 軸タイトルだけ設定（SeriesCountチェックはButton1Click側に移動）
Chart1->LeftAxis->Title->Caption = “Value1”;
Chart1->LeftAxis->Title->Visible = true;
Chart1->RightAxis->Title->Caption = “Value2”;
Chart1->RightAxis->Title->Visible = true;
Chart1->RightAxis->Visible = true;
Chart1->BottomAxis->Title->Caption = “Count”;
Chart1->BottomAxis->Visible = true;

}

// 構築後に親を設定する専用メソッド
void RecursiveWave::LinkParents(RecursiveWave* parent1, RecursiveWave* parent2) {
z_1 = parent1;
z_2 = parent2 ? parent2 : parent1; // z_2がなければz_1で代用

// Valueを親を使って計算
if (z_1) {
Value = z_1->Value * 1.7 – (z_2 ? z_2->Value * 0.97 : 0.0) + 1.0;
} else {
Value = 1.0; // 最上位の初期値
}
}

// 新規ヘルパー関数（トップダウンで親情報を伝播）
void RecursiveWave::SetupParents(RecursiveWave* node, RecursiveWave* p1, RecursiveWave* p2) {
if (!node) return;

node->LinkParents(p1, p2);
//node->LinkParents();

// 子に伝える（z_1 = node自身、z_2 = p1）
SetupParents(node->child1, node, p1); //再帰呼び出し
SetupParents(node->child2, node, p1); //再帰呼び出し
}

// ***** 順方向：逆方向選択版 コンストラクタ *****
RecursiveWave::RecursiveWave(int n, int dir): current_depth(n), Value(0.0), child1(NULL), child2(NULL),
z_1(NULL), z_2(NULL), direction(dir)
{
if (n >= max_depth) return;

if(direction == -1){
//逆方向
// 子を先に作成（再帰）
child1 = new RecursiveWave(n + 1, dir);
child2 = new RecursiveWave(n + 1, -dir);

// all_waves登録
all_waves1.push_back(this);
all_waves2.push_back(this);

}else if(direction == +1){
//順方向
// all_wavesに先に登録
all_waves1.push_back(this);
all_waves2.push_back(this);

// 子を先に作成（再帰）
child1 = new RecursiveWave(n + 1, dir);
child2 = new RecursiveWave(n + 1, -dir);

}

}

//**********************
//デストラクタ
RecursiveWave::~RecursiveWave()
{
if (child1 != NULL) {
delete child1; // 再帰的に子孫をすべて削除
child1 = NULL;
}

if (child2 != NULL) {
delete child1; // 再帰的に子孫をすべて削除
child1 = NULL;
}
// all_wavesから自分を除去（任意）
// （vectorからeraseするのは重いので、クリア時に一括処理する方が良い）
}

void RecursiveWave::Dump(TMemo *memo)
{
if (all_waves1.empty()) return;
memo->Hide();

AnsiString str;
int count = 0;

for (std::vector<RecursiveWave*>::const_iterator it = all_waves1.begin();
it != all_waves1.end(); ++it) {
RecursiveWave* neigh = *it;
if (neigh) {
memo->Lines->Add(str.sprintf(“depth= %d, Value= %g”, neigh->current_depth, neigh->Value) );
count++;
}
}

memo->Show();

/*memo->Lines->Add(“==============”);
for(int i=0; i<all_waves->size(); i++){
memo->Lines->Add(str.sprintf(“depth= %d, Value= %g”, i, (*all_waves)[i].Value ) );
}*/

}

void __fastcall TForm1::Button1Click(TObject *Sender)
{
// 既存のオブジェクトをクリア（メモリリーク防止）
for (size_t i = 0; i < all_waves1.size(); i++) {
delete all_waves1[i];
}
all_waves1.clear();
// 既存のオブジェクトをクリア（メモリリーク防止）
for (size_t i = 0; i < all_waves2.size(); i++) {
delete all_waves2[i];
}
all_waves2.clear();

//方向指定、初期化、生成
RecursiveWave *WaveTop;
if(RadioGroup1->ItemIndex == 1){
//順方向
WaveTop = new RecursiveWave(0, +1); // 引数不要！

}else if(RadioGroup1->ItemIndex == 0){
//逆方向
WaveTop = new RecursiveWave(0, -1); // 引数不要！
}

// 再帰的に親リンクを設定（トップダウン）
WaveTop->SetupParents(WaveTop, NULL, NULL); // 新規ヘルパー関数

WaveTop->Dump(Memo1);

ChartSetting();

// Seriesを安全に追加
Chart1->SeriesList->ClearValues(); // 念のため

// Series[0] と Series[1] が確実に存在することを保証
while (Chart1->SeriesCount() < 2) {
TLineSeries* s = new TLineSeries(Chart1);
Chart1->AddSeries(s);
}

int count = 0;
for (std::vector<RecursiveWave*>::const_iterator it = all_waves1.begin();
it != all_waves1.end(); ++it) {
RecursiveWave* neigh = *it;
if (neigh) {
Chart1->Series[0]->AddY(neigh->Value, “”, clBlue); // Value（左軸）
//Chart1->Series[0]->AddY(neigh->current_depth, “”, clBlue); // Depth（右軸） ← コメント解除OK
count++;
}
}
count = 0;
for (std::vector<RecursiveWave*>::const_iterator it = all_waves2.begin();
it != all_waves2.end(); ++it) {
RecursiveWave* neigh = *it;
if (neigh) {
Chart1->Series[1]->AddY(neigh->Value, “”, clRed); // Value（左軸）
//Chart1->Series[1]->AddY(neigh->current_depth, “”, clRed); // Depth（右軸） ← コメント解除OK
count++;
}
}

// 軸設定を再適用
Chart1->Series[0]->VertAxis = 0; // vaLeft
Chart1->Series[1]->VertAxis = 1; // vaRight

delete WaveTop;
all_waves1.clear();
all_waves2.clear();

}

//—————————————————————————

https://youtu.be/_tcyjQfWMFE?si=Mq_Zcw_aKWNyJOPD