1c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This file is part of Eigen, a lightweight C++ template library 2c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// for linear algebra. 3c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 4c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Copyright (C) 2008 Benoit Jacob <jacob.benoit.1@gmail.com> 5c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 6c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This Source Code Form is subject to the terms of the Mozilla 7c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Public License v. 2.0. If a copy of the MPL was not distributed 8c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// with this file, You can obtain one at http://mozilla.org/MPL/2.0/. 9c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 10c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include "mandelbrot.h" 11c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <iostream> 12c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include<QtGui/QPainter> 13c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include<QtGui/QImage> 14c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include<QtGui/QMouseEvent> 15c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include<QtCore/QTime> 16c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 17c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathvoid MandelbrotWidget::resizeEvent(QResizeEvent *) 18c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 19c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if(size < width() * height()) 20c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 21c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath std::cout << "reallocate buffer" << std::endl; 22c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath size = width() * height(); 23c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if(buffer) delete[]buffer; 24c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath buffer = new unsigned char[4*size]; 25c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 26c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 27c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 28c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<typename T> struct iters_before_test { enum { ret = 8 }; }; 29c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<> struct iters_before_test<double> { enum { ret = 16 }; }; 30c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 31c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<typename Real> void MandelbrotThread::render(int img_width, int img_height) 32c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 33c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath enum { packetSize = Eigen::internal::packet_traits<Real>::size }; // number of reals in a Packet 34c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef Eigen::Array<Real, packetSize, 1> Packet; // wrap a Packet as a vector 35c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 36c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath enum { iters_before_test = iters_before_test<Real>::ret }; 37c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath max_iter = (max_iter / iters_before_test) * iters_before_test; 38c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath const int alignedWidth = (img_width/packetSize)*packetSize; 39c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath unsigned char *const buffer = widget->buffer; 40c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath const double xradius = widget->xradius; 41c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath const double yradius = xradius * img_height / img_width; 42c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath const int threadcount = widget->threadcount; 43c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef Eigen::Array<Real, 2, 1> Vector2; 44c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Vector2 start(widget->center.x() - widget->xradius, widget->center.y() - yradius); 45c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Vector2 step(2*widget->xradius/img_width, 2*yradius/img_height); 46c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath total_iter = 0; 47c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 48c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int y = id; y < img_height; y += threadcount) 49c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 50c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath int pix = y * img_width; 51c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 52c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // for each pixel, we're going to do the iteration z := z^2 + c where z and c are complex numbers, 53c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // starting with z = c = complex coord of the pixel. pzi and pzr denote the real and imaginary parts of z. 54c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // pci and pcr denote the real and imaginary parts of c. 55c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 56c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Packet pzi_start, pci_start; 57c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int i = 0; i < packetSize; i++) pzi_start[i] = pci_start[i] = start.y() + y * step.y(); 58c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 59c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int x = 0; x < alignedWidth; x += packetSize, pix += packetSize) 60c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 61c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Packet pcr, pci = pci_start, pzr, pzi = pzi_start, pzr_buf; 62c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int i = 0; i < packetSize; i++) pzr[i] = pcr[i] = start.x() + (x+i) * step.x(); 63c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 64c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // do the iterations. Every iters_before_test iterations we check for divergence, 65c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // in which case we can stop iterating. 66c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath int j = 0; 67c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef Eigen::Matrix<int, packetSize, 1> Packeti; 68c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Packeti pix_iter = Packeti::Zero(), // number of iteration per pixel in the packet 69c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath pix_dont_diverge; // whether or not each pixel has already diverged 70c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath do 71c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 72c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int i = 0; i < iters_before_test/4; i++) // peel the inner loop by 4 73c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 74c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath# define ITERATE \ 75c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath pzr_buf = pzr; \ 76c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath pzr = pzr.square(); \ 77c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath pzr -= pzi.square(); \ 78c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath pzr += pcr; \ 79c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath pzi = (2*pzr_buf)*pzi; \ 80c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath pzi += pci; 81c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ITERATE ITERATE ITERATE ITERATE 82c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 83c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath pix_dont_diverge = ((pzr.square() + pzi.square()) 84c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath .eval() // temporary fix as what follows is not yet vectorized by Eigen 85c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath <= Packet::Constant(4)) 86c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // the 4 here is not a magic value, it's a math fact that if 87c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // the square modulus is >4 then divergence is inevitable. 88c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath .template cast<int>(); 89c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath pix_iter += iters_before_test * pix_dont_diverge; 90c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath j++; 91c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath total_iter += iters_before_test * packetSize; 92c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 93c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath while(j < max_iter/iters_before_test && pix_dont_diverge.any()); // any() is not yet vectorized by Eigen 94c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 95c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // compute pixel colors 96c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int i = 0; i < packetSize; i++) 97c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 98c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath buffer[4*(pix+i)] = 255*pix_iter[i]/max_iter; 99c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath buffer[4*(pix+i)+1] = 0; 100c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath buffer[4*(pix+i)+2] = 0; 101c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 102c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 103c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 104c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath // if the width is not a multiple of packetSize, fill the remainder in black 105c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int x = alignedWidth; x < img_width; x++, pix++) 106c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath buffer[4*pix] = buffer[4*pix+1] = buffer[4*pix+2] = 0; 107c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 108c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return; 109c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 110c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 111c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathvoid MandelbrotThread::run() 112c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 113c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath setTerminationEnabled(true); 114c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath double resolution = widget->xradius*2/widget->width(); 115c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath max_iter = 128; 116c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if(resolution < 1e-4f) max_iter += 128 * ( - 4 - std::log10(resolution)); 117c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath int img_width = widget->width()/widget->draft; 118c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath int img_height = widget->height()/widget->draft; 119c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath single_precision = resolution > 1e-7f; 120c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 121c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if(single_precision) 122c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath render<float>(img_width, img_height); 123c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath else 124c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath render<double>(img_width, img_height); 125c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 126c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 127c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathvoid MandelbrotWidget::paintEvent(QPaintEvent *) 128c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 129c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static float max_speed = 0; 130c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath long long total_iter = 0; 131c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 132c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath QTime time; 133c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath time.start(); 134c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int th = 0; th < threadcount; th++) 135c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath threads[th]->start(QThread::LowPriority); 136c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int th = 0; th < threadcount; th++) 137c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 138c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath threads[th]->wait(); 139c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath total_iter += threads[th]->total_iter; 140c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 141c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath int elapsed = time.elapsed(); 142c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 143c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if(draft == 1) 144c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 145c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath float speed = elapsed ? float(total_iter)*1000/elapsed : 0; 146c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath max_speed = std::max(max_speed, speed); 147c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath std::cout << threadcount << " threads, " 148c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath << elapsed << " ms, " 149c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath << speed << " iters/s (max " << max_speed << ")" << std::endl; 150c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath int packetSize = threads[0]->single_precision 151c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath ? int(Eigen::internal::packet_traits<float>::size) 152c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath : int(Eigen::internal::packet_traits<double>::size); 153c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath setWindowTitle(QString("resolution ")+QString::number(xradius*2/width(), 'e', 2) 154c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath +QString(", %1 iterations per pixel, ").arg(threads[0]->max_iter) 155c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath +(threads[0]->single_precision ? QString("single ") : QString("double ")) 156c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath +QString("precision, ") 157c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath +(packetSize==1 ? QString("no vectorization") 158c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath : QString("vectorized (%1 per packet)").arg(packetSize))); 159c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 160c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 161c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath QImage image(buffer, width()/draft, height()/draft, QImage::Format_RGB32); 162c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath QPainter painter(this); 163c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath painter.drawImage(QPoint(0, 0), image.scaled(width(), height())); 164c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 165c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if(draft>1) 166c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 167c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath draft /= 2; 168c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath setWindowTitle(QString("recomputing at 1/%1 resolution...").arg(draft)); 169c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath update(); 170c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 171c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 172c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 173c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathvoid MandelbrotWidget::mousePressEvent(QMouseEvent *event) 174c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 175c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if( event->buttons() & Qt::LeftButton ) 176c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 177c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath lastpos = event->pos(); 178c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath double yradius = xradius * height() / width(); 179c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath center = Eigen::Vector2d(center.x() + (event->pos().x() - width()/2) * xradius * 2 / width(), 180c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath center.y() + (event->pos().y() - height()/2) * yradius * 2 / height()); 181c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath draft = 16; 182c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int th = 0; th < threadcount; th++) 183c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath threads[th]->terminate(); 184c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath update(); 185c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 186c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 187c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 188c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathvoid MandelbrotWidget::mouseMoveEvent(QMouseEvent *event) 189c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 190c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath QPoint delta = event->pos() - lastpos; 191c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath lastpos = event->pos(); 192c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if( event->buttons() & Qt::LeftButton ) 193c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 194c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath double t = 1 + 5 * double(delta.y()) / height(); 195c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if(t < 0.5) t = 0.5; 196c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if(t > 2) t = 2; 197c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath xradius *= t; 198c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath draft = 16; 199c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int th = 0; th < threadcount; th++) 200c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath threads[th]->terminate(); 201c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath update(); 202c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 203c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 204c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 205c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathint main(int argc, char *argv[]) 206c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath{ 207c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath QApplication app(argc, argv); 208c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath MandelbrotWidget w; 209c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath w.show(); 210c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return app.exec(); 211c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 212c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 213c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include "mandelbrot.moc" 214