The purpose of the GMlib is to provide programmers of solvers or automated meshers in the field of scientific computing with an easy, fast and transparent way to port their codes on GPUs (Graphic Processing Units).
This library is based on the OpenCL language standard, thus taking advantage of almost every architectures supported by most platforms (Linux, macOS, Windows).
It is a simple loop parallelization scheme (known as kernels in the realm of GPU computing).
Provides the programer with pre defined mesh data structures.
Automatically vectorizes unstructured data like the ball of points or the edge shells.
It requires some knowledge on OpenCL programing, which akin to C and C .
Handles transparently the transfer and vectorization of mesh data structures.

Simply follow these steps:

• unarchive the ZIP file
• cd GMlib-master
• mkdir build
• cd build
• cmake ../
• sudo make install

• You first need to install CMake. Do not forget to choose "add cmake to the path for all users", from the install panel.
• Then you need a valid C compiler like the free Visual Studio Community 2019
• unarchive the ZIP file
• from the VS-2019 menu, open the x64 Native Tools Command Prompt
• cd GMlib-master
• mkdir build
• cd build
• cmake -DCMAKE_INSTALL_PREFIX=%HOMEPATH%\local -DCMAKE_BUILD_TYPE=Release ..\
• cmake --build . --target INSTALL

Optionally, you may download libMeshb to run the examples:

• you need to install the libMeshb from GitHub
• cd to /usr/local/GMlib/sample_meshes/
• uncompress them with lzip -d *.meshb.lz
• you may now enter /usr/local/GMlib/examples directory and run the various examples

The GMlib is written in ANSI C with some parts in OpenCL.
It is made of a single C file and a header file to be compiled and linked alongside the calling program.
It may be used in C and C programs (Fortran 77 and 90 APIs are under way).
Tested on Linux, macOS, Windows 7-10.

Here is a basic example that computes some triangles' barycenters on a GPU:

First the "C" part executed by the host CPU:

// Init the GMLIB with the first available GPU on the system
LibIdx = GmlInit(1);

// Create a vertex and a triangle datatype
VerIdx = GmlNewMeshData(LibIdx, GmlVertices,  NmbVer);
TriIdx = GmlNewMeshData(LibIdx, GmlTriangles, NmbTri);

// Fill the vertices with your mesh coordinates
for(i=0;i<NmbVer;i  )
   GmlSetDataLine(LibIdx, VerIdx, i, coords[i][0], coords[i][1], coords[i][2], VerRef[i]);

// Do the same with the elements
for(i=0;i<NmbTri;i  )
   GmlSetDataLine(LibIdx, TriIdx, i, TriVer[i][0], TriVer[i][1], TriVer[i][2], TriRef[i]);

// Create a raw datatype to store the calculated elements' centers
MidIdx = GmlNewSolutionData(LibIdx, GmlTriangles, 1, GmlFlt4, "TriMid");

// Compile the OpenCL source code with the two needed datatypes:
// the vertex coordinates (read) and the triangles centers (write)
CalMid = GmlCompileKernel( LibIdx,  TriangleCenter, "CalMid", GmlTriangles, 2,
                           VerIdx, GmlReadMode,  NULL,
                           MidIdx, GmlWriteMode, NULL );

// Launch the kernel on the GPU
GmlLaunchKernel(LibIdx, CalMid);

// Get the results back from the GPU and print it
for(i=0;i<NmbTri;i  )
{
   GmlGetDataLine(LibIdx, MidIdx, i, MidTab[i]);
   printf("triangle %d center = %g %g %g\n", i, MidTab[i][0], MidTab[i][1], MidTab[i][2]);
}

Then the "OpenCL" part executed by the GPU device:

TriMid = (TriCrd[0]   TriCrd[1]   TriCrd[2]) / 3.;