1 files changed, 145 insertions, 0 deletions
diff --git a/srgb-gamma.h b/srgb-gamma.h
new file mode 100644
index 0000000..3388497
--- /dev/null
+++ b/srgb-gamma.h
@@ -0,0 +1,145 @@
+/* See LICENSE file for copyright and license details. */
+
+static inline long double
+gamma_ullf(long double t)
+{
+	return GAMMA(long double, powl);
+}
+
+static inline double
+gamma_ulf(double t)
+{
+	return GAMMA(double, pow);
+}
+
+static inline float
+gamma_uf(float t)
+{
+	return GAMMA(float, powf);
+}
+
+static inline long double
+gamma_llf(long double t)
+{
+	return t >= 0 ? gamma_ullf(t) : -gamma_ullf(-t);
+}
+
+static inline double
+gamma_lf(double t)
+{
+	return t >= 0 ? gamma_ulf(t) : -gamma_ulf(-t);
+}
+
+static inline float
+gamma_f(float t)
+{
+	return t >= 0 ? gamma_uf(t) : -gamma_uf(-t);
+}
+
+static inline uint64_t
+gamma_u64(uint64_t t)
+{
+	return (uint64_t)(gamma_ullf((long double)t / UINT64_MAX) * UINT64_MAX);
+}
+
+static inline uint32_t
+gamma_u32(uint32_t t)
+{
+	return (uint32_t)(gamma_ulf((double)t / UINT32_MAX) * UINT32_MAX);
+}
+
+static inline uint16_t
+gamma_u16(uint32_t t)
+{
+	return (uint16_t)(gamma_ulf((double)t / UINT16_MAX) * UINT16_MAX);
+}
+
+void
+FUNCTION_NAME(struct libskrift_image *image, size_t x, size_t y, size_t width, size_t height)
+{
+	struct format_settings settings;
+	uint8_t *img = image->image;
+	size_t linesize, psize;
+
+	settings = libskrift_format_settings[image->format];
+	psize    = settings.spsize * (size_t)(4 - (settings.apos <= -2));
+	linesize = (size_t)image->width * psize;
+	img     += y * linesize + x * psize;
+	width   *= psize;
+
+	switch (settings.float_type) {
+	case 0:
+		switch (settings.spsize) {
+		case 1:
+			for (y = 0; y < height; y++, img += linesize) {
+				for (x = 0; x < width; x += psize) {
+					((uint8_t *)&img[x])[settings.rpos] = gamma_u8[((uint8_t *)&img[x])[settings.rpos]];
+					((uint8_t *)&img[x])[settings.gpos] = gamma_u8[((uint8_t *)&img[x])[settings.gpos]];
+					((uint8_t *)&img[x])[settings.bpos] = gamma_u8[((uint8_t *)&img[x])[settings.bpos]];
+				}
+			}
+			break;
+
+		case 2:
+			for (y = 0; y < height; y++, img += linesize) {
+				for (x = 0; x < width; x += psize) {
+					((uint16_t *)&img[x])[settings.rpos] = gamma_u16(((uint16_t *)&img[x])[settings.rpos]);
+					((uint16_t *)&img[x])[settings.gpos] = gamma_u16(((uint16_t *)&img[x])[settings.gpos]);
+					((uint16_t *)&img[x])[settings.bpos] = gamma_u16(((uint16_t *)&img[x])[settings.bpos]);
+				}
+			}
+			break;
+
+		case 4:
+			for (y = 0; y < height; y++, img += linesize) {
+				for (x = 0; x < width; x += psize) {
+					((uint32_t *)&img[x])[settings.rpos] = gamma_u32(((uint32_t *)&img[x])[settings.rpos]);
+					((uint32_t *)&img[x])[settings.gpos] = gamma_u32(((uint32_t *)&img[x])[settings.gpos]);
+					((uint32_t *)&img[x])[settings.bpos] = gamma_u32(((uint32_t *)&img[x])[settings.bpos]);
+				}
+			}
+			break;
+
+		default:
+			for (y = 0; y < height; y++, img += linesize) {
+				for (x = 0; x < width; x += psize) {
+					((uint64_t *)&img[x])[settings.rpos] = gamma_u64(((uint64_t *)&img[x])[settings.rpos]);
+					((uint64_t *)&img[x])[settings.gpos] = gamma_u64(((uint64_t *)&img[x])[settings.gpos]);
+					((uint64_t *)&img[x])[settings.bpos] = gamma_u64(((uint64_t *)&img[x])[settings.bpos]);
+				}
+			}
+			break;
+		}
+		break;
+
+	case 1:
+		for (y = 0; y < height; y++, img += linesize) {
+			for (x = 0; x < width; x += psize) {
+				((float *)&img[x])[settings.rpos] = gamma_f(((float *)&img[x])[settings.rpos]);
+				((float *)&img[x])[settings.gpos] = gamma_f(((float *)&img[x])[settings.gpos]);
+				((float *)&img[x])[settings.bpos] = gamma_f(((float *)&img[x])[settings.bpos]);
+			}
+		}
+		break;
+
+	case 2:
+		for (y = 0; y < height; y++, img += linesize) {
+			for (x = 0; x < width; x += psize) {
+				((double *)&img[x])[settings.rpos] = gamma_lf(((double *)&img[x])[settings.rpos]);
+				((double *)&img[x])[settings.gpos] = gamma_lf(((double *)&img[x])[settings.gpos]);
+				((double *)&img[x])[settings.bpos] = gamma_lf(((double *)&img[x])[settings.bpos]);
+			}
+		}
+		break;
+
+	default:
+		for (y = 0; y < height; y++, img += linesize) {
+			for (x = 0; x < width; x += psize) {
+				((long double *)&img[x])[settings.rpos] = gamma_llf(((long double *)&img[x])[settings.rpos]);
+				((long double *)&img[x])[settings.gpos] = gamma_llf(((long double *)&img[x])[settings.gpos]);
+				((long double *)&img[x])[settings.bpos] = gamma_llf(((long double *)&img[x])[settings.bpos]);
+			}
+		}
+		break;
+	}
+}