--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/mupdf-source/thirdparty/leptonica/src/fmorphgen.1.c	Mon Sep 15 11:43:07 2025 +0200
@@ -0,0 +1,273 @@
+/*====================================================================*
+ -  Copyright (C) 2001 Leptonica.  All rights reserved.
+ -
+ -  Redistribution and use in source and binary forms, with or without
+ -  modification, are permitted provided that the following conditions
+ -  are met:
+ -  1. Redistributions of source code must retain the above copyright
+ -     notice, this list of conditions and the following disclaimer.
+ -  2. Redistributions in binary form must reproduce the above
+ -     copyright notice, this list of conditions and the following
+ -     disclaimer in the documentation and/or other materials
+ -     provided with the distribution.
+ -
+ -  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+ -  ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+ -  LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+ -  A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL ANY
+ -  CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
+ -  EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
+ -  PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
+ -  PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY
+ -  OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
+ -  NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
+ -  SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+ *====================================================================*/
+
+/*!
+ *      Top-level fast binary morphology with auto-generated sels
+ *
+ *             PIX     *pixMorphDwa_1()
+ *             PIX     *pixFMorphopGen_1()
+ */
+
+#include <string.h>
+#include "allheaders.h"
+
+PIX *pixMorphDwa_1(PIX *pixd, PIX *pixs, l_int32 operation, char *selname);
+PIX *pixFMorphopGen_1(PIX *pixd, PIX *pixs, l_int32 operation, char *selname);
+l_int32 fmorphopgen_low_1(l_uint32 *datad, l_int32 w,
+                          l_int32 h, l_int32 wpld,
+                          l_uint32 *datas, l_int32 wpls,
+                          l_int32 index);
+
+static l_int32   NUM_SELS_GENERATED = 58;
+static char  SEL_NAMES[][80] = {
+                             "sel_2h",
+                             "sel_3h",
+                             "sel_4h",
+                             "sel_5h",
+                             "sel_6h",
+                             "sel_7h",
+                             "sel_8h",
+                             "sel_9h",
+                             "sel_10h",
+                             "sel_11h",
+                             "sel_12h",
+                             "sel_13h",
+                             "sel_14h",
+                             "sel_15h",
+                             "sel_20h",
+                             "sel_21h",
+                             "sel_25h",
+                             "sel_30h",
+                             "sel_31h",
+                             "sel_35h",
+                             "sel_40h",
+                             "sel_41h",
+                             "sel_45h",
+                             "sel_50h",
+                             "sel_51h",
+                             "sel_2v",
+                             "sel_3v",
+                             "sel_4v",
+                             "sel_5v",
+                             "sel_6v",
+                             "sel_7v",
+                             "sel_8v",
+                             "sel_9v",
+                             "sel_10v",
+                             "sel_11v",
+                             "sel_12v",
+                             "sel_13v",
+                             "sel_14v",
+                             "sel_15v",
+                             "sel_20v",
+                             "sel_21v",
+                             "sel_25v",
+                             "sel_30v",
+                             "sel_31v",
+                             "sel_35v",
+                             "sel_40v",
+                             "sel_41v",
+                             "sel_45v",
+                             "sel_50v",
+                             "sel_51v",
+                             "sel_2",
+                             "sel_3",
+                             "sel_4",
+                             "sel_5",
+                             "sel_2dp",
+                             "sel_2dm",
+                             "sel_5dp",
+                             "sel_5dm"};
+
+/*!
+ * \brief   pixMorphDwa_1()
+ *
+ * \param[in]    pixd         usual 3 choices: null, == pixs, != pixs
+ * \param[in]    pixs         1 bpp
+ * \param[in]    operation    L_MORPH_DILATE, L_MORPH_ERODE,
+ *                            L_MORPH_OPEN, L_MORPH_CLOSE
+ * \param[in]    sel          name
+ * \return    pixd
+ *
+ * <pre>
+ * Notes:
+ *      (1) This simply adds a border, calls the appropriate
+ *          pixFMorphopGen_*(), and removes the border.
+ *          See the notes for that function.
+ *      (2) The size of the border depends on the operation
+ *          and the boundary conditions.
+ * </pre>
+ */
+PIX *
+pixMorphDwa_1(PIX     *pixd,
+              PIX     *pixs,
+              l_int32  operation,
+              char    *selname)
+{
+l_int32  bordercolor, bordersize;
+PIX     *pixt1, *pixt2, *pixt3;
+
+    if (!pixs)
+        return (PIX *)ERROR_PTR("pixs not defined", __func__, pixd);
+    if (pixGetDepth(pixs) != 1)
+        return (PIX *)ERROR_PTR("pixs must be 1 bpp", __func__, pixd);
+
+        /* Set the border size */
+    bordercolor = getMorphBorderPixelColor(L_MORPH_ERODE, 1);
+    bordersize = 32;
+    if (bordercolor == 0 && operation == L_MORPH_CLOSE)
+        bordersize += 32;
+
+    pixt1 = pixAddBorder(pixs, bordersize, 0);
+    pixt2 = pixFMorphopGen_1(NULL, pixt1, operation, selname);
+    pixt3 = pixRemoveBorder(pixt2, bordersize);
+    pixDestroy(&pixt1);
+    pixDestroy(&pixt2);
+
+    if (!pixd)
+        return pixt3;
+
+    pixCopy(pixd, pixt3);
+    pixDestroy(&pixt3);
+    return pixd;
+}
+
+
+/*!
+ * \brief   pixFMorphopGen_1()
+ *
+ * \param[in]    pixd        usual 3 choices: null, == pixs, != pixs
+ * \param[in]    pixs        1 bpp
+ * \param[in]    operation   L_MORPH_DILATE, L_MORPH_ERODE,
+ *                           L_MORPH_OPEN, L_MORPH_CLOSE
+ * \param[in]    sel name
+ * \return     pixd
+ *
+ * <pre>
+ * Notes:
+ *      (1) This is a dwa operation, and the Sels must be limited in
+ *          size to not more than 31 pixels about the origin.
+ *      (2) A border of appropriate size (32 pixels, or 64 pixels
+ *          for safe closing with asymmetric b.c.) must be added before
+ *          this function is called.
+ *      (3) This handles all required setting of the border pixels
+ *          before erosion and dilation.
+ *      (4) The closing operation is safe; no pixels can be removed
+ *          near the boundary.
+ * </pre>
+ */
+PIX *
+pixFMorphopGen_1(PIX     *pixd,
+                 PIX     *pixs,
+                 l_int32  operation,
+                 char    *selname)
+{
+l_int32    i, index, found, w, h, wpls, wpld, bordercolor, erodeop, borderop;
+l_uint32  *datad, *datas, *datat;
+PIX       *pixt;
+
+    if (!pixs)
+        return (PIX *)ERROR_PTR("pixs not defined", __func__, pixd);
+    if (pixGetDepth(pixs) != 1)
+        return (PIX *)ERROR_PTR("pixs must be 1 bpp", __func__, pixd);
+
+        /* Get boundary colors to use */
+    bordercolor = getMorphBorderPixelColor(L_MORPH_ERODE, 1);
+    if (bordercolor == 1)
+        erodeop = PIX_SET;
+    else
+        erodeop = PIX_CLR;
+
+    found = FALSE;
+    for (i = 0; i < NUM_SELS_GENERATED; i++) {
+        if (strcmp(selname, SEL_NAMES[i]) == 0) {
+            found = TRUE;
+            index = 2 * i;
+            break;
+        }
+    }
+    if (found == FALSE)
+        return (PIX *)ERROR_PTR("sel index not found", __func__, pixd);
+
+    if (!pixd) {
+        if ((pixd = pixCreateTemplate(pixs)) == NULL)
+            return (PIX *)ERROR_PTR("pixd not made", __func__, NULL);
+    }
+    else  /* for in-place or pre-allocated */
+        pixResizeImageData(pixd, pixs);
+    wpls = pixGetWpl(pixs);
+    wpld = pixGetWpl(pixd);
+
+        /* The images must be surrounded, in advance, with a border of
+         * size 32 pixels (or 64, for closing), that we'll read from.
+         * Fabricate a "proper" image as the subimage within the 32
+         * pixel border, having the following parameters:  */
+    w = pixGetWidth(pixs) - 64;
+    h = pixGetHeight(pixs) - 64;
+    datas = pixGetData(pixs) + 32 * wpls + 1;
+    datad = pixGetData(pixd) + 32 * wpld + 1;
+
+    if (operation == L_MORPH_DILATE || operation == L_MORPH_ERODE) {
+        borderop = PIX_CLR;
+        if (operation == L_MORPH_ERODE) {
+            borderop = erodeop;
+            index++;
+        }
+        if (pixd == pixs) {  /* in-place; generate a temp image */
+            if ((pixt = pixCopy(NULL, pixs)) == NULL)
+                return (PIX *)ERROR_PTR("pixt not made", __func__, pixd);
+            datat = pixGetData(pixt) + 32 * wpls + 1;
+            pixSetOrClearBorder(pixt, 32, 32, 32, 32, borderop);
+            fmorphopgen_low_1(datad, w, h, wpld, datat, wpls, index);
+            pixDestroy(&pixt);
+        }
+        else { /* not in-place */
+            pixSetOrClearBorder(pixs, 32, 32, 32, 32, borderop);
+            fmorphopgen_low_1(datad, w, h, wpld, datas, wpls, index);
+        }
+    }
+    else {  /* opening or closing; generate a temp image */
+        if ((pixt = pixCreateTemplate(pixs)) == NULL)
+            return (PIX *)ERROR_PTR("pixt not made", __func__, pixd);
+        datat = pixGetData(pixt) + 32 * wpls + 1;
+        if (operation == L_MORPH_OPEN) {
+            pixSetOrClearBorder(pixs, 32, 32, 32, 32, erodeop);
+            fmorphopgen_low_1(datat, w, h, wpls, datas, wpls, index+1);
+            pixSetOrClearBorder(pixt, 32, 32, 32, 32, PIX_CLR);
+            fmorphopgen_low_1(datad, w, h, wpld, datat, wpls, index);
+        }
+        else {  /* closing */
+            pixSetOrClearBorder(pixs, 32, 32, 32, 32, PIX_CLR);
+            fmorphopgen_low_1(datat, w, h, wpls, datas, wpls, index);
+            pixSetOrClearBorder(pixt, 32, 32, 32, 32, erodeop);
+            fmorphopgen_low_1(datad, w, h, wpld, datat, wpls, index+1);
+        }
+        pixDestroy(&pixt);
+    }
+
+    return pixd;
+}
+