/* vim: set ts=8 sw=8 noexpandtab: */

//  qcms

//  Copyright (C) 2009 Mozilla Foundation

//  Copyright (C) 1998-2007 Marti Maria

//

// Permission is hereby granted, free of charge, to any person obtaining 

// a copy of this software and associated documentation files (the "Software"), 

// to deal in the Software without restriction, including without limitation 

// the rights to use, copy, modify, merge, publish, distribute, sublicense, 

// and/or sell copies of the Software, and to permit persons to whom the Software 

// is furnished to do so, subject to the following conditions:

//

// The above copyright notice and this permission notice shall be included in 

// all copies or substantial portions of the Software.

//

// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, 

// EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO 

// THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND 

// NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE 

// LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION 

// OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION 

// WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

#ifndef _QCMS_TRANSFORM_UTIL_H

#define _QCMS_TRANSFORM_UTIL_H

#include <stdlib.h>

#include <math.h>

#define CLU(table,x,y,z) table[(x*len + y*x_len + z*xy_len)*3]

//XXX: could use a bettername

typedef uint16_t uint16_fract_t;

float lut_interp_linear(double value, uint16_t *table, int length);

float lut_interp_linear_float(float value, float *table, int length);

uint16_t lut_interp_linear16(uint16_t input_value, uint16_t *table, int length);

static inline float lerp(float a, float b, float t)

{

        return a*(1.f-t) + b*t;

}

Unexecuted instantiation: chain.c:lerp

Unexecuted instantiation: transform.c:lerp

Unexecuted instantiation: transform_util.c:lerp

Unexecuted instantiation: Unified_cpp_gfx_tests_gtest1.cpp:lerp(float, float, float)

static inline unsigned char clamp_u8(float v)

{

  if (v > 255.)

    return 255;

  else if (v < 0)

    return 0;

  else

    return floorf(v+.5);

}

Unexecuted instantiation: chain.c:clamp_u8

Unexecuted instantiation: transform.c:clamp_u8

Unexecuted instantiation: transform_util.c:clamp_u8

Unexecuted instantiation: Unified_cpp_gfx_tests_gtest1.cpp:clamp_u8(float)

static inline float clamp_float(float a)

{

  /* One would naturally write this function as the following:

  if (a > 1.)

    return 1.;

  else if (a < 0)

    return 0;

  else

    return a;

  However, that version will let NaNs pass through which is undesirable

  for most consumers.

  */

  if (a > 1.)

    return 1.;

  else if (a >= 0)

    return a;

  else // a < 0 or a is NaN

    return 0;

}

Unexecuted instantiation: chain.c:clamp_float

Unexecuted instantiation: transform.c:clamp_float

Unexecuted instantiation: transform_util.c:clamp_float

Unexecuted instantiation: Unified_cpp_gfx_tests_gtest1.cpp:clamp_float(float)

static inline float u8Fixed8Number_to_float(uint16_t x)

{

  // 0x0000 = 0.

  // 0x0100 = 1.

  // 0xffff = 255  + 255/256

  return x/256.;

}

Unexecuted instantiation: chain.c:u8Fixed8Number_to_float

Unexecuted instantiation: transform.c:u8Fixed8Number_to_float

Unexecuted instantiation: transform_util.c:u8Fixed8Number_to_float

Unexecuted instantiation: Unified_cpp_gfx_tests_gtest1.cpp:u8Fixed8Number_to_float(unsigned short)

float *build_input_gamma_table(struct curveType *TRC);

struct matrix build_colorant_matrix(qcms_profile *p);

void build_output_lut(struct curveType *trc,

                      uint16_t **output_gamma_lut, size_t *output_gamma_lut_length);

struct matrix matrix_invert(struct matrix mat);

qcms_bool compute_precache(struct curveType *trc, uint8_t *output);

// Tested by GTest

#ifdef  __cplusplus

extern "C" {

#endif

uint16_fract_t lut_inverse_interp16(uint16_t Value, uint16_t LutTable[], int length);

#ifdef  __cplusplus

}

#endif

#endif