gumath 0.2.0dev5

data/ext/ruby_gumath/gumath/libgumath/Makefile.in

@@ -0,0 +1,170 @@
+
+# ==============================================================================
+#                          Unix Makefile for libxnd
+# ==============================================================================
+
+
+LIBSTATIC = @LIBSTATIC@
+LIBNAME = @LIBNAME@
+LIBSONAME = @LIBSONAME@
+LIBSHARED = @LIBSHARED@
+
+CC = @CC@
+LD = @LD@
+AR = @AR@
+RANLIB = @RANLIB@
+
+GM_INCLUDES = @CONFIGURE_INCLUDES@
+
+CONFIGURE_CFLAGS = @CONFIGURE_CFLAGS@
+GM_CFLAGS = $(strip -I.. -I$(GM_INCLUDES) $(CONFIGURE_CFLAGS) $(CFLAGS))
+GM_CFLAGS_SHARED = $(GM_CFLAGS) -fPIC
+
+CONFIGURE_LDFLAGS = @CONFIGURE_LDFLAGS@
+GM_LDFLAGS = $(strip $(CONFIGURE_LDFLAGS) $(LDFLAGS))
+
+
+default: $(LIBSTATIC) $(LIBSHARED)
+
+
+OBJS = apply.o func.o nploops.o tbl.o thread.o xndloops.o unary.o binary.o \
+       examples.o bfloat16.o graph.o quaternion.o pdist.o
+
+SHARED_OBJS = .objs/apply.o .objs/func.o .objs/nploops.o .objs/tbl.o .objs/thread.o .objs/xndloops.o \
+              .objs/unary.o .objs/binary.o .objs/examples.o .objs/bfloat16.o .objs/graph.o \
+              .objs/quaternion.o .objs/pdist.o
+
+
+$(LIBSTATIC): Makefile $(OBJS)
+	$(AR) rc $(LIBSTATIC) $(OBJS)
+	$(RANLIB) $(LIBSTATIC)
+
+$(LIBSHARED): Makefile $(SHARED_OBJS)
+	$(LD) $(GM_LDFLAGS) -o $(LIBSHARED) $(SHARED_OBJS)
+	ln -sf $(LIBSHARED) $(LIBNAME)
+	ln -sf $(LIBSHARED) $(LIBSONAME)
+
+
+apply.o:\
+Makefile apply.c gumath.h
+	$(CC) $(GM_CFLAGS) -c apply.c
+
+.objs/apply.o:\
+Makefile apply.c gumath.h
+	$(CC) $(GM_CFLAGS_SHARED) -c apply.c -o .objs/apply.o
+
+func.o:\
+Makefile func.c gumath.h
+	$(CC) $(GM_CFLAGS) -c func.c
+
+.objs/func.o:\
+Makefile func.c gumath.h
+	$(CC) $(GM_CFLAGS_SHARED) -c func.c -o .objs/func.o
+
+nploops.o:\
+Makefile nploops.c gumath.h
+	$(CC) $(GM_CFLAGS) -c nploops.c
+
+.objs/nploops.o:\
+Makefile nploops.c gumath.h
+	$(CC) $(GM_CFLAGS_SHARED) -c nploops.c -o .objs/nploops.o
+
+tbl.o:\
+Makefile tbl.c gumath.h
+	$(CC) $(GM_CFLAGS) -c tbl.c
+
+.objs/tbl.o:\
+Makefile tbl.c gumath.h
+	$(CC) $(GM_CFLAGS_SHARED) -c tbl.c -o .objs/tbl.o
+
+thread.o:\
+Makefile thread.c gumath.h
+	$(CC) $(GM_CFLAGS) -c thread.c
+
+.objs/thread.o:\
+Makefile thread.c gumath.h
+	$(CC) $(GM_CFLAGS_SHARED) -c thread.c -o .objs/thread.o
+
+xndloops.o:\
+Makefile xndloops.c gumath.h
+	$(CC) $(GM_CFLAGS) -c xndloops.c
+
+.objs/xndloops.o:\
+Makefile xndloops.c gumath.h
+	$(CC) $(GM_CFLAGS_SHARED) -c xndloops.c -o .objs/xndloops.o
+
+unary.o:\
+Makefile kernels/unary.c gumath.h
+	$(CC) -I. $(GM_CFLAGS) -Wno-absolute-value -c kernels/unary.c
+
+.objs/unary.o:\
+Makefile kernels/unary.c gumath.h
+	$(CC) -I. $(GM_CFLAGS_SHARED) -Wno-absolute-value -c kernels/unary.c -o .objs/unary.o
+
+binary.o:\
+Makefile kernels/binary.c gumath.h
+	$(CC) -I. $(GM_CFLAGS) -c kernels/binary.c
+
+.objs/binary.o:\
+Makefile kernels/binary.c gumath.h
+	$(CC) -I. $(GM_CFLAGS_SHARED) -c kernels/binary.c -o .objs/binary.o
+
+examples.o:\
+Makefile extending/examples.c gumath.h
+	$(CC) -I. $(GM_CFLAGS) -c extending/examples.c -o examples.o
+
+.objs/examples.o:\
+Makefile extending/examples.c gumath.h
+	$(CC) -I. $(GM_CFLAGS_SHARED) -c extending/examples.c -o .objs/examples.o
+
+bfloat16.o:\
+Makefile extending/bfloat16.c gumath.h
+	$(CC) -I. $(GM_CFLAGS) -c extending/bfloat16.c -o bfloat16.o
+
+.objs/bfloat16.o:\
+Makefile extending/bfloat16.c gumath.h
+	$(CC) -I. $(GM_CFLAGS_SHARED) -c extending/bfloat16.c -o .objs/bfloat16.o
+
+graph.o:\
+Makefile extending/graph.c gumath.h
+	$(CC) -I. $(GM_CFLAGS) -c extending/graph.c -o graph.o
+
+.objs/graph.o:\
+Makefile extending/graph.c gumath.h
+	$(CC) -I. $(GM_CFLAGS_SHARED) -c extending/graph.c -o .objs/graph.o
+
+quaternion.o:\
+Makefile extending/quaternion.c gumath.h
+	$(CC) -I. $(GM_CFLAGS) -c extending/quaternion.c -o quaternion.o
+
+.objs/quaternion.o:\
+Makefile extending/quaternion.c gumath.h
+	$(CC) -I. $(GM_CFLAGS_SHARED) -c extending/quaternion.c -o .objs/quaternion.o
+
+pdist.o:\
+Makefile extending/pdist.c gumath.h
+	$(CC) -I. $(GM_CFLAGS) -c extending/pdist.c -o pdist.o
+
+.objs/pdist.o:\
+Makefile extending/pdist.c gumath.h
+	$(CC) -I. $(GM_CFLAGS_SHARED) -c extending/pdist.c -o .objs/pdist.o
+
+
+# Coverage
+coverage:\
+Makefile clean runtest
+	./tests/runtest
+	for file in *.c; do gcov -l "$$file" > /dev/null 2>&1; done
+
+FORCE:
+
+clean: FORCE
+	rm -f *.o *.so *.gch *.gcda *.gcno *.gcov *.dyn *.dpi *.lock
+	rm -f $(LIBSTATIC) $(LIBSHARED) $(LIBSONAME) $(LIBNAME)
+	cd .objs && rm -f *.o *.so *.gch *.gcda *.gcno *.gcov *.dyn *.dpi *.lock
+
+distclean: clean
+	rm -f Makefile
+
+
+

data/ext/ruby_gumath/gumath/libgumath/Makefile.vc

@@ -0,0 +1,160 @@
+
+# ======================================================================
+#                  Visual C (nmake) Makefile for libgumath
+# ======================================================================
+
+LIBSTATIC = libgumath-0.2.0dev3.lib
+LIBIMPORT = libgumath-0.2.0dev3.dll.lib
+LIBSHARED = libgumath-0.2.0dev3.dll
+LIBNDTYPESIMPORT = libndtypes-0.2.0dev3.dll.lib
+LIBXNDIMPORT = libxnd-0.2.0dev3.dll.lib
+
+!ifndef LIBNDTYPESINCLUDE
+LIBNDTYPESINCLUDE = ..\ndtypes\libndtypes
+!endif
+
+!ifndef LIBNDTYPESDIR
+LIBNDTYPESDIR = ..\ndtypes\libndtypes
+!endif
+
+!ifndef LIBXNDINCLUDE
+LIBXNDINCLUDE = ..\xnd\libxnd
+!endif
+
+!ifndef LIBXNDDIR
+LIBXNDDIR = ..\xnd\libxnd
+!endif
+
+OPT = /MT /Ox /GS /EHsc /fp:strict
+OPT_SHARED = /DGM_EXPORT /DNDT_IMPORT /DXND_IMPORT /MD /Ox /GS /EHsc /fp:strict /Fo.objs^\
+
+COMMON_CFLAGS = /nologo /W4 /wd4200 /wd4201 /wd4204
+CFLAGS = $(COMMON_CFLAGS) $(OPT)
+CFLAGS_SHARED = $(COMMON_CFLAGS) $(OPT_SHARED)
+
+
+default: $(LIBSTATIC) $(LIBSHARED)
+	copy /y gumath.h ..\python\gumath
+	copy /y $(LIBSTATIC) ..\python\gumath
+	copy /y $(LIBIMPORT) ..\python\gumath
+	copy /y $(LIBSHARED) ..\python\gumath
+
+
+OBJS = apply.obj func.obj nploops.obj tbl.obj xndloops.obj \
+       unary.obj binary.obj examples.obj graph.obj pdist.obj
+
+SHARED_OBJS = .objs/apply.obj .objs/func.obj .objs/nploops.obj .objs/tbl.obj .objs/xndloops.obj \
+              .objs/unary.obj .objs/binary.obj .objs/examples.obj .objs/graph.obj .objs/pdist.obj
+
+
+$(LIBSTATIC):\
+Makefile $(OBJS)
+	-@if exist $@ del $(LIBSTATIC)
+	lib /nologo /out:$(LIBSTATIC) $(OBJS)
+
+$(LIBSHARED):\
+Makefile $(SHARED_OBJS)
+	-@if exist $@ del $(LIBSHARED)
+	link /nologo /DLL /MANIFEST /out:$(LIBSHARED) /implib:$(LIBIMPORT) $(SHARED_OBJS) \
+             "/LIBPATH:$(LIBNDTYPESDIR)" $(LIBNDTYPESIMPORT) \
+             "/LIBPATH:$(LIBXNDDIR)" $(LIBXNDIMPORT)
+	mt /nologo -manifest $(LIBSHARED).manifest -outputresource:$(LIBSHARED);2
+
+
+apply.obj:\
+Makefile apply.c gumath.h
+	$(CC) "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS) -c apply.c
+
+.objs\apply.obj:\
+Makefile apply.c gumath.h
+	$(CC) "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS_SHARED) -c apply.c
+
+func.obj:\
+Makefile func.c gumath.h
+	$(CC) "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS) -c func.c
+
+.objs\func.obj:\
+Makefile func.c gumath.h
+	$(CC) "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS_SHARED) -c func.c
+
+nploops.obj:\
+Makefile nploops.c gumath.h
+	$(CC) "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS) -c nploops.c
+
+.objs\nploops.obj:\
+Makefile nploops.c gumath.h
+	$(CC) "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS_SHARED) -c nploops.c
+
+tbl.obj:\
+Makefile tbl.c gumath.h
+	$(CC) "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS) -c tbl.c
+
+.objs\tbl.obj:\
+Makefile tbl.c gumath.h
+	$(CC) "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS_SHARED) -c tbl.c
+
+xndloops.obj:\
+Makefile xndloops.c gumath.h
+	$(CC) "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS) -c xndloops.c
+
+.objs\xndloops.obj:\
+Makefile xndloops.c gumath.h
+	$(CC) "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS_SHARED) -c xndloops.c
+
+unary.obj:\
+Makefile kernels\unary.c gumath.h
+	$(CC) -I. "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS) -c kernels\unary.c
+
+.objs\unary.obj:\
+Makefile kernels\unary.c gumath.h
+	$(CC) -I. "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS_SHARED) -c kernels\unary.c
+
+binary.obj:\
+Makefile kernels\binary.c gumath.h
+	$(CC) -I. "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS) -c kernels\binary.c
+
+.objs\binary.obj:\
+Makefile kernels\binary.c gumath.h
+	$(CC) -I. "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS_SHARED) -c kernels\binary.c
+
+examples.obj:\
+Makefile extending\examples.c gumath.h
+	$(CC) -I. "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS) -c extending\examples.c
+
+.objs\examples.obj:\
+Makefile extending\examples.c gumath.h
+	$(CC) -I. "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS_SHARED) -c extending\examples.c
+
+graph.obj:\
+Makefile extending\graph.c gumath.h
+	$(CC) -I. "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS) -c extending\graph.c
+
+.objs\graph.obj:\
+Makefile extending\graph.c gumath.h
+	$(CC) -I. "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS_SHARED) -c extending\graph.c
+
+pdist.obj:\
+Makefile extending\pdist.c gumath.h
+	$(CC) -I. "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS) -c extending\pdist.c
+
+.objs\pdist.obj:\
+Makefile extending\pdist.c gumath.h
+	$(CC) -I. "-I$(LIBNDTYPESINCLUDE)" "-I$(LIBXNDINCLUDE)" $(CFLAGS_SHARED) -c extending\pdist.c
+
+
+FORCE:
+
+clean: FORCE
+	del /q /f *.exe *.obj *.lib *.dll *.exp *.manifest 2>NUL
+	cd .objs && del /q /f *.obj 2>NUL
+	if exist "..\build\" rd /q /s "..\build\"
+	if exist "..\dist\" rd /q /s "..\dist\"
+	if exist "..\MANIFEST" del "..\MANIFEST"
+	if exist "..\record.txt" del "..\record.txt"
+	cd ..\python\gumath && del *.lib *.dll *.pyd gumath.h 2>NUL
+
+distclean: clean
+	del Makefile 2>NUL
+
+
+

data/ext/ruby_gumath/gumath/libgumath/apply.c

@@ -0,0 +1,201 @@
+/*
+* BSD 3-Clause License
+*
+* Copyright (c) 2017-2018, plures
+* 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.
+*
+* 3. Neither the name of the copyright holder nor the names of its
+*    contributors may be used to endorse or promote products derived from
+*    this software without specific prior written permission.
+*
+* 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 THE COPYRIGHT HOLDER OR 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.
+*/
+
+
+#include <stdlib.h>
+#include <stdint.h>
+#include <string.h>
+#include <inttypes.h>
+#include "ndtypes.h"
+#include "xnd.h"
+#include "gumath.h"
+
+
+static int
+sum_inner_dimensions(const xnd_t stack[], int nargs, int outer_dims)
+{
+    int sum = 0, n;
+    int i;
+
+    for (i = 0; i < nargs; i++) {
+        const ndt_t *t = stack[i].type;
+        n = t->ndim - outer_dims;
+        sum += n == 0 ? 1 : n;
+    }
+
+    return sum;
+}
+
+int
+gm_apply(const gm_kernel_t *kernel, xnd_t stack[], int outer_dims,
+         ndt_context_t *ctx)
+{
+    const int nargs = (int)kernel->set->sig->Function.nargs;
+
+    switch (kernel->flag) {
+    case NDT_ELEMWISE_1D: {
+        if (outer_dims == 0) {
+            ndt_err_format(ctx, NDT_RuntimeError,
+                "gm_xnd_map: optimized elementwise kernel called with "
+                "outer_dims==0");
+            return -1;
+        }
+
+        return gm_xnd_map(kernel->set->Opt, stack, nargs, outer_dims-1, ctx);
+    }
+
+    case NDT_C: {
+        return gm_xnd_map(kernel->set->C, stack, nargs, outer_dims, ctx);
+    }
+
+    case NDT_FORTRAN: {
+        return gm_xnd_map(kernel->set->Fortran, stack, nargs, outer_dims, ctx);
+    }
+
+    case NDT_XND: {
+        return gm_xnd_map(kernel->set->Xnd, stack, nargs, outer_dims, ctx);
+    }
+
+    case NDT_STRIDED: {
+        const int sum_inner = sum_inner_dimensions(stack, nargs, outer_dims);
+        const int dims_size = outer_dims + sum_inner;
+        const int steps_size = nargs * outer_dims + sum_inner;
+        ALLOCA(char *, args, nargs);
+        ALLOCA(intptr_t, dimensions, dims_size);
+        ALLOCA(intptr_t, steps, steps_size);
+
+        if (gm_np_convert_xnd(args, nargs,
+                              dimensions, dims_size,
+                              steps, steps_size,
+                              stack, outer_dims, ctx) < 0) {
+            return -1;
+        }
+
+        return gm_np_map(kernel->set->Strided, args, nargs,
+                         dimensions, steps, NULL, outer_dims);
+      }
+    }
+
+    /* NOT REACHED: tags should be exhaustive. */
+    ndt_internal_error("invalid tag");
+}
+
+static gm_kernel_t
+select_kernel(const ndt_apply_spec_t *spec, const gm_kernel_set_t *set,
+              ndt_context_t *ctx)
+{
+    gm_kernel_t kernel = {0U, NULL};
+
+    kernel.set = set;
+
+    if (set->Opt != NULL && (spec->flags&NDT_ELEMWISE_1D)) {
+        kernel.flag = NDT_ELEMWISE_1D;
+        return kernel;
+    }
+
+    if (set->C != NULL && (spec->flags&NDT_C)) {
+        kernel.flag = NDT_C;
+        return kernel;
+    }
+
+    if (set->Fortran != NULL && (spec->flags&NDT_FORTRAN)) {
+        kernel.flag = NDT_FORTRAN;
+        return kernel;
+    }
+
+    if (set->Strided != NULL && (spec->flags&NDT_STRIDED)) {
+        kernel.flag = NDT_STRIDED;
+        return kernel;
+    }
+
+    if (set->Xnd != NULL && (spec->flags&NDT_XND)) {
+        kernel.flag = NDT_XND;
+        return kernel;
+    }
+
+    kernel.set = NULL;
+    ndt_err_format(ctx, NDT_RuntimeError,
+        "could not find specialized kernel for '%s' input (available: %s, %s, %s, %s)",
+        ndt_apply_flags_as_string(spec),
+        set->C ? "C" : "_",
+        set->Fortran ? "Fortran" : "_",
+        set->Xnd ? "Xnd" : "_",
+        set->Strided ? "Strided" : "_");
+
+    return kernel;
+}
+
+/* Look up a multimethod by name and select a kernel. */
+gm_kernel_t
+gm_select(ndt_apply_spec_t *spec, const gm_tbl_t *tbl, const char *name,
+          const ndt_t *in_types[], int nin, const xnd_t args[],
+          ndt_context_t *ctx)
+{
+    gm_kernel_t empty_kernel = {0U, NULL};
+    const gm_func_t *f;
+    char *s;
+    int i;
+
+    f = gm_tbl_find(tbl, name, ctx);
+    if (f == NULL) {
+        return empty_kernel;
+    }
+
+    if (f->typecheck != NULL) {
+        const gm_kernel_set_t *set = f->typecheck(spec, f, in_types, nin, ctx);
+        if (set == NULL) {
+            return empty_kernel;
+        }
+        return select_kernel(spec, set, ctx);
+    }
+
+    for (i = 0; i < f->nkernels; i++) {
+        const gm_kernel_set_t *set = &f->kernels[i];
+        if (ndt_typecheck(spec, set->sig, in_types, nin, set->constraint, args,
+                          ctx) < 0) {
+            ndt_err_clear(ctx);
+            continue;
+        }
+        return select_kernel(spec, set, ctx);
+    }
+
+    s = ndt_list_as_string(in_types, nin, ctx);
+    if (s == NULL) {
+        return empty_kernel;
+    }
+
+    ndt_err_format(ctx, NDT_TypeError,
+        "could not find '%s' kernel for input types '%s'", name, s);
+    ndt_free(s);
+
+    return empty_kernel;
+}