//---------------------------------------------------------------------------------

//

//  Little Color Management System

//  Copyright (c) 1998-2024 Marti Maria Saguer

//

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

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

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

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

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

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

//

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

// all copies or substantial portions of the Software.

//

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

// EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO

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

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

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

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

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

//

//---------------------------------------------------------------------------------

//

#include "lcms2_internal.h"

// Transformations stuff

// -----------------------------------------------------------------------

#define DEFAULT_OBSERVER_ADAPTATION_STATE 1.0

// The Context0 observer adaptation state.

_cmsAdaptationStateChunkType _cmsAdaptationStateChunk = { DEFAULT_OBSERVER_ADAPTATION_STATE };

// Init and duplicate observer adaptation state

void _cmsAllocAdaptationStateChunk(struct _cmsContext_struct* ctx, 

                                   const struct _cmsContext_struct* src)

{

    static _cmsAdaptationStateChunkType AdaptationStateChunk = { DEFAULT_OBSERVER_ADAPTATION_STATE };

    void* from;

    if (src != NULL) {

        from = src ->chunks[AdaptationStateContext];       

    }

    else {

       from = &AdaptationStateChunk;

    }

    ctx ->chunks[AdaptationStateContext] = _cmsSubAllocDup(ctx ->MemPool, from, sizeof(_cmsAdaptationStateChunkType));     

}

// Sets adaptation state for absolute colorimetric intent in the given context.  Adaptation state applies on all 

// but cmsCreateExtendedTransformTHR().  Little CMS can handle incomplete adaptation states.

cmsFloat64Number CMSEXPORT cmsSetAdaptationStateTHR(cmsContext ContextID, cmsFloat64Number d)

{

    cmsFloat64Number prev;

    _cmsAdaptationStateChunkType* ptr = (_cmsAdaptationStateChunkType*) _cmsContextGetClientChunk(ContextID, AdaptationStateContext);

    // Get previous value for return

    prev = ptr ->AdaptationState;

    // Set the value if d is positive or zero

    if (d >= 0.0) {

        ptr ->AdaptationState = d;

    }

    // Always return previous value

    return prev;

}

// The adaptation state may be defaulted by this function. If you don't like it, use the extended transform routine

cmsFloat64Number CMSEXPORT cmsSetAdaptationState(cmsFloat64Number d)

{    

    return cmsSetAdaptationStateTHR(NULL, d);

}

// -----------------------------------------------------------------------

// Alarm codes for 16-bit transformations, because the fixed range of containers there are

// no values left to mark out of gamut. 

#define DEFAULT_ALARM_CODES_VALUE {0x7F00, 0x7F00, 0x7F00, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}

_cmsAlarmCodesChunkType _cmsAlarmCodesChunk = { DEFAULT_ALARM_CODES_VALUE };

// Sets the codes used to mark out-out-gamut on Proofing transforms for a given context. Values are meant to be 

// encoded in 16 bits.

void CMSEXPORT cmsSetAlarmCodesTHR(cmsContext ContextID, const cmsUInt16Number AlarmCodesP[cmsMAXCHANNELS])

{

    _cmsAlarmCodesChunkType* ContextAlarmCodes = (_cmsAlarmCodesChunkType*) _cmsContextGetClientChunk(ContextID, AlarmCodesContext);

    _cmsAssert(ContextAlarmCodes != NULL); // Can't happen

    memcpy(ContextAlarmCodes->AlarmCodes, AlarmCodesP, sizeof(ContextAlarmCodes->AlarmCodes));    

}

// Gets the current codes used to mark out-out-gamut on Proofing transforms for the given context.

// Values are meant to be encoded in 16 bits.

void CMSEXPORT cmsGetAlarmCodesTHR(cmsContext ContextID, cmsUInt16Number AlarmCodesP[cmsMAXCHANNELS])

{

    _cmsAlarmCodesChunkType* ContextAlarmCodes = (_cmsAlarmCodesChunkType*) _cmsContextGetClientChunk(ContextID, AlarmCodesContext);

    _cmsAssert(ContextAlarmCodes != NULL); // Can't happen

    memcpy(AlarmCodesP, ContextAlarmCodes->AlarmCodes, sizeof(ContextAlarmCodes->AlarmCodes));

}

void CMSEXPORT cmsSetAlarmCodes(const cmsUInt16Number NewAlarm[cmsMAXCHANNELS])

{

    _cmsAssert(NewAlarm != NULL);

    cmsSetAlarmCodesTHR(NULL, NewAlarm);

}

void CMSEXPORT cmsGetAlarmCodes(cmsUInt16Number OldAlarm[cmsMAXCHANNELS])

{ 

    _cmsAssert(OldAlarm != NULL);

    cmsGetAlarmCodesTHR(NULL, OldAlarm);

}

// Init and duplicate alarm codes

void _cmsAllocAlarmCodesChunk(struct _cmsContext_struct* ctx, 

                              const struct _cmsContext_struct* src)

{

    static _cmsAlarmCodesChunkType AlarmCodesChunk = { DEFAULT_ALARM_CODES_VALUE };

    void* from;

    if (src != NULL) {

        from = src ->chunks[AlarmCodesContext];       

    }

    else {

       from = &AlarmCodesChunk;

    }

    ctx ->chunks[AlarmCodesContext] = _cmsSubAllocDup(ctx ->MemPool, from, sizeof(_cmsAlarmCodesChunkType));     

}

// -----------------------------------------------------------------------

// Get rid of transform resources

void CMSEXPORT cmsDeleteTransform(cmsHTRANSFORM hTransform)

{

    _cmsTRANSFORM* p = (_cmsTRANSFORM*) hTransform;

    _cmsAssert(p != NULL);

    if (p -> GamutCheck)

        cmsPipelineFree(p -> GamutCheck);

    if (p -> Lut)

        cmsPipelineFree(p -> Lut);

    if (p ->InputColorant)

        cmsFreeNamedColorList(p ->InputColorant);

    if (p -> OutputColorant)

        cmsFreeNamedColorList(p ->OutputColorant);

    if (p ->Sequence)

        cmsFreeProfileSequenceDescription(p ->Sequence);

    if (p ->UserData)

        p ->FreeUserData(p ->ContextID, p ->UserData);

    _cmsFree(p ->ContextID, (void *) p);

}

static

cmsUInt32Number PixelSize(cmsUInt32Number Format)

{

    cmsUInt32Number fmt_bytes = T_BYTES(Format);

    // For double, the T_BYTES field is zero

    if (fmt_bytes == 0)

        return sizeof(cmsUInt64Number);

    // Otherwise, it is already correct for all formats

    return fmt_bytes;

}

// Apply transform.

void CMSEXPORT cmsDoTransform(cmsHTRANSFORM  Transform,

                              const void* InputBuffer,

                              void* OutputBuffer,

                              cmsUInt32Number Size)

{

    _cmsTRANSFORM* p = (_cmsTRANSFORM*) Transform;

    cmsStride stride;

    stride.BytesPerLineIn = 0;  // Not used

    stride.BytesPerLineOut = 0;

    stride.BytesPerPlaneIn = Size * PixelSize(p->InputFormat);

    stride.BytesPerPlaneOut = Size * PixelSize(p->OutputFormat);

    p -> xform(p, InputBuffer, OutputBuffer, Size, 1, &stride);

}

// This is a legacy stride for planar

void CMSEXPORT cmsDoTransformStride(cmsHTRANSFORM  Transform,

                              const void* InputBuffer,

                              void* OutputBuffer,

                              cmsUInt32Number Size, cmsUInt32Number Stride)

{

    _cmsTRANSFORM* p = (_cmsTRANSFORM*) Transform;

    cmsStride stride;

    stride.BytesPerLineIn = 0;  

    stride.BytesPerLineOut = 0;

    stride.BytesPerPlaneIn = Stride;

    stride.BytesPerPlaneOut = Stride;

    p -> xform(p, InputBuffer, OutputBuffer, Size, 1, &stride);

}

// This is the "fast" function for plugins

void CMSEXPORT cmsDoTransformLineStride(cmsHTRANSFORM  Transform,

                              const void* InputBuffer,

                              void* OutputBuffer,

                              cmsUInt32Number PixelsPerLine,

                              cmsUInt32Number LineCount,

                              cmsUInt32Number BytesPerLineIn,

                              cmsUInt32Number BytesPerLineOut,

                              cmsUInt32Number BytesPerPlaneIn,

                              cmsUInt32Number BytesPerPlaneOut)

{

    _cmsTRANSFORM* p = (_cmsTRANSFORM*) Transform;

    cmsStride stride;

    stride.BytesPerLineIn = BytesPerLineIn;

    stride.BytesPerLineOut = BytesPerLineOut;

    stride.BytesPerPlaneIn = BytesPerPlaneIn;

    stride.BytesPerPlaneOut = BytesPerPlaneOut;

    p->xform(p, InputBuffer, OutputBuffer, PixelsPerLine, LineCount, &stride);

}

// Transform routines ----------------------------------------------------------------------------------------------------------

// Float xform converts floats. Since there are no performance issues, one routine does all job, including gamut check.

// Note that because extended range, we can use a -1.0 value for out of gamut in this case.

static

void FloatXFORM(_cmsTRANSFORM* p,

                const void* in,

                void* out, 

                cmsUInt32Number PixelsPerLine,

                cmsUInt32Number LineCount,

                const cmsStride* Stride)

{

    cmsUInt8Number* accum;

    cmsUInt8Number* output;

    cmsFloat32Number fIn[cmsMAXCHANNELS], fOut[cmsMAXCHANNELS];

    cmsFloat32Number OutOfGamut;

    size_t i, j, c, strideIn, strideOut;

    _cmsHandleExtraChannels(p, in, out, PixelsPerLine, LineCount, Stride);

    strideIn = 0;

    strideOut = 0;

    memset(fIn, 0, sizeof(fIn));

    memset(fOut, 0, sizeof(fOut));

    for (i = 0; i < LineCount; i++) {

        accum = (cmsUInt8Number*)in + strideIn;

        output = (cmsUInt8Number*)out + strideOut;

        for (j = 0; j < PixelsPerLine; j++) {

            accum = p->FromInputFloat(p, fIn, accum, Stride->BytesPerPlaneIn);

            // Any gamut check to do?

            if (p->GamutCheck != NULL) {

                // Evaluate gamut marker.

                cmsPipelineEvalFloat(fIn, &OutOfGamut, p->GamutCheck);

                // Is current color out of gamut?

                if (OutOfGamut > 0.0) {

                    _cmsAlarmCodesChunkType* ContextAlarmCodes = (_cmsAlarmCodesChunkType*)_cmsContextGetClientChunk(p->ContextID, AlarmCodesContext);

                    // Certainly, out of gamut

                    for (c = 0; c < cmsMAXCHANNELS; c++)

                        fOut[c] = ContextAlarmCodes->AlarmCodes[c] / 65535.0F;

                }

                else {

                    // No, proceed normally

                    cmsPipelineEvalFloat(fIn, fOut, p->Lut);

                }

            }

            else {

                // No gamut check at all

                cmsPipelineEvalFloat(fIn, fOut, p->Lut);

            }

            output = p->ToOutputFloat(p, fOut, output, Stride->BytesPerPlaneOut);

        }

        strideIn += Stride->BytesPerLineIn;

        strideOut += Stride->BytesPerLineOut;

    }

}

static

void NullFloatXFORM(_cmsTRANSFORM* p,

                    const void* in,

                    void* out, 

                    cmsUInt32Number PixelsPerLine,

                    cmsUInt32Number LineCount,

                    const cmsStride* Stride)

{

    cmsUInt8Number* accum;

    cmsUInt8Number* output;

    cmsFloat32Number fIn[cmsMAXCHANNELS];

    size_t i, j, strideIn, strideOut;

    _cmsHandleExtraChannels(p, in, out, PixelsPerLine, LineCount, Stride);

    strideIn = 0;

    strideOut = 0;

    memset(fIn, 0, sizeof(fIn));

    for (i = 0; i < LineCount; i++) {

           accum = (cmsUInt8Number*) in + strideIn;

           output = (cmsUInt8Number*) out + strideOut;

           for (j = 0; j < PixelsPerLine; j++) {

                  accum = p->FromInputFloat(p, fIn, accum, Stride ->BytesPerPlaneIn);

                  output = p->ToOutputFloat(p, fIn, output, Stride->BytesPerPlaneOut);

           }

           strideIn += Stride->BytesPerLineIn;

           strideOut += Stride->BytesPerLineOut;

    }

}

// 16 bit precision -----------------------------------------------------------------------------------------------------------

// Null transformation, only applies formatters. No cache

static

void NullXFORM(_cmsTRANSFORM* p,

               const void* in,

               void* out,

               cmsUInt32Number PixelsPerLine,

               cmsUInt32Number LineCount,

               const cmsStride* Stride)

{

    cmsUInt8Number* accum;

    cmsUInt8Number* output;

    cmsUInt16Number wIn[cmsMAXCHANNELS];

    size_t i, j, strideIn, strideOut;

    _cmsHandleExtraChannels(p, in, out, PixelsPerLine, LineCount, Stride);

    strideIn = 0;

    strideOut = 0;

    memset(wIn, 0, sizeof(wIn));

    for (i = 0; i < LineCount; i++) {

        accum = (cmsUInt8Number*)in + strideIn;

        output = (cmsUInt8Number*)out + strideOut;

        for (j = 0; j < PixelsPerLine; j++) {

            accum = p->FromInput(p, wIn, accum, Stride->BytesPerPlaneIn);

            output = p->ToOutput(p, wIn, output, Stride->BytesPerPlaneOut);

        }

        strideIn += Stride->BytesPerLineIn;

        strideOut += Stride->BytesPerLineOut;

    }

}

// No gamut check, no cache, 16 bits

static

void PrecalculatedXFORM(_cmsTRANSFORM* p,

                        const void* in,

                        void* out, 

                        cmsUInt32Number PixelsPerLine,

                        cmsUInt32Number LineCount,

                        const cmsStride* Stride)

{

    CMSREGISTER cmsUInt8Number* accum;

    CMSREGISTER cmsUInt8Number* output;

    cmsUInt16Number wIn[cmsMAXCHANNELS], wOut[cmsMAXCHANNELS];

    size_t i, j, strideIn, strideOut;

    _cmsHandleExtraChannels(p, in, out, PixelsPerLine, LineCount, Stride);

    strideIn = 0;

    strideOut = 0;

    memset(wIn, 0, sizeof(wIn));

    memset(wOut, 0, sizeof(wOut));

    for (i = 0; i < LineCount; i++) {

        accum = (cmsUInt8Number*)in + strideIn;

        output = (cmsUInt8Number*)out + strideOut;

        for (j = 0; j < PixelsPerLine; j++) {

            accum = p->FromInput(p, wIn, accum, Stride->BytesPerPlaneIn);

            p->Lut->Eval16Fn(wIn, wOut, p->Lut->Data);

            output = p->ToOutput(p, wOut, output, Stride->BytesPerPlaneOut);

        }

        strideIn += Stride->BytesPerLineIn;

        strideOut += Stride->BytesPerLineOut;

    }

}

// Auxiliary: Handle precalculated gamut check. The retrieval of context may be alittle bit slow, but this function is not critical.

static

void TransformOnePixelWithGamutCheck(_cmsTRANSFORM* p,

                                     const cmsUInt16Number wIn[],

                                     cmsUInt16Number wOut[])

{

    cmsUInt16Number wOutOfGamut;

    p ->GamutCheck ->Eval16Fn(wIn, &wOutOfGamut, p ->GamutCheck ->Data);

    if (wOutOfGamut >= 1) {

        cmsUInt32Number i;

        _cmsAlarmCodesChunkType* ContextAlarmCodes = (_cmsAlarmCodesChunkType*) _cmsContextGetClientChunk(p->ContextID, AlarmCodesContext);        

        for (i=0; i < p ->Lut->OutputChannels; i++) {

            wOut[i] = ContextAlarmCodes ->AlarmCodes[i];

        }

    }

    else

        p ->Lut ->Eval16Fn(wIn, wOut, p -> Lut->Data);

}

// Gamut check, No cache, 16 bits.

static

void PrecalculatedXFORMGamutCheck(_cmsTRANSFORM* p,

                                  const void* in,

                                  void* out, 

                                  cmsUInt32Number PixelsPerLine,

                                  cmsUInt32Number LineCount,

                                  const cmsStride* Stride)

{

    cmsUInt8Number* accum;

    cmsUInt8Number* output;

    cmsUInt16Number wIn[cmsMAXCHANNELS], wOut[cmsMAXCHANNELS];

    size_t i, j, strideIn, strideOut;

    _cmsHandleExtraChannels(p, in, out, PixelsPerLine, LineCount, Stride);

    strideIn = 0;

    strideOut = 0;

    memset(wIn, 0, sizeof(wIn));

    memset(wOut, 0, sizeof(wOut));

    for (i = 0; i < LineCount; i++) {

        accum = (cmsUInt8Number*)in + strideIn;

        output = (cmsUInt8Number*)out + strideOut;

        for (j = 0; j < PixelsPerLine; j++) {

            accum = p->FromInput(p, wIn, accum, Stride->BytesPerPlaneIn);

            TransformOnePixelWithGamutCheck(p, wIn, wOut);

            output = p->ToOutput(p, wOut, output, Stride->BytesPerPlaneOut);

        }

        strideIn += Stride->BytesPerLineIn;

        strideOut += Stride->BytesPerLineOut;

    }

}

// No gamut check, Cache, 16 bits,

static

void CachedXFORM(_cmsTRANSFORM* p,

                 const void* in,

                 void* out,

                 cmsUInt32Number PixelsPerLine,

                 cmsUInt32Number LineCount,

                 const cmsStride* Stride)

{

    cmsUInt8Number* accum;

    cmsUInt8Number* output;

    cmsUInt16Number wIn[cmsMAXCHANNELS], wOut[cmsMAXCHANNELS];

    _cmsCACHE Cache;

    size_t i, j, strideIn, strideOut;

    _cmsHandleExtraChannels(p, in, out, PixelsPerLine, LineCount, Stride);

    // Empty buffers for quick memcmp

    memset(wIn, 0, sizeof(wIn));

    memset(wOut, 0, sizeof(wOut));

    // Get copy of zero cache

    memcpy(&Cache, &p->Cache, sizeof(Cache));

    strideIn = 0;

    strideOut = 0;

    for (i = 0; i < LineCount; i++) {

        accum = (cmsUInt8Number*)in + strideIn;

        output = (cmsUInt8Number*)out + strideOut;

        for (j = 0; j < PixelsPerLine; j++) {

            accum = p->FromInput(p, wIn, accum, Stride->BytesPerPlaneIn);

            if (memcmp(wIn, Cache.CacheIn, sizeof(Cache.CacheIn)) == 0) {

                memcpy(wOut, Cache.CacheOut, sizeof(Cache.CacheOut));

            }

            else {

                p->Lut->Eval16Fn(wIn, wOut, p->Lut->Data);

                memcpy(Cache.CacheIn, wIn, sizeof(Cache.CacheIn));

                memcpy(Cache.CacheOut, wOut, sizeof(Cache.CacheOut));

            }

            output = p->ToOutput(p, wOut, output, Stride->BytesPerPlaneOut);

        }

        strideIn += Stride->BytesPerLineIn;

        strideOut += Stride->BytesPerLineOut;

    }

}

// All those nice features together

static

void CachedXFORMGamutCheck(_cmsTRANSFORM* p,

                           const void* in,

                           void* out, 

                           cmsUInt32Number PixelsPerLine,

                           cmsUInt32Number LineCount,

                           const cmsStride* Stride)

{

    cmsUInt8Number* accum;

    cmsUInt8Number* output;

    cmsUInt16Number wIn[cmsMAXCHANNELS], wOut[cmsMAXCHANNELS];

    _cmsCACHE Cache;

    size_t i, j, strideIn, strideOut;

    _cmsHandleExtraChannels(p, in, out, PixelsPerLine, LineCount, Stride);

    // Empty buffers for quick memcmp

    memset(wIn, 0, sizeof(wIn));

    memset(wOut, 0, sizeof(wOut));

    // Get copy of zero cache

    memcpy(&Cache, &p->Cache, sizeof(Cache));

    strideIn = 0;

    strideOut = 0;

    for (i = 0; i < LineCount; i++) {

        accum = (cmsUInt8Number*)in + strideIn;

        output = (cmsUInt8Number*)out + strideOut;

        for (j = 0; j < PixelsPerLine; j++) {

            accum = p->FromInput(p, wIn, accum, Stride->BytesPerPlaneIn);

            if (memcmp(wIn, Cache.CacheIn, sizeof(Cache.CacheIn)) == 0) {

                memcpy(wOut, Cache.CacheOut, sizeof(Cache.CacheOut));

            }

            else {

                TransformOnePixelWithGamutCheck(p, wIn, wOut);

                memcpy(Cache.CacheIn, wIn, sizeof(Cache.CacheIn));

                memcpy(Cache.CacheOut, wOut, sizeof(Cache.CacheOut));

            }

            output = p->ToOutput(p, wOut, output, Stride->BytesPerPlaneOut);

        }

        strideIn += Stride->BytesPerLineIn;

        strideOut += Stride->BytesPerLineOut;

    }

}

// Transform plug-ins ----------------------------------------------------------------------------------------------------

// List of used-defined transform factories

typedef struct _cmsTransformCollection_st {

    _cmsTransform2Factory  Factory;

    cmsBool                OldXform;   // Factory returns xform function in the old style

    struct _cmsTransformCollection_st *Next;

} _cmsTransformCollection;

// The linked list head

_cmsTransformPluginChunkType _cmsTransformPluginChunk = { NULL };

// Duplicates the zone of memory used by the plug-in in the new context

static

void DupPluginTransformList(struct _cmsContext_struct* ctx, 

                                               const struct _cmsContext_struct* src)

{

   _cmsTransformPluginChunkType newHead = { NULL };

   _cmsTransformCollection*  entry;

   _cmsTransformCollection*  Anterior = NULL;

   _cmsTransformPluginChunkType* head = (_cmsTransformPluginChunkType*) src->chunks[TransformPlugin];

    // Walk the list copying all nodes

   for (entry = head->TransformCollection;

        entry != NULL;

        entry = entry ->Next) {

            _cmsTransformCollection *newEntry = ( _cmsTransformCollection *) _cmsSubAllocDup(ctx ->MemPool, entry, sizeof(_cmsTransformCollection));

            if (newEntry == NULL) 

                return;

            // We want to keep the linked list order, so this is a little bit tricky

            newEntry -> Next = NULL;

            if (Anterior)

                Anterior -> Next = newEntry;

            Anterior = newEntry;

            if (newHead.TransformCollection == NULL)

                newHead.TransformCollection = newEntry;

    }

  ctx ->chunks[TransformPlugin] = _cmsSubAllocDup(ctx->MemPool, &newHead, sizeof(_cmsTransformPluginChunkType));

}

// Allocates memory for transform plugin factory

void _cmsAllocTransformPluginChunk(struct _cmsContext_struct* ctx, 

                                        const struct _cmsContext_struct* src)

{

    if (src != NULL) {

        // Copy all linked list

        DupPluginTransformList(ctx, src);

    }

    else {

        static _cmsTransformPluginChunkType TransformPluginChunkType = { NULL };

        ctx ->chunks[TransformPlugin] = _cmsSubAllocDup(ctx ->MemPool, &TransformPluginChunkType, sizeof(_cmsTransformPluginChunkType));

    }

}

// Adaptor for old versions of plug-in

static

void _cmsTransform2toTransformAdaptor(struct _cmstransform_struct *CMMcargo,

                                      const void* InputBuffer,

                                      void* OutputBuffer,

                                      cmsUInt32Number PixelsPerLine,

                                      cmsUInt32Number LineCount,

                                      const cmsStride* Stride)

{

       size_t i, strideIn, strideOut;

       _cmsHandleExtraChannels(CMMcargo, InputBuffer, OutputBuffer, PixelsPerLine, LineCount, Stride);

       strideIn = 0;

       strideOut = 0;

       for (i = 0; i < LineCount; i++) {

              void *accum = (cmsUInt8Number*)InputBuffer + strideIn;

              void *output = (cmsUInt8Number*)OutputBuffer + strideOut;

              CMMcargo->OldXform(CMMcargo, accum, output, PixelsPerLine, Stride->BytesPerPlaneIn);

              strideIn += Stride->BytesPerLineIn;

              strideOut += Stride->BytesPerLineOut;

       }

}

// Register new ways to transform

cmsBool  _cmsRegisterTransformPlugin(cmsContext ContextID, cmsPluginBase* Data)

{

    cmsPluginTransform* Plugin = (cmsPluginTransform*) Data;

    _cmsTransformCollection* fl;

    _cmsTransformPluginChunkType* ctx = ( _cmsTransformPluginChunkType*) _cmsContextGetClientChunk(ContextID,TransformPlugin);

    if (Data == NULL) {

        // Free the chain. Memory is safely freed at exit

        ctx->TransformCollection = NULL;

        return TRUE;

    }

    // Factory callback is required

    if (Plugin->factories.xform == NULL) return FALSE;

    fl = (_cmsTransformCollection*) _cmsPluginMalloc(ContextID, sizeof(_cmsTransformCollection));

    if (fl == NULL) return FALSE;

    // Check for full xform plug-ins previous to 2.8, we would need an adapter in that case

    if (Plugin->base.ExpectedVersion < 2080) {

           fl->OldXform = TRUE;

    }

    else

           fl->OldXform = FALSE;

    // Copy the parameters

    fl->Factory = Plugin->factories.xform;

    // Keep linked list

    fl ->Next = ctx->TransformCollection;

    ctx->TransformCollection = fl;

    // All is ok

    return TRUE;

}

void CMSEXPORT _cmsSetTransformUserData(struct _cmstransform_struct *CMMcargo, void* ptr, _cmsFreeUserDataFn FreePrivateDataFn)

{

    _cmsAssert(CMMcargo != NULL);

    CMMcargo ->UserData = ptr;

    CMMcargo ->FreeUserData = FreePrivateDataFn;

}

// returns the pointer defined by the plug-in to store private data

void * CMSEXPORT _cmsGetTransformUserData(struct _cmstransform_struct *CMMcargo)

{

    _cmsAssert(CMMcargo != NULL);

    return CMMcargo ->UserData;

}

// returns the current formatters

void CMSEXPORT _cmsGetTransformFormatters16(struct _cmstransform_struct *CMMcargo, cmsFormatter16* FromInput, cmsFormatter16* ToOutput)

{

     _cmsAssert(CMMcargo != NULL);

     if (FromInput) *FromInput = CMMcargo ->FromInput;

     if (ToOutput)  *ToOutput  = CMMcargo ->ToOutput;

}

void CMSEXPORT _cmsGetTransformFormattersFloat(struct _cmstransform_struct *CMMcargo, cmsFormatterFloat* FromInput, cmsFormatterFloat* ToOutput)

{

     _cmsAssert(CMMcargo != NULL);

     if (FromInput) *FromInput = CMMcargo ->FromInputFloat;

     if (ToOutput)  *ToOutput  = CMMcargo ->ToOutputFloat;

}

// returns original flags

cmsUInt32Number CMSEXPORT _cmsGetTransformFlags(struct _cmstransform_struct* CMMcargo)

{

    _cmsAssert(CMMcargo != NULL);

    return CMMcargo->dwOriginalFlags;

}

// Returns the worker callback for parallelization plug-ins

_cmsTransform2Fn CMSEXPORT _cmsGetTransformWorker(struct _cmstransform_struct* CMMcargo)

{

    _cmsAssert(CMMcargo != NULL);

    return CMMcargo->Worker;

}

// This field holds maximum number of workers or -1 to auto 

cmsInt32Number CMSEXPORT _cmsGetTransformMaxWorkers(struct _cmstransform_struct* CMMcargo)

{

    _cmsAssert(CMMcargo != NULL);

    return CMMcargo->MaxWorkers;

}

// This field is actually unused and reserved

cmsUInt32Number CMSEXPORT _cmsGetTransformWorkerFlags(struct _cmstransform_struct* CMMcargo)

{

    _cmsAssert(CMMcargo != NULL);

    return CMMcargo->WorkerFlags;

}

// In the case there is a parallelization plug-in, let it to do its job

static

void ParalellizeIfSuitable(_cmsTRANSFORM* p)

{

    _cmsParallelizationPluginChunkType* ctx = (_cmsParallelizationPluginChunkType*)_cmsContextGetClientChunk(p->ContextID, ParallelizationPlugin);

    _cmsAssert(p != NULL);

    if (ctx != NULL && ctx->SchedulerFn != NULL) {

        p->Worker = p->xform;

        p->xform = ctx->SchedulerFn;

        p->MaxWorkers = ctx->MaxWorkers;

        p->WorkerFlags = ctx->WorkerFlags;

    }

}

/**

* An empty unroll to avoid a check with NULL on cmsDoTransform()

*/

static

cmsUInt8Number* UnrollNothing(CMSREGISTER _cmsTRANSFORM* info,

                              CMSREGISTER cmsUInt16Number wIn[],

                              CMSREGISTER cmsUInt8Number* accum,

                              CMSREGISTER cmsUInt32Number Stride)

{    

    return accum;

    cmsUNUSED_PARAMETER(info);

    cmsUNUSED_PARAMETER(wIn);

    cmsUNUSED_PARAMETER(Stride);

}

static

cmsUInt8Number* PackNothing(CMSREGISTER _cmsTRANSFORM* info,

                           CMSREGISTER cmsUInt16Number wOut[],

                           CMSREGISTER cmsUInt8Number* output,

                           CMSREGISTER cmsUInt32Number Stride)

{

    return output;

    cmsUNUSED_PARAMETER(info);

    cmsUNUSED_PARAMETER(wOut);

    cmsUNUSED_PARAMETER(Stride);

}

// Allocate transform struct and set it to defaults. Ask the optimization plug-in about if those formats are proper

// for separated transforms. If this is the case,

static

_cmsTRANSFORM* AllocEmptyTransform(cmsContext ContextID, cmsPipeline* lut,

                                               cmsUInt32Number Intent, cmsUInt32Number* InputFormat, cmsUInt32Number* OutputFormat, cmsUInt32Number* dwFlags)

{

     _cmsTransformPluginChunkType* ctx = ( _cmsTransformPluginChunkType*) _cmsContextGetClientChunk(ContextID, TransformPlugin);

     _cmsTransformCollection* Plugin;

       // Allocate needed memory

       _cmsTRANSFORM* p = (_cmsTRANSFORM*)_cmsMallocZero(ContextID, sizeof(_cmsTRANSFORM));

       if (!p) {

              cmsPipelineFree(lut);

              return NULL;

       }

       // Store the proposed pipeline

       p->Lut = lut;

       // Let's see if any plug-in want to do the transform by itself

       if (p->Lut != NULL) {

           if (!(*dwFlags & cmsFLAGS_NOOPTIMIZE))

           {

               for (Plugin = ctx->TransformCollection;

                   Plugin != NULL;

                   Plugin = Plugin->Next) {

                   if (Plugin->Factory(&p->xform, &p->UserData, &p->FreeUserData, &p->Lut, InputFormat, OutputFormat, dwFlags)) {

                       // Last plugin in the declaration order takes control. We just keep

                       // the original parameters as a logging. 

                       // Note that cmsFLAGS_CAN_CHANGE_FORMATTER is not set, so by default 

                       // an optimized transform is not reusable. The plug-in can, however, change

                       // the flags and make it suitable.

                       p->ContextID = ContextID;

                       p->InputFormat = *InputFormat;

                       p->OutputFormat = *OutputFormat;

                       p->dwOriginalFlags = *dwFlags;

                       // Fill the formatters just in case the optimized routine is interested.

                       // No error is thrown if the formatter doesn't exist. It is up to the optimization 

                       // factory to decide what to do in those cases.

                       p->FromInput = _cmsGetFormatter(ContextID, *InputFormat, cmsFormatterInput, CMS_PACK_FLAGS_16BITS).Fmt16;

                       p->ToOutput = _cmsGetFormatter(ContextID, *OutputFormat, cmsFormatterOutput, CMS_PACK_FLAGS_16BITS).Fmt16;

                       p->FromInputFloat = _cmsGetFormatter(ContextID, *InputFormat, cmsFormatterInput, CMS_PACK_FLAGS_FLOAT).FmtFloat;

                       p->ToOutputFloat = _cmsGetFormatter(ContextID, *OutputFormat, cmsFormatterOutput, CMS_PACK_FLAGS_FLOAT).FmtFloat;

                       // Save the day? (Ignore the warning)

                       if (Plugin->OldXform) {

                           p->OldXform = (_cmsTransformFn)(void*)p->xform;

                           p->xform = _cmsTransform2toTransformAdaptor;

                       }

                       ParalellizeIfSuitable(p);

                       return p;

                   }

               }

           }

           // Not suitable for the transform plug-in, let's check  the pipeline plug-in

           _cmsOptimizePipeline(ContextID, &p->Lut, Intent, InputFormat, OutputFormat, dwFlags);

       }

    // Check whatever this is a true floating point transform

    if (_cmsFormatterIsFloat(*InputFormat) || _cmsFormatterIsFloat(*OutputFormat)) {

        // Get formatter function always return a valid union, but the contents of this union may be NULL.

        p ->FromInputFloat = _cmsGetFormatter(ContextID, *InputFormat,  cmsFormatterInput, CMS_PACK_FLAGS_FLOAT).FmtFloat;

        p ->ToOutputFloat  = _cmsGetFormatter(ContextID, *OutputFormat, cmsFormatterOutput, CMS_PACK_FLAGS_FLOAT).FmtFloat;

        *dwFlags |= cmsFLAGS_CAN_CHANGE_FORMATTER;

        if (p ->FromInputFloat == NULL || p ->ToOutputFloat == NULL) {

            cmsSignalError(ContextID, cmsERROR_UNKNOWN_EXTENSION, "Unsupported raster format");

            cmsDeleteTransform(p);

            return NULL;

        }

        if (*dwFlags & cmsFLAGS_NULLTRANSFORM) {

            p ->xform = NullFloatXFORM;

        }

        else {

            // Float transforms don't use cache, always are non-NULL

            p ->xform = FloatXFORM;