#ifndef DT_DEVELOP_IMAGEOP_H

#define DT_DEVELOP_IMAGEOP_H

#include "common/darktable.h"

#include "control/settings.h"

#include "develop/pixelpipe.h"

#include "dtgtk/togglebutton.h"

#include <gmodule.h>

#include <gtk/gtk.h>

#include <sched.h>

struct dt_develop_t;

struct dt_dev_pixelpipe_t;

struct dt_dev_pixelpipe_iop_t;

struct dt_iop_roi_t;

struct dt_develop_blend_params_t;

struct dt_develop_tiling_t;

#define IOP_GROUP_BASIC 1

#define IOP_GROUP_COLOR 2

#define IOP_GROUP_CORRECT 4

#define IOP_GROUP_EFFECT 8

#define IOP_GROUP_TONE 16

#define IOP_SPECIAL_GROUP_ACTIVE_PIPE 32

#define IOP_SPECIAL_GROUP_USER_DEFINED 64

#define IOP_TAG_DISTORT 1

#define IOP_TAG_DECORATION 2

#define IOP_GROUP_ALL (IOP_GROUP_BASIC|IOP_GROUP_COLOR|IOP_GROUP_CORRECT|IOP_GROUP_EFFECT)

#define IOP_FLAGS_INCLUDE_IN_STYLES 1

#define IOP_FLAGS_SUPPORTS_BLENDING 2 // Does provide blending modes

#define IOP_FLAGS_DEPRECATED 4

#define IOP_FLAGS_BLEND_ONLY_LIGHTNESS 8 // Does only blend with L-channel in Lab space. Keeps a, b of original image.

#define IOP_FLAGS_ALLOW_TILING 16 // Does allow tile-wise processing (valid for CPU and GPU processing)

#define IOP_FLAGS_HIDDEN 32 // Hide the iop from userinterface

#define IOP_FLAGS_TILING_FULL_ROI 64 // Tiling code has to expect arbitrary roi's for this module (incl. flipping, mirroring etc.)

typedef struct dt_iop_params_t

{

int keep;

}

dt_iop_params_t;

typedef void dt_iop_gui_data_t;

typedef void dt_iop_data_t;

typedef void dt_iop_global_data_t;

struct dt_iop_module_so_t;

struct dt_iop_module_t;

typedef struct dt_iop_module_so_t

}

dt_iop_module_so_t;

typedef struct dt_iop_module_t

}

dt_iop_module_t;

void dt_iop_load_modules_so();

void dt_iop_unload_modules_so();

GList *dt_iop_load_modules(struct dt_develop_t *dev);

void dt_iop_cleanup_module(dt_iop_module_t *module);

void dt_iop_init_pipe(struct dt_iop_module_t *module, struct dt_dev_pixelpipe_t *pipe, struct dt_dev_pixelpipe_iop_t *piece);

gboolean dt_iop_is_hidden(dt_iop_module_t *module);

void dt_iop_gui_update(dt_iop_module_t *module);

void dt_iop_gui_reset(dt_iop_module_t *module);

void dt_iop_gui_set_expanded(dt_iop_module_t *module, gboolean expanded);

void dt_iop_commit_params(dt_iop_module_t *module, struct dt_iop_params_t *params, struct dt_develop_blend_params_t * blendop_params, struct dt_dev_pixelpipe_t *pipe, struct dt_dev_pixelpipe_iop_t *piece);

GtkWidget *dt_iop_gui_get_expander(dt_iop_module_t *module);

GtkWidget *dt_iop_gui_get_widget(dt_iop_module_t *module);

GtkWidget *dt_iop_gui_get_pluginui(dt_iop_module_t *module);

void dt_iop_request_focus(dt_iop_module_t *module);

void dt_iop_load_default_params(dt_iop_module_t *module);

void dt_iop_reload_defaults(dt_iop_module_t *module);

int dt_iop_breakpoint(struct dt_develop_t *dev, struct dt_dev_pixelpipe_t *pipe);

void dt_iop_nap(uint32_t usec);

typedef enum dt_iop_colorspace_type_t

}

dt_iop_colorspace_type_t;

dt_iop_colorspace_type_t dt_iop_module_colorspace(const dt_iop_module_t *module);

void dt_iop_flip_and_zoom_8( const uint8_t *in, int32_t iw, int32_t ih, uint8_t *out, int32_t ow, int32_t oh, const int32_t orientation, uint32_t *width, uint32_t *height);

void dt_iop_clip_and_zoom(float *out, const float *const in, const struct dt_iop_roi_t *const roi_out, const struct dt_iop_roi_t * const roi_in, const int32_t out_stride, const int32_t in_stride);

dt_iop_clip_and_zoom_demosaic_half_size(

const uint32_t filters);

dt_iop_clip_and_zoom_demosaic_half_size_f(

const float clip);

void dt_iop_clip_and_zoom_8(const uint8_t *i, int32_t ix, int32_t iy, int32_t iw, int32_t ih, int32_t ibw, int32_t ibh,

uint8_t *o, int32_t ox, int32_t oy, int32_t ow, int32_t oh, int32_t obw, int32_t obh);

void dt_iop_YCbCr_to_RGB(const float *yuv, float *rgb);

void dt_iop_RGB_to_YCbCr(const float *rgb, float *yuv);

dt_iop_module_t *get_colorout_module();

gchar *dt_iop_get_localized_name(const gchar * op);

void dt_iop_estimate_cubic(const float *const x, const float *const y, float *a);

static inline float dt_iop_eval_cubic(const float *const a, const float x)

{

// could be sse4.1 _mm_dot_ps

const float x4[4] = {x*x*x, x*x, x, 1.0f};

return a[3]*x4[3] + a[2]*x4[2] + a[1]*x4[1] + a[0]*x4[0];

}

static inline void dt_iop_estimate_exp(const float *const x, const float *const y, const int num, float *coeff)

{

// map every thing to y = y0*(x/x0)^g

// and fix (x0,y0) as the last point.

// assume (x,y) pairs are ordered by ascending x, so this is the last point:

float x0 = x[num-1], y0 = y[num-1];

float g = 0.0f;

int cnt = 0;

// solving for g yields

// g = log(y/y0)/log(x/x0)

//

// average that over the course of the other samples:

for(int k=0;k<num-1;k++)

{

const float yy = y[k]/y0, xx = x[k]/x0;

if(yy > 0.0f && xx > 0.0f)

{

const float gg = logf(y[k]/y0)/log(x[k]/x0);

g += gg;

cnt ++;

}

}

if(cnt)

g *= 1.0f/cnt;

else

g = 1.0f;

coeff[0] = 1.0f/x0;

coeff[1] = y0;

coeff[2] = g;

}

static inline float dt_iop_eval_exp(const float *const coeff, const float x)

{

return coeff[1] * powf(x*coeff[0], coeff[2]);

}

void dt_iop_connect_common_accels(dt_iop_module_t *module);

static inline void dt_iop_alpha_copy(const void *ivoid, void *ovoid, const int width, const int height)

{

#ifdef _OPENMP

#pragma omp parallel for schedule(static) default(none) shared(ovoid,ivoid)

for(int j=0; j<height; j++)

{

const float *in = ((const float *)ivoid)+4*width*j+3;

float *out = ((float *)ovoid)+4*width*j+3;

for(int i=0; i<width; i++)

{

*out = *in;

out += 4;

in += 4;

}

}

}