RenderingMandelbrotSetInCpp/main.cpp at main · Xue-Zhoujun-RhF6/RenderingMandelbrotSetInCpp · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
#define _USE_MATH_DEFINES
#include <iostream>
#include <Windows.h>
#include <fstream>
#include <vector>
#include <thread>
#include <easyx.h>
#include <graphics.h>
#define KEY_DOWN(VK_NONAME) ((GetAsyncKeyState(VK_NONAME) & 0x8000) ? 1:0)

using namespace std;
using ld = long double;
using ull = unsigned long long;

ull itPrec = 100; //迭代精度
ull cpPrec = 200; //计算每1^2方格的边长分值
ld vFs = 2.0; //判定发散值
ull mColorR = 114, mColorG = 51, mColorB = 4; //颜色值乘的
ull dColorR = 0, dColorG = 0, dColorB = 0; //颜色值加的
ull ctrlM = 2; //鼠标点击缩放倍数
ld sx1, sx2, sy1, sy2; //上次渲染的坐标
HWND jiaodian, tJiaodian; //绘图窗口的句柄与当前窗口的句柄
ull threads = 2;

void cp_q(ld x1, ld y1, ld x2, ld y2) {
	sx1 = x1, sx2 = x2, sy1 = y1, sy2 = y2;
	ld delta_c = 1.0 / cpPrec;
	for (ld re = x1; re < x2; re += delta_c) {
		for (ld im = y1; im < y2; im += delta_c) {
			ull tcl = 0;
			ld rec = 0, imc = 0, re2c = 0, im2c = 0, re_oldc = 0; //快速计算点
			for (ull i = 0; i < itPrec; i++) {
				rec = re2c - im2c + re;
				imc = 2 * re_oldc * imc + im;
				re_oldc = rec;
				re2c = rec * rec;
				im2c = imc * imc;
				if (sqrt(re2c + im2c) > vFs) {
					tcl = 200 * i;
					break;
				}
			}
			putpixel(((re - x1) * cpPrec), ((im - y1) * cpPrec), RGB(mColorR * tcl + dColorR, mColorG * tcl + dColorG, mColorB * tcl + dColorB));
		}
	}
}

void frm() {
	tJiaodian = GetForegroundWindow();
	POINT p;
	GetCursorPos(&p);
	ScreenToClient(tJiaodian, &p);
	ld x1, y1, x2, y2;
	x1 = 1.0 * p.x / cpPrec + sx1;
	y1 = 1.0 * p.y / cpPrec + sy1;
	if (KEY_DOWN(VK_LBUTTON) && tJiaodian == jiaodian && p.x>0 && p.y>0) {
		cpPrec = cpPrec * ctrlM;
		x2 = x1 + (sx2 - sx1) / ctrlM * 0.5;
		y2 = y1 + (sy2 - sy1) / ctrlM * 0.5;
		x1 = x1 - (sx2 - sx1) / ctrlM * 0.5;
		y1 = y1 - (sy2 - sy1) / ctrlM * 0.5;
		cout << x1 << " " << y1 << " " << x2 << " " << y2 << "\n精度1/" << cpPrec << "\n";
		cleardevice();
		cp_q(x1, y1, x2, y2);
		Sleep(100);
	}
	else if (KEY_DOWN(VK_RBUTTON) && tJiaodian == jiaodian && p.x > 0 && p.y > 0) {
		cpPrec = cpPrec / ctrlM;
		x2 = x1 + (sx2 - sx1) * ctrlM * 0.5;
		y2 = y1 + (sy2 - sy1) * ctrlM * 0.5;
		x1 = x1 - (sx2 - sx1) * ctrlM * 0.5;
		y1 = y1 - (sy2 - sy1) * ctrlM * 0.5;
		cout << x1 << " " << y1 << " " << x2 << " " << y2 << "\n精度1/" << cpPrec << "\n";
		cleardevice();
		cp_q(x1, y1, x2, y2);
		Sleep(100);
	}
    Sleep(50);
}

int main() {
	ios::sync_with_stdio(0);
	int tm1 = clock();
	initgraph(800, 800, EX_SHOWCONSOLE);
	jiaodian = GetForegroundWindow();
	cp_q(-2, -2, 2, 2);
	cout << clock() - tm1;
	while (1) frm();
	closegraph();
	return 0;
}