Brzina Pythona 3.5

U prosjeku Python je 5x sporiji nego C++ no moze biti cak do 15-20x, no to ne znaci da je los jezik posto u vrlo malo algoritama je potrebna velika brzina jer nisu toliko kompleksni. Meni crta Mandelbrota za 3.7s, sta je sasvim solidno. Ako ti je potrebna veca brzina onda predji na C++, i mozes onda bindat kod u Pythonu sa onim u C++.

Normalno da ti necu rijesiti, medutim, ako mozes staviti kod, stavi ga tako da te uputimo sto da promijenis.

Jedno pitanje, s cim crtas (tkinter, nesto drugo?) i crtas li piksel po piksel?

 Mislim da imas i Tkinter obradjen u 2. knjizi ali nisam 100% siguran

#Mandelbrotov set crtanje
from turtle import*
print('crtanje je kvadratno ')
#pocetak_x=float(input('Unesi početak intervala za realni dio '))
#kraj_x=float(input('Unesi kraj intervala za realni dio '))
#pocetak_y=float(input('Unesi početak intervala za imaginarni dio '))
#kraj_y=float(input('Unesi kraj intervala za imaginarni dio '))

dimenzija=float(input('Unesi dimenziju kvadratnog crtanja seta '))
gustoca=int(input('Unesi gustocu crtanja '))
provjera=int(input('Unesi broj provjera za svaki broj '))
def funkcija(x,y):
realni=x
imaginarni=y
c_realni=x
c_imaginarni=y
broj_provjera=1
rezultat=(realni**2 + imaginarni**2)**(1/2)
da_il_ne=0
if rezultat<2:
while broj_provjera<provjera and rezultat<2:
privremena=realni
privremena2=imaginarni
realni=privremena**2 - privremena2**2 + c_realni
imaginarni= 2*privremena*privremena2 + c_imaginarni
rezultat = ((realni)**2 + (imaginarni)**2)**(1/2)
broj_provjera= broj_provjera+1
else:
da_il_ne=0
if rezultat<2:
da_il_ne=1
else:
da_il_ne=0
c=[broj_provjera, da_il_ne]
return(c)
#neke provjere

#print('provjera za 0,i ',funkcija(0,1)[1])
#print('provjera za 0,0.3i ',funkcija(0,0.3)[1])

#crtanje, pomice crtez od sredista ekrana
#ht()
lt(90)
pu()
fd (100)
lt(90)
fd(100)
rt(180)
pd()

#razmak_x = (kraj_x-pocetak_x)/gustoca
#razmak_y = (kraj_y-pocetak_y)/gustoca
#neki_y=pocetak_y
#gustoca_y=(kraj_y-pocetak_y)*200
#gustoca_x=2

razmak_x = 2*dimenzija/gustoca
razmak_y= 2*dimenzija/gustoca
neki_x=-dimenzija
neki_y=dimenzija

for i in range (2*gustoca):
neki_x = -dimenzija
for j in range(2*gustoca):
if funkcija(neki_x,neki_y)[1]==1:
dot(2)
neki_x=neki_x+razmak_x
pu()
fd(1)
pd()
neki_y=neki_y-razmak_y
pu()
back(gustoca*2)
rt(90)
fd(1)
lt(90)
pd()

Naravno da zna. Je da je spor, al je zato svicarski nozic. Najbolji je.

Nego, daj formatiraj taj kod malo, imas na alatnoj traci format "Kod" pa onda napisi kod ovako kako sam ja u prijasnjem postu. Pogotovo jer je kod Pythona uvlacenje najbitnije.

EDIT: Ili jos bolje, stavi ga na pastebin.com pa ovdje samo link kopiraj.

Evo jedan savjet.

x * x

je jedno 10-ak puta brze od:

x ** 2

Tu mozes dobiti nesto brzine jer vidim da puno kvadriranja radis. Mislim cak da je za korijen math.sqrt() (koji se linka na C sqrt() funkciju) nesto brzi od **(1/2). Probaj.

Jos jedna stvar koji bi mogao napraviti da ubrzas kod je koristiti JIT kompajler, poput Numbe. Kratak primjer kako to izgleda u Pythonu:

from numba import jit

@jit
def fja(x):
    # radi neku tesku matematiku s x
    return rez

Ovo podebljano @jit je vazno, taj dekorator stavis iznad svake funkcije koju zelis da Numba kompajlira u nativne strojne instrukcije. To bi trebalo drasticno ubrzati stvari. Uz ono sto sam ti rekao vec prije sto se tice iscrtavanja, nikako ne piksel po piksel jer to trosi puno resursa. Bolje je unaprijed izracunati veci komad (liniju, stupac ili manju matricu) i onda to prikazati.

Takoder, opcenito savjet, pisi prozracniji kod. Ovo je sve natiskano i tesko citljivo. Stavljaj razmake oko znaka '=', osim u pozivu funkcije, dakle:

realni = x
# a ne
realni=x

Medutim, u pozivu funkcije:

fja(neki_parametar=x)
# a ne
fja(neki_parametar = x)

Za vise o tome kako pisati Python kod pogledaj PEP 8.

Neke druge savjete u principu vise nemam, nadam se da ce ti ovo pomoci.

To ti je paralelizacija algoritama i to je prilicno zeznuta stvar. Mislim da ti je to ipak prevelik zalogaj zasad (kazes da si jos srednja skola, znaci tek pocinjes).

No ako bas zelis i inzistiras: za paralelizam pogledaj Pythonovu dokumentaciju o dretvama, procesima i meduprocesnom komunikacijom (linkovi su za Python 3, ako koristis Python 2, lijevo gore ga odaberi). Potrazi na netu i MPI (Message Passing Interface), Python paket je mpi4py. Potrazi neke radove na temu paralelnog Mandelbrota, poput ovoga. Ako nista od ovoga ne kuzis, apsolutno je u redu, nemoj misliti da si glup jer stvarno nije trivijalno.