build_rgb_colorimetric_cache()

bool fl::colorimetric_response::build_rgb_colorimetric_cache ( const EmitterProfile & p, RgbColorimetricCache * cache )

inline

Definition at line 465 of file colorimetric_response.h.

                                                                                   {
    xyY_to_XYZ(p.xy_r[0], p.xy_r[1], p.lum_r, cache->P_R);
    xyY_to_XYZ(p.xy_g[0], p.xy_g[1], p.lum_g, cache->P_G);
    xyY_to_XYZ(p.xy_b[0], p.xy_b[1], p.lum_b, cache->P_B);
 
    const float P_RGB[3][3] = {
        { cache->P_R[0], cache->P_G[0], cache->P_B[0] },
        { cache->P_R[1], cache->P_G[1], cache->P_B[1] },
        { cache->P_R[2], cache->P_G[2], cache->P_B[2] },
    };
    const bool ok_rgb = invert3x3(P_RGB, cache->P_RGB_inv);
    const bool declares_source_space = (p.input_xy_w[1] > 1e-6f);
    cache->has_source_space = declares_source_space
                           && build_source_matrix(p.input_xy_r, p.input_xy_g,
                                                  p.input_xy_b, p.input_xy_w,
                                                  cache->M_src);
    if (cache->has_source_space && ok_rgb) {
        float X_w[3];
        xyY_to_XYZ(p.input_xy_w[0], p.input_xy_w[1], 1.0f, X_w);
        float t[3];
        matvec3(cache->P_RGB_inv, X_w, t);
        const float mt = fl::max(fl::max(fl::fabs(t[0]), fl::fabs(t[1])),
                                 fl::fabs(t[2]));
        if (mt > 1e-6f) {
            const float scale_k = 1.0f / mt;
            for (int i = 0; i < 3; ++i) {
                for (int j = 0; j < 3; ++j) {
                    cache->M_src[i][j] *= scale_k;
                }
            }
        }
    }
    if (!ok_rgb) {
        for (int i = 0; i < 3; ++i) {
            for (int j = 0; j < 3; ++j) {
                cache->P_RGB_inv[i][j] = 0.0f;
            }
        }
    }
    if (!cache->has_source_space) {
        for (int i = 0; i < 3; ++i) {
            for (int j = 0; j < 3; ++j) {
                cache->M_src[i][j] = 0.0f;
            }
        }
    }
    if (declares_source_space && !cache->has_source_space) {
        return false;
    }
    return ok_rgb;
}

References build_source_matrix(), fl::fabs(), fl::FL_NO_EXCEPT, invert3x3(), matvec3(), fl::max(), fl::t, and xyY_to_XYZ().

Here is the call graph for this function: