1
0

MandelWidget.cpp 25 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901902903904905906907908909910911912913914915916917918919920921922923924925926927928929930931932933934935936937938939940941942943944945946947948949950951952953954955956957958959960961962963964965966967968969970971972973974975976977978979980981982983984985
  1. #include "MandelWidget.h"
  2. #include <cmath>
  3. #include <sstream>
  4. using namespace mnd;
  5. #include <QPainter>
  6. #include <cstdio>
  7. Texture::Texture(const Bitmap<RGBColor>& bitmap, GLint param)
  8. {
  9. glGenTextures(1, &id);
  10. glBindTexture(GL_TEXTURE_2D, id);
  11. //int lineLength = (bitmap.width * 3 + 3) & ~3;
  12. /*std::unique_ptr<unsigned char[]> pixels = std::make_unique<unsigned char[]>(lineLength * bitmap.height);
  13. for (int i = 0; i < bitmap.width; i++) {
  14. for (int j = 0; j < bitmap.height; j++) {
  15. int index = i * 3 + j * lineLength;
  16. RGBColor c = bitmap.get(i, j);
  17. pixels[index] = c.r;
  18. pixels[index + 1] = c.g;
  19. pixels[index + 2] = c.b;
  20. }
  21. }*/
  22. glTexImage2D(GL_TEXTURE_2D, 0, GL_RGB, int(bitmap.width), int(bitmap.height), 0, GL_RGB, GL_UNSIGNED_BYTE, reinterpret_cast<char*> (bitmap.pixels.get()));
  23. glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_CLAMP_TO_EDGE);
  24. glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_CLAMP_TO_EDGE);
  25. glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, param);
  26. glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, param);
  27. }
  28. Texture::~Texture(void)
  29. {
  30. if (id != 0)
  31. glDeleteTextures(1, &id);
  32. }
  33. Texture::Texture(Texture&& other) :
  34. id{ other.id }
  35. {
  36. other.id = 0;
  37. }
  38. Texture& Texture::operator=(Texture&& other)
  39. {
  40. this->id = other.id;
  41. other.id = 0;
  42. return *this;
  43. }
  44. void Texture::bind(void) const
  45. {
  46. glBindTexture(GL_TEXTURE_2D, id);
  47. }
  48. void Texture::drawRect(float x, float y, float width, float height)
  49. {
  50. glColor3ub(255, 255, 255);
  51. glEnable(GL_TEXTURE_2D);
  52. bind();
  53. glBegin(GL_TRIANGLE_STRIP);
  54. glTexCoord2f(0, 0);
  55. glVertex2f(x, y);
  56. glTexCoord2f(1, 0);
  57. glVertex2f(x + width, y);
  58. glTexCoord2f(0, 1);
  59. glVertex2f(x, y + height);
  60. glTexCoord2f(1, 1);
  61. glVertex2f(x + width, y + height);
  62. glEnd();
  63. glDisable(GL_TEXTURE_2D);
  64. }
  65. CellImage::~CellImage(void)
  66. {
  67. }
  68. TextureClip::~TextureClip(void)
  69. {
  70. }
  71. void TextureClip::drawRect(float x, float y, float width, float height)
  72. {
  73. glColor3ub(255, 255, 255);
  74. glEnable(GL_TEXTURE_2D);
  75. glBindTexture(GL_TEXTURE_2D, texture->getId());
  76. glBegin(GL_TRIANGLE_STRIP);
  77. glTexCoord2f(tx, ty);
  78. glVertex2f(x, y);
  79. glTexCoord2f(tx + tw, ty);
  80. glVertex2f(x + width, y);
  81. glTexCoord2f(tx, ty + th);
  82. glVertex2f(x, y + height);
  83. glTexCoord2f(tx + tw, ty + th);
  84. glVertex2f(x + width, y + height);
  85. glEnd();
  86. glDisable(GL_TEXTURE_2D);
  87. }
  88. TextureClip TextureClip::clip(float x, float y, float w, float h)
  89. {
  90. float tx = this->tx + x * this->tw;
  91. float ty = this->ty + y * this->th;
  92. float tw = this->tw * w;
  93. float th = this->th * h;
  94. return TextureClip{ this->texture, tx, ty, tw, th };
  95. }
  96. std::shared_ptr<CellImage> TextureClip::clip(short i, short j)
  97. {
  98. return std::make_shared<TextureClip>(clip(i * 0.5f, j * 0.5f, 0.5f, 0.5f));
  99. }
  100. int TextureClip::getRecalcPriority() const
  101. {
  102. return int(1.0f / tw);
  103. }
  104. QuadImage::~QuadImage(void)
  105. {
  106. }
  107. void QuadImage::drawRect(float x, float y, float width, float height)
  108. {
  109. for (int i = 0; i < 2; i++) {
  110. for (int j = 0; j < 2; j++) {
  111. this->cells[i][j]->drawRect(x + i * 0.5f * width,
  112. y + j * 0.5f * height,
  113. width * 0.5f,
  114. height * 0.5f);
  115. }
  116. }
  117. }
  118. std::shared_ptr<CellImage> QuadImage::clip(short i, short j)
  119. {
  120. return cells[i][j];
  121. }
  122. int QuadImage::getRecalcPriority() const
  123. {
  124. return 1;
  125. }
  126. TexGrid::TexGrid(MandelView& owner, int level) :
  127. owner{ owner },
  128. level{ level },
  129. dpp{ owner.getDpp(level) }
  130. {
  131. }
  132. std::pair<GridIndex, GridIndex> TexGrid::getCellIndices(mnd::Real x, mnd::Real y)
  133. {
  134. return { GridIndex(mnd::floor(x / dpp / MandelView::chunkSize)), GridIndex(mnd::floor(y / dpp / MandelView::chunkSize)) };
  135. }
  136. std::pair<mnd::Real, mnd::Real> TexGrid::getPositions(GridIndex x, GridIndex y)
  137. {
  138. return { mnd::Real(x) * dpp * MandelView::chunkSize, mnd::Real(y) * dpp * MandelView::chunkSize };
  139. }
  140. GridElement* TexGrid::getCell(GridIndex i, GridIndex j)
  141. {
  142. auto cIt = cells.find({i, j});
  143. if (cIt != cells.end()) {
  144. return cIt->second.get();
  145. }
  146. else {
  147. return nullptr;
  148. }
  149. }
  150. void TexGrid::setCell(GridIndex i, GridIndex j, std::unique_ptr<GridElement> tex)
  151. {
  152. cells[{i, j}] = std::move(tex);
  153. }
  154. void TexGrid::clearCells(void)
  155. {
  156. cells.clear();
  157. }
  158. void TexGrid::clearUncleanCells(void)
  159. {
  160. for (auto it = cells.begin(); it != cells.end();) {
  161. if (it->second->img->getRecalcPriority() > 1)
  162. cells.erase(it++);
  163. else ++it;
  164. }
  165. }
  166. void Job::run(void)
  167. {
  168. auto [absX, absY] = grid->getPositions(i, j);
  169. mnd::Real gw = grid->dpp * MandelView::chunkSize;
  170. Bitmap<float> f(MandelView::chunkSize, MandelView::chunkSize);
  171. mnd::MandelInfo mi = owner.getMandelInfo();
  172. mi.view.x = absX;
  173. mi.view.y = absY;
  174. mi.view.width = mi.view.height = gw;
  175. mi.bWidth = mi.bHeight = MandelView::chunkSize;
  176. try {
  177. generator->generate(mi, f.pixels.get());
  178. auto* rgb = new Bitmap<RGBColor>(f.map<RGBColor>([&mi, this] (float i) {
  179. return i >= mi.maxIter ? RGBColor{ 0, 0, 0 } : gradient.get(i);
  180. }));
  181. emit done(level, i, j, calcState, rgb);
  182. }
  183. catch(std::exception& ex) {
  184. printf("wat: %s?!\n", ex.what()); fflush(stdout);
  185. exit(1);
  186. }
  187. catch(...) {
  188. printf("wat?!\n"); fflush(stdout);
  189. exit(1);
  190. }
  191. }
  192. Calcer::Calcer(mnd::MandelGenerator* generator, MandelWidget& owner) :
  193. jobsMutex{ QMutex::Recursive },
  194. generator{ generator },
  195. threadPool{ std::make_unique<QThreadPool>() },
  196. owner{ owner },
  197. gradient{ owner.getGradient() }
  198. {
  199. threadPool->setMaxThreadCount(1);
  200. }
  201. void Calcer::clearAll(void)
  202. {
  203. this->threadPool->clear();
  204. }
  205. void Calcer::calc(TexGrid& grid, int level, GridIndex i, GridIndex j, int priority)
  206. {
  207. jobsMutex.lock();
  208. if (jobs.find({ level, i, j }) == jobs.end()) {
  209. Job* job = new Job(generator, gradient, owner, &grid, level, i, j, calcState);
  210. connect(job, &Job::done, this, &Calcer::redirect);
  211. connect(job, &QObject::destroyed, this, [this, level, i, j] () { this->notFinished(level, i, j); });
  212. jobs.emplace(std::tuple{level, i, j}, job);
  213. threadPool->start(job, priority);
  214. }
  215. jobsMutex.unlock();
  216. }
  217. void Calcer::setCurrentLevel(int level)
  218. {
  219. if (this->currentLevel != level) {
  220. this->currentLevel = level;
  221. std::vector<QRunnable*> toCancel;
  222. jobsMutex.lock();
  223. for (auto&[tup, job] : jobs) {
  224. auto& [level, i, j] = tup;
  225. if(level != currentLevel) {
  226. toCancel.push_back(job);
  227. }
  228. }
  229. jobsMutex.unlock();
  230. for (auto* job : toCancel) {
  231. if (threadPool->tryTake(job)) {
  232. delete job;
  233. }
  234. }
  235. }
  236. }
  237. void Calcer::notFinished(int level, GridIndex i, GridIndex j)
  238. {
  239. jobsMutex.lock();
  240. jobs.erase({ level, i, j });
  241. jobsMutex.unlock();
  242. }
  243. void Calcer::redirect(int level, GridIndex i, GridIndex j, long calcState, Bitmap<RGBColor>* bmp)
  244. {
  245. jobsMutex.lock();
  246. jobs.erase({ level, i, j });
  247. jobsMutex.unlock();
  248. if (this->calcState == calcState) {
  249. emit done(level, i, j, bmp);
  250. }
  251. else {
  252. delete bmp;
  253. }
  254. }
  255. const int MandelView::chunkSize = 256;
  256. MandelView::MandelView(mnd::MandelGenerator* generator, MandelWidget& owner) :
  257. generator{ generator },
  258. calcer{ generator, owner },
  259. owner{ owner }
  260. {
  261. /*Bitmap<RGBColor> emp(8, 8);
  262. for(auto i = 0; i < emp.width; i++) {
  263. for(auto j = 0; j < emp.height; j++) {
  264. if((i + j) & 0x1) { // if i+j is odd
  265. emp.get(i, j) = RGBColor{ 255, 255, 255 };
  266. }
  267. else {
  268. emp.get(i, j) = RGBColor{ 120, 120, 120 };
  269. }
  270. }
  271. }*/
  272. Bitmap<RGBColor> emp(1, 1);
  273. emp.get(0, 0) = RGBColor{ 0, 0, 0 };
  274. empty = std::make_unique<Texture>(emp, GL_NEAREST);
  275. connect(&calcer, &Calcer::done, this, &MandelView::cellReady);
  276. }
  277. int MandelView::getLevel(mnd::Real dpp)
  278. {
  279. return int(mnd::log2(dpp / chunkSize));
  280. }
  281. mnd::Real MandelView::getDpp(int level)
  282. {
  283. return mnd::pow(mnd::Real(2), mnd::Real(level)) * chunkSize;
  284. }
  285. TexGrid& MandelView::getGrid(int level)
  286. {
  287. auto it = levels.find(level);
  288. if (it != levels.end()) {
  289. return it->second;
  290. }
  291. else {
  292. levels.insert(std::pair<int, TexGrid>{ level, TexGrid{ *this, level } });
  293. return levels.at(level);
  294. }
  295. }
  296. void MandelView::setGenerator(mnd::MandelGenerator* generator)
  297. {
  298. if (this->generator != generator) {
  299. this->generator = generator;
  300. calcer.setGenerator(generator);
  301. clearCells();
  302. emit redrawRequested();
  303. }
  304. }
  305. void MandelView::clearCells(void)
  306. {
  307. for(auto& [level, grid] : this->levels) {
  308. grid.clearCells();
  309. }
  310. }
  311. void MandelView::garbageCollect(int level, GridIndex /*i*/, GridIndex /*j*/)
  312. {
  313. for(auto& [l, grid] : levels) {
  314. int dist = ::abs(l - level);
  315. if (dist == 1) {
  316. grid.clearUncleanCells();
  317. }
  318. if (dist > 20) {
  319. grid.clearCells();
  320. }
  321. else if (dist > 10) {
  322. if (grid.countAllocatedCells() > 50)
  323. grid.clearCells();
  324. }
  325. else if (dist > 3) {
  326. if (grid.countAllocatedCells() > 150)
  327. grid.clearCells();
  328. }
  329. else if (dist > 0) {
  330. if (grid.countAllocatedCells() > 350)
  331. grid.clearCells();
  332. }
  333. else {
  334. if (grid.countAllocatedCells() > 2500)
  335. grid.clearCells();
  336. }
  337. }
  338. }
  339. GridElement* MandelView::searchAbove(int level, GridIndex i, GridIndex j, int recursionLevel)
  340. {
  341. auto& grid = getGrid(level);
  342. auto& gridAbove = getGrid(level + 1);
  343. GridIndex ai = (i < 0 ? (i - 1) : i) / 2;
  344. GridIndex aj = (j < 0 ? (j - 1) : j) / 2;
  345. GridElement* above = gridAbove.getCell(ai, aj);
  346. if (above == nullptr && recursionLevel > 0) {
  347. auto abFound = searchAbove(level + 1, ai, aj, recursionLevel - 1);
  348. if (abFound)
  349. above = abFound;
  350. }
  351. if (above != nullptr) {
  352. auto newElement = std::make_unique<GridElement>(
  353. false, above->img->clip(short(i & 1), short(j & 1))
  354. );
  355. GridElement* ret = newElement.get();
  356. grid.setCell(i, j, std::move(newElement));
  357. return ret;
  358. }
  359. else {
  360. return nullptr;
  361. }
  362. }
  363. GridElement* MandelView::searchUnder(int level, GridIndex i, GridIndex j, int recursionLevel)
  364. {
  365. if (recursionLevel == 0)
  366. return nullptr;
  367. auto& grid = getGrid(level);
  368. auto& gridUnder = getGrid(level - 1);
  369. GridIndex ai = i * 2;
  370. GridIndex aj = j * 2;
  371. GridElement* u00 = gridUnder.getCell(ai, aj);
  372. GridElement* u01 = gridUnder.getCell(ai, aj + 1);
  373. GridElement* u10 = gridUnder.getCell(ai + 1, aj);
  374. GridElement* u11 = gridUnder.getCell(ai + 1, aj + 1);
  375. /*if ( u00 == nullptr
  376. || u01 == nullptr
  377. || u10 == nullptr
  378. || u11 == nullptr) {
  379. auto abFound = searchUnder(level + 1, ai, aj, recursionLevel - 1);
  380. if (abFound)
  381. above = abFound;
  382. }*/
  383. if ( u00 != nullptr
  384. && u01 != nullptr
  385. && u10 != nullptr
  386. && u11 != nullptr) {
  387. auto newElement = std::make_unique<GridElement>(
  388. false, std::make_shared<QuadImage>(u00->img, u01->img, u10->img, u11->img)
  389. );
  390. GridElement* ret = newElement.get();
  391. grid.setCell(i, j, std::move(newElement));
  392. return ret;
  393. }
  394. else {
  395. return nullptr;
  396. }
  397. }
  398. void MandelView::paint(const mnd::MandelViewport& mvp)
  399. {
  400. mnd::Real dpp = mvp.width / width;
  401. int level = getLevel(dpp) - 1;
  402. auto& grid = getGrid(level);
  403. mnd::Real gw = getDpp(level) * chunkSize;
  404. auto [left, top] = grid.getCellIndices(mvp.x, mvp.y);
  405. auto [right, bottom] = grid.getCellIndices(mvp.right(), mvp.bottom());
  406. garbageCollect(level, (left + right) / 2, (top + bottom) / 2);
  407. emit calcer.setCurrentLevel(level);
  408. mnd::Real w = width * gw / mvp.width;
  409. auto [realXLeft, realYTop] = grid.getPositions(left, top);
  410. realXLeft = ((realXLeft - mvp.x) * mnd::Real(width)) / mvp.width;
  411. realYTop = ((realYTop - mvp.y) * mnd::Real(height)) / mvp.height;
  412. for(GridIndex i = left; i <= right; i++) {
  413. for(GridIndex j = top; j <= bottom; j++) {
  414. mnd::Real x = w * int(i - left) + realXLeft;
  415. mnd::Real y = w * int(j - top) + realYTop;
  416. GridElement* t = grid.getCell(i, j);
  417. if (t == nullptr) {
  418. auto under = searchUnder(level, i, j, 1);
  419. if (under) {
  420. t = under;
  421. }
  422. else {
  423. auto above = searchAbove(level, i, j, 3);
  424. if (above) {
  425. t = above;
  426. }
  427. }
  428. }
  429. if (t != nullptr) {
  430. t->img->drawRect(float(x), float(y), float(w), float(w));
  431. /*glBegin(GL_LINE_LOOP);
  432. glVertex2f(float(x), float(y));
  433. glVertex2f(float(x) + float(w), float(y));
  434. glVertex2f(float(x) + float(w), float(y) + float(w));
  435. glVertex2f(float(x), float(y) + float(w));
  436. glEnd();*/
  437. if (!t->enoughResolution) {
  438. calcer.calc(grid, level, i, j, t->img->getRecalcPriority());
  439. }
  440. }
  441. else {
  442. calcer.calc(grid, level, i, j, 1000);
  443. this->empty->drawRect(float(x), float(y), float(w), float(w));
  444. }
  445. }
  446. }
  447. }
  448. void MandelView::cellReady(int level, GridIndex i, GridIndex j, Bitmap<RGBColor>* bmp)
  449. {
  450. this->getGrid(level).setCell(i, j,
  451. std::make_unique<GridElement>(true, std::make_shared<TextureClip>(std::make_shared<Texture>(*bmp))));
  452. delete bmp;
  453. emit redrawRequested();
  454. }
  455. MandelWidget::MandelWidget(mnd::MandelContext& ctxt, mnd::MandelGenerator* generator, QWidget* parent) :
  456. QOpenGLWidget{ parent },
  457. mndContext{ ctxt },
  458. generator{ generator },
  459. gradient{ Gradient::defaultGradient() }
  460. {
  461. this->setContentsMargins(0, 0, 0, 0);
  462. this->setSizePolicy(QSizePolicy::Expanding,
  463. QSizePolicy::Expanding);
  464. qRegisterMetaType<GridIndex>("GridIndex");
  465. this->format().setSwapInterval(1);
  466. /*gradient = Gradient {
  467. {
  468. { RGBColor{ 0, 0, 0 }, 0 },
  469. { RGBColor{ 180, 20, 10 }, 30 },
  470. { RGBColor{ 210, 180, 15 }, 70 },
  471. { RGBColor{ 160, 220, 45 }, 170 },
  472. { RGBColor{ 50, 150, 170 }, 300 },
  473. }
  474. };*/
  475. }
  476. MandelWidget::~MandelWidget()
  477. {
  478. }
  479. void MandelWidget::setGradient(Gradient g)
  480. {
  481. this->gradient = std::move(g);
  482. if (mandelView) {
  483. mandelView->clearCells();
  484. mandelView->calcer.changeState();
  485. }
  486. emit update();
  487. }
  488. void MandelWidget::setSmoothColoring(bool sc)
  489. {
  490. if (sc != mandelInfo.smooth) {
  491. mandelInfo.smooth = sc;
  492. if (mandelView) {
  493. mandelView->clearCells();
  494. emit update();
  495. }
  496. }
  497. }
  498. void MandelWidget::setDisplayInfo(bool di)
  499. {
  500. if (di != this->displayInfo) {
  501. this->displayInfo = di;
  502. emit update();
  503. }
  504. }
  505. void MandelWidget::setMaxIterations(int maxIter)
  506. {
  507. if (mandelInfo.maxIter != maxIter) {
  508. mandelInfo.maxIter = maxIter;
  509. if (mandelView) {
  510. mandelView->clearCells();
  511. mandelView->calcer.clearAll();
  512. mandelView->calcer.changeState();
  513. }
  514. emit update();
  515. }
  516. }
  517. void MandelWidget::setJuliaPos(const mnd::Real& x, const mnd::Real& y)
  518. {
  519. mandelInfo.juliaX = x;
  520. mandelInfo.juliaY = y;
  521. if (mandelView)
  522. mandelView->calcer.changeState();
  523. emit update();
  524. }
  525. void MandelWidget::setGenerator(mnd::MandelGenerator* generator)
  526. {
  527. if (this->generator != generator) {
  528. this->generator = generator;
  529. if (mandelView)
  530. mandelView->setGenerator(generator);
  531. }
  532. }
  533. void MandelWidget::clearAll(void)
  534. {
  535. mandelView->clearCells();
  536. mandelView->calcer.clearAll();
  537. }
  538. void MandelWidget::initializeGL(void)
  539. {
  540. this->context()->functions()->glClearColor(0, 0, 0, 0);
  541. this->context()->makeCurrent(nullptr);
  542. glDisable(GL_DEPTH_TEST);
  543. // looks not even better
  544. glEnable(GL_FRAMEBUFFER_SRGB);
  545. //glShadeModel(GL_SMOOTH);
  546. mandelView = nullptr;
  547. requestRecalc();
  548. }
  549. void MandelWidget::paintGL(void)
  550. {
  551. if (mandelView == nullptr) {
  552. mandelView = std::make_unique<MandelView>(generator, *this);
  553. QObject::connect(mandelView.get(), &MandelView::redrawRequested, this, static_cast<void(QOpenGLWidget::*)(void)>(&QOpenGLWidget::update));
  554. }
  555. int width = this->width();
  556. int height = this->height();
  557. mandelView->width = width;
  558. mandelView->height = height;
  559. glViewport(0, 0, width, height);
  560. glMatrixMode(GL_PROJECTION);
  561. glLoadIdentity();
  562. #ifdef QT_OPENGL_ES_1
  563. glOrthof(0, width, height, 0, -1.0, 1.0);
  564. #else
  565. glOrtho(0, width, height, 0, -1.0, 1.0);
  566. #endif
  567. glMatrixMode(GL_MODELVIEW);
  568. glClear(GL_COLOR_BUFFER_BIT);
  569. glLoadIdentity();
  570. updateAnimations();
  571. mandelView->paint(this->currentViewport);
  572. if (rubberbanding)
  573. drawRubberband();
  574. if (displayInfo)
  575. drawInfo();
  576. if (selectingPoint)
  577. drawPoint();
  578. }
  579. void MandelWidget::updateAnimations(void)
  580. {
  581. if (mnd::abs(currentViewport.width / targetViewport.width - 1.0) < 0.1e-5
  582. && mnd::abs(currentViewport.height / targetViewport.height - 1.0) < 0.1e-5) {
  583. // animation finished
  584. currentViewport = targetViewport;
  585. }
  586. else {
  587. auto now = std::chrono::high_resolution_clock::now();
  588. auto millis = std::chrono::duration_cast<std::chrono::milliseconds>(now - lastAnimUpdate).count();
  589. const mnd::Real factor = mnd::Real(::pow(0.97, millis));
  590. const mnd::Real one(1.0);
  591. currentViewport.x = currentViewport.x * factor + targetViewport.x * (one - factor);
  592. currentViewport.y = currentViewport.y * factor + targetViewport.y * (one - factor);
  593. currentViewport.width = currentViewport.width * factor + targetViewport.width * (one - factor);
  594. currentViewport.height = currentViewport.height * factor + targetViewport.height * (one - factor);
  595. lastAnimUpdate = now;
  596. emit update();
  597. }
  598. }
  599. void MandelWidget::drawRubberband(void)
  600. {
  601. glColor3ub(10, 200, 10);
  602. glBegin(GL_LINE_LOOP);
  603. glVertex2d(rubberband.x(), rubberband.y());
  604. glVertex2d(rubberband.right(), rubberband.y());
  605. glVertex2d(rubberband.right(), rubberband.bottom());
  606. glVertex2d(rubberband.x(), rubberband.bottom());
  607. glEnd();
  608. glEnable(GL_BLEND);
  609. glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
  610. glColor4f(0.1f, 0.9f, 0.1f, 0.2f);
  611. glBegin(GL_TRIANGLE_FAN);
  612. glVertex2d(rubberband.x(), rubberband.y());
  613. glVertex2d(rubberband.right(), rubberband.y());
  614. glVertex2d(rubberband.right(), rubberband.bottom());
  615. glVertex2d(rubberband.x(), rubberband.bottom());
  616. glEnd();
  617. glDisable(GL_BLEND);
  618. }
  619. void MandelWidget::drawInfo(void)
  620. {
  621. const float DIST_FROM_BORDER = 15;
  622. float maxWidth = this->width() - 2 * DIST_FROM_BORDER;
  623. mnd::Real distPerPixel = currentViewport.width / this->width();
  624. float log10 = (mnd::convert<float>(mnd::log(distPerPixel)) + ::logf(maxWidth)) / ::logf(10);
  625. mnd::Real displayDist = mnd::pow(mnd::Real(10), ::floor(log10));
  626. float pixels = mnd::convert<float>(displayDist / distPerPixel);
  627. int factor = 1;
  628. for (int i = 9; i > 1; i--) {
  629. if (pixels * i < maxWidth) {
  630. factor *= i;
  631. pixels *= i;
  632. displayDist *= i;
  633. break;
  634. }
  635. }
  636. std::stringstream dis;
  637. if (::abs(log10) < 3) {
  638. dis << mnd::convert<float>(displayDist);
  639. }
  640. else {
  641. dis << factor << "e" << int(::floor(log10));
  642. }
  643. if (maxWidth > 400) {
  644. dis << "; per pixel: " << distPerPixel;
  645. }
  646. float lineY = this->height() - DIST_FROM_BORDER;
  647. float lineXEnd = DIST_FROM_BORDER + pixels;
  648. QPainter infoPainter{ this };
  649. infoPainter.setPen(Qt::white);
  650. infoPainter.setFont(QFont("Arial", 12));
  651. infoPainter.drawLine(QPointF{ DIST_FROM_BORDER, lineY }, QPointF{ lineXEnd, lineY });
  652. infoPainter.drawLine(QPointF{ DIST_FROM_BORDER, lineY }, QPointF{ DIST_FROM_BORDER, lineY - 5 });
  653. infoPainter.drawLine(QPointF{ lineXEnd, lineY }, QPointF{ lineXEnd, lineY - 5 });
  654. infoPainter.drawText(int(DIST_FROM_BORDER), int(lineY - 20), int(lineXEnd - DIST_FROM_BORDER), 20,
  655. Qt::AlignCenter, QString::fromStdString(dis.str()));
  656. infoPainter.end();
  657. }
  658. void MandelWidget::drawPoint(void)
  659. {
  660. glColor3ub(255, 255, 255);
  661. glBegin(GL_LINES);
  662. glVertex2f(0, pointY);
  663. glVertex2f(width(), pointY);
  664. glVertex2f(pointX, 0);
  665. glVertex2f(pointX, height());
  666. glEnd();
  667. }
  668. void MandelWidget::zoom(float scale, float x, float y)
  669. {
  670. targetViewport.zoom(scale, x, y);
  671. lastAnimUpdate = std::chrono::high_resolution_clock::now();
  672. //currentViewport.zoom(scale, x, y);
  673. requestRecalc();
  674. }
  675. void MandelWidget::setViewport(const mnd::MandelViewport& viewport)
  676. {
  677. targetViewport = viewport;
  678. targetViewport.adjustAspectRatio(this->width(), this->height());
  679. currentViewport = targetViewport;
  680. //lastAnimUpdate = std::chrono::high_resolution_clock::now();
  681. //currentViewport.zoom(scale, x, y);
  682. requestRecalc();
  683. }
  684. void MandelWidget::selectPoint(void)
  685. {
  686. this->selectingPoint = true;
  687. this->setMouseTracking(true);
  688. }
  689. void MandelWidget::stopSelectingPoint(void)
  690. {
  691. this->selectingPoint = false;
  692. this->setMouseTracking(false);
  693. }
  694. void MandelWidget::requestRecalc()
  695. {
  696. emit update();
  697. }
  698. void MandelWidget::resizeEvent(QResizeEvent* re)
  699. {
  700. QOpenGLWidget::resizeEvent(re);
  701. double aspect = double(geometry().width()) / geometry().height();
  702. currentViewport.height = currentViewport.width / aspect;
  703. targetViewport = currentViewport;
  704. if (mandelView.get() != nullptr) {
  705. mandelView->width = this->width();
  706. mandelView->height = this->height();
  707. }
  708. requestRecalc();
  709. }
  710. void MandelWidget::mousePressEvent(QMouseEvent* me)
  711. {
  712. QOpenGLWidget::mousePressEvent(me);
  713. if (me->button() == Qt::RightButton) {
  714. rubberbanding = true;
  715. rubberband.setCoords(me->x(), me->y(), me->x(), me->y());
  716. emit repaint();
  717. me->accept();
  718. }
  719. else if (me->button() == Qt::LeftButton) {
  720. dragging = true;
  721. dragX = me->x();
  722. dragY = me->y();
  723. me->accept();
  724. }
  725. }
  726. void MandelWidget::mouseMoveEvent(QMouseEvent* me)
  727. {
  728. QOpenGLWidget::mouseMoveEvent(me);
  729. if (rubberbanding) {
  730. QRectF& rect = rubberband;
  731. double aspect = double(geometry().width()) / geometry().height();
  732. rect.setBottomRight(QPoint(me->x(), me->y()));
  733. if (rect.width() > rect.height() * aspect)
  734. rect.setHeight(rect.width() / aspect);
  735. else
  736. rect.setWidth(rect.height() * aspect);
  737. emit repaint();
  738. }
  739. else if (selectingPoint) {
  740. pointX = me->x();
  741. pointY = me->y();
  742. emit repaint();
  743. }
  744. else if (dragging) {
  745. double deltaX = me->x() - dragX;
  746. double deltaY = me->y() - dragY;
  747. this->currentViewport.x -= deltaX * currentViewport.width / this->width();
  748. this->currentViewport.y -= deltaY * currentViewport.height / this->height();
  749. targetViewport = currentViewport;
  750. dragX = me->x(); dragY = me->y();
  751. emit repaint();
  752. }
  753. me->accept();
  754. }
  755. void MandelWidget::mouseReleaseEvent(QMouseEvent* me)
  756. {
  757. QOpenGLWidget::mouseReleaseEvent(me);
  758. if (rubberbanding) {
  759. QRect rect = rubberband.toRect();
  760. if(rect.width() != 0 && rect.height() != 0) {
  761. QRect full = this->geometry();
  762. targetViewport.x += mnd::Real(rect.left()) * targetViewport.width / full.width();
  763. targetViewport.y += mnd::Real(rect.top()) * targetViewport.height / full.height();
  764. targetViewport.width *= mnd::Real(rect.width()) / full.width();
  765. targetViewport.height *= mnd::Real(rect.height()) / full.height();
  766. targetViewport.normalize();
  767. currentViewport = targetViewport;
  768. }
  769. requestRecalc();
  770. rubberbanding = false;
  771. }
  772. else if (selectingPoint) {
  773. selectingPoint = false;
  774. this->setMouseTracking(false);
  775. mnd::Real x = currentViewport.x + currentViewport.width * mnd::convert<mnd::Real>(float(me->x()) / width());
  776. mnd::Real y = currentViewport.y + currentViewport.height * mnd::convert<mnd::Real>(float(me->y()) / height());
  777. emit pointSelected(x, y);
  778. emit repaint();
  779. }
  780. dragging = false;
  781. //requestRecalc();
  782. }
  783. void MandelWidget::wheelEvent(QWheelEvent* we)
  784. {
  785. QOpenGLWidget::wheelEvent(we);
  786. float x = float(we->x()) / this->width();
  787. float y = float(we->y()) / this->height();
  788. float scale = ::powf(0.9975f, we->angleDelta().y());
  789. zoom(scale, x, y);
  790. if (!we->pixelDelta().isNull())
  791. this->currentViewport = this->targetViewport;
  792. we->accept();
  793. }
  794. /*void MandelWidget::viewUpdated(Bitmap<RGBColor>* bitmap)
  795. {
  796. if (bitmap != nullptr) {
  797. delete bitmap;
  798. emit repaint();
  799. }
  800. }*/