/*
 * Copyright 2011-2013 Blender Foundation
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 * http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */

#ifndef __SVM_RAMP_H__
#define __SVM_RAMP_H__

CCL_NAMESPACE_BEGIN

/* NOTE: svm_ramp.h, svm_ramp_util.h and node_ramp_util.h must stay consistent */

ccl_device_inline float4 rgb_ramp_lookup(
    KernelGlobals *kg, int offset, float f, bool interpolate, bool extrapolate, int table_size)
{
  if ((f < 0.0f || f > 1.0f) && extrapolate) {
    float4 t0, dy;
    if (f < 0.0f) {
      t0 = fetch_node_float(kg, offset);
      dy = t0 - fetch_node_float(kg, offset + 1);
      f = -f;
    }
    else {
      t0 = fetch_node_float(kg, offset + table_size - 1);
      dy = t0 - fetch_node_float(kg, offset + table_size - 2);
      f = f - 1.0f;
    }
    return t0 + dy * f * (table_size - 1);
  }

  f = saturate(f) * (table_size - 1);

  /* clamp int as well in case of NaN */
  int i = clamp(float_to_int(f), 0, table_size - 1);
  float t = f - (float)i;

  float4 a = fetch_node_float(kg, offset + i);

  if (interpolate && t > 0.0f)
    a = (1.0f - t) * a + t * fetch_node_float(kg, offset + i + 1);

  return a;
}

ccl_device void svm_node_rgb_ramp(
    KernelGlobals *kg, ShaderData *sd, float *stack, uint4 node, int *offset)
{
  uint fac_offset, color_offset, alpha_offset;
  uint interpolate = node.z;

  svm_unpack_node_uchar3(node.y, &fac_offset, &color_offset, &alpha_offset);

  uint table_size = read_node(kg, offset).x;

  float fac = stack_load_float(stack, fac_offset);
  float4 color = rgb_ramp_lookup(kg, *offset, fac, interpolate, false, table_size);

  if (stack_valid(color_offset))
    stack_store_float3(stack, color_offset, float4_to_float3(color));
  if (stack_valid(alpha_offset))
    stack_store_float(stack, alpha_offset, color.w);

  *offset += table_size;
}

ccl_device void svm_node_curves(
    KernelGlobals *kg, ShaderData *sd, float *stack, uint4 node, int *offset)
{
  uint fac_offset, color_offset, out_offset;
  svm_unpack_node_uchar3(node.y, &fac_offset, &color_offset, &out_offset);

  uint table_size = read_node(kg, offset).x;

  float fac = stack_load_float(stack, fac_offset);
  float3 color = stack_load_float3(stack, color_offset);

  const float min_x = __int_as_float(node.z), max_x = __int_as_float(node.w);
  const float range_x = max_x - min_x;
  const float3 relpos = (color - make_float3(min_x, min_x, min_x)) / range_x;

  float r = rgb_ramp_lookup(kg, *offset, relpos.x, true, true, table_size).x;
  float g = rgb_ramp_lookup(kg, *offset, relpos.y, true, true, table_size).y;
  float b = rgb_ramp_lookup(kg, *offset, relpos.z, true, true, table_size).z;

  color = (1.0f - fac) * color + fac * make_float3(r, g, b);
  stack_store_float3(stack, out_offset, color);

  *offset += table_size;
}

CCL_NAMESPACE_END

#endif /* __SVM_RAMP_H__ */
