highperfocused/ffmpeg
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

vsrc_mandelbrot: simplify second bailout search by not recalculating values.

Browse Source 
Signed-off-by: Michael Niedermayer <michaelni@gmx.at>
This commit is contained in:
Michael Niedermayer 2011-11-19 15:00:26 +01:00
parent e749b3780c
commit da1344e0cd
1 changed files with 12 additions and 15 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
libavfilter/vsrc_mandelbrot.c
View File

@ -264,25 +264,22 @@ static void draw_mandelbrot(AVFilterContext *ctx, uint32_t *color, int linesize,
Z_Z2_C_ZYKLUS(t, zi, zr, zi, 0) Z_Z2_C_ZYKLUS(t, zi, zr, zi, 0)
i++; i++;
Z_Z2_C_ZYKLUS(zr, zi, t, zi, 1) Z_Z2_C_ZYKLUS(zr, zi, t, zi, 1)
if(zr*zr + zi*zi > mb->bailout) if(zr*zr + zi*zi > mb->bailout){
break;
}
i-= FFMIN(7, i); i-= FFMIN(7, i);
for(; i<mb->maxiter; i++){
zr= mb->zyklus[i][0]; zr= mb->zyklus[i][0];
zi= mb->zyklus[i][1]; zi= mb->zyklus[i][1];
for(; i<mb->maxiter; i++){
double t;
if(zr*zr + zi*zi > mb->bailout){ if(zr*zr + zi*zi > mb->bailout){
switch(mb->outer){ switch(mb->outer){
case ITERATION_COUNT: zr = i - (SQR(zr-cr)+SQR(zi-ci) > SQR(mb->bailout)); break; case ITERATION_COUNT: zr = i; break;
case NORMALIZED_ITERATION_COUNT: zr= i + log2(log(mb->bailout) / log(zr*zr + zi*zi)); break; case NORMALIZED_ITERATION_COUNT: zr= i + log2(log(mb->bailout) / log(zr*zr + zi*zi)); break;
} }
c= lrintf((sin(zr)+1)*127) + lrintf((sin(zr/1.234)+1)*127)*256*256 + lrintf((sin(zr/100)+1)*127)*256; c= lrintf((sin(zr)+1)*127) + lrintf((sin(zr/1.234)+1)*127)*256*256 + lrintf((sin(zr/100)+1)*127)*256;
break; break;
} }
Z_Z2_C_ZYKLUS(t, zi, zr, zi, 0) }
i++; break;
Z_Z2_C_ZYKLUS(zr, zi, t, zi, 1) }
} }
if(!c){ if(!c){
if(mb->inner==PERIOD){ if(mb->inner==PERIOD){