You cannot select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
101 lines
2.6 KiB
C
101 lines
2.6 KiB
C
#include <stdint.h>
|
|
#include <math.h>
|
|
#include <stdbool.h>
|
|
#include <stdio.h>
|
|
#include <stdlib.h>
|
|
|
|
#define PREFIX "img-"
|
|
|
|
#define VENTANA_ANCHO 1920
|
|
#define VENTANA_ALTO 1080
|
|
#define ANCHO VENTANA_ANCHO
|
|
#define ALTO VENTANA_ALTO
|
|
#define MAX_ITER 500
|
|
|
|
//long double xc = -1.740119020442992652289958266376L;
|
|
//long double yc = 0.028019209975637676263081582255L;
|
|
long double xc = 0.264226442163892496529931097626331393L, yc = 0.002572261418806922467217290773078275;
|
|
|
|
uint32_t computar_intensidad(long double x0, long double y0, long double zoom) {
|
|
//x0 = x0 * zoom - 1.7400623825793398724575; // Derecha: Achicar
|
|
//y0 = y0 * zoom + 0.0281753397792111; // Abajo: Achicar
|
|
x0 = x0 * zoom + xc; // Derecha: Achicar
|
|
y0 = y0 * zoom + yc; // Abajo: Achicar
|
|
//x0 = x0 * zoom - 1.78; // Derecha: Achicar
|
|
//y0 = y0 * zoom + 0.0; // Abajo: Achicar
|
|
long double x = 0;
|
|
long double y = 0;
|
|
int n = 0;
|
|
while(x*x + y*y <= 4 && n < MAX_ITER) {
|
|
long double xtemp = x*x - y*y + x0;
|
|
y = 2*x*y + y0;
|
|
x = xtemp;
|
|
n++;
|
|
}
|
|
|
|
if(n == MAX_ITER)
|
|
return 0;
|
|
|
|
//float h = (n + 1 - log(log2(sqrt(x*x+y*y)))) / MAX_ITER * 360;
|
|
//float h = n * 360.0 / MAX_ITER;
|
|
float h = (int)(n * 360.0 / 50) % 360;
|
|
float xxx = (1 - fabs(fmodf(h / 60.0, 2) - 1));
|
|
uint8_t xx = (xxx < 1 ? xxx : 1) * 255;
|
|
uint8_t r = 0, g = 0, b = 0;
|
|
|
|
if(h < 60)
|
|
r = 255, g = xx;
|
|
else if(h < 120)
|
|
g = 255, r = xx;
|
|
else if(h < 180)
|
|
g = 255, b = xx;
|
|
else if(h < 240)
|
|
b = 255, g = xx;
|
|
else if(h < 300)
|
|
b = 255, r = xx;
|
|
else
|
|
r = 255, b = xx;
|
|
|
|
return (r << 24) | (g << 16) | (b << 8);
|
|
}
|
|
|
|
|
|
int main() {
|
|
long double zoom = 0.179;
|
|
double factor = 0.99;
|
|
|
|
int i = 0;
|
|
|
|
while(1) {
|
|
if(zoom < 1e-20)
|
|
break;
|
|
|
|
zoom *= factor;
|
|
|
|
printf("%d\n", i);
|
|
|
|
char aux[100];
|
|
sprintf(aux, "%s%05d.ppm", PREFIX, i);
|
|
FILE *f = fopen(aux, "wt");
|
|
fprintf(f, "P3\n%d %d\n255\n", ANCHO, ALTO);
|
|
|
|
|
|
for(int vy = ALTO / 2; vy > - ALTO / 2; vy--)
|
|
for(int vx = - ANCHO / 2; vx < ANCHO / 2; vx++) {
|
|
uint32_t c = computar_intensidad(vx, vy, zoom);
|
|
fprintf(f, "%d %d %d\n", c >> 24, (c >> 16) & 0xFF, (c >> 8) & 0xFF);
|
|
}
|
|
fclose(f);
|
|
|
|
sprintf(aux, "convert %s%05d.ppm %s%05d.png; rm %s%05d.ppm", PREFIX, i, PREFIX, i, PREFIX, i);
|
|
system(aux);
|
|
|
|
i++;
|
|
}
|
|
|
|
system("ffmpeg -y -framerate 24 -i " PREFIX "%05d.png -c:v libx264 -crf 0 output.mp4");
|
|
|
|
return 0;
|
|
}
|
|
|