#pragma once

#ifndef _AFFINE_H_
#define _AFFINE_H_

#undef TNZAPI
#ifdef TNZ_IS_RASTERLIB
#define TNZAPI TNZ_EXPORT_API
#else
#define TNZAPI TNZ_IMPORT_API
#endif

typedef struct {        /* to be PRE-multiplied with a COLUMN vector, i.e., */
  double a11, a12, a13; /* x = a11*u + a12*v + a13 */
  double a21, a22, a23; /* y = a21*u + a22*v + a23 */
} AFFINE;

#define AFF_M_V_1(AFF, V1, V2) ((AFF).a11 * (V1) + (AFF).a12 * (V2) + (AFF).a13)
#define AFF_M_V_2(AFF, V1, V2) ((AFF).a21 * (V1) + (AFF).a22 * (V2) + (AFF).a23)

static const AFFINE Aff_I = {1.0, 0.0, 0.0, 0.0, 1.0, 0.0};

TNZAPI AFFINE aff_shift(double dx, double dy, AFFINE aff);
TNZAPI AFFINE aff_shift_post(double dx, double dy, AFFINE aff);
TNZAPI AFFINE aff_rot(double degrees, AFFINE aff);
TNZAPI AFFINE aff_scale(double sx, double sy, AFFINE aff);
TNZAPI AFFINE aff_shear(double shx, double shy, AFFINE aff);
TNZAPI AFFINE aff_place(double u, double v, double x, double y, AFFINE aff);
TNZAPI AFFINE aff_mult(AFFINE aff1, AFFINE aff2);
TNZAPI AFFINE aff_inv(AFFINE aff);
TNZAPI double aff_det(AFFINE aff);

#define AFF_FILL(AFF, A11, A12, A13, A21, A22, A23)                            \
  {                                                                            \
    (AFF).a11 = A11;                                                           \
    (AFF).a12 = A12;                                                           \
    (AFF).a13 = A13;                                                           \
    (AFF).a21 = A21;                                                           \
    (AFF).a22 = A22;                                                           \
    (AFF).a23 = A23;                                                           \
  }

#define AFF_EQ(A, B)                                                           \
  ((A).a11 == (B).a11 && (A).a12 == (B).a12 && (A).a13 == (B).a13 &&           \
   (A).a21 == (B).a21 && (A).a22 == (B).a22 && (A).a23 == (B).a23)

#define AFF_NE(A, B)                                                           \
  ((A).a11 != (B).a11 || (A).a12 != (B).a12 || (A).a13 != (B).a13 ||           \
   (A).a21 != (B).a21 || (A).a22 != (B).a22 || (A).a23 != (B).a23)

#define ROTATE(ang, cxin, cyin, cxout, cyout, aff)                             \
  aff_shift(cxout, cyout, aff_rot(ang, aff_shift(-cxin, -cyin, aff)))

#define ROTATE_POST(ang, cxin, cyin, cxout, cyout, aff)                        \
  aff_mult(aff, aff_shift(cxout, cyout,                                        \
                          aff_rot(ang, aff_shift(-cxin, -cyin, Aff_I))))

#define SCALE(sx, sy, cxin, cyin, cxout, cyout, aff)                           \
  aff_shift(cxout, cyout, aff_scale(sx, sy, aff_shift(-cxin, -cyin, aff)))

#define SCALE_POST(sx, sy, cxin, cyin, cxout, cyout, aff)                      \
  aff_mult(aff, aff_shift(cxout, cyout,                                        \
                          aff_scale(sx, sy, aff_shift(-cxin, -cyin, Aff_I))))

#endif