| Index: third_party/qcms/src/tests/qcms_test_output_trc.c
|
| diff --git a/third_party/qcms/src/tests/qcms_test_output_trc.c b/third_party/qcms/src/tests/qcms_test_output_trc.c
|
| index 9e25f8f165f082653902f5bca4795afad93728e8..049141b634b089e9fe89e4496870abfe0c20562e 100644
|
| --- a/third_party/qcms/src/tests/qcms_test_output_trc.c
|
| +++ b/third_party/qcms/src/tests/qcms_test_output_trc.c
|
| @@ -94,9 +94,12 @@ static int qcms_test_output_trc(size_t width,
|
| const char *out_path,
|
| const int force_software)
|
| {
|
| + uint16_t *gamma_table_out = NULL;
|
| + size_t output_size = 0;
|
| qcms_profile *profile;
|
| - uint16_t *gamma_table_out;
|
| - size_t output_size;
|
| + long time_stamp = (long)time(NULL);
|
| + char output_file_name[1024];
|
| + float scale_factor;
|
| size_t i;
|
|
|
| if (!in_path) {
|
| @@ -104,12 +107,9 @@ static int qcms_test_output_trc(size_t width,
|
| return EXIT_FAILURE;
|
| }
|
|
|
| - printf("Test qcms output gamma curve table integrity.\n");
|
| + printf("Test color profile gamma curves\n");
|
| + fflush(stdout);
|
|
|
| - // Create profiles and transforms, get table and then free resources to make sure none
|
| - // of the internal tables are initialized by previous calls.
|
| - gamma_table_out = NULL;
|
| - output_size = 0;
|
| if (get_output_gamma_table(in_path, &gamma_table_out, &output_size) != 0) {
|
| fprintf(stderr, "Unable to extract output gamma table\n");
|
| return EXIT_FAILURE;
|
| @@ -124,60 +124,115 @@ static int qcms_test_output_trc(size_t width,
|
| return EXIT_FAILURE;
|
| }
|
|
|
| - // Check only for red curve for now.
|
| if (profile->redTRC->type == PARAMETRIC_CURVE_TYPE) {
|
| + // Check the red TRC curve only for now.
|
| int type = - (int)(profile->redTRC->count + 1);
|
| uint16_t *gamma_table_in = NULL;
|
| size_t input_size = 0;
|
| - float scale_factor;
|
| FILE *output_file;
|
| - char output_file_name[1024];
|
| - long int time_stamp = (long int)time(NULL);
|
|
|
| printf("Detected parametric curve type = %d\n", profile->redTRC->count);
|
|
|
| - sprintf(output_file_name, "qcms-test-%ld-parametric-gamma-output-%s.csv", time_stamp, profile->description);
|
| + if (get_input_gamma_table(in_path, &gamma_table_in, &input_size) != 0) {
|
| + fprintf(stderr, "Failed to compute input gamma table\n");
|
| + qcms_profile_release(profile);
|
| + free(gamma_table_out);
|
| + return EXIT_FAILURE;
|
| + }
|
| +
|
| + // Write output to stdout and tables into a csv file.
|
| + sprintf(output_file_name, "qcms-test-%ld-parametric-gamma-output-%s.csv",
|
| + time_stamp, profile->description);
|
| printf("Writing output gamma tables to %s\n", output_file_name);
|
| + output_file = fopen(output_file_name, "w");
|
| +
|
| + printf("Parametric gamma values for profile %s description [%s]\n",
|
| + in_path, profile->description);
|
| + fprintf(output_file, "Parametric gamma values for profile %s description [%s]\n",
|
| + in_path, profile->description);
|
|
|
| printf("gamma = %.6f, a = %.6f, b = %.6f, c = %.6f, d = %.6f, e = %.6f, f = %.6f\n",
|
| profile->redTRC->parameter[0], profile->redTRC->parameter[1], profile->redTRC->parameter[2],
|
| profile->redTRC->parameter[3], profile->redTRC->parameter[4], profile->redTRC->parameter[5],
|
| profile->redTRC->parameter[6]);
|
|
|
| - // Write output to stdout and tables into a csv file.
|
| - output_file = fopen(output_file_name, "w");
|
| - fprintf(output_file, "Parametric gamma values for %s\n", profile->description);
|
| fprintf(output_file, "gamma, a, b, c, d, e, f\n");
|
| fprintf(output_file, "%.6f, %.6f, %.6f, %.6f, %.6f, %.6f, %.6f\n",
|
| profile->redTRC->parameter[0], profile->redTRC->parameter[1], profile->redTRC->parameter[2],
|
| profile->redTRC->parameter[3], profile->redTRC->parameter[4], profile->redTRC->parameter[5],
|
| profile->redTRC->parameter[6]);
|
|
|
| - get_input_gamma_table(in_path, &gamma_table_in, &input_size);
|
| - if (!gamma_table_in) {
|
| - fprintf(stderr, "Unable to compute input trc. Aborting\n");
|
| - fclose(output_file);
|
| + fprintf(output_file, "\nInput curve size: %zu", input_size);
|
| + fprintf(output_file, "\nOutput curve size: %zu", output_size);
|
| +
|
| + fprintf(output_file, "\n\nInput gamma, Output gamma, LCMS Output gamma, Output gamma error\n");
|
| + // Output gamma curve down-sample factor.
|
| + scale_factor = (float)(output_size - 1) / (input_size - 1);
|
| +
|
| + for (i = 0; i < input_size; ++i) {
|
| + float input = gamma_table_in[i * 4] * inverse65535;
|
| + size_t out_index = (size_t)floor(i * scale_factor + 0.5);
|
| + float output = gamma_table_out[out_index * 4] * inverse65535;
|
| + float x = out_index / (float)(output_size - 1);
|
| + float reference = clamp_float(evaluate_parametric_curve(type, profile->redTRC->parameter, x));
|
| + float difference = fabs(output - reference);
|
| +
|
| + fprintf(output_file, "%.6f, %.6f, %6f, %6f\n", input, output, reference, difference);
|
| + }
|
| +
|
| + fprintf(output_file, "\nNote: the output gamma curves are down-sampled by a factor of %zu / %zu\n",
|
| + output_size, input_size);
|
| +
|
| + fclose(output_file);
|
| + free(gamma_table_in);
|
| + } else {
|
| + uint16_t *gamma_table_in = NULL;
|
| + size_t input_size = 0;
|
| + FILE *output_file;
|
| +
|
| + if (get_input_gamma_table(in_path, &gamma_table_in, &input_size) != 0) {
|
| + fprintf(stderr, "Failed to compute input gamma table\n");
|
| qcms_profile_release(profile);
|
| free(gamma_table_out);
|
| return EXIT_FAILURE;
|
| }
|
|
|
| + // Write output to stdout and tables into a csv file.
|
| + sprintf(output_file_name, "qcms-test-%ld-gamma-output-%s.csv",
|
| + time_stamp, profile->description);
|
| + printf("Writing gamma tables to %s\n", output_file_name);
|
| + output_file = fopen(output_file_name, "w");
|
| +
|
| + printf("Gamma values for profile %s description [%s]\n",
|
| + in_path, profile->description);
|
| + fprintf(output_file, "Gamma values for profile %s description [%s]\n",
|
| + in_path, profile->description);
|
| +
|
| + if (profile->redTRC->count == 0) {
|
| + printf("Gamma LUT type 0: linear gamma\n");
|
| + fprintf(output_file, "Gamma LUT type 0: linear gamma\n");
|
| + } else if (profile->redTRC->count == 1) {
|
| + float gamma = profile->redTRC->data[0] / 256.0f;
|
| + printf("Gamma LUT type 1: gamma = %.6f\n", gamma);
|
| + fprintf(output_file, "Gamma LUT type 1: gamma = %.6f\n", gamma);
|
| + } else {
|
| + printf("Gamma LUT table size = %u\n", profile->redTRC->count);
|
| + fprintf(output_file, "Gamma LUT table size = %u\n", profile->redTRC->count);
|
| + }
|
| +
|
| fprintf(output_file, "\nInput curve size: %zu", input_size);
|
| fprintf(output_file, "\nOutput curve size: %zu", output_size);
|
|
|
| + fprintf(output_file, "\n\nInput gamma, Output gamma\n");
|
| + // Output gamma curve down-sample factor.
|
| scale_factor = (float)(output_size - 1) / (input_size - 1);
|
|
|
| - fprintf(output_file, "\n\nInput gamma, Output gamma, LCMS Output gamma, Output gamma error\n");
|
| -
|
| for (i = 0; i < input_size; ++i) {
|
| float input = gamma_table_in[i * 4] * inverse65535;
|
| size_t out_index = (size_t)floor(i * scale_factor + 0.5);
|
| - float p = out_index / (float)(output_size - 1);
|
| - float reference = clamp_float(evaluate_parametric_curve(type, profile->redTRC->parameter, p));
|
| - float actual = gamma_table_out[out_index * 4] * inverse65535;
|
| - float difference = fabs(actual - reference);
|
| + float output = gamma_table_out[out_index * 4] * inverse65535;
|
|
|
| - fprintf(output_file, "%.6f, %.6f, %6f, %6f\n", input, actual, reference, difference);
|
| + fprintf(output_file, "%.6f, %.6f\n", input, output);
|
| }
|
|
|
| fprintf(output_file, "\nNote: the output gamma curves are down-sampled by a factor of %zu / %zu\n",
|
|
|