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

//

//  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"

// Virtual (built-in) profiles

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

static

cmsBool SetTextTags(cmsHPROFILE hProfile, const wchar_t* Description)

{

    cmsMLU *DescriptionMLU, *CopyrightMLU;

    cmsBool  rc = FALSE;

    cmsContext ContextID = cmsGetProfileContextID(hProfile);

    DescriptionMLU  = cmsMLUalloc(ContextID, 1);

    CopyrightMLU    = cmsMLUalloc(ContextID, 1);

    if (DescriptionMLU == NULL || CopyrightMLU == NULL) goto Error;

    if (!cmsMLUsetWide(DescriptionMLU,  "en", "US", Description)) goto Error;

    if (!cmsMLUsetWide(CopyrightMLU,    "en", "US", L"No copyright, use freely")) goto Error;

    if (!cmsWriteTag(hProfile, cmsSigProfileDescriptionTag,  DescriptionMLU)) goto Error;

    if (!cmsWriteTag(hProfile, cmsSigCopyrightTag,           CopyrightMLU)) goto Error;

    rc = TRUE;

Error:

    if (DescriptionMLU)

        cmsMLUfree(DescriptionMLU);

    if (CopyrightMLU)

        cmsMLUfree(CopyrightMLU);

    return rc;

}

static

cmsBool  SetSeqDescTag(cmsHPROFILE hProfile, const char* Model)

{

    cmsBool  rc = FALSE;

    cmsContext ContextID = cmsGetProfileContextID(hProfile);

    cmsSEQ* Seq = cmsAllocProfileSequenceDescription(ContextID, 1);

    if (Seq == NULL) return FALSE;

    Seq->seq[0].deviceMfg = (cmsSignature) 0;

    Seq->seq[0].deviceModel = (cmsSignature) 0;

#ifdef CMS_DONT_USE_INT64

    Seq->seq[0].attributes[0] = 0;

    Seq->seq[0].attributes[1] = 0;

#else

    Seq->seq[0].attributes = 0;

#endif

    Seq->seq[0].technology = (cmsTechnologySignature) 0;

    cmsMLUsetASCII( Seq->seq[0].Manufacturer, cmsNoLanguage, cmsNoCountry, "Little CMS");

    cmsMLUsetASCII( Seq->seq[0].Model,        cmsNoLanguage, cmsNoCountry, Model);

    if (!_cmsWriteProfileSequence(hProfile, Seq)) goto Error;

    rc = TRUE;

Error:

    if (Seq)

        cmsFreeProfileSequenceDescription(Seq);

    return rc;

}

// This function creates a profile based on White point, primaries and

// transfer functions.

cmsHPROFILE CMSEXPORT cmsCreateRGBProfileTHR(cmsContext ContextID,

                                          const cmsCIExyY* WhitePoint,

                                          const cmsCIExyYTRIPLE* Primaries,

                                          cmsToneCurve* const TransferFunction[3])

{

    cmsHPROFILE hICC;

    cmsMAT3 MColorants;

    cmsCIEXYZTRIPLE Colorants;

    cmsCIExyY MaxWhite;

    cmsMAT3 CHAD;

    cmsCIEXYZ WhitePointXYZ;

    hICC = cmsCreateProfilePlaceholder(ContextID);

    if (!hICC)                          // can't allocate

        return NULL;

    cmsSetProfileVersion(hICC, 4.4);

    cmsSetDeviceClass(hICC,      cmsSigDisplayClass);

    cmsSetColorSpace(hICC,       cmsSigRgbData);

    cmsSetPCS(hICC,              cmsSigXYZData);

    cmsSetHeaderRenderingIntent(hICC,  INTENT_PERCEPTUAL);

    // Implement profile using following tags:

    //

    //  1 cmsSigProfileDescriptionTag

    //  2 cmsSigMediaWhitePointTag

    //  3 cmsSigRedColorantTag

    //  4 cmsSigGreenColorantTag

    //  5 cmsSigBlueColorantTag

    //  6 cmsSigRedTRCTag

    //  7 cmsSigGreenTRCTag

    //  8 cmsSigBlueTRCTag

    //  9 Chromatic adaptation Tag

    // This conforms a standard RGB DisplayProfile as says ICC, and then I add (As per addendum II)

    // 10 cmsSigChromaticityTag

    if (!SetTextTags(hICC, L"RGB built-in")) goto Error;

    if (WhitePoint) {

        if (!cmsWriteTag(hICC, cmsSigMediaWhitePointTag, cmsD50_XYZ())) goto Error;

        cmsxyY2XYZ(&WhitePointXYZ, WhitePoint);

        _cmsAdaptationMatrix(&CHAD, NULL, &WhitePointXYZ, cmsD50_XYZ());

        // This is a V4 tag, but many CMM does read and understand it no matter which version

        if (!cmsWriteTag(hICC, cmsSigChromaticAdaptationTag, (void*) &CHAD)) goto Error;

    }

    if (WhitePoint && Primaries) {

        MaxWhite.x =  WhitePoint -> x;

        MaxWhite.y =  WhitePoint -> y;

        MaxWhite.Y =  1.0;

        if (!_cmsBuildRGB2XYZtransferMatrix(&MColorants, &MaxWhite, Primaries)) goto Error;

        Colorants.Red.X   = MColorants.v[0].n[0];

        Colorants.Red.Y   = MColorants.v[1].n[0];

        Colorants.Red.Z   = MColorants.v[2].n[0];

        Colorants.Green.X = MColorants.v[0].n[1];

        Colorants.Green.Y = MColorants.v[1].n[1];

        Colorants.Green.Z = MColorants.v[2].n[1];

        Colorants.Blue.X  = MColorants.v[0].n[2];

        Colorants.Blue.Y  = MColorants.v[1].n[2];

        Colorants.Blue.Z  = MColorants.v[2].n[2];

        if (!cmsWriteTag(hICC, cmsSigRedColorantTag,   (void*) &Colorants.Red)) goto Error;

        if (!cmsWriteTag(hICC, cmsSigBlueColorantTag,  (void*) &Colorants.Blue)) goto Error;

        if (!cmsWriteTag(hICC, cmsSigGreenColorantTag, (void*) &Colorants.Green)) goto Error;

    }

    if (TransferFunction) {

        // Tries to minimize space. Thanks to Richard Hughes for this nice idea         

        if (!cmsWriteTag(hICC, cmsSigRedTRCTag,   (void*) TransferFunction[0])) goto Error;

        if (TransferFunction[1] == TransferFunction[0]) {

            if (!cmsLinkTag (hICC, cmsSigGreenTRCTag, cmsSigRedTRCTag)) goto Error;

        } else {

            if (!cmsWriteTag(hICC, cmsSigGreenTRCTag, (void*) TransferFunction[1])) goto Error;

        }

        if (TransferFunction[2] == TransferFunction[0]) {

            if (!cmsLinkTag (hICC, cmsSigBlueTRCTag, cmsSigRedTRCTag)) goto Error;

        } else {

            if (!cmsWriteTag(hICC, cmsSigBlueTRCTag, (void*) TransferFunction[2])) goto Error;

        }

    }

    if (Primaries) {

        if (!cmsWriteTag(hICC, cmsSigChromaticityTag, (void*) Primaries)) goto Error;

    }

    return hICC;

Error:

    if (hICC)

        cmsCloseProfile(hICC);

    return NULL;

}

cmsHPROFILE CMSEXPORT cmsCreateRGBProfile(const cmsCIExyY* WhitePoint,

                                          const cmsCIExyYTRIPLE* Primaries,

                                          cmsToneCurve* const TransferFunction[3])

{

    return cmsCreateRGBProfileTHR(NULL, WhitePoint, Primaries, TransferFunction);

}

// This function creates a profile based on White point and transfer function.

cmsHPROFILE CMSEXPORT cmsCreateGrayProfileTHR(cmsContext ContextID,

                                           const cmsCIExyY* WhitePoint,

                                           const cmsToneCurve* TransferFunction)

{

    cmsHPROFILE hICC;

    cmsCIEXYZ tmp;

    hICC = cmsCreateProfilePlaceholder(ContextID);

    if (!hICC)                          // can't allocate

        return NULL;

    cmsSetProfileVersion(hICC, 4.4);

    cmsSetDeviceClass(hICC,      cmsSigDisplayClass);

    cmsSetColorSpace(hICC,       cmsSigGrayData);

    cmsSetPCS(hICC,              cmsSigXYZData);

    cmsSetHeaderRenderingIntent(hICC,  INTENT_PERCEPTUAL);

    // Implement profile using following tags:

    //

    //  1 cmsSigProfileDescriptionTag

    //  2 cmsSigMediaWhitePointTag

    //  3 cmsSigGrayTRCTag

    // This conforms a standard Gray DisplayProfile

    // Fill-in the tags

    if (!SetTextTags(hICC, L"gray built-in")) goto Error;

    if (WhitePoint) {

        cmsxyY2XYZ(&tmp, WhitePoint);

        if (!cmsWriteTag(hICC, cmsSigMediaWhitePointTag, (void*) &tmp)) goto Error;

    }

    if (TransferFunction) {

        if (!cmsWriteTag(hICC, cmsSigGrayTRCTag, (void*) TransferFunction)) goto Error;

    }

    return hICC;

Error:

    if (hICC)

        cmsCloseProfile(hICC);

    return NULL;

}

cmsHPROFILE CMSEXPORT cmsCreateGrayProfile(const cmsCIExyY* WhitePoint,

                                                    const cmsToneCurve* TransferFunction)

{

    return cmsCreateGrayProfileTHR(NULL, WhitePoint, TransferFunction);

}

// This is a devicelink operating in the target colorspace with as many transfer functions as components

cmsHPROFILE CMSEXPORT cmsCreateLinearizationDeviceLinkTHR(cmsContext ContextID,

                                                          cmsColorSpaceSignature ColorSpace,

                                                          cmsToneCurve* const TransferFunctions[])

{

    cmsHPROFILE hICC;

    cmsPipeline* Pipeline;

    cmsInt32Number nChannels;

    hICC = cmsCreateProfilePlaceholder(ContextID);

    if (!hICC)

        return NULL;

    cmsSetProfileVersion(hICC, 4.4);

    cmsSetDeviceClass(hICC,      cmsSigLinkClass);

    cmsSetColorSpace(hICC,       ColorSpace);

    cmsSetPCS(hICC,              ColorSpace);

    cmsSetHeaderRenderingIntent(hICC,  INTENT_PERCEPTUAL);

    // Set up channels

    nChannels = cmsChannelsOfColorSpace(ColorSpace);

    // Creates a Pipeline with prelinearization step only

    Pipeline = cmsPipelineAlloc(ContextID, nChannels, nChannels);

    if (Pipeline == NULL) goto Error;

    // Copy tables to Pipeline

    if (!cmsPipelineInsertStage(Pipeline, cmsAT_BEGIN, cmsStageAllocToneCurves(ContextID, nChannels, TransferFunctions)))

        goto Error;

    // Create tags

    if (!SetTextTags(hICC, L"Linearization built-in")) goto Error;

    if (!cmsWriteTag(hICC, cmsSigAToB0Tag, (void*) Pipeline)) goto Error;

    if (!SetSeqDescTag(hICC, "Linearization built-in")) goto Error;

    // Pipeline is already on virtual profile

    cmsPipelineFree(Pipeline);

    // Ok, done

    return hICC;

Error:

    cmsPipelineFree(Pipeline);

    if (hICC)

        cmsCloseProfile(hICC);

    return NULL;

}

cmsHPROFILE CMSEXPORT cmsCreateLinearizationDeviceLink(cmsColorSpaceSignature ColorSpace,

                                                                 cmsToneCurve* const TransferFunctions[])

{

    return cmsCreateLinearizationDeviceLinkTHR(NULL, ColorSpace, TransferFunctions);

}

// Ink-limiting algorithm

//

//  Sum = C + M + Y + K

//  If Sum > InkLimit

//        Ratio= 1 - (Sum - InkLimit) / (C + M + Y)

//        if Ratio <0

//              Ratio=0

//        endif

//     Else

//         Ratio=1

//     endif

//

//     C = Ratio * C

//     M = Ratio * M

//     Y = Ratio * Y

//     K: Does not change

static

int InkLimitingSampler(CMSREGISTER const cmsUInt16Number In[], CMSREGISTER cmsUInt16Number Out[], CMSREGISTER void* Cargo)

{

    cmsFloat64Number InkLimit = *(cmsFloat64Number *) Cargo;

    cmsFloat64Number SumCMY, SumCMYK, Ratio;

    InkLimit = (InkLimit * 655.35);

    SumCMY   = (cmsFloat64Number) In[0]  + In[1] + In[2];

    SumCMYK  = SumCMY + In[3];

    if (SumCMYK > InkLimit) {

        Ratio = 1 - ((SumCMYK - InkLimit) / SumCMY);

        if (Ratio < 0)

            Ratio = 0;

    }

    else Ratio = 1;

    Out[0] = _cmsQuickSaturateWord(In[0] * Ratio);     // C

    Out[1] = _cmsQuickSaturateWord(In[1] * Ratio);     // M

    Out[2] = _cmsQuickSaturateWord(In[2] * Ratio);     // Y

    Out[3] = In[3];                                    // K (untouched)

    return TRUE;

}

// This is a devicelink operating in CMYK for ink-limiting

cmsHPROFILE CMSEXPORT cmsCreateInkLimitingDeviceLinkTHR(cmsContext ContextID,

                                                     cmsColorSpaceSignature ColorSpace,

                                                     cmsFloat64Number Limit)

{

    cmsHPROFILE hICC;

    cmsPipeline* LUT;

    cmsStage* CLUT;

    cmsInt32Number nChannels;

    if (ColorSpace != cmsSigCmykData) {

        cmsSignalError(ContextID, cmsERROR_COLORSPACE_CHECK, "InkLimiting: Only CMYK currently supported");

        return NULL;

    }

    if (Limit < 1.0 || Limit > 400) {

        cmsSignalError(ContextID, cmsERROR_RANGE, "InkLimiting: Limit should be between 1..400");

        if (Limit < 1) Limit = 1;

        if (Limit > 400) Limit = 400;

    }

    hICC = cmsCreateProfilePlaceholder(ContextID);

    if (!hICC)                          // can't allocate

        return NULL;

    cmsSetProfileVersion(hICC, 4.4);

    cmsSetDeviceClass(hICC,      cmsSigLinkClass);

    cmsSetColorSpace(hICC,       ColorSpace);

    cmsSetPCS(hICC,              ColorSpace);

    cmsSetHeaderRenderingIntent(hICC,  INTENT_PERCEPTUAL);

    // Creates a Pipeline with 3D grid only

    LUT = cmsPipelineAlloc(ContextID, 4, 4);

    if (LUT == NULL) goto Error;

    nChannels = cmsChannelsOf(ColorSpace);

    CLUT = cmsStageAllocCLut16bit(ContextID, 17, nChannels, nChannels, NULL);

    if (CLUT == NULL) goto Error;

    if (!cmsStageSampleCLut16bit(CLUT, InkLimitingSampler, (void*) &Limit, 0)) goto Error;

    if (!cmsPipelineInsertStage(LUT, cmsAT_BEGIN, _cmsStageAllocIdentityCurves(ContextID, nChannels)) ||

        !cmsPipelineInsertStage(LUT, cmsAT_END, CLUT) ||

        !cmsPipelineInsertStage(LUT, cmsAT_END, _cmsStageAllocIdentityCurves(ContextID, nChannels)))

        goto Error;

    // Create tags

    if (!SetTextTags(hICC, L"ink-limiting built-in")) goto Error;

    if (!cmsWriteTag(hICC, cmsSigAToB0Tag, (void*) LUT))  goto Error;

    if (!SetSeqDescTag(hICC, "ink-limiting built-in")) goto Error;

    // cmsPipeline is already on virtual profile

    cmsPipelineFree(LUT);

    // Ok, done

    return hICC;

Error:

    if (LUT != NULL)

        cmsPipelineFree(LUT);

    if (hICC != NULL)

        cmsCloseProfile(hICC);

    return NULL;

}

cmsHPROFILE CMSEXPORT cmsCreateInkLimitingDeviceLink(cmsColorSpaceSignature ColorSpace, cmsFloat64Number Limit)

{

    return cmsCreateInkLimitingDeviceLinkTHR(NULL, ColorSpace, Limit);

}

// Creates a fake Lab identity.

cmsHPROFILE CMSEXPORT cmsCreateLab2ProfileTHR(cmsContext ContextID, const cmsCIExyY* WhitePoint)

{

    cmsHPROFILE hProfile;

    cmsPipeline* LUT = NULL;

    hProfile = cmsCreateRGBProfileTHR(ContextID, WhitePoint == NULL ? cmsD50_xyY() : WhitePoint, NULL, NULL);

    if (hProfile == NULL) return NULL;

    cmsSetProfileVersion(hProfile, 2.1);

    cmsSetDeviceClass(hProfile, cmsSigAbstractClass);

    cmsSetColorSpace(hProfile,  cmsSigLabData);

    cmsSetPCS(hProfile,         cmsSigLabData);

    if (!SetTextTags(hProfile, L"Lab identity built-in")) return NULL;

    // An identity LUT is all we need

    LUT = cmsPipelineAlloc(ContextID, 3, 3);

    if (LUT == NULL) goto Error;

    if (!cmsPipelineInsertStage(LUT, cmsAT_BEGIN, _cmsStageAllocIdentityCLut(ContextID, 3)))

        goto Error;

    if (!cmsWriteTag(hProfile, cmsSigAToB0Tag, LUT)) goto Error;

    cmsPipelineFree(LUT);

    return hProfile;

Error:

    if (LUT != NULL)

        cmsPipelineFree(LUT);

    if (hProfile != NULL)

        cmsCloseProfile(hProfile);

    return NULL;

}

cmsHPROFILE CMSEXPORT cmsCreateLab2Profile(const cmsCIExyY* WhitePoint)

{

    return cmsCreateLab2ProfileTHR(NULL, WhitePoint);

}

// Creates a fake Lab V4 identity.

cmsHPROFILE CMSEXPORT cmsCreateLab4ProfileTHR(cmsContext ContextID, const cmsCIExyY* WhitePoint)

{

    cmsHPROFILE hProfile;

    cmsPipeline* LUT = NULL;

    cmsCIEXYZ xyz;

    if (WhitePoint == NULL)

        xyz = *cmsD50_XYZ();

    else

        cmsxyY2XYZ(&xyz, WhitePoint);

    hProfile = cmsCreateRGBProfileTHR(ContextID, NULL, NULL, NULL);

    if (hProfile == NULL) return NULL;

    cmsSetProfileVersion(hProfile, 4.4);

    cmsSetDeviceClass(hProfile, cmsSigAbstractClass);

    cmsSetColorSpace(hProfile,  cmsSigLabData);

    cmsSetPCS(hProfile,         cmsSigLabData);

    if (!cmsWriteTag(hProfile, cmsSigMediaWhitePointTag, &xyz)) goto Error;

    if (!SetTextTags(hProfile, L"Lab identity built-in")) goto Error;

    // An empty LUTs is all we need

    LUT = cmsPipelineAlloc(ContextID, 3, 3);

    if (LUT == NULL) goto Error;

    if (!cmsPipelineInsertStage(LUT, cmsAT_BEGIN, _cmsStageAllocIdentityCurves(ContextID, 3)))

        goto Error;

    if (!cmsWriteTag(hProfile, cmsSigAToB0Tag, LUT)) goto Error;

    cmsPipelineFree(LUT);

    return hProfile;

Error:

    if (LUT != NULL)

        cmsPipelineFree(LUT);

    if (hProfile != NULL)

        cmsCloseProfile(hProfile);

    return NULL;

}

cmsHPROFILE CMSEXPORT cmsCreateLab4Profile(const cmsCIExyY* WhitePoint)

{

    return cmsCreateLab4ProfileTHR(NULL, WhitePoint);

}

// Creates a fake XYZ identity

cmsHPROFILE CMSEXPORT cmsCreateXYZProfileTHR(cmsContext ContextID)

{

    cmsHPROFILE hProfile;

    cmsPipeline* LUT = NULL;

    hProfile = cmsCreateRGBProfileTHR(ContextID, cmsD50_xyY(), NULL, NULL);

    if (hProfile == NULL) return NULL;

    cmsSetProfileVersion(hProfile, 4.4);

    cmsSetDeviceClass(hProfile, cmsSigAbstractClass);

    cmsSetColorSpace(hProfile,  cmsSigXYZData);

    cmsSetPCS(hProfile,         cmsSigXYZData);

    if (!SetTextTags(hProfile, L"XYZ identity built-in")) goto Error;

    // An identity LUT is all we need

    LUT = cmsPipelineAlloc(ContextID, 3, 3);

    if (LUT == NULL) goto Error;

    if (!cmsPipelineInsertStage(LUT, cmsAT_BEGIN, _cmsStageAllocIdentityCurves(ContextID, 3)))

        goto Error;

    if (!cmsWriteTag(hProfile, cmsSigAToB0Tag, LUT)) goto Error;

    cmsPipelineFree(LUT);

    return hProfile;

Error:

    if (LUT != NULL)

        cmsPipelineFree(LUT);

    if (hProfile != NULL)

        cmsCloseProfile(hProfile);

    return NULL;

}

cmsHPROFILE CMSEXPORT cmsCreateXYZProfile(void)

{

    return cmsCreateXYZProfileTHR(NULL);

}

//sRGB Curves are defined by:

//

//If  R'sRGB,G'sRGB, B'sRGB < 0.04045

//

//    R =  R'sRGB / 12.92

//    G =  G'sRGB / 12.92

//    B =  B'sRGB / 12.92

//

//

//else if  R'sRGB,G'sRGB, B'sRGB >= 0.04045

//

//    R = ((R'sRGB + 0.055) / 1.055)^2.4

//    G = ((G'sRGB + 0.055) / 1.055)^2.4

//    B = ((B'sRGB + 0.055) / 1.055)^2.4

static

cmsToneCurve* Build_sRGBGamma(cmsContext ContextID)

{

    cmsFloat64Number Parameters[5];

    Parameters[0] = 2.4;

    Parameters[1] = 1. / 1.055;

    Parameters[2] = 0.055 / 1.055;

    Parameters[3] = 1. / 12.92;

    Parameters[4] = 0.04045;

    return cmsBuildParametricToneCurve(ContextID, 4, Parameters);

}

// Create the ICC virtual profile for sRGB space

cmsHPROFILE CMSEXPORT cmsCreate_sRGBProfileTHR(cmsContext ContextID)

{

       cmsCIExyY       D65 = { 0.3127, 0.3290, 1.0 };

       cmsCIExyYTRIPLE Rec709Primaries = {

                                   {0.6400, 0.3300, 1.0},

                                   {0.3000, 0.6000, 1.0},

                                   {0.1500, 0.0600, 1.0}

                                   };

       cmsToneCurve* Gamma22[3];

       cmsHPROFILE  hsRGB;

      // cmsWhitePointFromTemp(&D65, 6504);

       Gamma22[0] = Gamma22[1] = Gamma22[2] = Build_sRGBGamma(ContextID);

       if (Gamma22[0] == NULL) return NULL;

       hsRGB = cmsCreateRGBProfileTHR(ContextID, &D65, &Rec709Primaries, Gamma22);

       cmsFreeToneCurve(Gamma22[0]);

       if (hsRGB == NULL) return NULL;

       if (!SetTextTags(hsRGB, L"sRGB built-in")) {

           cmsCloseProfile(hsRGB);

           return NULL;

       }

       return hsRGB;

}

cmsHPROFILE CMSEXPORT cmsCreate_sRGBProfile(void)

{

    return cmsCreate_sRGBProfileTHR(NULL);

}

/**

* Oklab colorspace profile (experimental)

* 

* This virtual profile cannot be saved as an ICC file

*/

cmsHPROFILE CMSEXPORT cmsCreate_OkLabProfile(cmsContext ctx)

{

    cmsStage* XYZPCS = _cmsStageNormalizeFromXyzFloat(ctx);

    cmsStage* PCSXYZ = _cmsStageNormalizeToXyzFloat(ctx);

    const double M_D65_D50[] =

    {

       1.047886, 0.022919, -0.050216,

       0.029582, 0.990484, -0.017079,

      -0.009252, 0.015073,  0.751678

    };

    const double M_D50_D65[] =

    {

         0.955512609517083, -0.023073214184645,  0.063308961782107,

        -0.028324949364887,  1.009942432477107,  0.021054814890112,

         0.012328875695483, -0.020535835374141,  1.330713916450354

    };

    cmsStage* D65toD50 = cmsStageAllocMatrix(ctx, 3, 3, M_D65_D50, NULL);

    cmsStage* D50toD65 = cmsStageAllocMatrix(ctx, 3, 3, M_D50_D65, NULL);

    const double M_D65_LMS[] =

    {

        0.8189330101, 0.3618667424, -0.1288597137,

        0.0329845436, 0.9293118715,  0.0361456387,

        0.0482003018, 0.2643662691,  0.6338517070

    };

    const double M_LMS_D65[] =

    {

        1.227013851103521, -0.557799980651822,  0.281256148966468,

       -0.040580178423281,  1.112256869616830, -0.071676678665601,

       -0.076381284505707, -0.421481978418013,  1.586163220440795

    };

    cmsStage* D65toLMS = cmsStageAllocMatrix(ctx, 3, 3, M_D65_LMS, NULL);

    cmsStage* LMStoD65 = cmsStageAllocMatrix(ctx, 3, 3, M_LMS_D65, NULL);

    cmsToneCurve* CubeRoot = cmsBuildGamma(ctx, 1.0 / 3.0);

    cmsToneCurve* Cube     = cmsBuildGamma(ctx,  3.0);

    cmsToneCurve* Roots[3] = { CubeRoot, CubeRoot, CubeRoot };

    cmsToneCurve* Cubes[3] = { Cube, Cube, Cube };

    cmsStage* NonLinearityFw = cmsStageAllocToneCurves(ctx, 3, Roots);

    cmsStage* NonLinearityRv = cmsStageAllocToneCurves(ctx, 3, Cubes);

    const double M_LMSprime_OkLab[] =

    {

        0.2104542553,  0.7936177850, -0.0040720468,

        1.9779984951, -2.4285922050,  0.4505937099,

        0.0259040371,  0.7827717662, -0.8086757660

    };

    const double M_OkLab_LMSprime[] =

    {

        0.999999998450520,  0.396337792173768,  0.215803758060759,

        1.000000008881761, -0.105561342323656, -0.063854174771706,

        1.000000054672411, -0.089484182094966, -1.291485537864092

    };

    cmsStage* LMSprime_OkLab = cmsStageAllocMatrix(ctx, 3, 3, M_LMSprime_OkLab, NULL);

    cmsStage* OkLab_LMSprime = cmsStageAllocMatrix(ctx, 3, 3, M_OkLab_LMSprime, NULL);

    cmsPipeline* AToB = cmsPipelineAlloc(ctx, 3, 3);

    cmsPipeline* BToA = cmsPipelineAlloc(ctx, 3, 3);

    cmsHPROFILE hProfile = cmsCreateProfilePlaceholder(ctx);

    if (!hProfile)            // can't allocate

        goto error;

    cmsSetProfileVersion(hProfile, 4.4);

    cmsSetDeviceClass(hProfile, cmsSigColorSpaceClass);

    cmsSetColorSpace(hProfile, cmsSig3colorData);

    cmsSetPCS(hProfile, cmsSigXYZData);

    cmsSetHeaderRenderingIntent(hProfile, INTENT_RELATIVE_COLORIMETRIC);

    /**

    * Conversion PCS (XYZ/D50) to OkLab 

    */

    if (!cmsPipelineInsertStage(BToA, cmsAT_END, PCSXYZ)) goto error;

    if (!cmsPipelineInsertStage(BToA, cmsAT_END, D50toD65)) goto error;

    if (!cmsPipelineInsertStage(BToA, cmsAT_END, D65toLMS)) goto error;

    if (!cmsPipelineInsertStage(BToA, cmsAT_END, NonLinearityFw)) goto error;

    if (!cmsPipelineInsertStage(BToA, cmsAT_END, LMSprime_OkLab)) goto error;

    if (!cmsWriteTag(hProfile, cmsSigBToA0Tag, BToA)) goto error;

    if (!cmsPipelineInsertStage(AToB, cmsAT_END, OkLab_LMSprime)) goto error;

    if (!cmsPipelineInsertStage(AToB, cmsAT_END, NonLinearityRv)) goto error;

    if (!cmsPipelineInsertStage(AToB, cmsAT_END, LMStoD65)) goto error;

    if (!cmsPipelineInsertStage(AToB, cmsAT_END, D65toD50)) goto error;

    if (!cmsPipelineInsertStage(AToB, cmsAT_END, XYZPCS)) goto error;

    if (!cmsWriteTag(hProfile, cmsSigAToB0Tag, AToB)) goto error;

    cmsPipelineFree(BToA);

    cmsPipelineFree(AToB);

    cmsFreeToneCurve(CubeRoot);

    cmsFreeToneCurve(Cube);

    return hProfile;

error:

    cmsPipelineFree(BToA);

    cmsPipelineFree(AToB);

    cmsFreeToneCurve(CubeRoot);

    cmsFreeToneCurve(Cube);

    cmsCloseProfile(hProfile);

    return NULL;

}

typedef struct {

                cmsFloat64Number Brightness;

                cmsFloat64Number Contrast;

                cmsFloat64Number Hue;

                cmsFloat64Number Saturation;

                cmsBool          lAdjustWP;

                cmsCIEXYZ WPsrc, WPdest;

} BCHSWADJUSTS, *LPBCHSWADJUSTS;

static

int bchswSampler(CMSREGISTER const cmsUInt16Number In[], CMSREGISTER cmsUInt16Number Out[], CMSREGISTER void* Cargo)

{

    cmsCIELab LabIn, LabOut;

    cmsCIELCh LChIn, LChOut;

    cmsCIEXYZ XYZ;

    LPBCHSWADJUSTS bchsw = (LPBCHSWADJUSTS) Cargo;

    cmsLabEncoded2Float(&LabIn, In);

    cmsLab2LCh(&LChIn, &LabIn);

    // Do some adjusts on LCh

    LChOut.L = LChIn.L * bchsw ->Contrast + bchsw ->Brightness;

    LChOut.C = LChIn.C + bchsw -> Saturation;

    LChOut.h = LChIn.h + bchsw -> Hue;

    cmsLCh2Lab(&LabOut, &LChOut);

    // Move white point in Lab

    if (bchsw->lAdjustWP) {

           cmsLab2XYZ(&bchsw->WPsrc, &XYZ, &LabOut);

           cmsXYZ2Lab(&bchsw->WPdest, &LabOut, &XYZ);

    }

    // Back to encoded

    cmsFloat2LabEncoded(Out, &LabOut);

    return TRUE;

}

// Creates an abstract profile operating in Lab space for Brightness,

// contrast, Saturation and white point displacement

cmsHPROFILE CMSEXPORT cmsCreateBCHSWabstractProfileTHR(cmsContext ContextID,

                                                       cmsUInt32Number nLUTPoints,

                                                       cmsFloat64Number Bright,

                                                       cmsFloat64Number Contrast,

                                                       cmsFloat64Number Hue,

                                                       cmsFloat64Number Saturation,

                                                       cmsUInt32Number TempSrc,

                                                       cmsUInt32Number TempDest)

{

    cmsHPROFILE hICC;

    cmsPipeline* Pipeline;

    BCHSWADJUSTS bchsw;

    cmsCIExyY WhitePnt;

    cmsStage* CLUT;

    cmsUInt32Number Dimensions[MAX_INPUT_DIMENSIONS];

    cmsUInt32Number i;

    bchsw.Brightness = Bright;

    bchsw.Contrast   = Contrast;

    bchsw.Hue        = Hue;

    bchsw.Saturation = Saturation;

    if (TempSrc == TempDest) {

           bchsw.lAdjustWP = FALSE;

    }

    else {

           bchsw.lAdjustWP = TRUE;

           cmsWhitePointFromTemp(&WhitePnt, TempSrc);

           cmsxyY2XYZ(&bchsw.WPsrc, &WhitePnt);

           cmsWhitePointFromTemp(&WhitePnt, TempDest);

           cmsxyY2XYZ(&bchsw.WPdest, &WhitePnt);

    }

    hICC = cmsCreateProfilePlaceholder(ContextID);

    if (!hICC)                          // can't allocate

        return NULL;

    cmsSetDeviceClass(hICC,      cmsSigAbstractClass);

    cmsSetColorSpace(hICC,       cmsSigLabData);

    cmsSetPCS(hICC,              cmsSigLabData);

    cmsSetHeaderRenderingIntent(hICC,  INTENT_PERCEPTUAL);

    // Creates a Pipeline with 3D grid only

    Pipeline = cmsPipelineAlloc(ContextID, 3, 3);

    if (Pipeline == NULL) {

        cmsCloseProfile(hICC);

        return NULL;

    }

    for (i=0; i < MAX_INPUT_DIMENSIONS; i++) Dimensions[i] = nLUTPoints;

    CLUT = cmsStageAllocCLut16bitGranular(ContextID, Dimensions, 3, 3, NULL);

    if (CLUT == NULL) goto Error;