void wind_file_set_current_path(char *path)
{
w_int32_t len = wind_strlen(path);
if(curpath != W_NULL)
wind_free(curpath);
curpath = wind_malloc(len +1);
wind_memcpy(curpath,path,len+1);
}
static void pack_info_init(void)
{
w_int32_t i;
pack_info_s *pkinfo = &pack_info;
wind_memset(databuff,0,sizeof(databuff));
wind_memset(&imghead,0,sizeof(imghead));
if(pkinfo->file_cnt > 0)
{
for(i = 0;i < pkinfo->file_cnt;i ++)
wind_free(pkinfo->fileinfo[i].buff);
}
wind_memset(&pack_info,0,sizeof(pack_info));
}
static void release_file_buff(void)
{
w_int32_t i;
pack_info_s *pkinfo = &pack_info;
for(i = 0;i < pkinfo->file_cnt;i ++)
{
if(pkinfo->fileinfo[i].buff)
{
wind_free(pkinfo->fileinfo[i].buff);
pkinfo->fileinfo[i].buff = W_NULL;
}
}
}
w_err_t wind_full_path_release(char *path)
{
return wind_free(path);
}
static w_err_t pack_files(pack_info_s *info)
{
w_int32_t i;
w_uint8_t *buff;
char *img;
w_err_t err;
img_head_s *head;
w_encypt_ctx_s ctx;
w_uint32_t crc;
w_uint8_t key[] = ENCRYPT_KEY;
pack_info_s *pkinfo = &pack_info;
infile_info_s *finfo;
w_int32_t imglen = calc_img_lenth();
buff = wind_malloc(imglen);
WIND_ASSERT_RETURN(buff != W_NULL,W_ERR_MEM);
wind_memset(buff, 0, imglen);
WIND_ASSERT_RETURN(info->fileinfo[0].offset >= IMG_HEAD_LEN,W_ERR_FAIL);
for(i = 0;i < pkinfo->file_cnt;i ++)
{
finfo = &pkinfo->fileinfo[i];
wind_memcpy(&buff[finfo->offset],finfo->buff,finfo->flen);
}
if(ENCRYPT_TYPE)
{
wind_notice("encrypt bin file");
pkinfo->encrypt_type = ENCRYPT_TYPE;
wind_memcpy(pkinfo->keys,key,sizeof(key));
pkinfo->key_len = sizeof(key);
wind_encrypt_init(&ctx,pkinfo->keys,pkinfo->key_len);
wind_encrypt(&ctx,&buff[IMG_HEAD_LEN],imglen - IMG_HEAD_LEN);
}
wind_notice("calc bin file crc");
crc = wind_crc32(&buff[IMG_HEAD_LEN],imglen - IMG_HEAD_LEN,0xffffffff);
wind_notice("generate img file head info");
head = &imghead;
head->magic = IMG_MAGIC;
head->img_len = imglen;
head->head_len = IMG_HEAD_LEN;
head->head_ver = IMG_HEAD_VER;
head->hard_ver = pkinfo->hw_version;
head->soft_ver = pkinfo->sw_version;
head->bin_crc = crc;
head->bin_offset = IMG_HEAD_LEN;
head->encrypt_type = pkinfo->encrypt_type;
img = wind_strrchr(pkinfo->output_file, '\\');
if(img == W_NULL)
img = wind_strrchr(pkinfo->output_file, '/');
if(img == W_NULL)
img = pkinfo->output_file;
else
img ++;
wind_strncpy(head->img_name, img, 64);
wind_strncpy(head->board_name, pkinfo->board_name, 32);
wind_strncpy(head->cpu_name, pkinfo->cpu_name, 32);
wind_strncpy(head->arch_name, pkinfo->arch_name, 32);
boot_img_head_set(head,buff);
wind_notice("flush image file");
err = write_file(pkinfo->output_file,0, buff,imglen);
wind_free(buff);
return err;
}
void gravity(image_processor p, unsigned imageWidth, unsigned imageHeight, unsigned frames, unsigned nsSize, unsigned nsDepth) {
image* img = image_new(imageWidth, imageWidth);
fill_image(img, 0, 0, 0);
//Parameterization:
unsigned wind_count = 8; //Number of particle emmiters.
unsigned wind_size = 1000; //Max number of particles per emitter.
unsigned particlesPerFrame = 2000 / frames;
if(particlesPerFrame == 0) particlesPerFrame = 1;
unsigned particlesToDraw = 2;
float maxSpread = imageWidth / 500.0;
float particleNoiseVelocity = 1024.0 * imageWidth / (1000.0);
float particleEmitterMomentumScalar = 0.5;
float totalTime = TAU;
float wColorRevert = 0.01; //TODO these should be exponentiated by tStep
float wColorChange = 0.125;
float wAlpha = 1.0 / 8.0; //Alpha to use for drawing the wind.
unsigned tracer_count = 4; //Number of traced curves.
float tracerColorChange = 0.125;
float tracerColorReduction = 0.9;
////////////////////
//Particle Emitters:
wind* winds[wind_count];
for(unsigned i = 0; i < wind_count; i++) {
winds[i] = wind_new(wind_size);
//wind_randomize(winds[i], imageWidth, 0, 0, imageWidth, imageWidth, VMAX, 100);
}
float windDefaultColors[wind_count * 3];
float windColors[wind_count * 3];
for(unsigned i = 0; i < wind_count; i++) {
windDefaultColors[i * 3 + 0] = (i & 0x1) ? 3 : 0;
windDefaultColors[i * 3 + 1] = (i & 0x2) ? 3 : 0;
windDefaultColors[i * 3 + 2] = (i & 0x4) ? 3 : 0;
};
memcpy(windColors, windDefaultColors, sizeof(float) * wind_count * 3);
//////////////
//Wind Curves:
//Wind curve closures (Specific type storage).
ccl_1 pc1_cl[wind_count / 2];
ccl_2 pc2_cl[wind_count / 2];
//Parametric curve functions (Indirection layer to support generic curve types).
vec2(*parametricCurves[wind_count])(float, void*);
//Closures for (generic) parametric curve functions
void* cl[wind_count];
//Initialize wind closures.
for(unsigned i = 0; i < wind_count / 2; i++) {
randomize_ccl_1(pc1_cl + i, imageWidth / 2, imageWidth / 2, imageWidth / 4);
//Keep the amplitude low and the base high.
pc1_cl[i].base = (pc1_cl[i].base + imageWidth / 2) / 2;
pc1_cl[i].amplitude *= 0.5;
pc1_cl[i].frequency = 2 + pc1_cl[i].frequency / 4;
cl[i] = (void*) (pc1_cl + i);
parametricCurves[i] = parametric_curve_1;
randomize_ccl_2(pc2_cl + i, imageWidth / 2, imageWidth / 2, imageWidth / 2);
for(unsigned j = 0; j < 4; j++) {
pc2_cl[i].d[j] *= 1.0 / 32.0;
}
cl[i + wind_count / 2] = (void*) (pc2_cl + i);
parametricCurves[i + wind_count / 2] = parametric_curve_2;
}
/////////////////////
//Potential Function:
noise_sum* ncl = initialize_noise_sum_2d(nsSize, nsDepth);
noise_sum_scale_in(ncl, 1.0 / imageWidth);
centered_cl orbitcl = { imageWidth / 2.0, imageWidth / 2.0, imageWidth * imageWidth, (imageWidth / 8.0) * (imageWidth / 8.0) };
float(*potentialFunctions[POTENTIAL_COUNT])(float, float, void*) = {
noiseSumPotential,
distanceSquaredPotential
};
void* pcls[POTENTIAL_COUNT] = {ncl, &orbitcl};
float weights[POTENTIAL_COUNT] = {
80000000, 100000
};
poly_weighted_cl pcl = { potentialFunctions, pcls, weights, POTENTIAL_COUNT};
//////////
//Tracers:
//Offsets in time
float timeOffsets[tracer_count];
ccl_2 tracer_cl[tracer_count];
color tracer_color[tracer_count];
memset(tracer_cl, 0, tracer_count * sizeof(ccl_2));
memset(tracer_color, 0, tracer_count * sizeof(color));
for(unsigned i = 0; i < tracer_count; i++) {
tracer_cl[i].x0 = tracer_cl[i].y0 = -128;
tracer_cl[i].scale = 0;
timeOffsets[i] = 0;
//Start off offscreeen.
//randomize_ccl_2(tracer_cl + i, imageWidth / 2, imageWidth / 2, imageWidth / 2);
//randomize_color(tracer_color + i, -1, 2, 1);
}
/////////////////////
//Run the Simulation:
//Move along the curves, adding to the wind.
float tStep = totalTime / frames;
float drag = powf(0.75, tStep);
unsigned particleCount = 0; //Keep track of the total number of particles.
for(unsigned i = 0; i < frames; i++) {
float t = tStep * i;
//if(i % 10 == 0) {
printf("Frame %u / %u: Particles %u / %u\r", i, frames, particleCount, wind_size * wind_count);
fflush(stdout);
//}
particleCount = 0;
//Blur and darken the existing image.
/*
fill_image_a(img, 0, 0, 0, 1.0 / 32);
image_blur_inplace(img, (i / 4) % 1);
*/
if(i % 2 == 0) {
float m = 1.0 - 64.0 / imageWidth;
if(m < 0) m = 0;
//Darken one color channel at a time.
image_multiply_channel(img, (i / 2) % 3, m);
//fill_image_a(img, 0, 0, 0, 32.0 / imageWidth);
} else {
image_blur_fast_inplace(img, i / 2);
}
//"Gravity Vortex": bend light inward. (TODO move this to filters).
/*
for(unsigned i = 0; i < imageWidth * imageWidth / 128; i++) {
vec2 p = symmetricBall(imageWidth / 4);
unsigned x0 = imageWidth / 2 + (unsigned)p.x;
unsigned y0 = imageWidth / 2 + (unsigned)p.y;
unsigned x1 = imageWidth / 2 + (unsigned)(p.x * 1.125);
unsigned y1 = imageWidth / 2 + (unsigned)(p.y * 1.125);
for(unsigned c = 0; c < C; c++) {
float* p0 = image_pixel(img, x0, y0, c);
float* p1 = image_pixel(img, x1, y1, c);
*p0 = (*p0 + *p1 * 0.75);
*p1 *= 0.25;
}
}
*/
//image_blur_fast_inplace(img, i);
//image_blur_fast_inplace(img, i + 1);
if(i % 32 == 0) {
//Reset a curve
unsigned idx = uniformInt(0, tracer_count);
randomize_ccl_2(tracer_cl + idx, imageWidth / 2, imageWidth / 2, imageWidth / 2);
for(unsigned j = 0; j < 4; j++) {
tracer_cl[idx].d[j] *= 1.0 / 32.0;
}
//Randomize color.
//randomize_color(tracer_color + idx, -0.5, 1.5, 1);
//tracer_color[idx].a = uniformFloat(-.25, 1.5); //Negative alpha can cause some interesting effects (existing color is emphasized, new is subtracted).
//Randomize color, chosen to work well with additive coloring.
randomize_color(tracer_color + idx, -0.125, 0.5, 1);
timeOffsets[idx] = t;
}
/*
if(i % 64 == 0) {
randomize_ccl_2(&ccl, imageWidth / 2, imageWidth / 2, imageWidth);
for(unsigned i = 0; i < 4; i++) {
ccl.d[i] *= 1.0 / 16.0;
}
//draw_parametric_curve_uniform_space(img, draw_circle, parametric_curve_1, 0, TAU, 0.01, 8, 0.9, &drawcl_2, &cl);
}
*/
/*
float c2Speed = 1;
{
//color drawcl = {uniformFloat(0, 4), uniformFloat(0, 4), uniformFloat(0, 4), uniformFloat(1.0 / 16.0, 1)};
color drawcl = {uniformFloat(-1, 2), uniformFloat(-1, 2), uniformFloat(-1, 2), uniformFloat(1.0 / 16.0, 1)};
circle_cl drawcl_2 = {1, 0, 1, 0.5, 5.5};
draw_parametric_curve_uniform_time(img, draw_point, parametric_curve_2, i % 64 * c2Speed, (i % 64 + 1) * c2Speed, 0.01, &drawcl, &ccl);
//ccl.xS *= 0.95;
//ccl.yS *= 0.95;
}
*/
//Render a partial curve.
float tSpeed = 16.0;
for(unsigned j = 0; j < tracer_count; j ++) {
float renderSpacing = (1.0 / tracer_cl[j].scale) * (1.0);
if(tracer_cl[j].scale < 1) continue;
//Draw a partial curve
if(j < tracer_count / 2) {
float tt = (t - timeOffsets[j]) * tSpeed;
float tt2 = (t - timeOffsets[j] + tStep) * tSpeed;
//draw_parametric_curve_uniform_time(img, draw_point, parametric_curve_2, 0, uniformFloat(0, 200), 0.01, &drawcl, cl[j]);
draw_parametric_curve_uniform_time(img, draw_point_additive, parametric_curve_2, tt, tt2, renderSpacing * 0.1, &tracer_color[j], &tracer_cl[j]);
//color drawcl2 = {windColors[j * 3 + 1], windColors[j * 3 + 2], windColors[j * 3 + 0], 3};
//draw_parametric_curve_uniform_time(img, draw_point, parametricCurves[j], 0, tt, 0.01, &drawcl2, cl[j]);
} else {
//Flash and permute a full curve.
if(uniformInt(0, 2) == 0) {
continue;
}
float sd = 32.0;
float pd = TAU / 128.0;
tracer_cl[j].scale *= uniformFloat((sd - 1) / sd, (sd + 1) / sd); //Should shrink slightly over time.
tracer_cl[j].theta += uniformFloat(-pd, pd);
float dpMax = imageWidth / 1000.0;
vec2 dp = vScale(dpMax, symmetricUnitBall());
tracer_cl[j].x0 += dp.x;
tracer_cl[j].y0 += dp.y;
float dd = 128.0;
for(unsigned i = 0; i < PC2_SUMMANDS; i++) {
tracer_cl[j].d[i] *= uniformFloat((dd - 1) / dd, (dd + 1) / dd);
tracer_cl[j].f[i] *= uniformFloat((dd - 1) / dd, (dd + 1) / dd);
tracer_cl[j].p[i] += uniformFloat(-pd, pd);
}
//draw_parametric_curve_uniform_time(img, draw_point_additive, parametric_curve_2, 0, uniformFloat(0, 200), renderSpacing * (1 + uniformFloat(0, 1)), &tracer_color[j], &tracer_cl[j]);
//random_cl rcl = {3, 1, &tracer_color[j]};
//draw_parametric_curve_uniform_space(img, draw_random_additive, parametric_curve_2, 0, uniformFloat(0, 200), 0.01, 2, 0.9, &rcl, &tracer_cl[j]);
draw_parametric_curve_uniform_space(img, draw_point_additive, parametric_curve_2, 0, uniformFloat(0, 200), 0.01, 1, 0.9, &tracer_color[j], &tracer_cl[j]);
}
}
//Permute tracer colors.
for(unsigned i = 0; i < tracer_count; i++) {
for(unsigned c = 0; c < C; c++) {
tracer_color[i].c[c] += uniformFloat(-tracerColorChange, tracerColorChange);
tracer_color[i].c[c] *= uniformFloat(tracerColorReduction, 1);
}
//tracer_color[i].a += uniformFloat(-tracerColorChange, tracerColorChange);
//tracer_color[i].a *= tracerColorReduction;
}
for(unsigned i = 0; i < wind_count * 3; i++) {
windColors[i] = windColors[i] * (1 - wColorRevert) + windDefaultColors[i] * wColorRevert + uniformFloat(-wColorChange, wColorChange);
}
for(unsigned j = 0; j < wind_count; j++) {
//Randomly remove some particles.
wind_remove_rand(winds[j], 512);
//Apply a bit of drag.
wind_scale_velocity(winds[j], drag);
//Calculate the position and velocity of the emitter.
vec2 v = parametricCurves[j](t, cl[j]);
vec2 v2 = parametricCurves[j](t + DEPSILON, cl[j]);
vec2 dcdt = vScale((1.0 / DEPSILON), vMinus(v2, v));
for(unsigned k = 0; k < particlesPerFrame; k++) {
vec2 noise = vScale((particleNoiseVelocity), symmetricUnitBall());
vec2 dvdt = vPlus(vScale(particleEmitterMomentumScalar, dcdt), noise);
wind_append(winds[j], v.x, v.y, dvdt.x, dvdt.y, uniformFloat(1, 2));
}
//printf("Particles: %u\n", winds[j]->particles);
//wind_print(winds[j], stdout);
//Update each particle.
wind_update_bound(winds[j], tStep, sumWeightedPotential, &pcl, -(imageWidth / 2.0), -(imageHeight / 2.0), 3.0 * imageWidth / 2.0, 3.0 * imageWidth / 2.0);
//wind_update(winds[j], 1.0, sumWeightedPotential, &pcl);
//Draw the wind.
wind_draw_roffset(winds[j], img, windColors[j * 3 + 0], windColors[j * 3 + 1], windColors[j * 3 + 2], wAlpha, 0, 0, 1, particlesToDraw, maxSpread);
particleCount += winds[j]->particles;
//TODO want to change the noise over time so wind does not collect in low pressure areas.
//TODO render the noise function.
//TODO periodically reset wind emmiters for type 2 curves.
}
p(img, NULL); //Process the image (send it to the screen or to disk).
}
image_free(img);
for(unsigned i = 0; i < wind_count; i++) {
wind_free(winds[i]);
}
noise_sum_free(ncl);
}
