Files
Donovan Baarda 807520ca1d Update eigen to v5.0.1 and libigl to v2.6.0. (#11311)
* Update eigen from v3.3.7 to v5.0.1.

This updates eigen from v3.3.7 released on  December 11, 2018-12-11 to v5.0.1
released on 2025-11-11. There have be a large number of bug-fixes,
optimizations, and improvements between these releases. See the details at;

https://gitlab.com/libeigen/eigen/-/releases

It retains the previous custom minimal `CMakeLists.txt`, and adds a
README-OrcaSlicer.md that explains what version and parts of the upstream
eigen release have been included, and where the full release can be found.

* Update libigl from v2.0.0 (or older) to v2.6.0.

This updates libigl from what was probably v2.0.0 released on 2018-10-16 to
v2.6.0 released on 2025-05-15. It's possible the old version was even older
than that but there is no version indicators in the code and I ran out of
patience identifying missing changes and only went back as far as v2.0.0.

There have been a large number of bug-fixes, optimizations, and improvements
between these versions. See the following for details;

https://github.com/libigl/libigl/releases

I retained the minimal custom `CMakeLists.txt`, added `README.md` from the
libigl distribution which identifies the version, and added a
README-OrcaSlicer.md that details the version and parts that have been
included.

* Update libslic3r for libigl v2.6.0 changes.

This updates libslic3r for all changes moving to eigen v5.0.1 and libigl
v2.6.0. Despite the large number of updates to both dependencies, no changes
were required for the eigen update, and only one change was required for the
libigl update.

For libigl, `igl::Hit` was changed to a template taking the Scalar type to
use. Previously it was hard-coded to `float`, so to minimize possible impact
I've updated all places it is used from `igl::Hit` to `igl::Hit<float>`.

* Add compiler option `-DNOMINMAX` for libigl with MSVC.

MSVC by default defines `min(()` and `max()` macros that break
`std::numeric_limits<>::max()`. The upstream cmake that we don't include
adds `-DNOMINMAX` for the libigl module when compiling with MSVC, so we need
to add the same thing here.

* Fix src/libslic3r/TriangleMeshDeal.cpp for the unmodified upstream libigl.

This fixes `TriangleMeshDeal.cpp` to work with the unmodified upstream
libigl v2.6.0. loop.{h,cpp} implementation.

This file and feature was added in PR "BBS Port: Mesh Subdivision" (#12150)
which included changes to `loop.{h,cpp}` in the old version of libigl. This PR
avoids modifying the included dependencies, and uses the updated upstream
versions of those files without any modifications, which requires fixing
TriangleMeshDeal.cpp to work with them.

In particular, the modifications made to `loop.{h,cpp}` included changing the
return type from void to bool, adding additional validation checking of the
input meshes, and returning false if they failed validation. These added
checks looked unnecessary and would only have caught problems if the input
mesh was very corrupt.

To make `TriangleMeshDeal.cpp` work without this built-in checking
functionality, I removed checking/handling of any `false` return value.

There was also a hell of a lot of redundant copying and casting back and forth
between float and double, so I cleaned that up. The input and output meshs use
floats for the vertexes, and there would be no accuracy benefits from casting
to and from doubles for the simple weighted average operations done by
igl::loop(). So this just uses `Eigen:Map` to use the original input mesh
vertex data directly without requiring any copy or casting.

* Move eigen from included `deps_src` to externaly fetched `deps`.

This copys what PrusaSlicer did and moved it from an included dependency under
`deps_src` to an externaly fetched dependency under `deps`. This requires
updating some `CMakeList.txt` configs and removing the old and obsolete
`cmake/modules/FindEigen3.cmake`. The details of when this was done in
PrusaSlicer and the followup fixes are at;

* 21116995d7
* https://github.com/prusa3d/PrusaSlicer/issues/13608
* https://github.com/prusa3d/PrusaSlicer/pull/13609
* e3c277b9ee

For some reason I don't fully understand this also required fixing
`src/slic3r/GUI/GUI_App.cpp` by adding `#include <boost/nowide/cstdio.hpp>` to
fix an `error: ‘remove’ is not a member of ‘boost::nowide'`. The main thing I
don't understand is how it worked before. Note that this include is in the
PrusaSlicer version of this file, but it also significantly deviates from what
is currently in OrcaSlicer in many other ways.

* Whups... I missed adding the deps/Eigen/Eigen.cmake file...

* Tidy some whitespace indenting in CMakeLists.txt.

* Ugh... tabs indenting needing fixes.

* Change the include order of deps/Eigen.

It turns out that although Boost includes some references to Eigen, Eigen also
includes some references to Boost for supporting some of it's additional
numeric types.

I don't think it matters much since we are not using these features, but I
think technically its more correct to say Eigen depends on Boost than the
other way around, so I've re-ordered them.

* Add source for Eigen 5.0.1 download to flatpak yml config.

* Add explicit `DEPENDS dep_Boost to deps/Eigen.

I missed this before. This ensures we don't rely on include orders to make
sure Boost is installed before we configure Eigen.

* Add `DEPENDS dep_Boost dep_GMP dep_MPFR` to deps/Eigen.

It turns out Eigen can also use GMP and MPFR for multi-precision and
multi-precision-rounded numeric types if they are available.

Again, I don't think we are using these so it doesn't really matter, but it is
technically correct and ensures they are there if we ever do need them.

* Fix deps DEPENDENCY ordering for GMP, MPFR, Eigen, and CGAL.

I think this is finally correct. Apparently CGAL also optionally depends on
Eigen, so the correct dependency order from lowest to highest is GMP, MPFR, Eigen, and CGAL.

---------

Co-authored-by: Donovan Baarda <dbaarda@google.com>
Co-authored-by: Noisyfox <timemanager.rick@gmail.com>
2026-05-12 20:35:21 +08:00

551 lines
16 KiB
C++

////////////////////////////////////////////////////////////////////////////////
//
// WARNING
//
// THIS DOES NOT DEAL WITH VERTICALLY FLIPPED DATA CORRECTLY
//
////////////////////////////////////////////////////////////////////////////////
/* This file is derived from (actually an earlier version of)... */
/* The GIMP -- an image manipulation program
* Copyright (C) 1995 Spencer Kimball and Peter Mattis
*
* $Id: tga.cpp,v 1.1.2.5 2007-05-10 02:10:07 elif Exp $
* TrueVision Targa loading and saving file filter for the Gimp.
* Targa code Copyright (C) 1997 Raphael FRANCOIS and Gordon Matzigkeit
*
* The Targa reading and writing code was written from scratch by
* Raphael FRANCOIS <fraph@ibm.net> and Gordon Matzigkeit
* <gord@gnu.ai.mit.edu> based on the TrueVision TGA File Format
* Specification, Version 2.0:
*
* <URL:ftp://ftp.truevision.com/pub/TGA.File.Format.Spec/>
*
* It does not contain any code written for other TGA file loaders.
* Not even the RLE handling. ;)
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
*/
#include "tga.h"
#include "../../opengl2/glext.h"
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <assert.h>
static char error[256];
static unsigned int _verbose = 0;
static int totbytes = 0;
typedef struct {
unsigned char *statebuf;
int statelen;
int laststate;
} RLEstate;
IGL_INLINE static int
std_fread(RLEstate * /*rleInfo*/, unsigned char *buf, size_t datasize, size_t nelems, FILE *fp)
{
if (_verbose > 1) {
totbytes += nelems * datasize;
printf("TGA: std_fread %d (total %d)\n",
(int)(nelems * datasize), totbytes);
}
return fread(buf, datasize, nelems, fp);
}
#define MIN(a,b) (((a) < (b)) ? (a) : (b))
#define RLE_PACKETSIZE 0x80
/* Decode a bufferful of file. */
IGL_INLINE static int
rle_fread(RLEstate *rleInfo, unsigned char *vbuf, size_t datasize, size_t nelems, FILE *fp)
{
unsigned char *buf = vbuf;
int j, k;
int buflen, count, bytes, curbytes;
unsigned char *p;
/* Scale the buffer length. */
buflen = nelems * datasize;
j = 0;
curbytes = totbytes;
while (j < buflen) {
if (rleInfo->laststate < rleInfo->statelen) {
/* Copy bytes from our previously decoded buffer. */
bytes = MIN(buflen - j, rleInfo->statelen - rleInfo->laststate);
memcpy(buf + j, rleInfo->statebuf + rleInfo->laststate, bytes);
j += bytes;
rleInfo->laststate += bytes;
/* If we used up all of our state bytes, then reset them. */
if (rleInfo->laststate >= rleInfo->statelen) {
rleInfo->laststate = 0;
rleInfo->statelen = 0;
}
/* If we filled the buffer, then exit the loop. */
if (j >= buflen) break;
}
/* Decode the next packet. */
count = fgetc(fp);
if (count == EOF) {
if (_verbose) printf("TGA: hit EOF while looking for count\n");
return j / datasize;
}
/* Scale the byte length to the size of the data. */
bytes = ((count & ~RLE_PACKETSIZE) + 1) * datasize;
if (j + bytes <= buflen) {
/* We can copy directly into the image buffer. */
p = buf + j;
} else {
#ifdef PROFILE
printf("TGA: needed to use statebuf for %d bytes\n", buflen - j);
#endif
/* Allocate the state buffer if we haven't already. */
if (!rleInfo->statebuf) {
rleInfo->statebuf = (unsigned char *) malloc(RLE_PACKETSIZE * datasize);
}
p = rleInfo->statebuf;
}
if (count & RLE_PACKETSIZE) {
/* Fill the buffer with the next value. */
if (fread(p, datasize, 1, fp) != 1) {
if (_verbose) {
printf("TGA: EOF while reading %d/%d element RLE packet\n",
bytes, (int)datasize);
}
return j / datasize;
}
/* Optimized case for single-byte encoded data. */
if (datasize == 1) {
memset(p + 1, *p, bytes - 1);
} else {
for (k = datasize; k < bytes; k += datasize) {
memcpy(p + k, p, datasize);
}
}
} else {
/* Read in the buffer. */
if (fread(p, bytes, 1, fp) != 1) {
if (_verbose) {
printf("TGA: EOF while reading %d/%d element raw packet\n",
bytes, (int)datasize);
}
return j / datasize;
}
}
if (_verbose > 1) {
totbytes += bytes;
if (_verbose > 2) {
printf("TGA: %s packet %d/%d\n",
(count & RLE_PACKETSIZE) ? "RLE" : "raw",
bytes, totbytes);
}
}
/* We may need to copy bytes from the state buffer. */
if (p == rleInfo->statebuf) {
rleInfo->statelen = bytes;
} else {
j += bytes;
}
}
if (_verbose > 1) {
printf("TGA: rle_fread %d/%d (total %d)\n",
(int) ( nelems * datasize), totbytes - curbytes, totbytes);
}
return nelems;
}
IGL_INLINE igl::opengl::gliGenericImage *
igl::opengl::gliReadTGA(FILE *fp, char *name, int /*hflip*/, int vflip)
{
igl::opengl::TgaHeader tgaHeader;
igl::opengl::TgaFooter tgaFooter;
char horzrev, vertrev;
int width, height, bpp;
int start, end, dir;
int i, j, k;
int pelbytes, wbytes;
GLenum format;
int components;
RLEstate rleRec;
RLEstate *rleInfo;
int rle;
int index, colors, length;
GLubyte *cmap, *pixels, *data;
int (*myfread)(RLEstate *rleInfo, unsigned char*, size_t, size_t, FILE*);
igl::opengl::gliGenericImage *genericImage;
/* Check the footer. */
if (fseek(fp, 0L - sizeof(tgaFooter), SEEK_END)
|| fread(&tgaFooter, sizeof(tgaFooter), 1, fp) != 1) {
sprintf(error, "TGA: Cannot read footer from \"%s\"", name);
if (_verbose) printf("%s\n", error);
return NULL;
}
/* Check the signature. */
if (memcmp(tgaFooter.signature, TGA_SIGNATURE,
sizeof(tgaFooter.signature)) == 0) {
if (_verbose) printf("TGA: found New TGA\n");
} else {
if (_verbose) printf("TGA: found Original TGA\n");
}
if (fseek(fp, 0, SEEK_SET) ||
fread(&tgaHeader, sizeof(tgaHeader), 1, fp) != 1) {
sprintf(error, "TGA: Cannot read header from \"%s\"", name);
if (_verbose) printf("%s\n", error);
return NULL;
}
if (_verbose && tgaHeader.idLength) {
char *idString = (char*) malloc(tgaHeader.idLength);
if (fread(idString, tgaHeader.idLength, 1, fp) != 1) {
sprintf(error, "TGA: Cannot read ID field in \"%s\"", name);
printf("%s\n", error);
} else {
printf("TGA: ID field: \"%*s\"\n", tgaHeader.idLength, idString);
}
free(idString);
} else {
/* Skip the image ID field. */
if (tgaHeader.idLength && fseek(fp, tgaHeader.idLength, SEEK_CUR)) {
sprintf(error, "TGA: Cannot skip ID field in \"%s\"", name);
if (_verbose) printf("%s\n", error);
return NULL;
}
}
/* Reassemble the multi-byte values correctly, regardless of
host endianness. */
width = (tgaHeader.widthHi << 8) | tgaHeader.widthLo;
height = (tgaHeader.heightHi << 8) | tgaHeader.heightLo;
bpp = tgaHeader.bpp;
if (_verbose) {
printf("TGA: width=%d, height=%d, bpp=%d\n", width, height, bpp);
}
horzrev = tgaHeader.descriptor & TGA_DESC_HORIZONTAL;
vertrev = tgaHeader.descriptor & TGA_DESC_VERTICAL;
//vertrev=0;
// // JASON - we can force this stuff if we want
// if( hflip )
// horzrev = 1;
if( vflip )
vertrev = 1;
if (_verbose && horzrev) printf("TGA: horizontal reversed\n");
if (_verbose && vertrev) printf("TGA: vertical reversed\n");
rle = 0;
switch (tgaHeader.imageType) {
case TGA_TYPE_MAPPED_RLE:
rle = 1;
if (_verbose) printf("TGA: run-length encoded\n");
case TGA_TYPE_MAPPED:
/* Test for alpha channel. */
format = GL_COLOR_INDEX;
components = 1;
if (_verbose) {
printf("TGA: %d bit indexed image (%d bit palette)\n",
tgaHeader.colorMapSize, bpp);
}
break;
case TGA_TYPE_GRAY_RLE:
rle = 1;
if (_verbose) printf("TGA: run-length encoded\n");
case TGA_TYPE_GRAY:
format = GL_LUMINANCE;
components = 1;
if (_verbose) printf("TGA: %d bit grayscale image\n", bpp);
break;
case TGA_TYPE_COLOR_RLE:
rle = 1;
if (_verbose) printf("TGA: run-length encoded\n");
case TGA_TYPE_COLOR:
/* Test for alpha channel. */
if (bpp == 32) {
format = GL_BGRA_EXT;
components = 4;
if (_verbose) {
printf("TGA: %d bit color image with alpha channel\n", bpp);
}
} else {
format = GL_BGR_EXT;
components = 3;
if (_verbose) printf("TGA: %d bit color image\n", bpp);
}
break;
default:
sprintf(error,
"TGA: unrecognized image type %d\n", tgaHeader.imageType);
if (_verbose) printf("%s\n", error);
return NULL;
}
if ((format == GL_BGRA_EXT && bpp != 32) ||
(format == GL_BGR_EXT && bpp != 24) ||
((format == GL_LUMINANCE || format == GL_COLOR_INDEX) && bpp != 8)) {
/* FIXME: We haven't implemented bit-packed fields yet. */
fprintf(stderr, "bpp %d, format %x\n", bpp, (unsigned int)format);
sprintf(error, "TGA: channel sizes other than 8 bits are unimplemented");
if (_verbose) printf("%s\n", error);
return NULL;
}
/* Check that we have a color map only when we need it. */
if (format == GL_COLOR_INDEX) {
if (tgaHeader.colorMapType != 1) {
sprintf(error, "TGA: indexed image has invalid color map type %d\n",
tgaHeader.colorMapType);
if (_verbose) printf("%s\n", error);
return NULL;
}
} else if (tgaHeader.colorMapType != 0) {
sprintf(error, "TGA: non-indexed image has invalid color map type %d\n",
tgaHeader.colorMapType);
if (_verbose) printf("%s\n", error);
return NULL;
}
if (tgaHeader.colorMapType == 1) {
/* We need to read in the colormap. */
index = (tgaHeader.colorMapIndexHi << 8) | tgaHeader.colorMapIndexLo;
length = (tgaHeader.colorMapLengthHi << 8) | tgaHeader.colorMapLengthLo;
if (_verbose) {
printf("TGA: reading color map (%d + %d) * (%d / 8)\n",
index, length, tgaHeader.colorMapSize);
}
if (length == 0) {
sprintf(error, "TGA: invalid color map length %d", length);
if (_verbose) printf("%s\n", error);
return NULL;
}
if (tgaHeader.colorMapSize != 24) {
/* We haven't implemented bit-packed fields yet. */
sprintf(error, "TGA: channel sizes other than 8 bits are unimplemented");
if (_verbose) printf("%s\n", error);
return NULL;
}
pelbytes = tgaHeader.colorMapSize / 8;
colors = length + index;
cmap = (GLubyte*)malloc (colors * pelbytes);
/* Zero the entries up to the beginning of the map. */
memset(cmap, 0, index * pelbytes);
/* Read in the rest of the colormap. */
if (fread(cmap, pelbytes, length, fp) != (size_t) length) {
sprintf(error, "TGA: error reading colormap (ftell == %ld)\n",
ftell (fp));
if (_verbose) printf("%s\n", error);
free(cmap);
return NULL;
}
if (pelbytes >= 3) {
/* Rearrange the colors from BGR to RGB. */
int tmp;
for (j = index; j < length * pelbytes; j += pelbytes) {
tmp = cmap[j];
cmap[j] = cmap[j + 2];
cmap[j + 2] = tmp;
}
}
} else {
colors = 0;
cmap = NULL;
}
/* Allocate the data. */
pelbytes = bpp / 8;
pixels = (unsigned char *) malloc (width * height * pelbytes);
if (rle) {
rleRec.statebuf = 0;
rleRec.statelen = 0;
rleRec.laststate = 0;
rleInfo = &rleRec;
myfread = rle_fread;
} else {
rleInfo = NULL;
myfread = std_fread;
}
wbytes = width * pelbytes;
if (vertrev) {
start = 0;
end = height;
dir = 1;
} else {
/* We need to reverse loading order of rows. */
start = height-1;
end = -1;
dir = -1;
}
for (i = start; i != end; i += dir) {
data = pixels + i*wbytes;
/* Suck in the data one row at a time. */
if (myfread(rleInfo, data, pelbytes, width, fp) != width) {
/* Probably premature end of file. */
if (_verbose) {
printf ("TGA: error reading (ftell == %ld, width=%d)\n",
ftell(fp), width);
}
return NULL;
}
if (horzrev) {
/* We need to mirror row horizontally. */
for (j = 0; j < width/2; j++) {
GLubyte tmp;
for (k = 0; k < pelbytes; k++) {
tmp = data[j*pelbytes+k];
data[j*pelbytes+k] = data[(width-j-1)*pelbytes+k];
data[(width-j-1)*pelbytes+k] = tmp;
}
}
}
}
if (rle) {
free(rleInfo->statebuf);
}
if (fgetc (fp) != EOF) {
if (_verbose) printf ("TGA: too much input data, ignoring extra...\n");
}
genericImage = (igl::opengl::gliGenericImage*) malloc(sizeof(igl::opengl::gliGenericImage));
genericImage->width = width;
genericImage->height = height;
genericImage->format = format;
genericImage->components = components;
genericImage->cmapEntries = colors;
genericImage->cmapFormat = GL_BGR_EXT; // XXX fix me
genericImage->cmap = cmap;
genericImage->pixels = pixels;
return genericImage;
}
IGL_INLINE int igl::opengl::gli_verbose(int new_verbose)
{
_verbose = new_verbose;
return _verbose;
}
// added 10/2005, Denis Zorin
// a very simple TGA output, supporting
// uncompressed luminance RGB and RGBA
// G22.2270 students: this is C (no C++)
// so this is not the style I would encourage
// you to use; I used it for consistency
// with the rest of the code in this file
// fixed header values for the subset of TGA we use for writing
unsigned char TGAHeaderColor[12] =
{ 0,// 0 ID length = no id
0,// 1 color map type = no color map
2,// 2 image type = uncompressed true color
0, 0, 0, 0, 0,// color map spec = empty
0, 0, // x origin of image
0, 0 // y origin of image
};
unsigned char TGAHeaderBW[12] =
{ 0,// 0 ID length = no id
0,// 1 color map type = no color map
3,// 3 image type = uncompressed black and white
0, 0, 0, 0, 0,// color map spec = empty
0, 0, // x origin of image
0, 0 // y origin of image
};
// this makes sure that
// image size is written in correct format
// and byte order (least first)
IGL_INLINE void write16bit(int n, FILE* fp) {
unsigned char bytes[] = { static_cast<unsigned char>(n % 256), static_cast<unsigned char>(n / 256) };
fwrite(bytes, 2, sizeof(unsigned char),fp);
}
IGL_INLINE void igl::opengl::writeTGA( igl::opengl::gliGenericImage* image, FILE *fp) {
assert(!image->cmap); // we do not deal with color map images
if(image->components == 3 || image->components == 4)
fwrite(TGAHeaderColor, 12, sizeof(unsigned char),fp);
else {
if(image->components == 1 )
fwrite(TGAHeaderBW, 12, sizeof(unsigned char),fp);
else { fprintf(stderr,"Supported component number: 1,3 or 4\n"); exit(1); }
}
write16bit(image->width,fp);
write16bit(image->height,fp);
switch (image->components ) {
case 1:
putc(8,fp);
break;
case 3:
putc(24,fp);
break;
case 4:
putc(32,fp);
break;
default: fprintf(stderr,"Supported component number: 1,3 or 4\n"); exit(1);
};
if(image-> components == 4)
putc(0x04,fp); // bottom left image (0x00) + 8 bit alpha (0x4)
else
putc(0x00,fp);
fwrite(image->pixels, image->height*image->width*image->components, sizeof(char),fp);
}