IccMpeBasic.cpp

Go to the documentation of this file.

/** @file
    File:       IccMpeBasic.cpp
 
    Contains:   Implementation of Basic Multi Processing Elements
 
    Version:    V1
 
    Copyright:  � see ICC Software License
*/
 
/*
 * The ICC Software License, Version 0.2
 *
 *
 * Copyright (c) 2003-2012 The International Color Consortium. All rights 
 * reserved.
 *
 * Redistribution and use in source and binary forms, with or without
 * modification, are permitted provided that the following conditions
 * are met:
 *
 * 1. Redistributions of source code must retain the above copyright
 *    notice, this list of conditions and the following disclaimer. 
 *
 * 2. Redistributions in binary form must reproduce the above copyright
 *    notice, this list of conditions and the following disclaimer in
 *    the documentation and/or other materials provided with the
 *    distribution.
 *
 * 3. In the absence of prior written permission, the names "ICC" and "The
 *    International Color Consortium" must not be used to imply that the
 *    ICC organization endorses or promotes products derived from this
 *    software.
 *
 *
 * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESSED OR IMPLIED
 * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
 * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
 * DISCLAIMED.  IN NO EVENT SHALL THE INTERNATIONAL COLOR CONSORTIUM OR
 * ITS CONTRIBUTING MEMBERS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
 * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
 * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF
 * USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
 * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
 * OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT
 * OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
 * SUCH DAMAGE.
 * ====================================================================
 *
 * This software consists of voluntary contributions made by many
 * individuals on behalf of the The International Color Consortium. 
 *
 *
 * Membership in the ICC is encouraged when this software is used for
 * commercial purposes. 
 *
 *  
 * For more information on The International Color Consortium, please
 * see <http://www.color.org/>.
 *  
 * 
 */
 
////////////////////////////////////////////////////////////////////// 
// HISTORY:
//
// -Initial implementation by Max Derhak 1-30-2006
//
//////////////////////////////////////////////////////////////////////
 
#ifdef WIN32
#pragma warning( disable: 4786) //disable warning in <list.h>
#endif
 
#include <stdio.h>
#include <math.h>
#include <string.h>
#include <stdlib.h>
#include "IccMpeBasic.h"
#include "IccIO.h"
#include <map>
#include "IccUtil.h"
#include "IccCAM.h"
#include "IccMpeCalc.h"
 
#ifdef USEREFICCMAXNAMESPACE
namespace refIccMAX {
#endif
 
/**
 ******************************************************************************
 * Name: CIccFormulaCurveSegment::CIccFormulaCurveSegment
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccFormulaCurveSegment::CIccFormulaCurveSegment(icFloatNumber start, icFloatNumber end)
{
  m_nReserved = 0;
  m_nReserved2 = 0;
  m_startPoint = start;
  m_endPoint = end;
 
  m_nFunctionType = 0;
  m_nShortcutType = 0;
  m_nParameters = 0;
  m_params = NULL;
}
 
/**
 ******************************************************************************
 * Name: CIccFormulaCurveSegment::CIccFormulaCurveSegment
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccFormulaCurveSegment::CIccFormulaCurveSegment(const CIccFormulaCurveSegment &seg)
{
  m_nReserved = seg.m_nReserved;
  m_nReserved2 = seg.m_nReserved2;
  m_startPoint = seg.m_startPoint;
  m_endPoint = seg.m_endPoint;
 
  m_nFunctionType = seg.m_nFunctionType;
  m_nShortcutType = seg.m_nShortcutType;
  m_nParameters = seg.m_nParameters;
 
  if (seg.m_params) {
    m_params = (icFloatNumber*)malloc(m_nParameters*sizeof(icFloatNumber));
    if (m_params)
      memcpy(m_params, seg.m_params, m_nParameters*sizeof(icFloatNumber));
  }
  else
    m_params = NULL;
}
 
/**
 ******************************************************************************
 * Name: &CIccFormulaCurveSegment::operator=
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccFormulaCurveSegment &CIccFormulaCurveSegment::operator=(const CIccFormulaCurveSegment &seg)
{
  if (m_params)
    free(m_params);
 
  m_nReserved = seg.m_nReserved;
  m_nReserved2 = seg.m_nReserved2;
  m_startPoint = seg.m_startPoint;
  m_endPoint = seg.m_endPoint;
 
  m_nFunctionType = seg.m_nFunctionType;
  m_nShortcutType = seg.m_nShortcutType;
  m_nParameters = seg.m_nParameters;
  if (seg.m_params) {
    m_params = (icFloatNumber*)malloc(m_nParameters*sizeof(icFloatNumber));
    if (m_params)
      memcpy(m_params, seg.m_params, m_nParameters*sizeof(icFloatNumber));
  }
  else
    m_params = NULL;
 
  return (*this);
}
 
/**
 ******************************************************************************
 * Name: CIccFormulaCurveSegment::~CIccFormulaCurveSegment
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccFormulaCurveSegment::~CIccFormulaCurveSegment()
{
  if (m_params) {
    free(m_params);
  }
}
 
/**
 ******************************************************************************
 * Name: CIccFormulaCurveSegment::Describe
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccFormulaCurveSegment::Describe(std::string &sDescription, int nVerboseness)
{
  icChar buf[256];
 
  sDescription += "Segment [";
  if (m_startPoint==icMinFloat32Number)
    sDescription += "-Infinity, ";
  else {
  sprintf(buf, "%.8f, ", m_startPoint);
  sDescription += buf;
  }
  if (m_endPoint==icMaxFloat32Number)
    sDescription += "+Infinity";
  else {
  sprintf(buf, "%.8f", m_endPoint);
  sDescription += buf;
  }
  sprintf(buf, "]\nFunctionType: %04Xh\n", m_nFunctionType);
  sDescription += buf;
 
  switch(m_nFunctionType) {
  case 0x0000:
    if (icIsNear(m_params[1], 0.0) && icIsNear(m_params[2], 0.0))
      sprintf(buf, "Y = %.8f\n\n", m_params[3]);
    else if (icIsNear(m_params[0], 1.0) && icIsNear(m_params[1], 1.0) &&
             icIsNear(m_params[2], 0.0) && icIsNear(m_params[3], 0.0))
      sprintf(buf, "Y = X\n\n");
    else if (icIsNear(m_params[0], 1.0) && icIsNear(m_params[2], 0.0))
      sprintf(buf, "Y = %.8f * X + %.8f\n\n", 
              m_params[1], m_params[3]);
    else
      sprintf(buf, "Y = (%.8f * X + %.8f)^%.4f + %.8f\n\n", 
              m_params[1], m_params[2], m_params[0], m_params[3]);
    sDescription += buf;
    return;
 
  case 0x0001:
    sprintf(buf, "Y = %.8f * log (%.8f * (X ^ %.8f)  + %.8f) + %.8f\n\n",
            m_params[1], m_params[2], m_params[0], m_params[3], m_params[4]);
    sDescription += buf;
    return;
 
  case 0x0002:
    sprintf(buf, "Y = %.8f * (%.8f ^ (%.8f * X + %.8f)) + %.8f\n\n",
            m_params[0], m_params[1], m_params[2], m_params[3], m_params[4]);
    sDescription += buf;
    return;
 
  case 0x0003:
    if (icIsNear(m_params[1],0.0) && icIsNear(m_params[2], 0.0))
      sprintf(buf, "Y = %.8f\n\n", m_params[3]);
    else if (icIsNear(m_params[0], 1.0) && icIsNear(m_params[1], 1.0) &&
             icIsNear(m_params[2], 1.0) && icIsNear(m_params[3], 0.0) &&
             icIsNear(m_params[4], 0.0))
      sprintf(buf, "Y = X\n\n");
    else if (icIsNear(m_params[0], 1.0) && icIsNear(m_params[1], 1.0) &&
             icIsNear(m_params[3], 0.0))
      sprintf(buf, "Y = %.8f * X + %.8f\n\n", 
              m_params[2], m_params[3]);
    else if (icIsNear(m_params[0], 1.0) && icIsNear(m_params[2], 1.0) &&
             icIsNear(m_params[3], 0.0))
      sprintf(buf, "Y = %.8f * X + %.8f\n\n", 
      m_params[1], m_params[3]);
    else
      sprintf(buf, "Y = %8f * (%.8f * X + %.8f)^%.4f + %.8f\n\n", 
              m_params[1], m_params[2], m_params[3], m_params[0], m_params[4]);
    sDescription += buf;
    return;
 
  case 0x0004:
    if (!icIsNear(m_params[0], 1.0))
      sprintf(buf, "Y = %.8f * ln(%.8f * X^%8.f - %.8f) + %.8f\n\n",
        m_params[1], m_params[4], m_params[0], m_params[2], m_params[3]);
    else
      sprintf(buf, "Y = %.8f * ln(%.8f * X - %.8f) + %.8f\n\n",
        m_params[1], m_params[4], m_params[2], m_params[3]);
 
    sDescription += buf;
    return;
 
  case 0x0005:
    if (!icIsNear(m_params[0], 1.0))
      sprintf(buf, "Y = %.8f * exp((%.8f * X^%.8f - %.8f) / %.8f) + %.8f\n\n",
        m_params[5], m_params[4], m_params[0], m_params[3], m_params[1], m_params[2]);
    else
      sprintf(buf, "Y = %.8f * exp((%.8f * X - %.8f) / %.8f) + %.8f\n\n",
        m_params[5], m_params[4], m_params[3], m_params[1], m_params[2]);
 
    sDescription += buf;
    return;
 
  case 0x0006:
    if (!icIsNear(m_params[1], 1.0))
      sprintf(buf, "Y = %.8f * (max[(%.8f * X^%.8f - %.8f), 0] / (%.8f - %.8f * X^%.8f))^%.8f\n\n",
        m_params[5], m_params[6], m_params[1], m_params[2], m_params[3], m_params[4], m_params[1], m_params[0]);
    else
      sprintf(buf, "Y = %.8f * (max[(%.8f * X - %.8f), 0] / (%.8f - %.8f * X))^%.8f\n\n",
        m_params[5], m_params[6], m_params[2], m_params[3], m_params[4], m_params[0]);
    sDescription += buf;
    return;
 
  case 0x0007:
    if (!icIsNear(m_params[1], 1.0))
      sprintf(buf, "Y = %.8f * ((%.8f + %.8f * X^%.8f) / (1 + %.8f * X^%.8f))^%.8f\n\n",
        m_params[5], m_params[2], m_params[3], m_params[1], m_params[4], m_params[1], m_params[0]);
    else
      sprintf(buf, "Y = %.8f * ((%.8f + %.8f * X) / (1 + %.8f * X))^%.8f\n\n",
        m_params[5], m_params[2], m_params[3], m_params[4], m_params[0]);
    sDescription += buf;
    return;
 
  default:
    int i;
    sprintf(buf, "Unknown Function with %d parameters:\n\n", m_nParameters);
    sDescription += buf;
 
    for (i=0; i<m_nParameters; i++) {
      sprintf(buf, "Param[%d] = %.8lf\n\n", i, m_params[i]);
      sDescription += buf;
    }
  }
}
 
 
/**
 ******************************************************************************
 * Name: CIccFormulaCurveSegment::SetFunction
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccFormulaCurveSegment::SetFunction(icUInt16Number functionType, icUInt8Number num_parameters, icFloatNumber *parameters)
{
  if (m_params)
    free(m_params);
 
  if (num_parameters) {
    m_params = (icFloatNumber*)malloc(num_parameters * sizeof(icFloatNumber));
    if (m_params)
      memcpy(m_params, parameters, num_parameters * sizeof(icFloatNumber));
  }
  else
    m_params = NULL;
 
  m_nFunctionType = functionType;
  m_nParameters = num_parameters;
}
 
/**
 ******************************************************************************
 * Name: CIccFormulaCurveSegment::Read
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccFormulaCurveSegment::Read(icUInt32Number size, CIccIO *pIO)
{
  icCurveSegSignature sig;
 
  icUInt32Number headerSize = sizeof(icTagTypeSignature) + 
    sizeof(icUInt32Number) + 
    sizeof(icUInt16Number) + 
    sizeof(icUInt16Number);
 
  if (headerSize > size)
    return false;
 
  if (!pIO) {
    return false;
  }
 
  if (!pIO->Read32(&sig))
    return false;
 
  if (!pIO->Read32(&m_nReserved))
    return false;
 
  if (!pIO->Read16(&m_nFunctionType))
    return false;
 
  if (!pIO->Read16(&m_nReserved2))
    return false;
 
  if (m_params) {
    free(m_params);
  }
 
  switch(m_nFunctionType) {
    case 0x0000:
      m_nParameters = 4;
      break;
 
    case 0x0001:
    case 0x0002:
    case 0x0003:
    case 0x0004:
      m_nParameters = 5;
      break;
 
    case 0x0005:
    case 0x0007:
      m_nParameters = 6;
      break;
 
    case 0x0006:
      m_nParameters = 7;
      break;
 
    default:
      return false;
  }
 
  if (m_nParameters) {
 
    m_params = (icFloatNumber*)malloc(m_nParameters * sizeof(icFloatNumber));
    if (!m_params)
      return false;
 
    if (pIO->ReadFloat32Float(m_params, m_nParameters)!= m_nParameters) {
      return false;
    }
  }
  else
    m_params = NULL;
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccFormulaCurveSegment::Write
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccFormulaCurveSegment::Write(CIccIO *pIO)
{
  icCurveSegSignature sig = GetType();
 
  if (!pIO)
    return false;
 
  if (!pIO->Write32(&sig))
    return false;
 
  if (!pIO->Write32(&m_nReserved))
    return false;
 
  if (!pIO->Write16(&m_nFunctionType))
    return false;
 
  if (!pIO->Write16(&m_nReserved2))
    return false;
 
  switch(m_nFunctionType) {
    case 0x0000:
      if (m_nParameters!=4)
        return false;
      break;
 
    case 0x0001:
    case 0x0002:
    case 0x0003:
    case 0x0004:
      if (m_nParameters!=5)
        return false;
      break;
 
    case 0x0005:
    case 0x0007:
      if (m_nParameters != 6)
        return false;
      break;
 
    case 0x0006:
      if (m_nParameters != 7)
        return false;
      break;
  }
 
  if (m_nParameters) {
    if (pIO->WriteFloat32Float(m_params, m_nParameters)!=m_nParameters)
      return false;
  }
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccFormulaCurveSegment::Begin
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccFormulaCurveSegment::Begin(CIccCurveSegment *pPrevSeg = NULL)
{
  switch (m_nFunctionType) {
  case 0x0000:
    if (!m_params || m_nParameters<4)
      return false;
 
    if (icIsNear(m_params[0], 1.0)) { //don't apply gamma
      if (icIsNear(m_params[2], 0.0) && icIsNear(m_params[3], 0.0))
        m_nShortcutType = 1;
      else if (icIsNear(m_params[2], 0.0))
        m_nShortcutType = 2;
      else if (icIsNear(m_params[3], 0.0))
        m_nShortcutType = 3;
      else
        m_nShortcutType = 4;
    }
    else {
      m_nShortcutType = 0;
    }
 
    return true;
 
  case 0x0001:
    if (!m_params || m_nParameters<5)
      return false;
 
    return true;
 
  case 0x0002:
  case 0x0003:
    if (!m_params || m_nParameters < 5)
      return false;
 
  case 0x0004:
    if (!m_params || m_nParameters < 5)
      return false;
 
    if (icIsNear(m_params[0], 1.0)) { //don't apply gamma
      m_nShortcutType = 1;
    }
    else {
      m_nShortcutType = 0;
    }
    return true;
 
  case 0x0005:
    if (!m_params || m_nParameters < 6)
      return false;
 
    if (icIsNear(m_params[0], 1.0)) { //don't apply gamma
      m_nShortcutType = 1;
    }
    else {
      m_nShortcutType = 0;
    }
    return true;
 
  case 0x0006:
    if (!m_params || m_nParameters < 7)
      return false;
 
    if (icIsNear(m_params[1], 1.0)) { //don't apply gamma
      m_nShortcutType = 1;
    }
    else {
      m_nShortcutType = 0;
    }
    return true;
 
  case 0x0007:
    if (!m_params || m_nParameters < 6)
      return false;
 
    if (icIsNear(m_params[1], 1.0)) { //don't apply gamma
      m_nShortcutType = 1;
    }
    else {
      m_nShortcutType = 0;
    }
    return true;
 
 
  default:
    return false;
  }
 
  return true;
}
 
static icFloatNumber clipPow(double v, double g)
{
  if (v <= 0)
    return 0;
  return (icFloatNumber)pow(v, g);
}
 
/**
 ******************************************************************************
 * Name: CIccFormulaCurveSegment::Apply
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
icFloatNumber CIccFormulaCurveSegment::Apply(icFloatNumber v) const
{
  switch (m_nFunctionType) {
  case 0x0000:
    //Y = (a * X + b) ^ g  + c        : g a b c
    switch (m_nShortcutType) {
    case 0:
    default:
      return (icFloatNumber)(pow(m_params[1] * v + m_params[2], m_params[0]) + m_params[3]);
    case 1:
      return (m_params[1] * v);
    case 2:
      return (m_params[1] * v + m_params[3]);
    case 3:
      return (m_params[1] * v + m_params[2]);
    case 4:
      return (m_params[1] * v + m_params[2] + m_params[3]);
    }
 
  case 0x0001:
    // Y = a * log (b * X^g + c) + d      : g a b c d
    return (icFloatNumber)(m_params[1] * log10(m_params[2] * pow(v, m_params[0]) + m_params[3]) + m_params[4]);
 
  case 0x0002:
    //Y = a * b^(c*X+d) + e               : a b c d e
    return (icFloatNumber)(m_params[0] * pow(m_params[1], m_params[2] * v + m_params[3]) + m_params[4]);
 
  case 0x0003:
    //Y = a * (b * X + c) ^ g  + d        : g a b c d
    return (icFloatNumber)(m_params[1] * clipPow(m_params[2] * v + m_params[3], m_params[0]) + m_params[4]);
 
  case 0x0004:
    //Y = a * ln(d * X^g - b) + c         : g a b c d
    if (m_nShortcutType != 1)
      return (icFloatNumber)(m_params[1] * log(m_params[4] * pow(v, m_params[0]) - m_params[2]) + m_params[3]);
    else
      return (icFloatNumber)(m_params[1] * log(m_params[4] * v - m_params[2]) + m_params[3]);
 
  case 0x0005:
    //Y = e * exp((d * X^g - c) / a) + b  : g a b c d e
    if (m_nShortcutType != 1)
      return (icFloatNumber)(m_params[5] * exp((m_params[4] * pow(v, m_params[0]) - m_params[3]) / m_params[1]) + m_params[2]);
    else
      return (icFloatNumber)(m_params[5] * exp((m_params[4] * v - m_params[3]) / m_params[1]) + m_params[2]);
 
  case 0x0006:
    //Y = d * (max(e * X^g - a, 0)/(b - c * X^g))^w  : w g a b c d e
    if (m_nShortcutType!=1) {
      return (icFloatNumber)(m_params[5] * pow(icMax((icFloatNumber)(m_params[6] * pow(v, m_params[1]) - m_params[2]), 0.0f) /
                                               (m_params[3] - m_params[4] * pow(v, m_params[1])), m_params[0]));
    }
    else { //m_nShortcutType == 1
      return (icFloatNumber)(m_params[5] * pow(icMax((icFloatNumber)(m_params[6] * v - m_params[2]), 0.0f) /
                                               (m_params[3] - m_params[4] * v), m_params[0]));
    }
 
  case 0x0007:
    //Y = d * ((a + b * X^g)/(1 + c * X^g)) ^ w     : w g a b c d
    if (m_nShortcutType != 1) {
      return (icFloatNumber)(m_params[5] * pow((m_params[2] + m_params[3] * pow(v, m_params[1])) /
                                               (1.0 + m_params[4] * pow(v, m_params[1])), m_params[0]));
    }
    else { //m_nShortcutType == 1
      return (icFloatNumber)(m_params[5] * pow((m_params[2] + m_params[3] * v) /
                                               (1.0 + m_params[4] * v), m_params[0]));
    }
  }
 
  //Shouldn't get here!
  return v;
}
 
/**
 ******************************************************************************
 * Name: CIccFormulaCurveSegment::Validate
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
icValidateStatus CIccFormulaCurveSegment::Validate(std::string sigPath, std::string &sReport, const CIccTagMultiProcessElement* pMPE/*=NULL*/, const CIccProfile *pProfile/*=NULL*/) const
{
  CIccInfo Info;
  std::string sSigPathName = Info.GetSigPathName(sigPath);
 
  icValidateStatus rv = icValidateOK;
  if (m_nReserved || m_nReserved2) {
    sReport += icMsgValidateWarning;
    sReport += sSigPathName;
    sReport += " formula curve has non zero reserved data.\n";
    rv = icMaxStatus(rv, icValidateWarning);
  }
 
  switch (m_nFunctionType) {
  case 0x0000:
    if (!m_params || m_nParameters<4) {
      sReport += icMsgValidateCriticalError;
      sReport += sSigPathName;
      sReport += " formula curve has Invalid formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateCriticalError);
    }
    else if (m_nParameters > 4) {
      sReport += icMsgValidateWarning;
      sReport += sSigPathName;
      sReport += " formula curve has too many formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateWarning);
    }
    break;
 
  case 0x0001:
    if (!m_params || m_nParameters<5) {
      sReport += icMsgValidateCriticalError;
      sReport += sSigPathName;
      sReport += " formula curve has Invalid formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateCriticalError);
    }
    else if (m_nParameters > 5) {
      sReport += icMsgValidateWarning;
      sReport += sSigPathName;
      sReport += " formula curve has too many formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateWarning);
    }
    break;
 
  case 0x0002:
    if (!m_params || m_nParameters<5) {
      sReport += icMsgValidateCriticalError;
      sReport += sSigPathName;
      sReport += " formula curve has Invalid formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateCriticalError);
    }
    else if (m_nParameters > 5) {
      sReport += icMsgValidateWarning;
      sReport += sSigPathName;
      sReport += " formula curve has too many formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateWarning);
    }
    break;
 
  case 0x0003:
    if (!m_params || m_nParameters<5) {
      sReport += icMsgValidateCriticalError;
      sReport += sSigPathName;
      sReport += " formula curve has Invalid formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateCriticalError);
    }
    else if (m_nParameters > 5) {
      sReport += icMsgValidateWarning;
      sReport += sSigPathName;
      sReport += " formula curve has too many formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateWarning);
    }
    break;
 
  case 0x0004:
    if (!m_params || m_nParameters < 5) {
      sReport += icMsgValidateCriticalError;
      sReport += sSigPathName;
      sReport += " formula curve has Invalid formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateCriticalError);
    }
    else if (m_nParameters > 5) {
      sReport += icMsgValidateWarning;
      sReport += sSigPathName;
      sReport += " formula curve has too many formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateWarning);
    }
    break;
 
  case 0x0005:
    if (!m_params || m_nParameters < 6) {
      sReport += icMsgValidateCriticalError;
      sReport += sSigPathName;
      sReport += " formula curve has Invalid formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateCriticalError);
    }
    else if (m_nParameters > 6) {
      sReport += icMsgValidateWarning;
      sReport += sSigPathName;
      sReport += " formula curve has too many formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateWarning);
    }
    break;
 
  case 0x0006:
    if (!m_params || m_nParameters < 6) {
      sReport += icMsgValidateCriticalError;
      sReport += sSigPathName;
      sReport += " formula curve has Invalid formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateCriticalError);
    }
    else if (m_nParameters > 6) {
      sReport += icMsgValidateWarning;
      sReport += sSigPathName;
      sReport += " formula curve has too many formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateWarning);
    }
    break;
 
  case 0x0007:
    if (!m_params || m_nParameters < 6) {
      sReport += icMsgValidateCriticalError;
      sReport += sSigPathName;
      sReport += " formula curve has Invalid formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateCriticalError);
    }
    else if (m_nParameters > 6) {
      sReport += icMsgValidateWarning;
      sReport += sSigPathName;
      sReport += " formula curve has too many formulaCurveSegment parameters.\n";
      rv = icMaxStatus(rv, icValidateWarning);
    }
    break;
 
  default:
    {
      icChar buf[128];
      sReport += icMsgValidateCriticalError;
      sReport += sSigPathName;
      sprintf(buf, " formula curve uses unknown formulaCurveSegment function type %d\n", m_nFunctionType);
      sReport += buf;
      rv = icMaxStatus(rv, icValidateCriticalError);
    }
  }
 
  return rv;
}
 
/**
 ******************************************************************************
 * Name: CIccSampledCurveSegment::CIccSampledCurveSegment
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccSampledCurveSegment::CIccSampledCurveSegment(icFloatNumber start, icFloatNumber end)
{
  m_nReserved = 0;
  m_startPoint = start;
  m_endPoint = end;
  m_range = end - start;
  m_nCount = 0;
  m_pSamples = 0;
  m_last = 0;
}
 
/**
 ******************************************************************************
 * Name: CIccSampledCurveSegment::CIccSampledCurveSegment
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccSampledCurveSegment::CIccSampledCurveSegment(const CIccSampledCurveSegment &curve)
{
  m_nReserved = curve.m_nReserved;
  m_startPoint = curve.m_startPoint;
  m_endPoint = curve.m_endPoint;
  m_range = m_endPoint - m_startPoint;
  m_nCount = curve.m_nCount;
  m_last = curve.m_last;
 
  if (m_nCount) {
    m_pSamples = (icFloatNumber*)malloc(m_nCount * sizeof(icFloatNumber));
    if (m_pSamples)
      memcpy(m_pSamples, curve.m_pSamples, m_nCount * sizeof(icFloatNumber));
    else
      m_nCount = 0;
  }
  else{
    m_pSamples = NULL;
  }
}
 
/**
 ******************************************************************************
 * Name: &CIccSampledCurveSegment::operator=
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccSampledCurveSegment &CIccSampledCurveSegment::operator=(const CIccSampledCurveSegment &curve)
{
  if (m_pSamples)
    free(m_pSamples);
 
  m_nReserved = curve.m_nReserved;
  m_startPoint = curve.m_startPoint;
  m_endPoint = curve.m_endPoint;
  m_nCount = curve.m_nCount;
 
  if (m_nCount) {
    m_pSamples = (icFloatNumber*)malloc(m_nCount * sizeof(icFloatNumber));
    if (m_pSamples)
      memcpy(m_pSamples, curve.m_pSamples, m_nCount * sizeof(icFloatNumber));
    else
      m_nCount = 0;
  }
  else {
    m_pSamples = NULL;
  }
  return (*this);
}
 
/**
 ******************************************************************************
 * Name: CIccSampledCurveSegment::~CIccSampledCurveSegment
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccSampledCurveSegment::~CIccSampledCurveSegment()
{
  if (m_pSamples)
    free(m_pSamples);
}
 
/**
 ******************************************************************************
 * Name: CIccSampledCurveSegment::SetSize
 * 
 * Purpose: 
 *  Sets size of sampled lookup table.  Previous data (if exists) is lost.
 * 
 * Args: 
 *  nCount = number of elements in lut (must be >= 2).  Note: the m_pSample[0] is
 *    initialized from the the previous segment.  It is not saved as part
 *    of Write(), or loaded as part of Read().  It will be initialized during
 *    the call to Begin(),  The actual count of elements written to the file
 *    will be nCount-1
 *  bZeroAlloc = flag to decide if memory should be set to zero.
 * 
 * Return: 
 *  true if allocation successful.
 ******************************************************************************/
bool CIccSampledCurveSegment::SetSize(icUInt32Number nCount, bool bZeroAlloc/*=true*/)
{
  if (!nCount) {
    if (m_pSamples)
      free(m_pSamples);
    m_pSamples = NULL;
    m_nCount = nCount;
    return true;
  }
 
  if (m_pSamples) {
    free(m_pSamples);
  }
 
  if (bZeroAlloc)
    m_pSamples = (icFloatNumber*)calloc(nCount, sizeof(icFloatNumber));
  else
    m_pSamples = (icFloatNumber*)malloc(nCount * sizeof(icFloatNumber));
 
  if (m_pSamples)
    m_nCount = nCount;
  else
    m_nCount = 0;
 
  return (m_pSamples != NULL);
}
 
/**
 ******************************************************************************
 * Name: CIccSampledCurveSegment::Describe
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccSampledCurveSegment::Describe(std::string &sDescription, int nVerboseness)
{
  icChar buf[128];
 
  if (m_nCount<2) {
    sDescription += "Empty Segment [";
    if (m_startPoint==icMinFloat32Number)
      sDescription += "-Infinity, ";
    else {
    sprintf(buf, "%.8f, ", m_startPoint);
    sDescription += buf;
    }
    if (m_endPoint==icMaxFloat32Number)
      sDescription += "+Infinity";
    else {
    sprintf(buf, "%.8f", m_endPoint);
    sDescription += buf;
    }
 
    sprintf(buf, "]\n");
    sDescription += buf;
  }
  else {
    sDescription += "Sampled Segment [";
    if (m_startPoint==icMinFloat32Number)
      sDescription += "-Infinity, ";
    else {
    sprintf(buf, "%.8f, ", m_startPoint);
    sDescription += buf;
    }
    if (m_endPoint==icMaxFloat32Number)
      sDescription += "+Infinity";
    else {
    sprintf(buf, "%.8f", m_endPoint);
    sDescription += buf;
    }
    sprintf(buf, "]\n");
    sDescription += buf;
    sDescription += "IN  OUT\n";
 
    icUInt32Number i;
 
    icFloatNumber range = m_endPoint - m_startPoint;
    icFloatNumber last = (icFloatNumber)(m_nCount-1);
 
    for (i=1; i<m_nCount; i++) {
      sprintf(buf, "%.8f %.8f\n", m_startPoint + (icFloatNumber)i*range/last, m_pSamples[i]);
      sDescription += buf;
    }
  }
  sDescription += "\n";
}
 
/**
 ******************************************************************************
 * Name: CIccSampledCurveSegment::Read
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccSampledCurveSegment::Read(icUInt32Number size, CIccIO *pIO)
{
  icCurveSegSignature sig;
 
  icUInt32Number headerSize = sizeof(icTagTypeSignature) + 
    sizeof(icUInt32Number) + 
    sizeof(icUInt32Number);
 
  if (headerSize > size)
    return false;
 
  if (!pIO) {
    return false;
  }
 
  if (!pIO->Read32(&sig))
    return false;
 
  if (!pIO->Read32(&m_nReserved))
    return false;
 
  if (!pIO->Read32(&m_nCount))
    return false;
 
  if (m_nCount > size - headerSize || m_nCount * sizeof(icFloat32Number) > size - headerSize)
    return false;
 
  //Reserve room for first point who's value comes from previous segment
  m_nCount++;
 
  if (!SetSize(m_nCount, false))
    return false;
 
  if (m_nCount) {
    if (pIO->ReadFloat32Float(m_pSamples+1, m_nCount-1)!=(icInt32Number)(m_nCount-1))
      return false;
  }
 
  //Initialize first point with zero.  Properly initialized during Begin()
  m_pSamples[0] = 0;
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccSampledCurveSegment::Write
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccSampledCurveSegment::Write(CIccIO *pIO)
{
  icCurveSegSignature sig = GetType();
 
  if (!pIO)
    return false;
 
  if (!pIO->Write32(&sig))
    return false;
 
  if (!pIO->Write32(&m_nReserved))
    return false;
 
  icUInt32Number nCount;
 
  if (m_nCount)
    nCount = m_nCount -1;
  else
    nCount = 0;
 
  if (!pIO->Write32(&nCount))
    return false;
 
  //First point in samples is ONLY for interpolation (not saved)
  if (nCount) {
    if (pIO->WriteFloat32Float(m_pSamples+1, nCount)!=(icInt32Number)nCount)
      return false;
  }
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccSampledCurveSegment::Begin
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccSampledCurveSegment::Begin(CIccCurveSegment *pPrevSeg = NULL)
{
  if (m_nCount<2)
    return false;
 
  m_range = m_endPoint - m_startPoint;
  m_last = (icFloatNumber)(m_nCount - 1);
 
  if (m_endPoint-m_startPoint == 0.0)
    return false;
 
  if (!pPrevSeg)
    return false;
 
  //Set up interpolation from Application of last segment
  m_pSamples[0] = pPrevSeg->Apply(m_startPoint);
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccSampledCurveSegment::Apply
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
icFloatNumber CIccSampledCurveSegment::Apply(icFloatNumber v) const
{
  if (v<m_startPoint)
    v=m_startPoint;
  else if (v>m_endPoint)
    v=m_endPoint;
 
  icFloatNumber pos = (v-m_startPoint)/m_range * m_last;
  icUInt32Number index = (icUInt32Number) pos;
  icFloatNumber remainder = pos - (icFloatNumber)index;
 
  if (remainder==0.0)
    return m_pSamples[index];
 
  return (icFloatNumber)((1.0-remainder)*m_pSamples[index] + remainder*m_pSamples[index+1]);
}
 
/**
 ******************************************************************************
 * Name: CIccSampledCurveSegment::Validate
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
icValidateStatus CIccSampledCurveSegment::Validate(std::string sigPath, std::string &sReport, const CIccTagMultiProcessElement* pMPE/*=NULL*/, const CIccProfile *pProfile/*=NULL*/) const
{
  CIccInfo Info;
  std::string sSigPathName = Info.GetSigPathName(sigPath);
 
  icValidateStatus rv = icValidateOK;
  if (m_nReserved) {
    sReport += icMsgValidateWarning;
    sReport += sSigPathName;
    sReport += " sampled curve has non zero reserved data.\n";
    rv = icValidateWarning;
  }
 
  if (m_nCount<1) {
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " sampled curve has too few sample points.\n";
    rv = icMaxStatus(rv, icValidateCriticalError);
  }
  else if (m_endPoint-m_startPoint == 0.0) {
    sReport += icMsgValidateWarning;
    sReport += sSigPathName;
    sReport += " sampled curve has a range of zero.\n";
    rv = icMaxStatus(rv, icValidateWarning);
  }
 
  return rv;
}
 
 
 
/**
******************************************************************************
* Name: CIccSingleSampledCurve::CIccSingleSampledCurve
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
CIccSingleSampledCurve::CIccSingleSampledCurve(icFloatNumber first, icFloatNumber last)
{
  m_nReserved = 0;
  m_nCount = 0;
  m_pSamples = 0;
  m_last = 0;
 
  m_storageType = icValueTypeFloat32;
  m_extensionType = icClipSingleSampledCurve;
 
  if (first < last) {
    m_firstEntry = first;
    m_lastEntry = last;
  }
  else {
    m_firstEntry = last;
    m_lastEntry = first;
  }
 
  m_range = m_lastEntry - m_firstEntry;
  m_loIntercept = 0;
  m_loSlope = 0;
  m_hiIntercept = 1.0;
  m_hiSlope = 0;
}
 
/**
******************************************************************************
* Name: CIccSingleSampledCurve::CIccSingleSampledCurve
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
CIccSingleSampledCurve::CIccSingleSampledCurve(const CIccSingleSampledCurve &curve)
{
  m_nReserved = curve.m_nReserved;
  m_nCount = curve.m_nCount;
 
  m_storageType = curve.m_storageType;
  m_extensionType = curve.m_extensionType;
 
  if (m_nCount) {
    m_pSamples = (icFloatNumber*)malloc(m_nCount * sizeof(icFloatNumber));
    if (m_pSamples)
      memcpy(m_pSamples, curve.m_pSamples, m_nCount * sizeof(icFloatNumber));
    else
      m_nCount = 0;
  }
  else{
    m_pSamples = NULL;
  }
 
  m_firstEntry = curve.m_firstEntry;
  m_lastEntry = curve.m_lastEntry;
  m_range = m_lastEntry - m_firstEntry;
 
  m_loIntercept = curve.m_loIntercept;
  m_loSlope = curve.m_loSlope;
  m_hiIntercept = curve.m_hiIntercept;
  m_hiSlope = curve.m_hiSlope;
}
 
/**
******************************************************************************
* Name: &CIccSingleSampledCurve::operator=
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
CIccSingleSampledCurve &CIccSingleSampledCurve::operator=(const CIccSingleSampledCurve &curve)
{
  if (m_pSamples)
    free(m_pSamples);
 
  m_nReserved = curve.m_nReserved;
  m_nCount = curve.m_nCount;
 
  m_storageType = curve.m_storageType;
  m_extensionType = curve.m_extensionType;
 
  if (m_nCount) {
    m_pSamples = (icFloatNumber*)malloc(m_nCount * sizeof(icFloatNumber));
    if (m_pSamples)
      memcpy(m_pSamples, curve.m_pSamples, m_nCount * sizeof(icFloatNumber));
    else
      m_nCount = 0;
  }
  else {
    m_pSamples = NULL;
  }
 
  m_firstEntry = curve.m_firstEntry;
  m_lastEntry = curve.m_lastEntry;
 
  m_loIntercept = curve.m_loIntercept;
  m_loSlope = curve.m_loSlope;
  m_hiIntercept = curve.m_hiIntercept;
  m_hiSlope = curve.m_hiSlope;
  return (*this);
}
 
/**
******************************************************************************
* Name: CIccSingleSampledCurve::~CIccSingleSampledCurve
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
CIccSingleSampledCurve::~CIccSingleSampledCurve()
{
  if (m_pSamples)
    free(m_pSamples);
}
 
/**
******************************************************************************
* Name: CIccSingleSampledCurve::SetRange
* 
* Purpose: 
*  Sets range of sampled lookup table.
* 
* Args: 
*  first = first entry input position
*  last = last entry input position
******************************************************************************/
void CIccSingleSampledCurve::SetRange(icFloatNumber first/* =0.0f */, icFloatNumber last/* =1.0f */)
{
  if (first < last) {
    m_firstEntry = first;
    m_lastEntry = last;
  }
  else {
    m_firstEntry = last;
    m_lastEntry = first;
  }
}
 
 
/**
******************************************************************************
* Name: CIccSingleSampledCurve::SetExtensionType
*
* Purpose:
*  Sets extension type of the single sampled curve
*
* Args:
*  nExtensionType is type of extension to use
*
* Return:
*  true if valid extension type, false otherwise
******************************************************************************/
bool CIccSingleSampledCurve::SetExtensionType(icUInt16Number nExtensionType)
{
  m_extensionType = nExtensionType;
 
  switch (nExtensionType) {
    case icClipSingleSampledCurve:
    case icExtendSingleSampledCurve:
      return true;
  }
  return false;
}
 
/**
******************************************************************************
* Name: CIccSingleSampledCurve::SetStorageType
*
* Purpose:
*  Sets storate type of the data stored in the single sampled curve
*
* Args:
*  nStorageType is type of data to use
*
* Return:
*  true if valid storage type, false otherwise
******************************************************************************/
bool CIccSingleSampledCurve::SetStorageType(icUInt16Number nStorateType)
{
  m_storageType = nStorateType;
 
  switch (nStorateType) {
    case icValueTypeUInt8:
    case icValueTypeUInt16:
    case icValueTypeFloat16:
    case icValueTypeFloat32:
      return true;
  }
  return false;
}
 
 
/**
******************************************************************************
* Name: CIccSingleSampledCurve::SetSize
* 
* Purpose: 
*  Sets size of sampled lookup table.  Previous data (if exists) is lost.
* 
* Args: 
*  nCount = number of elements in lut (must be >= 2).  
*  bZeroAlloc = flag to decide if memory should be set to zero.
* 
* Return: 
*  true if allocation successful.
******************************************************************************/
bool CIccSingleSampledCurve::SetSize(icUInt32Number nCount, bool bZeroAlloc/*=true*/)
{
  if (!nCount) {
    if (m_pSamples)
      free(m_pSamples);
    m_pSamples = NULL;
    m_nCount = nCount;
    return true;
  }
 
  if (m_pSamples) {
    free(m_pSamples);
  }
 
  if (bZeroAlloc)
    m_pSamples = (icFloatNumber*)calloc(nCount, sizeof(icFloatNumber));
  else
    m_pSamples = (icFloatNumber*)malloc(nCount * sizeof(icFloatNumber));
 
  if (m_pSamples)
    m_nCount = nCount;
  else
    m_nCount = 0;
 
  return (m_pSamples != NULL);
}
 
/**
******************************************************************************
* Name: CIccSingleSampledCurve::Describe
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
void CIccSingleSampledCurve::Describe(std::string &sDescription, int nVerboseness)
{
  icChar buf[128];
 
  if (m_nCount<2) {
    sDescription += "Empty Single Sampled Curve [";
    sprintf(buf, "%.8f, ", m_firstEntry);
    sDescription += buf;
 
    sprintf(buf, "%.8f", m_lastEntry);
    sDescription += buf;
 
    sprintf(buf, "]\n");
    sDescription += buf;
  }
  else {
    sDescription += "Single Sampled Curve [";
 
    sprintf(buf, "%.8f, ", m_firstEntry);
    sDescription += buf;
 
    sprintf(buf, "%.8f", m_lastEntry);
    sDescription += buf;
 
    sprintf(buf, "]\n");
    sDescription += buf;
 
    switch(m_storageType) {
      case icValueTypeUInt8:
        sDescription +="UInt8Number Array\n";
        break;
      case icValueTypeUInt16:
        sDescription +="UInt16Number Array\n";
        break;
      case icValueTypeFloat16:
        sDescription +="Float16Number Array\n";
        break;
      case icValueTypeFloat32:
        sDescription +="Float32Number Array\n";
        break;
      default:
        sDescription +="Unknown Data Type\n";
        break;
    }
 
    switch(m_extensionType) {
      case icClipSingleSampledCurve:
        sDescription += "Clip End Points\n";
        break;
      case icExtendSingleSampledCurve:
        sDescription += "Extend End Points\n";
        break;
      default:
        sDescription += "Unknown extension handling\n";
        break;
    }
    sDescription += "IN  OUT\n";
 
    icUInt32Number i;
 
    icFloatNumber range = m_lastEntry - m_firstEntry;
    icFloatNumber last = (icFloatNumber)(m_nCount-1);
 
    for (i=0; i<m_nCount; i++) {
      sprintf(buf, "%.8f %.8f\n", m_firstEntry + (icFloatNumber)i*range/last, m_pSamples[i]);
      sDescription += buf;
    }
  }
  sDescription += "\n";
}
 
/**
******************************************************************************
* Name: CIccSingleSampledCurve::Read
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
bool CIccSingleSampledCurve::Read(icUInt32Number size, CIccIO *pIO)
{
  icCurveSegSignature sig;
 
  icUInt32Number headerSize = sizeof(icTagTypeSignature) + 
    sizeof(icUInt32Number) + 
    sizeof(icUInt32Number) +
    sizeof(icFloat32Number) +
    sizeof(icFloat32Number) + 
    sizeof(icUInt16Number) + 
    sizeof(icUInt16Number);
 
  if (headerSize > size)
    return false;
 
  if (!pIO) {
    return false;
  }
 
  if (!pIO->Read32(&sig) || sig!=GetType())
    return false;
 
  if (!pIO->Read32(&m_nReserved))
    return false;
 
  if (!pIO->Read32(&m_nCount))
    return false;
 
  if (!pIO->ReadFloat32Float(&m_firstEntry))
    return false;
 
  if (!pIO->ReadFloat32Float(&m_lastEntry))
    return false;
 
  if (!pIO->Read16(&m_extensionType))
    return false;
 
  if (!pIO->Read16(&m_storageType))
    return false;
 
  if (!SetSize(m_nCount, false))
    return false;
 
  if (m_nCount > size - headerSize)
    return false;
 
  if (m_nCount) {
    switch(m_storageType) {
      case icValueTypeUInt8:
        if (m_nCount * sizeof(icUInt8Number) > size - headerSize)
          return false;
 
        if (pIO->ReadUInt8Float(m_pSamples, m_nCount)!=(icInt32Number)(m_nCount))
          return false;
        break;
      case icValueTypeUInt16:
        if (m_nCount * 2 > size - headerSize)
          return false;
 
        if (pIO->ReadUInt16Float(m_pSamples, m_nCount)!=(icInt32Number)(m_nCount))
          return false;
        break;
      case icValueTypeFloat16:
        if (m_nCount * 2 > size - headerSize)
          return false;
 
        if (pIO->ReadFloat16Float(m_pSamples, m_nCount)!=(icInt32Number)(m_nCount))
          return false;
        break;
      case icValueTypeFloat32:
        if (m_nCount * sizeof(icFloat32Number) > size - headerSize)
          return false;
 
        if (pIO->ReadFloat32Float(m_pSamples, m_nCount)!=(icInt32Number)(m_nCount))
          return false;
        break;
      default:
        return false;
    }
  }
 
  return true;
}
 
/**
******************************************************************************
* Name: CIccSingleSampledCurve::Write
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
bool CIccSingleSampledCurve::Write(CIccIO *pIO)
{
  icCurveElemSignature sig = GetType();
 
  if (!pIO)
    return false;
 
  if (!pIO->Write32(&sig))
    return false;
 
  if (!pIO->Write32(&m_nReserved))
    return false;
 
  if (!pIO->Write32(&m_nCount))
    return false;
 
  if (!pIO->WriteFloat32Float(&m_firstEntry))
    return false;
 
  if (!pIO->WriteFloat32Float(&m_lastEntry))
    return false;
 
  if (!pIO->Write16(&m_extensionType))
    return false;
 
  if (!pIO->Write16(&m_storageType))
    return false;
 
  if (m_nCount) {
 
    switch(m_storageType) {
      case icValueTypeUInt8:
        if (pIO->WriteUInt8Float(m_pSamples, m_nCount)!=(icInt32Number)(m_nCount))
          return false;
        break;
      case icValueTypeUInt16:
        if (pIO->WriteUInt16Float(m_pSamples, m_nCount)!=(icInt32Number)(m_nCount))
          return false;
        break;
      case icValueTypeFloat16:
        if (pIO->WriteFloat16Float(m_pSamples, m_nCount)!=(icInt32Number)(m_nCount))
          return false;
        break;
      case icValueTypeFloat32:
        if (pIO->WriteFloat32Float(m_pSamples, m_nCount)!=(icInt32Number)(m_nCount))
          return false;
        break;
      default:
        return false;
    }
  }
 
  return true;
}
 
/**
******************************************************************************
* Name: CIccSingleSampledCurve::Begin
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
bool CIccSingleSampledCurve::Begin(icElemInterp nInterp, CIccTagMultiProcessElement *pMPE)
{
  if (m_nCount<2)
    return false;
 
  m_range = m_lastEntry - m_firstEntry;
  m_last = (icFloatNumber)(m_nCount - 1);
  icFloatNumber stepSize = m_range / m_last;
 
  if (m_range == 0.0)
    return false;
 
  switch(m_extensionType) {
    case icClipSingleSampledCurve:
      m_loSlope=0.0;
      m_loIntercept = m_pSamples[0];
 
      m_hiSlope=0.0;
      m_hiIntercept = m_pSamples[m_nCount-1];
      break;
 
    case icExtendSingleSampledCurve:
      m_loSlope = (m_pSamples[1]-m_pSamples[0])/stepSize;
      m_loIntercept = m_pSamples[0] - m_loSlope*m_firstEntry;
 
      m_hiSlope = (m_pSamples[m_nCount-1] - m_pSamples[m_nCount-2])/stepSize;
      m_hiIntercept = m_pSamples[m_nCount-1] - m_hiSlope*m_lastEntry;
      break;
 
    default:
      return false;
  }
 
  return true;
}
 
/**
******************************************************************************
* Name: CIccSingleSampledCurve::Apply
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
icFloatNumber CIccSingleSampledCurve::Apply(icFloatNumber v) const
{
  if (v<m_firstEntry) {
    return m_loSlope * v + m_loIntercept;;
  }
  else if (v>m_lastEntry) {
    return m_hiSlope * v + m_hiIntercept;
  }
 
  icFloatNumber pos = (v-m_firstEntry)/m_range * m_last;
  icUInt32Number index = (icUInt32Number) pos;
  icFloatNumber remainder = pos - (icFloatNumber)index;
 
  if (remainder==0.0)
    return m_pSamples[index];
 
  return (icFloatNumber)((1.0-remainder)*m_pSamples[index] + remainder*m_pSamples[index+1]);
}
 
/**
******************************************************************************
* Name: CIccSingleSampledCurve::Validate
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
icValidateStatus CIccSingleSampledCurve::Validate(std::string sigPath, std::string &sReport, const CIccTagMultiProcessElement* pMPE/*=NULL*/, const CIccProfile *pProfile/*=NULL*/) const
{
  CIccInfo Info;
  std::string sSigPathName = Info.GetSigPathName(sigPath);
 
  icValidateStatus rv = icValidateOK;
  if (m_nReserved) {
    sReport += icMsgValidateWarning;
    sReport += sSigPathName;
    sReport += " single sampled curve has non zero reserved data.\n";
    rv = icValidateWarning;
  }
 
  if (m_extensionType > icMaxSingleSampledCurveType) {
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " single sampled curve has unknown extension type\n";
    rv = icMaxStatus(rv, icValidateCriticalError);
  }
 
  if (m_storageType > icMaxValueType) {
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " single sampled curve uses unknown value type\n";
    rv = icMaxStatus(rv, icValidateCriticalError);
  }
 
  if (m_nCount<2) {
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " single sampled curve has too few sample points.\n";
    rv = icMaxStatus(rv, icValidateCriticalError);
  }
  
  if (m_lastEntry-m_firstEntry <= 0.0) {
    sReport += icMsgValidateWarning;
    sReport += sSigPathName;
    sReport += " single sampled curve has an invalid sample range.\r\n";
    rv = icMaxStatus(rv, icValidateWarning);
  }
 
  return rv;
}
 
 
/**
******************************************************************************
* Name: CIccSampledCalculatorCurve::CIccSampledCalculatorCurve
*
* Purpose:
*
* Args:
*
* Return:
******************************************************************************/
CIccSampledCalculatorCurve::CIccSampledCalculatorCurve(icFloatNumber first, icFloatNumber last)
{
  m_nReserved = 0;
  m_nReserved2 = 0;
 
  m_pCalc = 0;
 
  m_nCount = 0;
  m_pSamples = 0;
 
  m_extensionType = icClipSingleSampledCurve;
 
  if (first < last) {
    m_firstEntry = first;
    m_lastEntry = last;
  }
  else {
    m_firstEntry = last;
    m_lastEntry = first;
  }
 
  m_loIntercept = 0;
  m_loSlope = 0;
  m_hiIntercept = 1.0;
  m_hiSlope = 0;
}
 
/**
******************************************************************************
* Name: CIccSampledCalculatorCurve::CIccSampledCalculatorCurve
*
* Purpose:
*
* Args:
*
* Return:
******************************************************************************/
CIccSampledCalculatorCurve::CIccSampledCalculatorCurve(const CIccSampledCalculatorCurve &curve)
{
  m_nReserved = curve.m_nReserved;
  m_nReserved2 = curve.m_nReserved2;
 
  m_nCount = curve.m_nCount;
 
  m_extensionType = curve.m_extensionType;
 
  m_nDesiredSize = curve.m_nDesiredSize;
 
  if (curve.m_pCalc)
    m_pCalc = curve.m_pCalc->NewCopy();
 
  if (m_nCount) {
    m_pSamples = (icFloatNumber*)malloc(m_nCount * sizeof(icFloatNumber));
    if (m_pSamples)
      memcpy(m_pSamples, curve.m_pSamples, m_nCount * sizeof(icFloatNumber));
    else
      m_nCount = 0;
  }
  else {
    m_pSamples = NULL;
  }
 
  m_firstEntry = curve.m_firstEntry;
  m_lastEntry = curve.m_lastEntry;
 
  m_loIntercept = curve.m_loIntercept;
  m_loSlope = curve.m_loSlope;
  m_hiIntercept = curve.m_hiIntercept;
  m_hiSlope = curve.m_hiSlope;
}
 
/**
******************************************************************************
* Name: &CIccSampledCalculatorCurve::operator=
*
* Purpose:
*
* Args:
*
* Return:
******************************************************************************/
CIccSampledCalculatorCurve &CIccSampledCalculatorCurve::operator=(const CIccSampledCalculatorCurve &curve)
{
  if (m_pCalc)
    delete m_pCalc;
 
  if (m_pSamples)
    free(m_pSamples);
 
  m_nReserved = curve.m_nReserved;
  m_nReserved2 = curve.m_nReserved2;
 
  m_nCount = curve.m_nCount;
 
  m_nDesiredSize = curve.m_nDesiredSize;
 
  m_extensionType = curve.m_extensionType;
 
  if (curve.m_pCalc)
    m_pCalc = curve.m_pCalc->NewCopy();
 
  if (m_nCount) {
    m_pSamples = (icFloatNumber*)malloc(m_nCount * sizeof(icFloatNumber));
    if (m_pSamples)
      memcpy(m_pSamples, curve.m_pSamples, m_nCount * sizeof(icFloatNumber));
    else
      m_nCount = 0;
  }
  else {
    m_pSamples = NULL;
  }
 
  m_firstEntry = curve.m_firstEntry;
  m_lastEntry = curve.m_lastEntry;
 
  m_loIntercept = curve.m_loIntercept;
  m_loSlope = curve.m_loSlope;
  m_hiIntercept = curve.m_hiIntercept;
  m_hiSlope = curve.m_hiSlope;
  return (*this);
}
 
/**
******************************************************************************
* Name: CIccSampledCalculatorCurve::~CIccSampledCalculatorCurve
*
* Purpose:
*
* Args:
*
* Return:
******************************************************************************/
CIccSampledCalculatorCurve::~CIccSampledCalculatorCurve()
{
  if (m_pSamples)
    free(m_pSamples);
 
  if (m_pCalc)
    delete m_pCalc;
}
 
/**
******************************************************************************
* Name: CIccSampledCalculatorCurve::SetRange
*
* Purpose:
*  Sets range of sampled lookup table.
*
* Args:
*  first = first entry input position
*  last = last entry input position
******************************************************************************/
void CIccSampledCalculatorCurve::SetRange(icFloatNumber first/* =0.0f */, icFloatNumber last/* =1.0f */)
{
  if (first < last) {
    m_firstEntry = first;
    m_lastEntry = last;
  }
  else {
    m_firstEntry = last;
    m_lastEntry = first;
  }
}
 
 
/**
******************************************************************************
* Name: CIccSampledCalculatorCurve::SetExtensionType
*
* Purpose:
*  Sets extension type of the single sampled curve
*
* Args:
*  nExtensionType is type of extension to use
*
* Return:
*  true if valid extension type, false otherwise
******************************************************************************/
bool CIccSampledCalculatorCurve::SetExtensionType(icUInt16Number nExtensionType)
{
  m_extensionType = nExtensionType;
 
  switch (nExtensionType) {
  case icClipSingleSampledCurve:
  case icExtendSingleSampledCurve:
    return true;
  }
  return false;
}
 
 
/**
******************************************************************************
* Name: CIccSampledCalculatorCurve::SetCalculator
* Purpose:
*  Sets calculator associated with the curve
*
* Args:
*  nSize = desired number of elements in lut (must be >= 2).
*
*
* Return:
*  true if calculator element is valid for curve, false otherwise
******************************************************************************/
bool CIccSampledCalculatorCurve::SetCalculator(CIccMpeCalculator *pCalc)
{
  if (m_pCalc)
    delete m_pCalc;
 
  m_pCalc = pCalc;
 
  return true;
}
 
 
/**
******************************************************************************
* Name: CIccSampledCalculatorCurve::SetRecommendedSize
*
* Purpose:
*  Sets size of sampled lookup table.  Previous data (if exists) is lost.
*
* Args:
*  nSize = desired number of elements in lut (must be >= 2).
*
* Return:
*  true if amount is valid.
******************************************************************************/
bool CIccSampledCalculatorCurve::SetRecommendedSize(icUInt32Number nSize)
{
  if (nSize < 2) {
    return false;
  }
 
  m_nDesiredSize = nSize;
 
  return true;
}
 
 
/**
******************************************************************************
* Name: CIccSampledCalculatorCurve::SetSize
*
* Purpose:
*  Sets size of sampled lookup table.  Previous data (if exists) is lost.
*
* Args:
*  nCount = number of elements in lut (must be >= 2).
*  bZeroAlloc = flag to decide if memory should be set to zero.
*
* Return:
*  true if allocation successful.
******************************************************************************/
bool CIccSampledCalculatorCurve::SetSize(icUInt32Number nCount, bool bZeroAlloc/*=true*/)
{
  if (nCount<2) {
    if (m_pSamples)
      free(m_pSamples);
    m_pSamples = NULL;
    m_nCount = nCount;
    return nCount == 0;
  }
 
  if (m_pSamples) {
    free(m_pSamples);
  }
 
  if (bZeroAlloc)
    m_pSamples = (icFloatNumber*)calloc(nCount, sizeof(icFloatNumber));
  else
    m_pSamples = (icFloatNumber*)malloc(nCount * sizeof(icFloatNumber));
 
  if (m_pSamples)
    m_nCount = nCount;
  else
    m_nCount = 0;
 
  return (m_pSamples != NULL);
}
 
/**
******************************************************************************
* Name: CIccSampledCalculatorCurve::Describe
*
* Purpose:
*
* Args:
*
* Return:
******************************************************************************/
void CIccSampledCalculatorCurve::Describe(std::string &sDescription, int nVerboseness/*=*100*/)
{
  icChar buf[128];
 
  if (!m_pCalc) {
    sDescription += "Empty Sampled Calculator Curve [";
    sprintf(buf, "%.8f, ", m_firstEntry);
    sDescription += buf;
 
    sprintf(buf, "%.8f,", m_lastEntry);
    sDescription += buf;
 
    sprintf(buf, "%d", m_nDesiredSize);
    sDescription += buf;
 
    sprintf(buf, "]\r\n");
    sDescription += buf;
  }
  else {
    sDescription += "Sampled Calculator Curve [";
 
    sprintf(buf, "%.8f, ", m_firstEntry);
    sDescription += buf;
 
    sprintf(buf, "%.8f,", m_lastEntry);
    sDescription += buf;
 
    sprintf(buf, "%d", m_nDesiredSize);
    sDescription += buf;
 
    sprintf(buf, "]\r\n");
    sDescription += buf;
 
    switch (m_extensionType) {
    case icClipSingleSampledCurve:
      sDescription += "Clip End Points\n";
      break;
    case icExtendSingleSampledCurve:
      sDescription += "Extend End Points\n";
      break;
    default:
      sDescription += "Unknown extension handling\n";
      break;
    }
    m_pCalc->Describe(sDescription, nVerboseness);
  }
  sDescription += "\r\n";
}
 
/**
******************************************************************************
* Name: CIccSampledCalculatorCurve::Read
*
* Purpose:
*
* Args:
*
* Return:
******************************************************************************/
bool CIccSampledCalculatorCurve::Read(icUInt32Number size, CIccIO *pIO)
{
  icCurveSegSignature sig;
 
  icUInt32Number headerSize = sizeof(icTagTypeSignature) +
    sizeof(icUInt32Number) +
    sizeof(icUInt32Number) +
    sizeof(icFloat32Number) +
    sizeof(icFloat32Number) +
    sizeof(icUInt16Number) +
    sizeof(icUInt16Number);
 
  if (headerSize >= size)
    return false;
 
  if (!pIO) {
    return false;
  }
 
  if (!pIO->Read32(&sig) || sig != GetType())
    return false;
 
  if (!pIO->Read32(&m_nReserved))
    return false;
 
  if (!pIO->Read32(&m_nDesiredSize))
    return false;
 
  if (!pIO->ReadFloat32Float(&m_firstEntry))
    return false;
 
  if (!pIO->ReadFloat32Float(&m_lastEntry))
    return false;
 
  if (!pIO->Read16(&m_extensionType))
    return false;
 
  if (!pIO->Read16(&m_nReserved2))
    return false;
 
  if (!SetCalculator((CIccMpeCalculator*)CIccMultiProcessElement::Create(icSigCalculatorElemType)) || !m_pCalc)
    return false;
 
  if (!m_pCalc->Read(size - headerSize, pIO))
    return false;
 
  return true;
}
 
/**
******************************************************************************
* Name: CIccSampledCalculatorCurve::Write
*
* Purpose:
*
* Args:
*
* Return:
******************************************************************************/
bool CIccSampledCalculatorCurve::Write(CIccIO *pIO)
{
  icCurveElemSignature sig = GetType();
 
  if (!pIO)
    return false;
 
  if (!pIO->Write32(&sig))
    return false;
 
  if (!pIO->Write32(&m_nReserved))
    return false;
 
  if (!pIO->Write32(&m_nDesiredSize))
    return false;
 
  if (!pIO->WriteFloat32Float(&m_firstEntry))
    return false;
 
  if (!pIO->WriteFloat32Float(&m_lastEntry))
    return false;
 
  if (!pIO->Write16(&m_extensionType))
    return false;
 
  if (!pIO->Write16(&m_nReserved2))
    return false;
 
  if (m_pCalc && m_pCalc->Write(pIO)) {
    return true;
  }
 
  return false;
}
 
/**
******************************************************************************
* Name: CIccSampledCalculatorCurve::Begin
*
* Purpose:
*
* Args:
*
* Return:
******************************************************************************/
bool CIccSampledCalculatorCurve::Begin(icElemInterp nInterp, CIccTagMultiProcessElement *pMPE)
{
  if (m_nDesiredSize < 2 || !m_pCalc || !m_pCalc->Begin(nInterp, pMPE))
    return false;
  
  icUInt32Number nSize = m_nDesiredSize;
  if (nSize < 2)
    nSize = 2;
  if (nSize > ICC_MAXCALCCURVESIZE)
    nSize = ICC_MAXCALCCURVESIZE;
 
  SetSize(nSize);
 
  m_range = m_lastEntry - m_firstEntry;
 
  if (m_range == 0.0)
    return false;
 
  m_last = (icFloatNumber)(m_nCount - 1);
  icFloatNumber stepSize = m_range / m_last;
 
  //Use calculator element to populate lookup table
  CIccApplyMpe *pApply = m_pCalc->GetNewApply(NULL);
  for (icUInt32Number i = 0; i < m_nCount; i++) {
    icFloatNumber src, dst;
    src = (icFloatNumber)i / m_last * m_range + m_firstEntry;
 
    m_pCalc->Apply(pApply, &dst, &src);
    m_pSamples[i] = dst;
  }
 
  if (pApply)
    delete pApply;
 
  switch (m_extensionType) {
  case icClipSingleSampledCurve:
    m_loSlope = 0.0;
    m_loIntercept = m_pSamples[0];
 
    m_hiSlope = 0.0;
    m_hiIntercept = m_pSamples[m_nCount - 1];
    break;
 
  case icExtendSingleSampledCurve:
    m_loSlope = (m_pSamples[1] - m_pSamples[0]) / stepSize;
    m_loIntercept = m_pSamples[0] - m_loSlope*m_firstEntry;
 
    m_hiSlope = (m_pSamples[m_nCount - 1] - m_pSamples[m_nCount - 2]) / stepSize;
    m_hiIntercept = m_pSamples[m_nCount - 1] - m_hiSlope*m_lastEntry;
    break;
 
  default:
    return false;
  }
 
  return true;
}
 
/**
******************************************************************************
* Name: CIccSampledCalculatorCurve::Apply
*
* Purpose:
*
* Args:
*
* Return:
******************************************************************************/
icFloatNumber CIccSampledCalculatorCurve::Apply(icFloatNumber v) const
{
  if (v < m_firstEntry) {
    return m_loSlope * v + m_loIntercept;;
  }
  else if (v > m_lastEntry) {
    return m_hiSlope * v + m_hiIntercept;
  }
 
  icFloatNumber pos = (v - m_firstEntry) / m_range * m_last;
  icUInt32Number index = (icUInt32Number)pos;
  icFloatNumber remainder = pos - (icFloatNumber)index;
 
  if (remainder == 0.0)
    return m_pSamples[index];
 
  return (icFloatNumber)((1.0 - remainder)*m_pSamples[index] + remainder*m_pSamples[index + 1]);
}
 
/**
******************************************************************************
* Name: CIccSampledCalculatorCurve::Validate
*
* Purpose:
*
* Args:
*
* Return:
******************************************************************************/
icValidateStatus CIccSampledCalculatorCurve::Validate(std::string sigPath, std::string &sReport, const CIccTagMultiProcessElement* pMPE/*=NULL*/, const CIccProfile *pProfile/*=NULL*/) const
{
  CIccInfo Info;
  std::string sSigPathName = Info.GetSigPathName(sigPath);
 
  icValidateStatus rv = icValidateOK;
  if (m_nReserved || m_nReserved2) {
    sReport += icMsgValidateWarning;
    sReport += sSigPathName;
    sReport += " sampled calculator curve has non zero reserved data.\r\n";
    rv = icValidateWarning;
  }
 
  if (m_extensionType > icMaxSingleSampledCurveType) {
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " sampled calculator curve has unknown extension type\r\n";
    rv = icMaxStatus(rv, icValidateCriticalError);
  }
 
  if (m_nDesiredSize < 2) {
    sReport += icMsgValidateWarning;
    sReport += sSigPathName;
    sReport += " sampled calculator curve has too few desired sample points.\r\n";
    rv = icMaxStatus(rv, icValidateWarning);
  }
 
  if (m_lastEntry - m_firstEntry <= 0.0) {
    sReport += icMsgValidateWarning;
    sReport += sSigPathName;
    sReport += " sampled calculator curve has an invalid sample range.\n";
    rv = icMaxStatus(rv, icValidateWarning);
  }
 
  if (pProfile && pProfile->m_Header.version < icVersionNumberV5_1) {
    sReport += icMsgValidateWarning;
    sReport += sSigPathName;
    sReport += " sampled calculator curve is not supported by version of profile.\n";
    rv = icMaxStatus(rv, icValidateWarning);
  }
 
  if (!m_pCalc) {
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " sampled calculator curve has no calculator element\r\n";
    rv = icMaxStatus(rv, icValidateCriticalError);
  }
  else {
    icValidateStatus stat = m_pCalc->Validate(sSigPathName, sReport, pMPE, pProfile);
 
    rv = icMaxStatus(rv, stat);
  }
 
  return rv;
}
 
 
/**
 ******************************************************************************
 * Name: CIccCurveSegment::Create
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccCurveSegment* CIccCurveSegment::Create(icCurveSegSignature sig, icFloatNumber start, icFloatNumber end)
{
  switch(sig) {
  case icSigFormulaCurveSeg:
    return new CIccFormulaCurveSegment(start, end);
  case icSigSampledCurveSeg:
    return new CIccSampledCurveSegment(start, end);
  default:
    return NULL;
  }
 
}
 
/**
 ******************************************************************************
 * Name: CIccSegmentedCurve::CIccSegmentedCurve
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccSegmentedCurve::CIccSegmentedCurve()
{
  m_list = new CIccCurveSegmentList();
  m_nReserved1 = 0;
  m_nReserved2 = 0;
 
}
 
 
/**
 ******************************************************************************
 * Name: CIccSegmentedCurve::CIccSegmentedCurve
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccSegmentedCurve::CIccSegmentedCurve(const CIccSegmentedCurve &curve)
{
  CIccCurveSegmentList::iterator i;
 
  m_list = new CIccCurveSegmentList();
 
  for (i=curve.m_list->begin(); i!=curve.m_list->end(); i++) {
    m_list->push_back((*i)->NewCopy());
  }
  m_nReserved1 = curve.m_nReserved1;
  m_nReserved2 = curve.m_nReserved2;
}
 
 
/**
 ******************************************************************************
 * Name: &CIccSegmentedCurve::operator=
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccSegmentedCurve &CIccSegmentedCurve::operator=(const CIccSegmentedCurve &curve)
{
  Reset();
 
  CIccCurveSegmentList::iterator i;
 
  for (i=curve.m_list->begin(); i!=curve.m_list->end(); i++) {
    m_list->push_back((*i)->NewCopy());
  }
  m_nReserved1 = curve.m_nReserved1;
  m_nReserved2 = curve.m_nReserved2;
 
  return (*this);
}
 
 
/**
 ******************************************************************************
 * Name: CIccSegmentedCurve::~CIccSegmentedCurve
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccSegmentedCurve::~CIccSegmentedCurve()
{
  Reset();
  delete m_list;
}
 
 
/**
 ******************************************************************************
 * Name: CIccSegmentedCurve::Describe
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccSegmentedCurve::Describe(std::string &sDescription, int nVerboseness)
{
  CIccCurveSegmentList::iterator i;
 
  sDescription += "BEGIN_CURVE\n";
  for (i=m_list->begin(); i!=m_list->end(); i++) {
    (*i)->Describe(sDescription, nVerboseness);
  }
}
 
 
 
/**
 ******************************************************************************
 * Name: CIccSegmentedCurve::Read
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccSegmentedCurve::Read(icUInt32Number size, CIccIO *pIO)
{
  icCurveElemSignature sig;
 
  icUInt32Number startPos = pIO->Tell();
  
  icUInt32Number headerSize = sizeof(icCurveElemSignature) + 
    sizeof(icUInt32Number) + 
    sizeof(icUInt16Number) + 
    sizeof(icUInt16Number);
 
  if (headerSize > size)
    return false;
 
  if (!pIO) {
    return false;
  }
 
  if (!pIO->Read32(&sig))
    return false;
 
  if (sig!=GetType())
    return false;
 
  if (!pIO->Read32(&m_nReserved1))
    return false;
 
  icUInt16Number nSegments;
 
  if (!pIO->Read16(&nSegments))
    return false;
 
  if (!pIO->Read16(&m_nReserved2))
    return false;
 
  Reset();
 
  icUInt32Number pos = pIO->Tell();
  icCurveSegSignature segSig;
  CIccCurveSegment *pSeg;
 
  if (nSegments==1) {
    if (sizeof(segSig) > size - headerSize)
      return false;
 
    if (!pIO->Read32(&segSig))
      return false;
    pSeg = CIccCurveSegment::Create(segSig, icMinFloat32Number, icMaxFloat32Number);
    if (!pSeg)
      return false;
 
    pIO->Seek(pos, icSeekSet);
 
    if (!pSeg->Read(size-(pos-startPos), pIO)) {
      delete pSeg;
      return false;
    }
 
    m_list->push_back(pSeg);
  }
  else if (nSegments) {
    if (nSegments > size - headerSize)
      return false;
 
    icFloatNumber *breakpoints=(icFloatNumber*)calloc(nSegments-1, sizeof(icFloatNumber));
 
    if (!breakpoints)
      return false;
 
    if (pIO->ReadFloat32Float(breakpoints, nSegments-1)!=nSegments-1) {
      free(breakpoints);
      return false;
    }
 
    int i;
    for (i=0; i<nSegments; i++) {
      pos = pIO->Tell();
      if (!pIO->Read32(&segSig)) {
        free(breakpoints);
        return false;
      }
      if (pIO->Seek(pos, icSeekSet)!=(icInt32Number)pos)
        return false;;
 
      if (!i)
        pSeg = CIccCurveSegment::Create(segSig, icMinFloat32Number, breakpoints[i]);
      else if (i==nSegments-1) 
        pSeg = CIccCurveSegment::Create(segSig, breakpoints[i-1], icMaxFloat32Number);
      else
        pSeg = CIccCurveSegment::Create(segSig, breakpoints[i-1], breakpoints[i]);
      
      if (!pSeg) {
        free(breakpoints);
        return false;
      }
 
      if (!pSeg->Read(size-(pos-startPos), pIO)) {
        delete pSeg;
        free(breakpoints);
        return false;
      }
 
      m_list->push_back(pSeg);
    }
 
    free(breakpoints);
  }
 
  return true;
}
 
 
/**
 ******************************************************************************
 * Name: CIccSegmentedCurve::Write
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccSegmentedCurve::Write(CIccIO *pIO)
{
  icCurveElemSignature sig = GetType();
 
  if (!pIO)
    return false;
 
  if (!pIO->Write32(&sig))
    return false;
 
  if (!pIO->Write32(&m_nReserved1))
    return false;
 
  icUInt16Number nSegments = (icUInt16Number)m_list->size();
 
  if (!pIO->Write16(&nSegments))
    return false;
 
  if (!pIO->Write16(&m_nReserved2))
    return false;
 
  CIccCurveSegmentList::iterator i;
  if (nSegments>1) {
    icFloatNumber breakpoint;
    
    i=m_list->begin();
    for (i++; i!=m_list->end(); i++) {
      breakpoint = (*i)->StartPoint();
      if (!pIO->WriteFloat32Float(&breakpoint))
        return false;
    }
  }
  for (i=m_list->begin(); i!=m_list->end(); i++) {
    if (!(*i)->Write(pIO))
      return false;
  }
  
  return true;
}
 
 
/**
 ******************************************************************************
 * Name: CIccSegmentedCurve::Reset
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccSegmentedCurve::Reset()
{
  CIccCurveSegmentList::iterator i;
 
  for (i=m_list->begin(); i!=m_list->end(); i++) {
    delete (*i);
  }
  m_list->clear();
}
 
 
/**
 ******************************************************************************
 * Name: CIccSegmentedCurve::Insert
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccSegmentedCurve::Insert(CIccCurveSegment *pCurveSegment)
{
  CIccCurveSegmentList::reverse_iterator last = m_list->rbegin();
 
  if (last!=m_list->rend()) {
    if (pCurveSegment->StartPoint() == (*last)->EndPoint()) {
      m_list->push_back(pCurveSegment);
      return true;
    }
  }
  else {
    m_list->push_back(pCurveSegment);
    return true;
  }
 
  return false;
}
 
 
/**
 ******************************************************************************
 * Name: CIccSegmentedCurve::Begin
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccSegmentedCurve::Begin(icElemInterp nInterp, CIccTagMultiProcessElement *pMPE)
{
  if (m_list->size()==0)
    return false;
 
 CIccCurveSegmentList::iterator i;
 CIccCurveSegment *pLast = NULL;
 
  for (i=m_list->begin(); i!=m_list->end(); i++) {
    if (!(*i)->Begin(pLast))
      return false;
    pLast = *i;
  }
 
  return true;
}
 
 
/**
 ******************************************************************************
 * Name: CIccSegmentedCurve::Apply
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
icFloatNumber CIccSegmentedCurve::Apply(icFloatNumber v) const
{
 CIccCurveSegmentList::iterator i;
 
  for (i=m_list->begin(); i!=m_list->end(); i++) {
    if (v <= (*i)->EndPoint())
      return (*i)->Apply(v);
  }
  return v;
}
 
/**
 ******************************************************************************
 * Name: CIccSegmentedCurve::Validate
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
icValidateStatus CIccSegmentedCurve::Validate(std::string sigPath, std::string &sReport, const CIccTagMultiProcessElement* pMPE/*=NULL*/, const CIccProfile *pProfile/*=NULL*/) const
{
  CIccInfo Info;
  std::string sSigPathName = Info.GetSigPathName(sigPath);
 
  icValidateStatus rv = icValidateOK;
  if (m_nReserved1 || m_nReserved2) {
    sReport += icMsgValidateWarning;
    sReport += sSigPathName;
    sReport += " Segmented curve has non zero reserved data.\n";
    rv = icValidateWarning;
  }
 
  if (m_list->size()==0) {
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " Has Empty CurveSegment!\n";
    return icMaxStatus(rv, icValidateCriticalError);
  }
 
  CIccCurveSegmentList::iterator i;
 
  for (i=m_list->begin(); i!=m_list->end(); i++) {
    rv = icMaxStatus(rv, (*i)->Validate(sigPath+icGetSigPath(GetType()), sReport, pMPE, pProfile));
  }
 
  return rv;
}
 
 
/**
 ******************************************************************************
 * Name: CIccCurveSetCurve::Create
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccCurveSetCurve* CIccCurveSetCurve::Create(icCurveElemSignature sig)
{
  switch (sig) {
  case icSigSegmentedCurve:
    return new CIccSegmentedCurve();
  case icSigSingleSampledCurve:
    return new CIccSingleSampledCurve();
  case icSigSampledCalculatorCurve:
    return new CIccSampledCalculatorCurve();
  default:
    return NULL;
  }
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCurveSet::CIccMpeCurveSet
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeCurveSet::CIccMpeCurveSet(int nSize/*=0*/)
{
  m_nReserved = 0;
  if (nSize) {
    m_nInputChannels = m_nOutputChannels = nSize;
    m_curve = (icCurveSetCurvePtr*)calloc(nSize, sizeof(icCurveSetCurvePtr));
    m_position = (icPositionNumber*)calloc(nSize, sizeof(icPositionNumber));
  }
  else {
    m_nInputChannels = m_nOutputChannels = 0;
    m_curve = NULL;
    m_position = NULL;
  }
}
 
typedef std::map<icCurveSetCurvePtr, icCurveSetCurvePtr> icCurveMap;
 
/**
 ******************************************************************************
 * Name: CIccMpeCurveSet::CIccMpeCurveSet
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeCurveSet::CIccMpeCurveSet(const CIccMpeCurveSet &curveSet)
{
  m_nReserved = curveSet.m_nReserved;
 
  if (curveSet.m_nInputChannels) {
    int i;
 
    m_nInputChannels = m_nOutputChannels = curveSet.m_nInputChannels;
    m_curve = (icCurveSetCurvePtr*)calloc(m_nInputChannels, sizeof(icCurveSetCurvePtr));
    m_position = (icPositionNumber*)calloc(m_nInputChannels, sizeof(icPositionNumber));
 
    icCurveMap map;
    for (i=0; i<m_nInputChannels; i++) {
      CIccCurveSetCurve *ptr = curveSet.m_curve[i];
      if (ptr) {
        if (!map[ptr]) {
          m_curve[i] = ptr->NewCopy();
          map[ptr] = m_curve[i];
        }
        else
          m_curve[i] = map[ptr];
      }
    }
  }
  else {
    m_nInputChannels = m_nOutputChannels = 0;
    m_curve = NULL;
  }
}
 
/**
 ******************************************************************************
 * Name: &CIccMpeCurveSet::operator=
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeCurveSet &CIccMpeCurveSet::operator=(const CIccMpeCurveSet &curveSet)
{
  m_nReserved = m_nReserved;
 
  if (m_curve) {
    free(m_curve);
  }
 
  if (curveSet.m_nInputChannels) {
    int i;
 
    m_nInputChannels = m_nOutputChannels = curveSet.m_nInputChannels;
    m_curve = (icCurveSetCurvePtr*)calloc(m_nInputChannels, sizeof(icCurveSetCurvePtr));
    m_position = (icPositionNumber*)calloc(m_nInputChannels, sizeof(icPositionNumber));
 
    icCurveMap map;
    for (i=0; i<m_nInputChannels; i++) {
      CIccCurveSetCurve *ptr = curveSet.m_curve[i];
      if (ptr) {
        if (!map[ptr]) {
          m_curve[i] = ptr->NewCopy();
          map[ptr] = m_curve[i];
        }
        else
          m_curve[i] = map[ptr];
      }
    }
  }
  else {
    m_nInputChannels = m_nOutputChannels = 0;
    m_curve = NULL;
  }
 
  return *this;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCurveSet::~CIccMpeCurveSet
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeCurveSet::~CIccMpeCurveSet()
{
  SetSize(0);
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCurveSet::SetSize
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccMpeCurveSet::SetSize(int nNewSize)
{
  if (m_curve) {
    icCurveMap map;
    int i;
 
    for (i=0; i<m_nInputChannels; i++) {
      if (!map[m_curve[i]]) {
        map[m_curve[i]] = (CIccCurveSetCurve*)1;
        delete m_curve[i];
      }
    }
    free(m_curve);
  }
  if (m_position) {
    free(m_position);
  }
 
  if (nNewSize) {
    m_curve = (icCurveSetCurvePtr*)calloc(nNewSize, sizeof(icCurveSetCurvePtr));
 
    if (!m_curve) {
      m_position = NULL;
      m_nInputChannels = m_nOutputChannels = 0;
      return false;
    }
 
    m_position = (icPositionNumber*)calloc(nNewSize, sizeof(icPositionNumber));
 
    if (!m_position) {
      free(m_curve);
      m_curve = NULL;
 
      m_nInputChannels = m_nOutputChannels = 0;
      return false;
    }
    m_nInputChannels = m_nOutputChannels = nNewSize;
  }
  else {
    m_curve = NULL;
    m_nInputChannels = m_nOutputChannels = 0;
  }
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCurveSet::SetCurve
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccMpeCurveSet::SetCurve(int nIndex, icCurveSetCurvePtr newCurve)
{
  if (nIndex<0 || nIndex>m_nInputChannels)
    return false;
 
  if (m_curve) {
    int i;
 
    for (i=0; i<m_nInputChannels; i++)
      if (i!=nIndex && m_curve[i]==m_curve[nIndex])
        break;
 
    if (i==m_nInputChannels && m_curve[nIndex]) {
      delete m_curve[nIndex];
    }
  }
  m_curve[nIndex] = newCurve;
  
  return true;
}
 
 
/**
 ******************************************************************************
 * Name: CIccMpeCurveSet::Describe
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccMpeCurveSet::Describe(std::string &sDescription, int nVerboseness)
{
  if (m_curve) {
    icChar buf[81];
    int i;
 
    sprintf(buf, "BEGIN_CURVE_SET %d\n", m_nInputChannels);
    sDescription += buf;
 
    for (i=0; i<m_nInputChannels; i++) {
      sprintf(buf, "Curve %d of %d\n", i+1, m_nInputChannels);
      sDescription += buf;
      if (m_curve[i]) {
        m_curve[i]->Describe(sDescription, nVerboseness);
      }
    }
    sDescription += "END_CURVE_SET\n";
  }
}
 
typedef std::map<icUInt32Number, CIccCurveSetCurve*> icCurveOffsetMap;
typedef std::map<CIccCurveSetCurve*, icPositionNumber> icCurvePtrMap;
 
/**
 ******************************************************************************
 * Name: CIccMpeCurveSet::Read
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccMpeCurveSet::Read(icUInt32Number size, CIccIO *pIO)
{
  icElemTypeSignature sig;
 
  icUInt32Number startPos = pIO->Tell();
  
  icUInt32Number headerSize = sizeof(icTagTypeSignature) + 
    sizeof(icUInt32Number) + 
    sizeof(icUInt16Number) + 
    sizeof(icUInt16Number);
 
  if (headerSize > size)
    return false;
 
  if (!pIO) {
    return false;
  }
 
  icUInt16Number nInputChannels, nOutputChannels;
 
  if (!pIO->Read32(&sig))
    return false;
 
  if (!pIO->Read32(&m_nReserved))
    return false;
 
  if (!pIO->Read16(&nInputChannels))
    return false;
 
  if (!pIO->Read16(&nOutputChannels))
    return false;
 
  if (nInputChannels != nOutputChannels)
    return false;
 
  if (nInputChannels > size - headerSize || nInputChannels * 2*sizeof(icUInt32Number) > size - headerSize)
    return false;
 
  if (!SetSize(nInputChannels))
    return false;
 
  if (m_curve) {
    int i;
 
    if (headerSize + m_nInputChannels*2*sizeof(icUInt32Number) > size)
      return false;
 
    for (i=0; i<m_nInputChannels; i++) {
      if (!pIO->Read32(&m_position[i].offset)) {
        return false;
      }
      if (!pIO->Read32(&m_position[i].size)) {
        return false;
      }
    }
 
    icCurveOffsetMap map;
    icCurveElemSignature curveSig;
    for (i=0; i<m_nInputChannels; i++) {
      if (!map[m_position[i].offset]) {
        icUInt32Number pos;
        if (!m_position[i].offset || !m_position[i].size) {
          return false;
        }
 
        pos = startPos + m_position[i].offset;
        if (pIO->Seek(pos, icSeekSet)!=(icInt32Number)pos) {
          return false;
        }
        
        if (!pIO->Read32(&curveSig)) {
          return false;
        }
        m_curve[i] = CIccCurveSetCurve::Create(curveSig);
 
        if (!m_curve[i]) {
          return false;
        }
 
        if (pIO->Seek(pos, icSeekSet)!=(icInt32Number)pos) {
          return false;
        }
      
        if (!m_curve[i]->Read(m_position[i].size, pIO)) {
          return false;
        }
 
        map[m_position[i].offset] = m_curve[i];
      }
      else {
        m_curve[i] = map[m_position[i].offset];
      }
    }
  }
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCurveSet::Write
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccMpeCurveSet::Write(CIccIO *pIO)
{
  icElemTypeSignature sig = GetType();
 
  if (!pIO)
    return false;
 
  icUInt32Number elemStart = pIO->Tell();
 
  if (!pIO->Write32(&sig))
    return false;
 
  if (!pIO->Write32(&m_nReserved))
    return false;
 
  if (!pIO->Write16(&m_nInputChannels))
    return false;
 
  if (!pIO->Write16(&m_nInputChannels))
    return false;
 
  if (m_curve && m_nInputChannels) {
    int i;
    icCurvePtrMap map;
    icUInt32Number start, end;
    icUInt32Number zeros[2] = { 0, 0};
    icPositionNumber position;
 
    icUInt32Number startTable = pIO->Tell();
 
    //First write empty position table
    for (i=0; i<m_nInputChannels; i++) {
      if (pIO->Write32(&zeros[0], 2)!=2)
        return false;
    }
 
    //Now Write curves
    for (i=0; i<m_nInputChannels; i++) {
      if (m_curve[i]) {
        if (map.find(m_curve[i])==map.end()) {
          start = pIO->Tell();
          m_curve[i]->Write(pIO);
          end = pIO->Tell();
          pIO->Align32();
          position.offset = start - elemStart;
          position.size = end - start;
          map[m_curve[i]] = position;
        }
        m_position[i] = map[m_curve[i]];
      }
    }
    end = pIO->Tell();
    
    //Back fill position table
    pIO->Seek(startTable, icSeekSet);
    for (i=0; i<m_nInputChannels; i++) {
      if (!pIO->Write32(&m_position[i].offset))
        return false;
      if (!pIO->Write32(&m_position[i].size))
        return false;
    }
 
    pIO->Seek(end, icSeekSet);
  }
  
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCurveSet::Begin
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccMpeCurveSet::Begin(icElemInterp nInterp, CIccTagMultiProcessElement *pMPE)
{
  if (!m_curve)
    return false;
 
  icCurveMap map;
 
  int i;
  for (i=0; i<m_nInputChannels; i++) {
    if (!m_curve[i])
      return false;
 
    if (!map[m_curve[i]]) {
      if (!m_curve[i]->Begin(nInterp, pMPE)) {
        return false;
      }
      //ensure that Begin is called only once for each curve
      map[m_curve[i]] = (CIccCurveSetCurve*)1;
    }
  }
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCurveSet::Apply
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccMpeCurveSet::Apply(CIccApplyMpe *pApply, icFloatNumber *pDestPixel, const icFloatNumber *pSrcPixel) const
{
  int i;
  for (i=0; i<m_nInputChannels; i++) {
    *pDestPixel++ = m_curve[i]->Apply(*pSrcPixel++);
  }
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCurveSet::Validate
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
icValidateStatus CIccMpeCurveSet::Validate(std::string sigPath, std::string &sReport, const CIccTagMultiProcessElement* pMPE/*=NULL*/, const CIccProfile *pProfile/*=NULL*/) const
{
  std::string mpeSigPath = sigPath + icGetSigPath(GetType());
  icValidateStatus rv = CIccMultiProcessElement::Validate(sigPath, sReport, pMPE, pProfile);
 
  bool empty=false;
  if (m_curve) {
    int i;
    for (i=0; !empty && i<m_nInputChannels; i++) {
      if (!m_curve[i]) {
        empty = true;
      }
      else {
        rv = icMaxStatus(rv, m_curve[i]->Validate(mpeSigPath, sReport, pMPE, pProfile));
      }
    }
  }
  else
    empty = true;
  
  if (empty) {
    CIccInfo Info;
    std::string sSigPathName = Info.GetSigPathName(mpeSigPath);
 
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " -  Has Empty Curve Element(s)!\n";
    return icValidateCriticalError;
  }
 
  return rv;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeTintArray::CIccMpeTintArray
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeTintArray::CIccMpeTintArray(int nVectorSize/*=0*/)
{
  m_nReserved = 0;
  m_nInputChannels = 1;
  m_nOutputChannels = nVectorSize;
  m_Array = NULL;
}
 
typedef std::map<icCurveSetCurvePtr, icCurveSetCurvePtr> icCurveMap;
 
/**
 ******************************************************************************
 * Name: CIccMpeTintArray::CIccMpeTintArray
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeTintArray::CIccMpeTintArray(const CIccMpeTintArray &tintArray)
{
  m_nReserved = tintArray.m_nReserved;
 
  m_nInputChannels = tintArray.m_nInputChannels;
  m_nOutputChannels = tintArray.m_nOutputChannels;
 
  if (tintArray.m_Array)
    m_Array = (CIccTagNumArray*)tintArray.m_Array->NewCopy();
 
}
 
/**
 ******************************************************************************
 * Name: &CIccMpeTintArray::operator=
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeTintArray &CIccMpeTintArray::operator=(const CIccMpeTintArray &tintArray)
{
  m_nReserved = m_nReserved;
 
  if (m_Array) {
    delete m_Array;
  }
 
  m_nInputChannels = tintArray.m_nInputChannels;
  m_nOutputChannels = tintArray.m_nOutputChannels;
 
  if (tintArray.m_Array)
    m_Array = (CIccTagNumArray*)tintArray.m_Array->NewCopy();
 
  return *this;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeTintArray::~CIccMpeTintArray
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeTintArray::~CIccMpeTintArray()
{
  if (m_Array)
    delete m_Array;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeTintArray::SetVectorSize
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccMpeTintArray::SetVectorSize(int nVectorSize)
{
  m_nOutputChannels = nVectorSize;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeTintArray::SetCurve
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccMpeTintArray::SetArray(CIccTagNumArray *pArray)
{
  if (m_Array)
    delete m_Array;
 
  m_Array = pArray;
}
 
 
/**
 ******************************************************************************
 * Name: CIccMpeTintArray::Describe
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccMpeTintArray::Describe(std::string &sDescription, int nVerboseness)
{
  if (m_Array) {
    icChar buf[81];
 
    sprintf(buf, "BEGIN_TINT_ARRAY %d\n", m_nOutputChannels);
    sDescription += buf;
 
    m_Array->Describe(sDescription, nVerboseness);
 
    sDescription += "END_TINT_ARRAY\n";
  }
}
 
/**
 ******************************************************************************
 * Name: CIccMpeTintArray::Read
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccMpeTintArray::Read(icUInt32Number size, CIccIO *pIO)
{
  if (m_Array)
    delete m_Array;
  m_Array = NULL;
 
  icElemTypeSignature sig;
 
  icUInt32Number startPos = pIO->Tell();
  
  icUInt32Number headerSize = sizeof(icElemTypeSignature) + 
    sizeof(icUInt32Number) + 
    sizeof(icUInt16Number) + 
    sizeof(icUInt16Number);
 
  if (headerSize +sizeof(icTagTypeSignature) > size)
    return false;
 
  if (!pIO) {
    return false;
  }
 
  icUInt16Number nInputChannels, nOutputChannels;
 
  if (!pIO->Read32(&sig))
    return false;
 
  if (!pIO->Read32(&m_nReserved))
    return false;
 
  if (!pIO->Read16(&nInputChannels))
    return false;
 
  if (!pIO->Read16(&nOutputChannels))
    return false;
 
  if (nInputChannels != 1 || !nOutputChannels)
    return false;
 
  m_nInputChannels = nInputChannels;
  m_nOutputChannels = nOutputChannels;
 
  icUInt32Number arrayPos = pIO->Tell();
 
  icTagTypeSignature tagType;
  if (!pIO->Read32(&tagType))
    return false;
 
  CIccTag *pTag = CIccTag::Create(tagType);
  if (!pTag)
    return false;
 
  if (!pTag->IsNumArrayType()) {
    delete pTag;
  }
 
  m_Array = (CIccTagNumArray*)pTag;
  pIO->Seek(arrayPos, icSeekSet);
  if (!m_Array->Read(size-headerSize, pIO)) {
    return false;
  }
 
  icUInt32Number nVals = m_Array->GetNumValues();
  if (nVals/m_nOutputChannels <2 || (nVals%m_nOutputChannels) != 0) {
    return false;
  }
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeTintArray::Write
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccMpeTintArray::Write(CIccIO *pIO)
{
  if (!m_Array)
    return false;
 
  icElemTypeSignature sig = GetType();
 
  if (!pIO)
    return false;
 
  if (!pIO->Write32(&sig))
    return false;
 
  if (!pIO->Write32(&m_nReserved))
    return false;
 
  if (!pIO->Write16(&m_nInputChannels))
    return false;
 
  if (!pIO->Write16(&m_nOutputChannels))
    return false;
 
  return m_Array->Write(pIO);
}
 
/**
 ******************************************************************************
 * Name: CIccMpeTintArray::Begin
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccMpeTintArray::Begin(icElemInterp nInterp, CIccTagMultiProcessElement *pMPE)
{
  if (!m_Array)
    return false;
 
  icUInt32Number nVals = m_Array->GetNumValues();
 
  if (nVals/m_nOutputChannels<2 || nVals % m_nOutputChannels != 0)
    return false;
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeTintArray::Apply
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccMpeTintArray::Apply(CIccApplyMpe *pApply, icFloatNumber *pDestPixel, const icFloatNumber *pSrcPixel) const
{
  if (m_Array) {
    m_Array->Interpolate(pDestPixel, *pSrcPixel, m_nOutputChannels);
  }
  else
    memset(pDestPixel, 0, m_nOutputChannels*sizeof(icFloatNumber));
}
 
/**
 ******************************************************************************
 * Name: CIccMpeTintArray::Validate
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
icValidateStatus CIccMpeTintArray::Validate(std::string sigPath, std::string &sReport, const CIccTagMultiProcessElement* pMPE/*=NULL*/, const CIccProfile *pProfile/*=NULL*/) const
{
  std::string mpeSigPath = sigPath + icGetSigPath(GetType());
  icValidateStatus rv = CIccMultiProcessElement::Validate(sigPath, sReport, pMPE, pProfile);
 
  if (m_nInputChannels!=1) {
    CIccInfo Info;
    std::string sSigPathName = Info.GetSigPathName(mpeSigPath);
 
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " -  Bad number of input channels!\n";
    return icValidateCriticalError;
  }
 
  if (!m_nOutputChannels) {
    CIccInfo Info;
    std::string sSigPathName = Info.GetSigPathName(mpeSigPath);
 
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " -  Bad number of output channels!\n";
    return icValidateCriticalError;
  }
 
  if (!m_Array) {
    CIccInfo Info;
    std::string sSigPathName = Info.GetSigPathName(mpeSigPath);
 
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " -  Has no tint values(s)!\n";
    return icValidateCriticalError;
  }
  else {
    icUInt32Number nVals = m_Array->GetNumValues();
    bool bBad = false;
    
    /*
     * NOTE: nVals is unsigned as is m_nOutputChannels
     * so the result will never be < 0
     */
    /*
    if (nVals/m_nOutputChannels<0) {
      CIccInfo Info;
      std::string sSigPathName = Info.GetSigPathName(mpeSigPath);
 
      sReport += icMsgValidateCriticalError;
      sReport += sSigPathName;
      sReport += " -  Needs two or more tint steps!\n";
      bBad = true;
    }
    */
    if ((nVals % m_nOutputChannels)!=0) {
      CIccInfo Info;
      std::string sSigPathName = Info.GetSigPathName(mpeSigPath);
 
      sReport += icMsgValidateCriticalError;
      sReport += sSigPathName;
      sReport += " -  Array size must be multiple of output channels!\n";
      bBad = true;
    }
 
    if (bBad) {
      return icValidateCriticalError;
    }
  }
 
  return rv;
}
 
CIccToneMapFunc::CIccToneMapFunc()
{
  m_nFunctionType = 0;
  m_nParameters = 0;
  m_params = NULL;
  m_nReserved = 0;
  m_nReserved2 = 0;
}
 
 
CIccToneMapFunc::~CIccToneMapFunc()
{
  if (m_params)
    free(m_params);
}
 
CIccToneMapFunc& CIccToneMapFunc::operator=(const CIccToneMapFunc& toneMapFunc)
{
  if (&toneMapFunc == this)
    return *this;
 
  m_nFunctionType = toneMapFunc.m_nFunctionType;
  m_nParameters = toneMapFunc.m_nParameters;
  
  if (m_params)
    free(m_params);
  if (toneMapFunc.m_nParameters && toneMapFunc.m_params) {
    m_params = (icFloatNumber*)malloc(m_nParameters * sizeof(icFloatNumber));
    if (m_params)
      memcpy(m_params, toneMapFunc.m_params, m_nParameters * sizeof(icFloatNumber));
  }
  else
    m_params = NULL;
 
  m_nReserved = toneMapFunc.m_nReserved;
  m_nReserved2 = toneMapFunc.m_nReserved2;
 
  return *this;
}
 
CIccToneMapFunc* CIccToneMapFunc::NewCopy() const
{
  CIccToneMapFunc* rv = new CIccToneMapFunc();
 
  if (rv)
    *rv = *this;
 
  return rv;
}
 
int CIccToneMapFunc::NumArgs() const
{
  if (m_nFunctionType == 0)
    return 3;
  return 0;
}
 
bool CIccToneMapFunc::SetFunction(icUInt16Number nFunc, icUInt8Number nParams, icFloatNumber* pParams)
{
  m_nFunctionType = nFunc;
 
  if (m_params)
    free(m_params);
 
  int nArgs = NumArgs();
  m_nParameters = (icUInt8Number)icIntMin(nParams, NumArgs());
  if (nArgs && pParams) {
    m_params = (icFloatNumber*)calloc(nArgs, sizeof(icFloatNumber));
    if (m_params)
      memcpy(m_params, pParams, icIntMin(nArgs, nParams) * sizeof(icFloatNumber));
  }
 
  return nParams == NumArgs();
}
 
void CIccToneMapFunc::Describe(std::string& sDescription, int nVerboseness)
{
  icChar buf[128];
 
  sprintf(buf, "ToneFunctionType: %04Xh\n", m_nFunctionType);
  sDescription += buf;
 
  switch (m_nFunctionType) {
  case 0x0000:
    sprintf(buf, "Y = %.8f * M * ( X + %.8f) + %.8f\n\n", m_params[0], m_params[1], m_params[2]);
    sDescription += buf;
    return;
 
  default:
    int i;
    sprintf(buf, "Unknown Function with %d parameters:\n\n", m_nParameters);
    sDescription += buf;
 
    for (i = 0; i < m_nParameters; i++) {
      sprintf(buf, "Param[%d] = %.8lf\n\n", i, m_params[i]);
      sDescription += buf;
    }
 
  }
}
 
bool CIccToneMapFunc::Read(icUInt32Number size, CIccIO* pIO)
{
  icCurveSegSignature sig;
 
  icUInt32Number headerSize = sizeof(icTagTypeSignature) +
    sizeof(icUInt32Number) +
    sizeof(icUInt16Number) +
    sizeof(icUInt16Number);
 
  if (headerSize > size)
    return false;
 
  if (!pIO) {
    return false;
  }
 
  if (!pIO->Read32(&sig))
    return false;
 
  if (sig != icSigToneMapFunction)
    return false;
 
  if (!pIO->Read32(&m_nReserved))
    return false;
 
  if (!pIO->Read16(&m_nFunctionType))
    return false;
 
  if (!pIO->Read16(&m_nReserved2))
    return false;
 
  if (m_params) {
    free(m_params);
  }
 
  m_nParameters = (icUInt8Number)((size - headerSize)/sizeof(icFloatNumber));
 
  if (m_nParameters) {
 
    m_params = (icFloatNumber*)malloc(m_nParameters * sizeof(icFloatNumber));
    if (!m_params)
      return false;
 
    if (pIO->ReadFloat32Float(m_params, m_nParameters) != m_nParameters) {
      return false;
    }
  }
  else
    m_params = NULL;
 
  return true;
 
}
 
bool CIccToneMapFunc::Write(CIccIO* pIO)
{
  icToneFunctionSignature sig = GetType();
 
  if (!pIO)
    return false;
 
  if (!pIO->Write32(&sig))
    return false;
 
  if (!pIO->Write32(&m_nReserved))
    return false;
 
  if (!pIO->Write16(&m_nFunctionType))
    return false;
 
  if (!pIO->Write16(&m_nReserved2))
    return false;
 
 
  if (m_nParameters != NumArgs())
    return false;
 
  if (m_nParameters) {
    if (pIO->WriteFloat32Float(m_params, m_nParameters) != m_nParameters)
      return false;
  }
 
  return true;
 
}
 
bool CIccToneMapFunc::Begin()
{
  if (m_nFunctionType != 0x0000)
    return false;
 
  if (m_nParameters < NumArgs() || !m_params)
    return false;
 
  return true;
}
 
icFloatNumber CIccToneMapFunc::Apply(icFloatNumber lumValue, icFloatNumber pixelValue) const
{
  if (!m_nFunctionType && m_params) {
    return m_params[0] * lumValue * (pixelValue + m_params[1]) + m_params[2];
  }
 
  return 0;
}
 
icValidateStatus CIccToneMapFunc::Validate(std::string& sReport, int nVerboseness) const
{
  CIccInfo Info;
 
  icValidateStatus rv = icValidateOK;
 
  if (m_nReserved || m_nReserved2) {
    sReport += icMsgValidateWarning;
    sReport += " formula curve has non zero reserved data.\n";
    rv = icMaxStatus(rv, icValidateWarning);
  }
 
  switch (m_nFunctionType) {
    case 0x0000:
      if (!m_params || m_nParameters < 3) {
        sReport += icMsgValidateCriticalError;
        sReport += " Tone mapping function has invalid parameters.\n";
        rv = icMaxStatus(rv, icValidateCriticalError);
      }
      else if (m_nParameters > 3) {
        sReport += icMsgValidateWarning;
        sReport += " tone mapping function has too many formulaCurveSegment parameters.\n";
        rv = icMaxStatus(rv, icValidateWarning);
      }
      break;
 
    default:
    {
      icChar buf[128];
      sReport += icMsgValidateCriticalError;
      sprintf(buf, " tone mapping function uses unknown function type %d\n", m_nFunctionType);
      sReport += buf;
      rv = icMaxStatus(rv, icValidateCriticalError);
    }
  }
 
  return rv;
}
 
 
 
 
/**
 ******************************************************************************
 * Name: CIccMpeToneMap::CIccMpeToneMap
 *
 * Purpose:
 *
 * Args:
 *
 * Return:
 ******************************************************************************/
CIccMpeToneMap::CIccMpeToneMap(icUInt16Number nOutputChannels/*=0*/)
{
  m_nReserved = 0;
  m_nInputChannels = nOutputChannels + 1;
  m_nOutputChannels = nOutputChannels;
  m_pLumCurve = NULL;
  m_pToneFuncs = NULL;
 
  SetNumOutputChannels(nOutputChannels);
 
}
 
/**
 ******************************************************************************
 * Name: CIccMpeToneMap::CIccMpeToneMap
 *
 * Purpose:
 *
 * Args:
 *
 * Return:
 ******************************************************************************/
CIccMpeToneMap::CIccMpeToneMap(const CIccMpeToneMap& toneMap)
{
  m_nReserved = toneMap.m_nReserved;
 
  m_nInputChannels = toneMap.m_nInputChannels;
  m_nOutputChannels = toneMap.m_nOutputChannels;
 
  if (toneMap.m_pLumCurve)
    m_pLumCurve = toneMap.m_pLumCurve->NewCopy();
 
  m_nFunc = toneMap.m_nFunc;
  if (toneMap.m_pToneFuncs)
    m_pToneFuncs = toneMap.CopyToneFuncs();
  else
    m_pToneFuncs = NULL;
}
 
/**
 ******************************************************************************
 * Name: &CIccMpeToneMap::operator=
 *
 * Purpose:
 *
 * Args:
 *
 * Return:
 ******************************************************************************/
CIccMpeToneMap& CIccMpeToneMap::operator=(const CIccMpeToneMap& toneMap)
{
  if (&toneMap == this)
    return *this;
 
  m_nReserved = toneMap.m_nReserved;
 
  m_nInputChannels = toneMap.m_nInputChannels;
  m_nOutputChannels = toneMap.m_nOutputChannels;
 
  if (m_pLumCurve)
    delete m_pLumCurve;
 
  if (toneMap.m_pLumCurve)
    m_pLumCurve = toneMap.m_pLumCurve->NewCopy();
  else
    m_pLumCurve = NULL;
 
  ClearToneFuncs();
 
  m_nFunc = toneMap.m_nFunc;
  if (toneMap.m_pToneFuncs)
    m_pToneFuncs = toneMap.CopyToneFuncs();
  else
    m_pToneFuncs = NULL;
 
  return *this;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeToneMap::~CIccMpeToneMap
 *
 * Purpose:
 *
 * Args:
 *
 * Return:
 ******************************************************************************/
CIccMpeToneMap::~CIccMpeToneMap()
{
  if (m_pLumCurve)
    delete m_pLumCurve;
 
  ClearToneFuncs();
}
 
/**
 ******************************************************************************
 * Name: CIccMpeToneMap::CopyToneFuncs
 *
 * Purpose:
 *
 * Args:
 *
 * Return:
 ******************************************************************************/
CIccToneMapFunc** CIccMpeToneMap::CopyToneFuncs() const
{
  if (!m_pToneFuncs || m_nOutputChannels==0)
    return NULL;
 
  CIccToneMapFunc** rv = (CIccToneMapFunc**)calloc(m_nOutputChannels, sizeof(CIccToneMapFunc*));
  if (!rv)
    return NULL;
 
  rv[0] = m_pToneFuncs[0] ? m_pToneFuncs[0]->NewCopy() : NULL;
 
  int j;
  for (int i = 1; i < m_nOutputChannels; i++) {
    for (j = 0; j < i; j++)
      if (m_pToneFuncs[j] == m_pToneFuncs[i])
        break;
    if (j < i)
      rv[i] = rv[j];
    else
      rv[i] = m_pToneFuncs[i] ? m_pToneFuncs[i]->NewCopy() : NULL;
  }
 
  return rv;
}
 
 
/**
 ******************************************************************************
 * Name: CIccMpeToneMap::ClearToneFuncs
 *
 * Purpose:
 *
 * Args:
 *
 * Return:
 ******************************************************************************/
void CIccMpeToneMap::ClearToneFuncs()
{
  if (m_pToneFuncs) {
 
    if (m_nOutputChannels) {
      if (m_pToneFuncs[0])
        delete m_pToneFuncs[0];
 
      int j;
      for (int i = 1; i < m_nOutputChannels; i++) {
        for (j = 0; j < i; j++)
          if (m_pToneFuncs[j] == m_pToneFuncs[i])
            break;
        if (j == i && m_pToneFuncs[i])
          delete m_pToneFuncs[i];
      }
    }
 
    free(m_pToneFuncs);
    m_pToneFuncs = NULL;
  }
}
 
/**
 ******************************************************************************
 * Name: CIccMpeToneMap::SetLumCurve
 *
 * Purpose:
 *
 * Args:
 *
 * Return:
 ******************************************************************************/
void CIccMpeToneMap::SetLumCurve(CIccCurveSetCurve *pLumCurve)
{
  if (pLumCurve == m_pLumCurve)
    return;
 
  if (m_pLumCurve)
    delete m_pLumCurve;
 
  m_pLumCurve = pLumCurve;
}
 
 
/**
 ******************************************************************************
 * Name: CIccMpeToneMap::SetNumOutputChannels
 *
 * Purpose:
 *
 * Args:
 *
 * Return:
 ******************************************************************************/
void CIccMpeToneMap::SetNumOutputChannels(icUInt16Number nVectorSize)
{
  ClearToneFuncs();
  m_nOutputChannels = nVectorSize;
  if (m_nOutputChannels) {
    m_pToneFuncs = (CIccToneMapFunc**)calloc(m_nOutputChannels, sizeof(CIccToneMapFunc*));
  }
  m_nFunc = 0;
}
 
 
/**
 ******************************************************************************
 * Name: CIccMpeToneMap::PushBackFunc
 *
 * Purpose:
 *
 * Args:
 *
 * Return:
 ******************************************************************************/
bool CIccMpeToneMap::Insert(CIccToneMapFunc* pToneMapFunc)
{
  if (!m_pToneFuncs || m_nFunc >= m_nOutputChannels)
    return false;
 
  m_pToneFuncs[m_nFunc++] = pToneMapFunc;
 
  return true;
}
 
 
/**
 ******************************************************************************
 * Name: CIccMpeToneMap::Describe
 *
 * Purpose:
 *
 * Args:
 *
 * Return:
 ******************************************************************************/
void CIccMpeToneMap::Describe(std::string& sDescription, int nVerboseness)
{
  char buf[256];
 
  sprintf(buf, "BEGIN_TONE_MAP %d\n", m_nOutputChannels);
  sDescription += buf;
 
  if (m_pLumCurve) {
    sDescription += "LUM_FUNCTION:\n";
    m_pLumCurve->Describe(sDescription, nVerboseness);
  }
 
  if (m_pToneFuncs) {
    for (int i = 0; i < m_nOutputChannels; i++) {
      sprintf(buf, "TONE_FUNCTION_%d:\n", i);
      sDescription += buf;
      if (m_pToneFuncs[i])
        m_pToneFuncs[i]->Describe(sDescription, nVerboseness);
    }
  }
 
  sDescription += "END_TINT_ARRAY\n";
}
 
/**
 ******************************************************************************
 * Name: CIccMpeToneMap::Read
 *
 * Purpose:
 *
 * Args:
 *
 * Return:
 ******************************************************************************/
bool CIccMpeToneMap::Read(icUInt32Number size, CIccIO* pIO)
{
  if (m_pLumCurve)
    delete m_pLumCurve;
  m_pLumCurve = NULL;
 
  ClearToneFuncs();
 
  icElemTypeSignature sig;
 
  icUInt32Number startPos = pIO->Tell();
 
  icUInt32Number headerSize = sizeof(icElemTypeSignature) +
    sizeof(icUInt32Number) +
    sizeof(icUInt16Number) +
    sizeof(icUInt16Number) +
    sizeof(icUInt32Number)*2;
 
  if (headerSize + sizeof(icTagTypeSignature) > size)
    return false;
 
  if (!pIO) {
    return false;
  }
 
  icUInt16Number nInputChannels, nOutputChannels;
 
  if (!pIO->Read32(&sig))
    return false;
 
  if (!pIO->Read32(&m_nReserved))
    return false;
 
  if (!pIO->Read16(&nInputChannels))
    return false;
 
  if (!pIO->Read16(&nOutputChannels))
    return false;
 
  if (nInputChannels != nOutputChannels+1 || !nOutputChannels)
    return false;
 
  m_nInputChannels = nInputChannels;
  SetNumOutputChannels(nOutputChannels);
 
  //setup to read luminance curve based on position
  icPositionNumber lumPos;
  if (!pIO->Read32(&lumPos.offset) ||
      !pIO->Read32(&lumPos.size))
    return false;
 
  icUInt32Number curPos = pIO->Tell();
 
  //We need signature of curve type to construct and read it
  icCurveElemSignature curveSig;
  if (lumPos.size < sizeof(curveSig))
    return false;
 
  pIO->Seek(startPos + lumPos.offset, icSeekSet);
  if (!pIO->Read32(&curveSig))
    return false;
  pIO->Seek(startPos + lumPos.offset, icSeekSet);
  
  CIccCurveSetCurve* pCurve = CIccCurveSetCurve::Create(curveSig);
 
  if (!pCurve) {
    return false;
  }
 
  if (!pCurve->Read(lumPos.size, pIO)) {
    return false;
  }
 
  SetLumCurve(pCurve);
 
  //setup positions to get the output channel mapping functions
  pIO->Seek(curPos, icSeekSet);
 
  headerSize += nOutputChannels + sizeof(icUInt32Number)*2;
 
  if (headerSize + sizeof(icTagTypeSignature) > size)
    return false;
 
  if (!icValidTagPos(lumPos, headerSize, size))
    return false;
 
  icPositionNumber* funcPos = new icPositionNumber[m_nOutputChannels];
  if (!funcPos)
    return false;
 
  //read output tone mapping function positions
  int j;
  for (int i = 0; i < m_nOutputChannels; i++) {
    if (!pIO->Read32(&funcPos[i].offset) ||
      !pIO->Read32(&funcPos[i].size) ||
      !icValidTagPos(funcPos[i], headerSize, size)) {
      delete[] funcPos;
      return false;
    }
 
    if (!icValidOverlap(lumPos, funcPos[i], false)) {
      delete[] funcPos;
      return false;
    }
    for (j = 0; j < i; j++) {
      if (!icValidOverlap(funcPos[j], funcPos[i], true)) {
        delete[] funcPos;
        return false;
      }
    }
  }
 
  //Read tone mapping functions making copies as we go
  m_pToneFuncs[0] = NewToneMapFunc();
  if (!m_pToneFuncs[0]) {
    delete[] funcPos;
    return false;
  }
  pIO->Seek(startPos + funcPos[0].offset, icSeekSet);
  if (!m_pToneFuncs[0]->Read(funcPos[0].size, pIO)) {
    delete[] funcPos;
    return false;
  }
 
  for (int i = 1; i < m_nOutputChannels; i++) {
    for (j = 0; j < i; i++) {
      if (funcPos[j].offset == funcPos[i].offset)
        break;
    }
    if (j != i)
      m_pToneFuncs[i] = m_pToneFuncs[j];
    else {
      m_pToneFuncs[i] = NewToneMapFunc();
      if (!m_pToneFuncs[i]) {
        delete[] funcPos;
        return false;
      }
      pIO->Seek(startPos + funcPos[i].offset, icSeekSet);
      if (!m_pToneFuncs[i]->Read(funcPos[i].size, pIO)) {
        delete[] funcPos;
        return false;
      }
    }
  }
 
  //keep track of functions read in
  m_nFunc = m_nOutputChannels;
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeToneMap::Write
 *
 * Purpose:
 *
 * Args:
 *
 * Return:
 ******************************************************************************/
bool CIccMpeToneMap::Write(CIccIO* pIO)
{
  if (!m_pLumCurve || !m_pToneFuncs || !m_nOutputChannels)
    return false;
 
  icElemTypeSignature sig = GetType();
 
  if (!pIO)
    return false;
 
  icUInt32Number nTagStartPos = pIO->Tell();
 
  if (!pIO->Write32(&sig))
    return false;
 
  if (!pIO->Write32(&m_nReserved))
    return false;
 
  if (!pIO->Write16(&m_nInputChannels))
    return false;
 
  if (!pIO->Write16(&m_nOutputChannels))
    return false;
 
  icUInt32Number lumOffset = pIO->Tell();
 
  //Reserve position entry for Luminance Curve
  icUInt32Number zero = 0;
  if (!pIO->Write32(&zero) ||
      !pIO->Write32(&zero))
    return false;
 
  //Reserve Position entry for mapping functions
  for (int i = 0; i < m_nOutputChannels; i++) {
    if (!pIO->Write32(&zero) ||
      !pIO->Write32(&zero))
      return false;
  }
 
  //write out luminance curve
  icPositionNumber lumPos;
  lumPos.offset = pIO->Tell()- nTagStartPos;
 
  if (!m_pLumCurve->Write(pIO))
    return false;
 
  lumPos.size = pIO->Tell() - (lumPos.offset + nTagStartPos);
 
  //Keep track of tone function positions
  icPositionNumber* funcPos = new icPositionNumber[m_nOutputChannels];
  if (!funcPos)
    return false;
 
  //write out first tone function
  int j;
  funcPos[0].offset = pIO->Tell() - nTagStartPos;
  if (!m_pToneFuncs[0]->Write(pIO)) {
    delete[] funcPos;
    return false;
  }
  funcPos[0].size = pIO->Tell() - (funcPos[0].offset + nTagStartPos);
 
  //write out additional non-copied tone functions
  for (int i = 1; i < m_nOutputChannels; i++) {
    for (j = 0; j < i; j++)
      if (m_pToneFuncs[j] == m_pToneFuncs[i])
        break;
    if (j < i) {
      funcPos[i] = funcPos[j];
    }
    else {
      funcPos[i].offset = pIO->Tell() - nTagStartPos;
      if (!m_pToneFuncs[i]->Write(pIO)) {
        delete[] funcPos;
        return false;
      }
      funcPos[i].size = pIO->Tell() - (funcPos[i].offset + nTagStartPos);
    }
  }
 
  //Everything but positions is written so make sure we end on 32 bit boundary
  pIO->Align32();
  icUInt32Number endOffset = pIO->Tell();
 
  //write out luminance curve position
  pIO->Seek(lumOffset, icSeekSet);
 
  if (!pIO->Write32(&lumPos.offset) ||
      !pIO->Write32(&lumPos.size)) {
    delete[] funcPos;
    return false;
  }
 
  //write out tone mapping function positions
  for (int i = 0; i < m_nOutputChannels; i++) {
    if (!pIO->Write32(&funcPos[i].offset) ||
      !pIO->Write32(&funcPos[i].size)) {
      delete[] funcPos;
      return false;
    }
  }
  delete[] funcPos;
 
  //Go back to end and we are done
  pIO->Seek(endOffset, icSeekSet);
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeToneMap::Begin
 *
 * Purpose:
 *
 * Args:
 *
 * Return:
 ******************************************************************************/
bool CIccMpeToneMap::Begin(icElemInterp nInterp, CIccTagMultiProcessElement* pMPE)
{
  if (!m_pLumCurve || !m_pToneFuncs)
    return false;
 
  if (!m_pLumCurve->Begin(nInterp, pMPE))
    return false;
 
  if (!m_pToneFuncs[0] || !m_pToneFuncs[0]->Begin())
    return false;
 
  int j;
  for (int i=1; i<m_nOutputChannels; i++) {
    for (j = 0; j < i; j++)
      if (m_pToneFuncs[j] == m_pToneFuncs[i])
        break;
    if (j==i && (!m_pToneFuncs[i] || !m_pToneFuncs[i]->Begin()))
      return false;
  }
  return true;
}
 
 
/**
 ******************************************************************************
 * Name: CIccMpeToneMap::Apply
 *
 * Purpose:
 *
 * Args:
 *
 * Return:
 ******************************************************************************/
void CIccMpeToneMap::Apply(CIccApplyMpe* pApply, icFloatNumber* pDestPixel, const icFloatNumber* pSrcPixel) const
{
  icFloatNumber lum = m_pLumCurve->Apply(pSrcPixel[m_nOutputChannels]);
  for (int i = 0; i < m_nOutputChannels; i++) {
    pDestPixel[i] = m_pToneFuncs[i]->Apply(lum, pSrcPixel[i]);
  }
}
 
/**
 ******************************************************************************
 * Name: CIccMpeToneMap::Validate
 *
 * Purpose:
 *
 * Args:
 *
 * Return:
 ******************************************************************************/
icValidateStatus CIccMpeToneMap::Validate(std::string sigPath, std::string& sReport, const CIccTagMultiProcessElement* pMPE/*=NULL*/, const CIccProfile* pProfile/*=NULL*/) const
{
  std::string mpeSigPath = sigPath + icGetSigPath(GetType());
  icValidateStatus rv = CIccMultiProcessElement::Validate(sigPath, sReport, pMPE, pProfile);
 
  if (m_nInputChannels != m_nOutputChannels + 1) {
    CIccInfo Info;
    std::string sSigPathName = Info.GetSigPathName(mpeSigPath);
 
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " -  Bad number of input channels!\n";
    return icValidateCriticalError;
  }
 
  if (!m_nOutputChannels) {
    CIccInfo Info;
    std::string sSigPathName = Info.GetSigPathName(mpeSigPath);
 
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " -  Bad number of output channels!\n";
    return icValidateCriticalError;
  }
 
  if (!m_pLumCurve) {
    CIccInfo Info;
    std::string sSigPathName = Info.GetSigPathName(mpeSigPath);
 
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " -  Has no luminance mapping curve!\n";
    return icValidateCriticalError;
  }
  else {
    rv = icMaxStatus(rv, m_pLumCurve->Validate(mpeSigPath, sReport, pMPE, pProfile));
  }
 
  if (!m_pToneFuncs) {
    CIccInfo Info;
    std::string sSigPathName = Info.GetSigPathName(mpeSigPath);
 
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " -  Has no tone mapping functions!\n";
    return icValidateCriticalError;
  }
  else {
    for (int i = 0; i < m_nOutputChannels; i++) {
      std::string funcReport;
      char buf[20];
      sprintf(buf, "#%d", i);
      if (!m_pToneFuncs[i]) {
        CIccInfo Info;
        std::string sSigPathName = Info.GetSigPathName(mpeSigPath);
 
        sReport += icMsgValidateCriticalError;
        sReport += sSigPathName;
        sReport += std::string(" -  Missing Tone function #") + buf + "\n";
 
        return icValidateCriticalError;
      }
 
      icValidateStatus crv = m_pToneFuncs[i]->Validate(funcReport);
      if (crv != icValidateOK) {
        CIccInfo Info;
        std::string sSigPathName = Info.GetSigPathName(mpeSigPath);
 
        sReport += icMsgValidateInformation;
        sReport += sSigPathName;
        sReport += std::string(" - Tone Function #") + buf + " status:\n" + funcReport;
 
        rv = icMaxStatus(rv, crv);
      }
    }
  }
 
  return rv;
}
 
 
 
/**
 ******************************************************************************
 * Name: CIccMpeMatrix::CIccMpeMatrix
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeMatrix::CIccMpeMatrix()
{
  m_nReserved = 0;
  m_nInputChannels = m_nOutputChannels = 0;
  m_size = 0;
  m_pMatrix = NULL;
  m_pConstants = NULL;
}
 
 
/**
 ******************************************************************************
 * Name: CIccMpeMatrix::CIccMpeMatrix
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeMatrix::CIccMpeMatrix(const CIccMpeMatrix &matrix)
{
  m_nReserved = matrix.m_nReserved;
 
  m_nInputChannels = matrix.m_nInputChannels;
  m_nOutputChannels = matrix.m_nOutputChannels;
 
  m_size = matrix.m_size;
  if(matrix.m_pMatrix) {
    int num = m_size * sizeof(icFloatNumber);
    m_pMatrix = (icFloatNumber*)malloc(num);
    memcpy(m_pMatrix, matrix.m_pMatrix, num);
  }
  else
    m_pMatrix = NULL;
 
  if (matrix.m_pConstants) {
    int num = m_nOutputChannels*sizeof(icFloatNumber);
    m_pConstants = (icFloatNumber*)malloc(num);
    memcpy(m_pConstants, matrix.m_pConstants, num);
  }
  else
    m_pConstants = NULL;
 
  m_bApplyConstants = true;
}
 
/**
 ******************************************************************************
 * Name: &CIccMpeMatrix::operator=
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeMatrix &CIccMpeMatrix::operator=(const CIccMpeMatrix &matrix)
{
  m_nReserved = matrix.m_nReserved;
 
  m_nInputChannels = matrix.m_nInputChannels;
  m_nOutputChannels = matrix.m_nOutputChannels;
 
  if (m_pMatrix)
    free(m_pMatrix);
 
  m_size = matrix.m_size;
  if (matrix.m_pMatrix) {
    int num = m_size * sizeof(icFloatNumber);
    m_pMatrix = (icFloatNumber*)malloc(num);
    memcpy(m_pMatrix, matrix.m_pMatrix, num);
  }
  else
    m_pMatrix = NULL;
 
  if (m_pConstants)
    free(m_pConstants);
 
  if (matrix.m_pConstants) {
    int num = m_nOutputChannels*sizeof(icFloatNumber);
    m_pConstants = (icFloatNumber*)malloc(num);
    memcpy(m_pConstants, matrix.m_pConstants, num);
  }
  else
    m_pConstants = NULL;
 
  m_bApplyConstants = matrix.m_bApplyConstants;
 
  return *this;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeMatrix::~CIccMpeMatrix
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeMatrix::~CIccMpeMatrix()
{
  if (m_pMatrix)
    free(m_pMatrix);
 
  if (m_pConstants)
    free(m_pConstants);
}
 
/**
 ******************************************************************************
 * Name: CIccMpeMatrix::SetSize
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccMpeMatrix::SetSize(icUInt16Number nInputChannels, icUInt16Number nOutputChannels, bool bUseConstants)
{
  if (m_pMatrix) {
    free(m_pMatrix);
    m_pMatrix = NULL;
  }
  if (m_pConstants) {
    free(m_pConstants);
    m_pConstants = NULL;
  }
 
  m_size = (icUInt32Number)nInputChannels * nOutputChannels;
 
  if (m_size) {
    m_pMatrix = (icFloatNumber*)calloc(m_size, sizeof(icFloatNumber));
    if (!m_pMatrix)
      return false;
  }
 
  if (bUseConstants) {
    m_pConstants = (icFloatNumber*)calloc(nOutputChannels, sizeof(icFloatNumber));
 
    if (!m_pConstants)
      return false;
  }
 
  m_nInputChannels = nInputChannels;
  m_nOutputChannels = nOutputChannels;
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeMatrix::Describe
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccMpeMatrix::Describe(std::string &sDescription, int nVerboseness)
{
  icChar buf[81];
  int i, j;
  icFloatNumber *data = m_pMatrix;
 
  sprintf(buf, "BEGIN_ELEM_MATRIX %d %d\n", m_nInputChannels, m_nOutputChannels);
  sDescription += buf;
 
  for (j=0; j<m_nOutputChannels; j++) {
    if (data) {
      for (i=0; i<m_nInputChannels; i++) {
        if (i)
          sDescription += " ";
        sprintf(buf, "%12.8lf", data[i]);
        sDescription += buf;
      }
      if (m_pConstants) {
        sprintf(buf, "  +  %12.8lf\n", m_pConstants[j]);
        sDescription += buf;
      }
      data += i;
    }
    else {
      sprintf(buf, "ZeroRow  +  %12.8lf\n", m_pConstants[j]);
      sDescription += buf;
    }
  }
  sDescription += "END_ELEM_MATRIX\n";
}
 
/**
 ******************************************************************************
 * Name: CIccMpeMatrix::Read
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccMpeMatrix::Read(icUInt32Number size, CIccIO *pIO)
{
  icElemTypeSignature sig;
  
  icUInt32Number headerSize = sizeof(icElemTypeSignature) + 
    sizeof(icUInt32Number) + 
    sizeof(icUInt16Number) + 
    sizeof(icUInt16Number);
 
  if (headerSize > size)
    return false;
 
  icUInt32Number dataSize = size - headerSize;
 
  if (!pIO) {
    return false;
  }
 
  icUInt16Number nInputChannels, nOutputChannels;
 
  if (!pIO->Read32(&sig))
    return false;
 
  if (!pIO->Read32(&m_nReserved))
    return false;
 
  if (!pIO->Read16(&nInputChannels))
    return false;
 
  if (!pIO->Read16(&nOutputChannels))
    return false;
 
  if (dataSize >= (icUInt32Number)nInputChannels * nOutputChannels * sizeof(icFloatNumber) &&
      dataSize < ((icUInt32Number)nInputChannels+1) * nOutputChannels * sizeof(icFloatNumber)) {
    //Matrix with no constants
    if (!SetSize(nInputChannels, nOutputChannels, false))
      return false;
 
    if (!m_pMatrix)
      return false;
 
    //Read Matrix data
    if (pIO->ReadFloat32Float(m_pMatrix, m_size) != (icInt32Number)m_size)
      return false;
  }
  else if (dataSize < (icUInt32Number)nInputChannels * nOutputChannels *sizeof(icFloatNumber) &&
           dataSize >= (icUInt32Number)nOutputChannels * sizeof(icFloatNumber)) {
    //Constants with no matrix
    if (!SetSize(0, nOutputChannels))
      return false;
 
    m_nInputChannels = nInputChannels;
 
    //Read Constant data
    if (pIO->ReadFloat32Float(m_pConstants, m_nOutputChannels)!=m_nOutputChannels)
      return false;
  }
  else {
    if ((icUInt32Number)nInputChannels*nOutputChannels > dataSize ||
        ((icUInt32Number)nInputChannels*nOutputChannels + nOutputChannels) > dataSize ||
        ((icUInt32Number)nInputChannels*nOutputChannels + nOutputChannels) * sizeof(icFloat32Number) > dataSize)
      return false;
 
    //Matrix with constants
    if (!SetSize(nInputChannels, nOutputChannels))
      return false;
 
    if (!m_pMatrix) 
      return false;
 
    if ((m_size + nOutputChannels)*sizeof(icFloat32Number) > dataSize)
      return false;
 
    //Read Matrix data
    if (pIO->ReadFloat32Float(m_pMatrix, m_size)!=(icInt32Number)m_size)
      return false;
 
    //Read Constant data
    if (pIO->ReadFloat32Float(m_pConstants, m_nOutputChannels)!=m_nOutputChannels)
      return false;
  }
  
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeMatrix::Write
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccMpeMatrix::Write(CIccIO *pIO)
{
  icElemTypeSignature sig = GetType();
 
  if (!pIO)
    return false;
 
  if (!pIO->Write32(&sig))
    return false;
 
  if (!pIO->Write32(&m_nReserved))
    return false;
 
  if (!pIO->Write16(&m_nInputChannels))
    return false;
 
  if (!pIO->Write16(&m_nOutputChannels))
    return false;
 
  if (m_pMatrix) {
    if (pIO->WriteFloat32Float(m_pMatrix, m_size)!=(icInt32Number)m_size)
      return false;
  }
 
  //Write Constant data
  if (m_pConstants) {
    if (pIO->WriteFloat32Float(m_pConstants, m_nOutputChannels)!=m_nOutputChannels)
      return false;
  }
  else {
    //Write zero constants because spec says it they have to be there
    icFloat32Number zero = 0;
    int i;
    for (i = 0; i < m_nOutputChannels; i++) {
      if (!pIO->WriteFloat32Float(&zero, 1))
        return false;
    }
  }
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeMatrix::Begin
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccMpeMatrix::Begin(icElemInterp nInterp, CIccTagMultiProcessElement *pMPE)
{
  m_bApplyConstants = false;
  if (m_pConstants) {
    int i;
    for (i = 0; i < m_nOutputChannels; i++) {
      if (icNotZero(m_pConstants[i])) {
        m_bApplyConstants = true;
        break;
      }
    }
  }
 
  if (m_nInputChannels==3 && m_nOutputChannels==3)
    m_type = ic3x3Matrix;
  else if (m_nInputChannels==3 && m_nOutputChannels==4)
    m_type = ic3x4Matrix;
  else if (m_nInputChannels==4 && m_nOutputChannels==3)
    m_type = ic4x3Matrix;
  else if (m_nInputChannels==4 && m_nOutputChannels==4)
    m_type = ic4x4Matrix;
  else
    m_type = icOtherMatrix;
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeMatrix::Apply
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccMpeMatrix::Apply(CIccApplyMpe *pApply, icFloatNumber *dstPixel, const icFloatNumber *srcPixel) const
{
  icFloatNumber *data = m_pMatrix;
  if (data) {
    if (m_bApplyConstants) {
      switch (m_type) {
        case ic3x3Matrix:
          *dstPixel++ = data[ 0]*srcPixel[0] + data[ 1]*srcPixel[1] + data[ 2]*srcPixel[2] + m_pConstants[0];
          *dstPixel++ = data[ 3]*srcPixel[0] + data[ 4]*srcPixel[1] + data[ 5]*srcPixel[2] + m_pConstants[1];
          *dstPixel   = data[ 6]*srcPixel[0] + data[ 7]*srcPixel[1] + data[ 8]*srcPixel[2] + m_pConstants[2];
          break;
 
        case ic3x4Matrix:
          *dstPixel++ = data[ 0]*srcPixel[0] + data[ 1]*srcPixel[1] + data[ 2]*srcPixel[2] + m_pConstants[0];
          *dstPixel++ = data[ 3]*srcPixel[0] + data[ 4]*srcPixel[1] + data[ 5]*srcPixel[2] + m_pConstants[1];
          *dstPixel++ = data[ 6]*srcPixel[0] + data[ 7]*srcPixel[1] + data[ 8]*srcPixel[2] + m_pConstants[2];
          *dstPixel   = data[ 9]*srcPixel[0] + data[10]*srcPixel[1] + data[11]*srcPixel[2] + m_pConstants[3];
          break;
 
        case ic4x3Matrix:
          *dstPixel++ = data[ 0]*srcPixel[0] + data[ 1]*srcPixel[1] + data[ 2]*srcPixel[2] + data[ 3]*srcPixel[3] + m_pConstants[0];
          *dstPixel++ = data[ 4]*srcPixel[0] + data[ 5]*srcPixel[1] + data[ 6]*srcPixel[2] + data[ 7]*srcPixel[3] + m_pConstants[1];
          *dstPixel   = data[ 8]*srcPixel[0] + data[ 9]*srcPixel[1] + data[10]*srcPixel[2] + data[11]*srcPixel[3] + m_pConstants[2];
          break;
 
        case ic4x4Matrix:
          *dstPixel++ = data[ 0]*srcPixel[0] + data[ 1]*srcPixel[1] + data[ 2]*srcPixel[2] + data[ 3]*srcPixel[3] + m_pConstants[0];
          *dstPixel++ = data[ 4]*srcPixel[0] + data[ 5]*srcPixel[1] + data[ 6]*srcPixel[2] + data[ 7]*srcPixel[3] + m_pConstants[1];
          *dstPixel++ = data[ 8]*srcPixel[0] + data[ 9]*srcPixel[1] + data[10]*srcPixel[2] + data[11]*srcPixel[3] + m_pConstants[2];
          *dstPixel   = data[12]*srcPixel[0] + data[13]*srcPixel[1] + data[14]*srcPixel[2] + data[15]*srcPixel[3] + m_pConstants[3];
          break;
      
        case icOtherMatrix:
        default:
          {
            int i, j;
 
            for (j=0; j<m_nOutputChannels; j++) {
              *dstPixel = m_pConstants[j];
 
              for (i=0; i<m_nInputChannels; i++) {
                *dstPixel += data[i]*srcPixel[i];
              }
 
              data += i;
              dstPixel++;
            }
          }
          break;
      }
    }
    else {
      switch (m_type) {
        case ic3x3Matrix:
          *dstPixel++ = data[0] * srcPixel[0] + data[1] * srcPixel[1] + data[2] * srcPixel[2];
          *dstPixel++ = data[3] * srcPixel[0] + data[4] * srcPixel[1] + data[5] * srcPixel[2];
          *dstPixel = data[6] * srcPixel[0] + data[7] * srcPixel[1] + data[8] * srcPixel[2];
          break;
 
        case ic3x4Matrix:
          *dstPixel++ = data[0] * srcPixel[0] + data[1] * srcPixel[1] + data[2] * srcPixel[2];
          *dstPixel++ = data[3] * srcPixel[0] + data[4] * srcPixel[1] + data[5] * srcPixel[2];
          *dstPixel++ = data[6] * srcPixel[0] + data[7] * srcPixel[1] + data[8] * srcPixel[2];
          *dstPixel = data[9] * srcPixel[0] + data[10] * srcPixel[1] + data[11] * srcPixel[2];
          break;
 
        case ic4x3Matrix:
          *dstPixel++ = data[0] * srcPixel[0] + data[1] * srcPixel[1] + data[2] * srcPixel[2] + data[3] * srcPixel[3];
          *dstPixel++ = data[4] * srcPixel[0] + data[5] * srcPixel[1] + data[6] * srcPixel[2] + data[7] * srcPixel[3];
          *dstPixel = data[8] * srcPixel[0] + data[9] * srcPixel[1] + data[10] * srcPixel[2] + data[11] * srcPixel[3];
          break;
 
        case ic4x4Matrix:
          *dstPixel++ = data[0] * srcPixel[0] + data[1] * srcPixel[1] + data[2] * srcPixel[2] + data[3] * srcPixel[3];
          *dstPixel++ = data[4] * srcPixel[0] + data[5] * srcPixel[1] + data[6] * srcPixel[2] + data[7] * srcPixel[3];
          *dstPixel++ = data[8] * srcPixel[0] + data[9] * srcPixel[1] + data[10] * srcPixel[2] + data[11] * srcPixel[3];
          *dstPixel = data[12] * srcPixel[0] + data[13] * srcPixel[1] + data[14] * srcPixel[2] + data[15] * srcPixel[3];
          break;
 
        case icOtherMatrix:
        default:
          {
            int i, j;
 
            for (j = 0; j < m_nOutputChannels; j++) {
              *dstPixel = 0.0f;
 
              for (i = 0; i < m_nInputChannels; i++) {
                *dstPixel += data[i] * srcPixel[i];
              }
 
              data += i;
              dstPixel++;
            }
          }
          break;
      }
    }
  }
  else if (m_bApplyConstants) {
    memcpy(dstPixel, m_pConstants, m_nOutputChannels*sizeof(icFloatNumber));
  }
  else {
    memset(dstPixel, 0, m_nOutputChannels * sizeof(icFloatNumber));
  }
}
 
/**
 ******************************************************************************
 * Name: CIccMpeMatrix::Validate
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
icValidateStatus CIccMpeMatrix::Validate(std::string sigPath, std::string &sReport, const CIccTagMultiProcessElement* pMPE/*=NULL*/, const CIccProfile *pProfile/*=NULL*/) const
{
  std::string mpeSigPath = sigPath + icGetSigPath(GetType());
  icValidateStatus rv = CIccMultiProcessElement::Validate(sigPath, sReport, pMPE, pProfile);
  
  if (!m_pConstants) {
    CIccInfo Info;
    std::string sSigPathName = Info.GetSigPathName(mpeSigPath);
 
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " - Has Empty Matrix Constant data!\n";
    return icValidateCriticalError;
  }
 
  return rv;
}
 
 
static icFloatNumber NoClip(icFloatNumber v)
{
  return v;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCLUT::CIccMpeCLUT
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeCLUT::CIccMpeCLUT()
{
  m_pCLUT = NULL;
  m_nInputChannels = 0;
  m_nOutputChannels = 0;
  m_interpType  = ic1dInterp;
  m_nReserved = 0;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCLUT::CIccMpeCLUT
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeCLUT::CIccMpeCLUT(const CIccMpeCLUT &clut)
{
  if (clut.m_pCLUT)
    m_pCLUT = new CIccCLUT(*clut.m_pCLUT);
  else
    m_pCLUT = NULL;
 
  m_nReserved = clut.m_nReserved;
  m_nInputChannels = clut.m_nInputChannels;
  m_nOutputChannels = clut.m_nOutputChannels;
  m_interpType = ic1dInterp;
}
 
/**
 ******************************************************************************
 * Name: &CIccMpeCLUT::operator=
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeCLUT &CIccMpeCLUT::operator=(const CIccMpeCLUT &clut)
{
  if (m_pCLUT)
    delete m_pCLUT;
 
  if (clut.m_pCLUT)
    m_pCLUT = new CIccCLUT(*clut.m_pCLUT);
  else
    m_pCLUT = NULL;
 
  m_nReserved = clut.m_nReserved;
  m_nInputChannels = clut.m_nInputChannels;
  m_nOutputChannels = clut.m_nOutputChannels;
 
  return *this;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCLUT::~CIccMpeCLUT
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
CIccMpeCLUT::~CIccMpeCLUT()
{
  if (m_pCLUT)
    delete m_pCLUT;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCLUT::SetCLUT
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccMpeCLUT::SetCLUT(CIccCLUT *pCLUT)
{
  if (m_pCLUT)
    delete m_pCLUT;
 
  m_pCLUT = pCLUT;
  if (pCLUT) {
    pCLUT->SetClipFunc(NoClip);
    m_nInputChannels = pCLUT->GetInputDim();
    m_nOutputChannels = pCLUT->GetOutputChannels();
  }
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCLUT::Describe
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccMpeCLUT::Describe(std::string &sDescription, int nVerboseness)
{
  if (m_pCLUT) {
    m_pCLUT->DumpLut(sDescription, "ELEM_CLUT", icSigUnknownData, icSigUnknownData, nVerboseness);
  }
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCLUT::Read
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccMpeCLUT::Read(icUInt32Number size, CIccIO *pIO)
{
  icTagTypeSignature sig;
  
  icUInt32Number headerSize = sizeof(icTagTypeSignature) + 
    sizeof(icUInt32Number) + 
    sizeof(icUInt16Number) + 
    sizeof(icUInt16Number) + 
    16 * sizeof(icUInt8Number);
 
  if (headerSize > size)
    return false;
 
  icUInt32Number dataSize = size - headerSize;
 
  if (!pIO) {
    return false;
  }
 
  if (!pIO->Read32(&sig))
    return false;
 
  if (!pIO->Read32(&m_nReserved))
    return false;
 
  if (!pIO->Read16(&m_nInputChannels))
    return false;
 
  if (!pIO->Read16(&m_nOutputChannels))
    return false;
 
  icUInt8Number gridPoints[16];
 
  if (pIO->Read8(gridPoints, 16)!=16) {
    return false;
  }
 
  icUInt32Number nPoints = (icUInt32Number)m_nInputChannels * m_nOutputChannels;
 
  if (m_nInputChannels > 16 || nPoints > dataSize || nPoints * sizeof (icFloat32Number) > dataSize)
    return false;
 
  m_pCLUT = new CIccCLUT((icUInt8Number)m_nInputChannels, (icUInt16Number)m_nOutputChannels, 4);
 
  if (!m_pCLUT)
    return false;
 
  m_pCLUT->SetClipFunc(NoClip);
 
  m_pCLUT->Init(gridPoints);
 
  icFloatNumber *pData = m_pCLUT->GetData(0);
 
  if (!pData)
    return false;
 
  nPoints = m_pCLUT->NumPoints()*m_nOutputChannels;
 
  if (pIO->ReadFloat32Float(pData,nPoints)!= nPoints)
    return false;
  
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCLUT::Write
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccMpeCLUT::Write(CIccIO *pIO)
{
  icElemTypeSignature sig = GetType();
 
  if (!pIO)
    return false;
 
  if (!pIO->Write32(&sig))
    return false;
 
  if (!pIO->Write32(&m_nReserved))
    return false;
 
  if (!pIO->Write16(&m_nInputChannels))
    return false;
 
  if (!pIO->Write16(&m_nOutputChannels))
    return false;
 
  if (m_pCLUT) {
    icUInt8Number gridPoints[16];
    int i;
 
    for (i=0; i<16; i++)
      gridPoints[i] = m_pCLUT->GridPoint(i);
 
    if (pIO->Write8(gridPoints, 16)!=16)
      return false;
 
    icFloatNumber *pData = m_pCLUT->GetData(0);
    icInt32Number nPoints = m_pCLUT->NumPoints()*m_nOutputChannels;
 
    if (pIO->WriteFloat32Float(pData, nPoints) != nPoints) 
      return false;
  }
 
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCLUT::Begin
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
bool CIccMpeCLUT::Begin(icElemInterp nInterp, CIccTagMultiProcessElement *pMPE)
{
  if (!m_pCLUT)
    return false;
 
  m_pCLUT->Begin();
 
  switch (m_nInputChannels) {
  case 1:
    m_interpType = ic1dInterp;
    break;
  case 2:
    m_interpType = ic2dInterp;
    break;
  case 3:
    if (nInterp==icElemInterpTetra)
      m_interpType = ic3dInterpTetra;
    else
      m_interpType = ic3dInterp;
    break;
  case 4:
    m_interpType = ic4dInterp;
    break;
  case 5:
    m_interpType = ic5dInterp;
    break;
  case 6:
    m_interpType = ic6dInterp;
    break;
  default:
    m_interpType = icNdInterp;
    break;
  }
  return true;
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCLUT::Apply
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
void CIccMpeCLUT::Apply(CIccApplyMpe *pApply, icFloatNumber *dstPixel, const icFloatNumber *srcPixel) const
{
  const CIccCLUT *pCLUT = m_pCLUT;
 
  switch(m_interpType) {
  case ic1dInterp:
    pCLUT->Interp1d(dstPixel, srcPixel);
    break;
  case ic2dInterp:
    pCLUT->Interp2d(dstPixel, srcPixel);
    break;
  case ic3dInterpTetra:
    pCLUT->Interp3dTetra(dstPixel, srcPixel);
    break;
  case ic3dInterp:
    pCLUT->Interp3d(dstPixel, srcPixel);
    break;
  case ic4dInterp:
    pCLUT->Interp4d(dstPixel, srcPixel);
    break;
  case ic5dInterp:
    pCLUT->Interp5d(dstPixel, srcPixel);
    break;
  case ic6dInterp:
    pCLUT->Interp6d(dstPixel, srcPixel);
    break;
  case icNdInterp:
    CIccApplyMpeCLUT* pApplyCLUT = (CIccApplyMpeCLUT*)pApply;
    pCLUT->InterpND(dstPixel, srcPixel, pApplyCLUT->m_pApply);
    break;
  }
}
 
/**
 ******************************************************************************
 * Name: CIccMpeCLUT::Validate
 * 
 * Purpose: 
 * 
 * Args: 
 * 
 * Return: 
 ******************************************************************************/
icValidateStatus CIccMpeCLUT::Validate(std::string sigPath, std::string &sReport, const CIccTagMultiProcessElement* pMPE/*=NULL*/, const CIccProfile *pProfile/*=NULL*/) const
{
  std::string mpeSigPath = sigPath + icGetSigPath(GetType());
  icValidateStatus rv = CIccMultiProcessElement::Validate(sigPath, sReport, pMPE, pProfile);
 
  if (!m_pCLUT) {
    CIccInfo Info;
    std::string sSigPathName = Info.GetSigPathName(mpeSigPath);
 
    sReport += icMsgValidateCriticalError;
    sReport += sSigPathName;
    sReport += " - Has No CLUT!\n";
    return icValidateCriticalError;
  }
 
  return rv;
}
 
 
/**
******************************************************************************
* Name: CIccMpeCLUT::Begin
*
* Purpose:
*
* Args:
*
* Return:
******************************************************************************/
CIccApplyMpe* CIccMpeCLUT::GetNewApply(CIccApplyTagMpe* pApplyTag)
{
  if (!m_pCLUT) {
    return NULL;
  }
  
  CIccApplyCLUT* pApply = m_pCLUT->GetNewApply();
  if (!pApply)
    return NULL;
 
  CIccApplyMpeCLUT * rv = new CIccApplyMpeCLUT(this, pApply);
  if (!rv)
    delete pApply;
 
  return rv;
}
 
/**
**************************************************************************
* Name: CIccApplyMpeCLUT::CIccApplyMpeCLUT
*
* Purpose:
*  Constructor
**************************************************************************
*/
CIccApplyMpeCLUT::CIccApplyMpeCLUT(CIccMultiProcessElement* pElem, CIccApplyCLUT *pApply) : CIccApplyMpe(pElem)
{
  m_pApply = pApply;
}
 
 
/**
**************************************************************************
* Name: CIccApplyMpeCLUT::~CIccApplyMpeCLUT
*
* Purpose:
*  Destructor
**************************************************************************
*/
CIccApplyMpeCLUT::~CIccApplyMpeCLUT()
{
  if (m_pApply)
    delete m_pApply;
}
 
 
/**
******************************************************************************
* Name: CIccMpeExtCLUT::CIccMpeExtCLUT
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
CIccMpeExtCLUT::CIccMpeExtCLUT()
{
  m_pCLUT = NULL;
  m_nInputChannels = 0;
  m_nOutputChannels = 0;
 
  m_nReserved = 0;
  m_nReserved2 = 0;
 
  m_storageType = icValueTypeFloat32;
}
 
/**
******************************************************************************
* Name: CIccMpeExtCLUT::CIccMpeExtCLUT
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
CIccMpeExtCLUT::CIccMpeExtCLUT(const CIccMpeExtCLUT &clut) : CIccMpeCLUT(clut)
{
  m_nReserved2 = clut.m_nReserved2;
  m_storageType = clut.m_storageType;
}
 
/**
******************************************************************************
* Name: &CIccMpeExtCLUT::operator=
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
CIccMpeExtCLUT &CIccMpeExtCLUT::operator=(const CIccMpeExtCLUT &clut)
{
  CIccMpeCLUT::operator=(clut);
 
  m_nReserved2 = clut.m_nReserved2;
  m_storageType = clut.m_storageType;
 
  return *this;
}
 
/**
******************************************************************************
* Name: CIccMpeExtCLUT::Describe
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
void CIccMpeExtCLUT::Describe(std::string &sDescription, int nVerboseness)
{
  if (m_pCLUT) {
    char desc[256];
    sprintf(desc, "EXT_ELEM_CLUT(%d)", m_storageType);
 
    m_pCLUT->DumpLut(sDescription, desc, icSigUnknownData, icSigUnknownData, nVerboseness);
  }
}
 
/**
******************************************************************************
* Name: CIccMpeExtCLUT::SetStorageType
*
* Purpose:
*  Sets storage type of the data stored in the CLUT
*
* Args:
*  nStorageType is type of data to use
*
* Return:
*  true if valid storage type, false otherwise
******************************************************************************/
bool CIccMpeExtCLUT::SetStorageType(icUInt16Number nStorateType)
{
  m_storageType = nStorateType;
 
  switch (nStorateType) {
  case icValueTypeUInt8:
  case icValueTypeUInt16:
  case icValueTypeFloat16:
  case icValueTypeFloat32:
    return true;
  }
  return false;
}
 
 
/**
******************************************************************************
* Name: CIccMpeExtCLUT::Read
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
bool CIccMpeExtCLUT::Read(icUInt32Number size, CIccIO *pIO)
{
  icTagTypeSignature sig;
 
  icUInt32Number headerSize = sizeof(icTagTypeSignature) + 
    sizeof(icUInt32Number) + 
    sizeof(icUInt16Number) + 
    sizeof(icUInt16Number) + 
    sizeof(icUInt16Number) +
    sizeof(icUInt16Number) +
    16 * sizeof(icUInt8Number);
 
  if (headerSize > size)
    return false;
 
  icUInt32Number dataSize = size - headerSize;
 
  if (!pIO) {
    return false;
  }
 
  if (!pIO->Read32(&sig))
    return false;
 
  if (!pIO->Read32(&m_nReserved))
    return false;
 
  if (!pIO->Read16(&m_nInputChannels))
    return false;
 
  if (!pIO->Read16(&m_nOutputChannels))
    return false;
 
  if (!pIO->Read16(&m_storageType))
    return false;
 
  if (!pIO->Read16(&m_nReserved2))
    return false;
 
  icUInt8Number gridPoints[16];
 
  if (pIO->Read8(gridPoints, 16)!=16) {
    return false;
  }
 
  icUInt32Number nPoints = (icUInt32Number)m_nInputChannels * m_nOutputChannels;
 
  if (m_nInputChannels > 16 || nPoints > dataSize)
    return false;
  
  m_pCLUT = new CIccCLUT((icUInt8Number)m_nInputChannels, (icUInt16Number)m_nOutputChannels, 4);
 
  if (!m_pCLUT)
    return false;
 
  m_pCLUT->SetClipFunc(NoClip);
 
  if (!m_pCLUT->Init(gridPoints, dataSize, icGetStorageTypeBytes(m_storageType)))
    return false;
 
  icFloatNumber *pData = m_pCLUT->GetData(0);
 
  if (!pData)
    return false;
 
  nPoints = m_pCLUT->NumPoints()*m_nOutputChannels;
 
  if (nPoints > dataSize)
    return false;
 
  switch(m_storageType) {
    case icValueTypeUInt8:
      if (pIO->ReadUInt8Float(pData,nPoints)!= nPoints)
        return false;
      break;
 
    case icValueTypeUInt16:
      if (nPoints * 2 > dataSize)
        return false;
 
      if (pIO->ReadUInt16Float(pData,nPoints)!= nPoints)
        return false;
      break;
 
    case icValueTypeFloat16:
      if (nPoints * 2 > dataSize)
        return false;
 
      if (pIO->ReadFloat16Float(pData,nPoints)!= nPoints)
        return false;
      break;
 
    case icValueTypeFloat32:
      if (pIO->ReadFloat32Float(pData,nPoints)!= nPoints)
        return false;
      break;
 
    default:
      return false;
  }
  return true;
}
 
/**
******************************************************************************
* Name: CIccMpeExtCLUT::Write
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/
bool CIccMpeExtCLUT::Write(CIccIO *pIO)
{
  icElemTypeSignature sig = GetType();
 
  if (!pIO)
    return false;
 
  if (!pIO->Write32(&sig))
    return false;
 
  if (!pIO->Write32(&m_nReserved))
    return false;
 
  if (!pIO->Write16(&m_nInputChannels))
    return false;
 
  if (!pIO->Write16(&m_nOutputChannels))
    return false;
 
  if (!pIO->Write16(&m_storageType))
    return false;
 
  if (!pIO->Write16(&m_nReserved2))
    return false;
 
  if (m_pCLUT) {
    icUInt8Number gridPoints[16];
    int i;
 
    for (i=0; i<16; i++)
      gridPoints[i] = m_pCLUT->GridPoint(i);
 
    if (pIO->Write8(gridPoints, 16)!=16)
      return false;
 
    icFloatNumber *pData = m_pCLUT->GetData(0);
    icInt32Number nPoints = m_pCLUT->NumPoints()*m_nOutputChannels;
 
    switch(m_storageType) {
    case icValueTypeUInt8:
      if (pIO->WriteUInt8Float(pData,nPoints)!= nPoints)
        return false;
      break;
 
    case icValueTypeUInt16:
      if (pIO->WriteUInt16Float(pData,nPoints)!= nPoints)
        return false;
      break;
 
    case icValueTypeFloat16:
      if (pIO->WriteFloat16Float(pData,nPoints)!= nPoints)
        return false;
      break;
 
    case icValueTypeFloat32:
      if (pIO->WriteFloat32Float(pData,nPoints)!= nPoints)
        return false;
      break;
 
    default:
      return false;
    }
  }
 
  return true;
}
 
 
/**
******************************************************************************
* Name: CIccMpeExtCLUT::Validate
* 
* Purpose: 
* 
* Args: 
* 
* Return: 
******************************************************************************/