mandelbrot.h

/*
 * Dynamic MPI Mandelbrot algorithm
 * Copyright (C) 2015  Martin Ohmann <martin@mohmann.de>
 *
 * 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/>.
 */
#ifndef _MO_MANDELBROT_H
#define _MO_MANDELBROT_H

#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <getopt.h>
#include <mpi.h>

/* 
 * default values
 */
#define MO_N           2                    /* size of problem space (x, y from -MO_N to MO_N) */
#define MO_THRESHOLD   4                    /* divergence threshold (usually 4) */
#define MO_SIZE        1024                 /* height & width of the canvas in pixels */
#define MO_FILENAME    "image"              /* default filename of resulting bitmap */
#define MO_MAXITER     2000                 /* default maximum iterations */
#define MO_BLOCKSIZE   1                    /* default # of blocks (rows) to be assigned for 
                                               computation to a slave each time */
#define MO_COLORMIN    0x000000             /* default min color (0x000000 = black) */
#define MO_COLORMAX    0xffffff             /* default max color (0xffffff = white) */
#define MO_COLORMASK   0xffffff             /* default color mask (all tones) */
#define MO_PROGRESS    0                    /* show (1) or hide (0) progress  */
#define MO_PWIDTH      50                   /* progress bar width */
#define MO_PUPDATE     20                   /* update progress bar MO_UPDATE times */
#define MO_ZOOMSTEPS   1                    /* default number of zoom steps to take */                                               
#define MO_ZOOMFACTOR  0.75                 /* default zoom factor per step */

/*
 * communication flags
 */
#define MO_CALC        1                    /* "calc" message flag (master to slave) */
#define MO_DATA        2                    /* "data" message flag (slave to master) */
#define MO_STOP        3                    /* "stop" message flag (master to slave) */

/*
 * function marcos
 */
#define eprintf(...) { \
    fprintf(stderr, "Error: "); \
    fprintf(stderr, __VA_ARGS__); \
}

#define finalize_exit(error_code) { \
    MPI_Finalize(); \
    exit((error_code)); \
}

/* 
 * structdef for complex numbers 
 */
typedef struct _mo_complex
{
    double re;                  /* real part */
    double im;                  /* imaginary part */
} mo_complex_t;

/* 
 * stuctdef for core options 
 */
typedef struct _mo_opts
{
    int width, height;          /* image width/height */
    char *filename;             /* filename of resulting bitmap */
    double max_re, min_re;      /* real value range */
    double max_im, min_im;      /* imaginary value range */
    int max_iterations;         /* iteration maximum */
    int blocksize;              /* # of rows to be send to slave at once */
    long min_color, max_color;  /* color ranges */
    long color_mask;            /* color mask */
    int show_progress;          /* if 1, show progress */
    int zoom_steps;             /* number of zoom steps */
    double zoom_factor;         /* zoom factor per step */
    double axis_length;         /* axis length */
    double x_offset, y_offset;  /* x & y offset */
} mo_opts_t;

/*
 * structdef for scaling
 */
typedef struct _mo_scale
{
    double color, im, re;       /* scalings of color, imaginary and 
                                   real part for computation */
} mo_scale_t;

/*
 * structdef for bmp header
 * derived from http://cpansearch.perl.org/src/DHUNT/PDL-Planet-0.05/libimage/bmp.c
 */
typedef struct _mo_bmp_header
{
    char type[2];               /* "BM" */
    int fsize;                  /* Size of file in bytes */
    int reserved;               /* set to 0 */
    int offbits;                /* Byte offset to actual bitmap data (= 54) */
    int hsize;                  /* Size of BITMAPINFOHEADER, in bytes (= 40) */
    int width;                  /* Width of image, in pixels */
    int height;                 /* Height of images, in pixels */
    short planes;               /* Number of planes in target device (set to 1) */
    short bit_count;            /* Bits per pixel (24 in this case) */
    int compression;            /* Type of compression (0 if no compression) */
    int size_image;             /* Image size, in bytes (0 if no compression) */
    int x_pels_per_meter;       /* Resolution in pixels/meter of display device */
    int y_pels_per_meter;       /* Resolution in pixels/meter of display device */
    int clr_used;               /* Number of colors in the color table (if 0, use 
                                   maximum allowed by bit_count) */
    int clr_important;          /* Number of important colors. If 0, all colors 
                                   are important */
} mo_bmp_header_t;

typedef struct COLOR 
{
    char r,g,b;
} COLOR_t;


/* 
 * method prototypes 
 */
static int parse_args(int, char **, mo_opts_t *, int, int);
static void print_params(mo_opts_t *);
static void print_usage(char **);
static int master_proc(int, mo_opts_t *, int);
static int slave_proc(int, mo_opts_t *);
static long mandelbrot(int, int, mo_scale_t *, mo_opts_t *);
static inline void print_progress(int, int);
static int write_bitmap(const char *, int, int, char *);
static COLOR_t HSVtoRGB(double, double, double);


#endif /* _MO_MANDELBROT_H */