CIE: Style fix

Fallout from 3788e861e545d43d5911677250c0f55d5f7ccedf

CIE: Add an SSE2 version of "RGBA float" to "CIE L float"

On an Intel i7 Haswell, it now takes 0.056s to convert a 15 megapixel
buffer from "RGBA float" to "CIE L float" instead of the earlier
0.107s.

https://bugzilla.gnome.org/show_bug.cgi?id=795686

CIE: Add an SSE2 version of "RGBA float" to "CIE Lab alpha float"

On an Intel i7 Haswell, it now takes 0.13s to convert a 15 megapixel
buffer from "RGBA float" to "CIE Lab alpha float" instead of the
earlier 0.27s.

SSEx doesn't have integer multiplication or division operations, and
using bit shifts to implement integer divisions by powers of 2 seems to
introduce errors. Therefore, it was problematic to use the cube root
approximation from Hacker's Delight, which uses quite a few integer
divisions to make the initial guess. Instead, Halley's method of
approximating the cube root seems more SSEx friendly because the
initial guess requires only one integer division, which we can manage
by jumping through a relatively small number of hoops.

The scalar version of Halley's method seems to have originated from
http://metamerist.com/cbrt/cbrt.htm but that's not accessible anymore.
At present there's a copy in CubeRoot.cpp in the Skia sources that's
licensed under a BSD-style license. There's some discussion on the
implementation at http://www.voidcn.com/article/p-gpwztojr-wt.html.

Note that Darktable also has an SSE2 version of the same algorithm,
but uses only a single iteration of Halley's method, which is too
coarse.

Here's some more discussion on the cube root approximation algorithms:
https://bugzilla.gnome.org/show_bug.cgi?id=791837

https://bugzilla.gnome.org/show_bug.cgi?id=795686

build: Add scaffolding for SSE3

It's currently unused, but is useful to have it around when comparing
different vectorization strategies.

https://bugzilla.gnome.org/show_bug.cgi?id=795686

CIE: Make the XYZ to LAB code look similar to the RGB to LAB code

Such superficial differences make the code harder to scan when looking
for bugs and optimization opportunities.

babl-internal: print debug output to stderr

extensions/two-table: remove unused variable

type-u8: fix u8 <> double conversions for chroma type

postrelease version bump

release 0.1.46

NEWS: update

tools: add html variant of path matrix

docs: remove short-cut coverage, which is too expensive to generate

extensions: remove many inline statements

When the functions are not reused by other fast path definitions, they are only
ever used as function pointers and inline is not neccesary.

extensions/double: adds some more paths

extensions: rename some extensions

extensions: add more u32 short cut utility conversions

extensions: u32<->float conversions in simple extension

They were using 'float' rather than 'double' data types for intermediate values
the lack of precision made the conversions be discarded by babl.

extensions: add missing 32bit.c

extensions: add some u32 fast paths

float-half: half fast paths for premultiplied alpha

extensions: add a double extension

This will permit slightly faster than the generic fast paths for some possibly
desired intermediate conversions. This also adds alpha stripping - which fills
in the last babl fast path currently reported missing on every launch of GIMP.

extensions: simplify gggl/conv_rgbF_rgbaF slightly

extensions/cairo: optimize R'G'B'A u8 to cairo-ARGB32

A combination of SIMD paths via float was beating the naive implementation in
the cairo extension, now the cairo extension comes out on top again.

extension: add some Y' u8 fast paths

GIMP makes use of "Y' u8" to "R'G'B' u8" and we had no paths covering it,
making it go through floating point. This commit adds code paths for expanding
single grayscale u8 to R'G'B and R'G'B'A as well as R'aG'aB'aA u8.

extensions/gegl-fixups: swap increments of src and dst in rgba8_rgb8

They were reversed, causing memory corruptions fixing bugs 794973 and 794972.

extensions: add rgba8_rgb8 conversion that copies 24/32bit chunks

Slightly speeding up RGB u8 (without alpha) painting in GIMP, the rgba to rgb
and reverse conversions are used together with the very fast SIMD conversions
that operate on chunks of 4 components.

post release version bump

release 0.1.44

babl: fix minor leak reported in bug #790924

Added mingw template pkgbuild.

Fix windows dirs_sep. Cleanup.

Remove check for linux header on win32

Remove os_win32

Temporarily disable checking for headers

PKGBUILD points again to gnome upstream

Fix : libdir for autotools and meson are not the same…

babl: fix icc matrix profile fast paths

For core members, babl relies on shared data layout between fish and
conversion to be able to do dispatch, and have source/target match -
though only for simple fishes.

The ICC fast path conversions also used the user data - but directly
instead of throug the argument, changing this makes the matrix and
luts work again - and let the conversions be used.

meson: Add missing single quote

Fixes build since ad0eac19c7a1f0046979ad8923f6ac1b7c2d7696 ("release 0.1.42")

Signed-off-by: Jan Vesely <jano.vesely@gmail.com>

post release version bump

release 0.1.42

babl: initialize the bpp of fishes loaded from cache

The cache loading code manually constructs path fishes, which did not
get their bpp initialized, this resulted in conversions of long runs to
fail, not crash - but stop processing after MAX_BUFFER_SIZE pixels.

post release version bump

NEWS: update

release 0.1.40

babl: add a test iteration multiplier

On my system many of the conversions got measured to costs of 3 4 and 5, these
makes a big room for measurement errors, to increase reliability the tests are
now run 4 times - making it a ranking of contenders with fewer collisions.

This will slightly increase incurred delayed on search for conversions.

babl: make pixel-count instrumentation opt-in

Checking if an integer is 0/1 is lower overhead than incrementing a long stored
kt a memory location. As a side effect, the entries in
~/.cache/babl/babl-fishes will not be sorted by most used to least used unless
the BABL_INSTRUMENT environment variable is set to a value.

babl: set user_data ptr for quick dispatch of 1 step path

babl: indentation and code cleanup

babl: use memcpy dpsatch path for all identity transforms

babl: reduce amount of instrumentation

The integer counter for number of processings overflows in some scenarios, thus
cannot be relied on for more than being informative.

We might ultimately want to also drop the instrumentation.

babl: store pointer to where user data resides instead of copy

This will make the scenario of changing user data work, we do it
unconditionally to avoid introducing a branch.

babl: refactor code to get rid of more branches in babl_process

In the best case scenario now, babl_process increments instrumentation
counters and directly calls the relevant registered fast path function,
with no additional call stack frames in-between.

babl: use babl_conversion_process where appropriate

Instead of overloading babl_process() with two tasks, both fishes and
conversions, use the conversion specific API for conversions directly.

babl: simplify branching in babl_process dispatch

babl: factor out some branches from path fish dispatch

babl: refactor conversion dispatch to do less branches at runtime

babl.map: update symbols for meson build

babl: add babl_process_rows api

Useful when fetching buffers smaller than a tile or partial tiles in APIs like
gegl_buffer_get.

sse2-float: use same TRC for "leftover" samples

The SSE2 versions of the gamma <-> linear TRC functions can produce
slightly different results than the scalar ones.  Consequently,
last commit didn't resolve all discrepancies between the main and
leftover samples (and didn't address the gamma -> linear
direction).

Use the same SSE2 TRC functions for converting the leftover
samples, rather than the scalar functions, to fix that.

sse2-float: use babl_linear_to_gamma_2_2f() (instead of double)

Use the single-precision babl_linear_to_gamma_2_2f() function,
instead of the double-precision version, when converting "leftover"
samples, so that the result is consistent with the vectorized
samples.

util, sse2-float: make sure 1 maps to 1 in linear -> gamma conversions

Add a small offset to the result of the linear -> gamma
conversions, such that an input value of 1.0 maps to an output
value of 1.0, exactly.

babl: shuffle some dispatch functions to be inlineable

babl: strip out some internal length returning code

babl: prepare babl_fish_process for direct calling

palette: expect palettes/formats to be gamma corrected in 8bit

A correctness correction related to bug #763581 the formats wanted
by GIMP and expected are gamma corrected not linear, I expect this
fix either does not much for GIMP or fixes an unknown 8bit linear
vs gamma quantization problem that has been lurking in INDEXED mode.

babl: improve babl_trc_formula_srgb sRGB detection constants

For bug #791816. Use constants as passed from GIMP when its built in sRGB ICC
profile is passed.

CIE: Use a faster cbrtf implementation

This is the approximate cube root of an IEEE float implementation from
Hacker's Delight. The elimination of all conditional branches probably
makes it a better candidate for future SIMD accelerated code paths.

On an Intel i7 Haswell, it now takes 0.27s to convert a 15 megapixel
buffer from "RGBA float" to "CIE Lab alpha float" instead of the
earlier 0.35s. A "Y float" to "CIE L float" conversion takes 0.085s
instead of 0.102s.

Original code: http://www.hackersdelight.org/hdcodetxt/acbrt.c.txt
Permissions: http://www.hackersdelight.org/permissions.htm

https://bugzilla.gnome.org/show_bug.cgi?id=791837

extensions-gggl: clamp when converting float to 16bit

Improve the quality of a conversion that was picked in error for bug #790658

babl: re-enable earlier bailing for big path errors

This is one way of dealing with 790658 - that in general is a good idea as
well.

meson.build: babl's license is LGPL3+

Simpler library handling. PKGBUILD+release, lto.

fix check version

Fix library name

Added PKGBUILD for archlinux, based on AUR.

Remove useless join_paths arg

use find_program. Use env for multiple commands as it allows to split into string list.

use env instead of bash

Use formatting instead of weird path reconstruction

Dispatch the custom_target in two custom_targets for cleaner dependency declaration.

Fix push_web, add index_html generation. TODO: simplify this someday

Add linker script for symbols export

some small cleanup, export script for scp

add docs. fix install libs. Seems finished.~

add extensions, tools

Add tests

Fix build, links.

Added basic meson files

babl: make conversion to model double more generic

Iterate over the components of the model, and fetch the relevant
component from the source format or use a filler value if not found.

Some changes by Debarshi Ray.

https://bugzilla.gnome.org/show_bug.cgi?id=790860

CIE: refactor away some intermediate variables/assignments

CIE: add CIE XYZ and CIE XYZ alpha models

Even more useful now with custom RGB primaries - and being able to use babl for
conversions when needed is convenient. This is with reference paths - which
means fast paths for double precision float - single precision would have to be
added for it to be used and higher throughput in most common scenarios.

tests: not all UNIX-like OSes have a libpthread.

In particular Android systems don't need to link with -lpthread
(actually the link would fail with "ld: cannot find -lpthread").
Use the $(THREAD_LIB) variable which is set correctly during configure
since commit c02af82.

CIE: Simplify code

Exploit the fact that fy is not needed to calculate the a and b
components. This is consistent with the other monochrome conversions.

configure.ac: post release version bump

configure.ac: release 0.1.38

NEWS: update

babl_pow_24f(): add missing float suffix

CIE: Add "CIE L float"

Some of these conversions will be leveraged by gegl:shadows-highlights
which needs to go from "Y float" or "YaA float" to "CIE L float".

The conversion from "RGBA float" was added to aid "YaA float" to
"CIE L float" fishes. They go via:
  "YaA float"     to "RaGaBaA float"
  "RaGaBaA float" to "RGBA float"
  "RGBA float"    to "Y float"
  "Y float"       to "CIE L float"

A direct conversion from "YaA float" to "Y float" in simple C is
hindered by the need to check every pixel's alpha value to avoid
dividing by zero. The pipeline stalls make it lose out to the look-up
table and SIMD based RGB conversions to unassociated alpha.

However, we can trivially cut out the third step and still reduce some
memory traffic.

https://bugzilla.gnome.org/show_bug.cgi?id=790111

babl: fix memory leaks on error paths

The database-focused functions of babl-cache are prone to
memory leaks when errors occur.

Signed-off-by: Tobias Stoeckmann <tobias@stoeckmann.org>

extensions: protect another conversion against /0.0