Tri algoritmy pro stredovani
============================
system setup: 
	BH (M_BH=3.5e+6 Msun), c=0.2
	+ Sfericka porucha (M_sphe=1.0e-3 M_BH, r0=2.0*a(S2), alpha=1.4)
	+ Ring (M_ring=5.0e-2 M_BH, a_ring=10.0*a(S2), z_ring=0.0)
IVD: 
	e: 0..0.975, 45 hodnot
	omega: 0°..359°, 360 hodnot
	CELKEM: 14400 HODNOT

Testovaci PC: CPU 3.066 GHz, 1 GB RAM

1] Pevne nastaveny maximalni krok hmax=Perioda/Nsteps, Nsteps=5000, vystredovat po 2*Nsteps krocich; ukladani hodnot
do poli a nasledne hromadne stredovani
Doba vypoctu dat pro isokontury: 931.66user 8.49system 15:43.07elapsed 99%CPU 
ex: ring5e-2_sphere1e-3
bug: pro ruzne e potrebuji ruzne pocty kroku => pocitam zbytecne dlouho
NAME: TYP_1
ID: r5e-2

2] Predpocitavani poctu kroku pro 1 neporusenou periodu (poustim z apocentra a hledam apocentrum porovnavanim 3
nasledujicich r), vystredovat po tomto poctu kroku; ukladani hodnot do poli a nasledne hromadne stredovani
Doba vypoctu dat pro isokontury: 768.55user 2.89system 12:59.36elapsed 98%CPU 
ex: integrace_02_ring5e-2_sphere1e-3
bug: pro e=0 se spatne urci apocentrum => spatne se vypocte pocet kroku => spatne se vystreduje potencial
NAME: TYP_2
ID: ir5e-2

3] Zpatky na stromy aneb pocitej (pricitej) hodnoty potencialu, dokud nepocitas dele nez 1 neporusenou periodu; scita
se okamzite, zadne ukladani do poli, streduje se nakonec
Doba vypoctu dat pro isokontury: 253.91user 0.74system 4:16.73elapsed 99%CPU 
ex: iso_integrace_konec_po_periode
bug: zbytecny vypocet eliptickych integralu v derivs, odstraneno v TYP_3a
NAME: TYP_3
ID: konec

3a] TYP_3 po odstraneni zbytecnych vypoctu eliptickych integralu v dervis
Doba vypoctu dat pro isokontury: 129.79user 0.48system 2:10.33elapsed 99%CPU 
ex: iso_integrace_konec_po_periode_bez_ElInt
bug: zatim nezjisten
NAME: TYP_3a
ID: konec_a

Shrnuti:
Nejlepsi je algoritmus TYP_3*, nebot pri pohledu na obrazek timing.png bystry ctenar zjisti, ze casy stredovani tvori
temer delta-funkci, kdezto v pripade algoritmu TYP_1 jsou casy rozlozeny v intervalu o delce zhruba 1/4 neporusene
periody. TYP_2 dava celkem slusne vysledky pro nenulove ekcentricity. Vzhledem k rychlosti je nejlepsi TYP_3a, protoze
je nejlepsi ve stredovani a 2x rychlejsi nez TYP_3.