summaryrefslogtreecommitdiff
path: root/flam3-animate.c
diff options
context:
space:
mode:
Diffstat (limited to 'flam3-animate.c')
-rw-r--r--flam3-animate.c342
1 files changed, 342 insertions, 0 deletions
diff --git a/flam3-animate.c b/flam3-animate.c
new file mode 100644
index 0000000..3c32628
--- /dev/null
+++ b/flam3-animate.c
@@ -0,0 +1,342 @@
+/*
+ FLAM3 - cosmic recursive fractal flames
+ Copyright (C) 1992-2009 Spotworks LLC
+
+ This program is free software; you can redistribute it and/or modify
+ it under the terms of the GNU General Public License as published by
+ the Free Software Foundation; either version 3 of the License, or
+ (at your option) any later version.
+
+ This program is distributed in the hope that it will be useful,
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ GNU General Public License for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with this program. If not, see <http://www.gnu.org/licenses/>.
+*/
+
+
+#include "private.h"
+#include "img.h"
+#include "isaacs.h"
+
+int main(int argc, char **argv) {
+ char *ai, *fname;
+ char flamename[256];
+ char *prefix = args("prefix", "");
+ int first_frame = argi("begin", 0);
+ int last_frame = argi("end", 0);
+ int frame_time = argi("time", 0);
+ int dtime = argi("dtime", 1);
+ int do_fields = argi("fields", 0);
+ int write_genome = argi("write_genome",0);
+ double qs = argf("qs", 1.0);
+ double ss = argf("ss", 1.0);
+ char *format = getenv("format");
+ int verbose = argi("verbose", 1);
+ int transparency = argi("transparency", 0);
+ int bits = argi("bits", 33);
+ int sub_batch_size = argi("sub_batch_size", 10000);
+ int bpc = argi("bpc",8);
+ int earlyclip = argi("earlyclip",0);
+ int ftime, channels;
+ unsigned char *image;
+ flam3_genome *cps,center_cp;
+ int i, ncps = 0;
+ char *inf = getenv("in");
+ double pixel_aspect = argf("pixel_aspect", 1.0);
+ int num_threads = argi("nthreads",0);
+ FILE *in,*fp,*genfp;
+ flam3_frame f;
+ flam3_img_comments fpc;
+ stat_struct stats,stats2;
+
+ char badval_string[64];
+ char numiter_string[64];
+ char rtime_string[64];
+
+#ifdef WIN32
+
+ char *slashloc;
+ char exepath[256];
+ char palpath[256];
+ memset(exepath,0,256);
+ memset(palpath,0,256);
+ slashloc = strrchr(argv[0],'\\');
+ if (NULL==slashloc) {
+ sprintf(palpath,"flam3_palettes=flam3-palettes.xml");
+ } else {
+ strncpy(exepath,argv[0],slashloc-argv[0]+1);
+ sprintf(palpath,"flam3_palettes=%sflam3-palettes.xml",exepath);
+ }
+ putenv(palpath);
+
+#endif
+
+
+ memset(&center_cp,0, sizeof(flam3_genome));
+
+ if (1 != argc) {
+ docstring();
+ exit(0);
+ }
+
+ /* Init random number generators */
+ flam3_init_frame(&f);
+ flam3_srandom();
+
+ /* Set the number of threads */
+ if (num_threads==0) {
+ num_threads = flam3_count_nthreads();
+ if (verbose > 1)
+ fprintf(stderr,"Automatically detected %d core(s)...\n",num_threads);
+ } else{
+ if (verbose)
+ fprintf(stderr,"Manually specified %d thread(s)...\n",num_threads);
+ }
+
+ if (getenv("frame")) {
+ if (getenv("time")) {
+ fprintf(stderr, "cannot specify both time and frame.\n");
+ exit(1);
+ }
+ if (getenv("begin") || getenv("end")) {
+ fprintf(stderr, "cannot specify both frame and begin or end.\n");
+ exit(1);
+ }
+ first_frame = last_frame = atoi(getenv("frame"));
+ }
+
+ if (getenv("time")) {
+ if (getenv("begin") || getenv("end")) {
+ fprintf(stderr, "cannot specify both time and begin or end.\n");
+ exit(1);
+ }
+ first_frame = last_frame = frame_time;
+ }
+
+ if (NULL == format) format = "png";
+ if (strcmp(format, "jpg") &&
+ strcmp(format, "ppm") &&
+ strcmp(format, "png")) {
+ fprintf(stderr, "format must be either jpg, ppm, or png, not %s.\n", format);
+ exit(1);
+ }
+
+
+ if (pixel_aspect <= 0.0) {
+ fprintf(stderr, "pixel aspect ratio must be positive, not %g.\n",
+ pixel_aspect);
+ exit(1);
+ }
+
+ if (inf)
+ in = fopen(inf, "rb");
+ else
+ in = stdin;
+ if (NULL == in) {
+ perror(inf);
+ exit(1);
+ }
+
+ cps = flam3_parse_from_file(in, inf, flam3_defaults_on, &ncps);
+ if (inf)
+ fclose(in);
+ if (NULL == cps) {
+ fprintf(stderr," error reading genomes.\n");
+ exit(1);
+ }
+ if (0 == ncps) {
+ fprintf(stderr, "error: no genomes.\n");
+ exit(1);
+ }
+ for (i = 0; i < ncps; i++) {
+ cps[i].sample_density *= qs;
+ cps[i].height = (int)(cps[i].height * ss);
+ cps[i].width = (int)(cps[i].width * ss);
+ cps[i].pixels_per_unit *= ss;
+ if (cps[i].height<=0 || cps[i].width<=0) {
+ fprintf(stderr,"output image has dimension <=0, aborting.\n");
+ exit(1);
+ }
+ if (i > 0 && cps[i].time <= cps[i-1].time) {
+ fprintf(stderr, "error: control points must be sorted by time, but %g <= %g, index %d.\n",
+ cps[i].time, cps[i-1].time, i);
+ exit(1);
+ }
+ if ((cps[i].width != cps[0].width) ||
+ (cps[i].height != cps[0].height)) {
+ fprintf(stderr, "warning: flame %d at time %g size mismatch. "
+ "(%d,%d) should be (%d,%d).\n",
+ i, cps[i].time,
+ cps[i].width, cps[i].height,
+ cps[0].width, cps[0].height);
+ cps[i].width = cps[0].width;
+ cps[i].height = cps[0].height;
+ }
+ }
+ if (!getenv("time") && !getenv("frame")) {
+ if (!getenv("begin")) {
+ first_frame = (int) cps[0].time;
+ }
+ if (!getenv("end")) {
+ last_frame = (int) cps[ncps-1].time - 1;
+ if (last_frame < first_frame) last_frame = first_frame;
+ }
+ }
+ channels = strcmp(format, "png") ? 3 : 4;
+
+ /* Check for 16-bit-per-channel processing */
+ if ( (16 == bpc) && (strcmp(format,"png") != 0)) {
+ fprintf(stderr,"Support for 16 bpc images is only present for the png format.\n");
+ exit(1);
+ } else if (bpc != 8 && bpc != 16) {
+ fprintf(stderr,"Unexpected bpc specified (%d)\n",bpc);
+ exit(1);
+ }
+
+// f.temporal_filter_radius = argf("blur", 0.5);
+ f.pixel_aspect_ratio = pixel_aspect;
+ f.genomes = cps;
+ f.ngenomes = ncps;
+ f.verbose = verbose;
+ f.bits = bits;
+ f.progress = 0;
+ f.earlyclip = earlyclip;
+ f.nthreads = num_threads;
+ f.sub_batch_size = sub_batch_size;
+
+ if (16==bpc)
+ f.bytes_per_channel = 2;
+ else
+ f.bytes_per_channel = 1;
+
+
+ image = (void *) calloc((size_t)channels *
+ (size_t)cps[0].width *
+ (size_t)cps[0].height * f.bytes_per_channel, sizeof(char));
+
+ if (dtime < 1) {
+ fprintf(stderr, "dtime must be positive, not %d.\n", dtime);
+ exit(1);
+ }
+
+ for (ftime = first_frame; ftime <= last_frame; ftime += dtime) {
+ f.time = (double) ftime;
+
+ if (verbose && ((last_frame-first_frame)/dtime) >= 1) {
+ fprintf(stderr, "time = %d/%d/%d\n", ftime, last_frame, dtime);
+ }
+
+ if (do_fields) {
+
+ if (flam3_render(&f, image, flam3_field_even, channels, transparency,&stats)) {
+ fprintf(stderr,"error rendering image: aborting.\n");
+ exit(1);
+ }
+ f.time += 0.5;
+ if (flam3_render(&f, image, flam3_field_odd, channels, transparency,&stats2)) {
+ fprintf(stderr,"error rendering image: aborting.\n");
+ exit(1);
+ }
+
+ stats.badvals+=stats2.badvals;
+ stats.render_seconds+=stats2.render_seconds;
+ stats.num_iters+=stats2.num_iters;
+ } else {
+ if (flam3_render(&f, image, flam3_field_both, channels, transparency,&stats)) {
+ fprintf(stderr,"error rendering image: aborting.\n");
+ exit(1);
+ }
+ }
+
+ if (getenv("out"))
+ fname = getenv("out");
+ else {
+ fname = malloc(strlen(prefix) + 20);
+ sprintf(fname, "%s%05d.%s", prefix, ftime, format);
+ }
+
+ if (verbose)
+ fprintf(stderr, "writing %s...", fname);
+
+ if (write_genome) {
+ sprintf(flamename,"%s.flam3",fname);
+
+ /* get center genome */
+ flam3_interpolate(f.genomes, f.ngenomes, f.time, 0, &center_cp);
+
+ /* write it out */
+ genfp = fopen(flamename,"w");
+ if (NULL == genfp) {
+ perror(flamename);
+ exit(1);
+ }
+
+ flam3_print(genfp, &center_cp, NULL, flam3_print_edits);
+
+ fclose(genfp);
+ }
+
+ fp = fopen(fname, "wb");
+ if (NULL == fp) {
+ perror(fname);
+ exit(1);
+ }
+
+ /* Get center cp for embedding in png file */
+ flam3_interpolate(f.genomes, f.ngenomes, f.time, 0, &center_cp);
+
+ /* Convert to string */
+ fpc.genome = flam3_print_to_string(&center_cp);
+ sprintf(badval_string, "%g",stats.badvals/(double)stats.num_iters);
+ fpc.badvals = badval_string;
+ sprintf(numiter_string,"%g",(double)stats.num_iters);
+ fpc.numiters = numiter_string;
+ sprintf(rtime_string,"%d",stats.render_seconds);
+ fpc.rtime = rtime_string;
+
+ if (!strcmp(format, "png")) {
+
+ write_png(fp, image, cps[0].width, cps[0].height, &fpc, f.bytes_per_channel);
+
+ } else if (!strcmp(format, "jpg")) {
+
+ write_jpeg(fp, image, cps[0].width, cps[0].height, &fpc);
+
+ } else {
+ int size = 3 * cps[0].width * cps[0].height;
+ fprintf(fp, "P6\n");
+ fprintf(fp, "%d %d\n255\n", cps[0].width, cps[0].height);
+ if (size != fwrite(image, 1, size, fp)) {
+ perror(fname);
+ }
+ }
+
+ /* Free string */
+ free(fpc.genome);
+
+ clear_cp(&center_cp,0);
+
+ fclose(fp);
+
+ if (!getenv("out"))
+ free(fname);
+
+ }
+
+ for (i = 0; i < ncps; i++) {
+ xmlFreeDoc(cps[i].edits);
+ clear_cp(&cps[i],0);
+ }
+ free(cps);
+
+ free(image);
+
+ if (verbose)
+ fprintf(stderr, "done.\n");
+
+ fflush(stderr);
+ return 0;
+}