--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/Meerwald-dir/coeff.c	Fri Dec 20 13:08:59 2024 +0100
@@ -0,0 +1,72 @@
+#include "wm.h"
+#include "coeff.h"
+
+double **alloc_coeffs(int cols, int rows) {
+  double **p;
+  int i;
+
+  p = (double **)malloc(rows * sizeof(double *));
+  if (!p) {
+#ifdef DEBUG
+   fprintf(stderr, "alloc_coeffs(): malloc() failed\n");
+   exit(1);
+#else
+   return NULL;
+#endif
+  }
+  p[0] = malloc(rows * cols * sizeof(double));
+  if (!p[0]) {
+#ifdef DEBUG
+    fprintf(stderr, "alloc_coeffs(): malloc() failed\n");
+    exit(1);
+#else
+    free(p);
+    return NULL;
+#endif
+  }
+  for (i = 1; i < rows; i++) {
+    p[i] = &(p[0][i * cols]);
+  }
+
+  return p;
+}
+
+void free_coeffs(double **coeffs) {
+  free(coeffs[0]);
+  free(coeffs);
+}
+
+double **alloc_coeffs_8x8() {
+  return alloc_coeffs(8, 8);
+}
+
+void print_coeffs_8x8(double **coeffs) {
+ int i, j;
+
+  for (i = 0; i < 8; i++) {
+    for (j = 0; j < 8; j++)
+      fprintf(stderr, "%8.2f ", coeffs[i][j]);
+    fprintf(stderr, "\n");
+  }
+}
+
+void print_coeffs(double **coeffs, int c, int r, int w, int h) {
+  int i, j;
+  double *p;
+
+#ifdef DEBUG
+  if (!coeffs) {
+    fprintf(stderr, "print_coeffs(): NULL pixels\n");
+  }
+  if (w <= 0 || h <= 0 || c < 0 || r < 0) {
+    fprintf(stderr, "print_coeffs(): block dimension out of range\n");
+  }
+#endif
+
+  for (j = r; j < r + h; j++) {
+    p = &coeffs[j][c];
+    for (i = 0; i < w; i++)
+      fprintf(stderr, "%8.2f ", *(p++));
+    fprintf(stderr, "\n");
+  }
+}