libgap: export some matrix APIs

src/libgap-api.c

@@ -13,6 +13,7 @@
 #include <signal.h>
 
 #include "libgap-api.h"
+#include "libgap_intern.h"
 
 #include "ariths.h"
 #include "bool.h"
@@ -27,6 +28,7 @@
 #include "integer.h"
 #include "lists.h"
 #include "macfloat.h"
+#include "modules.h"
 #include "opers.h"
 #include "plist.h"
 #include "precord.h"
@@ -332,6 +334,59 @@ Obj GAP_NewPlist(Int capacity)
     return NEW_PLIST(T_PLIST_EMPTY, capacity);
 }
 
+
+////
+//// matrix obj
+////
+
+static Obj IsMatrixObjFilt;
+static Obj NrRowsAttr;
+static Obj NrColsAttr;
+
+// Returns 1 if <obj> is a GAP matrix obj, 0 if not.
+int GAP_IsMatrixObj(Obj obj)
+{
+    return obj && DoFilter(IsMatrixObjFilt, obj) == True;
+}
+
+// Returns the number of rows of the given GAP matrix obj.
+// If <mat> is not a GAP matrix obj, an error may be raised.
+UInt GAP_NrRows(Obj mat)
+{
+    Obj nrows = CALL_1ARGS(NrRowsAttr, mat);
+    return UInt_ObjInt(nrows);
+}
+
+// Returns the number of columns of the given GAP matrix obj.
+// If <mat> is not a GAP matrix obj, an error may be raised.
+UInt GAP_NrCols(Obj mat)
+{
+    Obj ncols = CALL_1ARGS(NrColsAttr, mat);
+    return UInt_ObjInt(ncols);
+}
+
+// Assign <val> at position <pos> into the GAP matrix obj <mat>.
+// If <val> is zero, then this unbinds the list entry.
+// If <mat> is not a GAP matrix obj, an error may be raised.
+void GAP_AssMat(Obj mat, UInt row, UInt col, Obj val)
+{
+    Obj r = ObjInt_UInt(row);
+    Obj c = ObjInt_UInt(col);
+    ASS_MAT(mat, r, c, val);
+}
+
+// Returns the element at the <row>, <col> in the GAP matrix obj <mat>.
+// Returns 0 if <row> or <col> are out of bounds, i.e., if either
+// is zero, or larger than the number of rows respectively columns of the list.
+// If <mat> is not a GAP matrix obj, an error may be raised.
+Obj GAP_ElmMat(Obj mat, UInt row, UInt col)
+{
+    Obj r = ObjInt_UInt(row);
+    Obj c = ObjInt_UInt(col);
+    return ELM_MAT(mat, r, c);
+}
+
+
 ////
 //// records
 ////
@@ -465,3 +520,34 @@ void GAP_Error_Postjmp_Returning_(void)
         EnterStackCount = -EnterStackCount;
     }
 }
+
+
+/****************************************************************************
+**
+*F  InitKernel( <module> )  . . . . . . . . initialise kernel data structures
+*/
+static Int InitKernel(StructInitInfo * module)
+{
+    ImportFuncFromLibrary("IsMatrixObj", &IsMatrixObjFilt);
+    ImportFuncFromLibrary("NrRows", &NrRowsAttr);
+    ImportFuncFromLibrary("NrCols", &NrColsAttr);
+
+    return 0;
+}
+
+/****************************************************************************
+**
+*F  InitInfoLibGapApi() . . . . . . . . . . . . . . . table of init functions
+*/
+static StructInitInfo module = {
+    // init struct using C99 designated initializers; for a full list of
+    // fields, please refer to the definition of StructInitInfo
+    .type = MODULE_BUILTIN,
+    .name = "libgap",
+    .initKernel = InitKernel,
+};
+
+StructInitInfo * InitInfoLibGapApi ( void )
+{
+    return &module;
+}

src/libgap-api.h

@@ -325,7 +325,7 @@ UInt GAP_LenList(Obj list);
 // If <list> is not a GAP list, an error may be raised.
 void GAP_AssList(Obj list, UInt pos, Obj val);
 
-// Returns the element at the position <pos> in the list <list>.
+// Returns the element at the position <pos> in the GAP list <list>.
 // Returns 0 if there is no entry at the given position.
 // Also returns 0 if <pos> is out of bounds, i.e., if <pos> is zero,
 // or larger than the length of the list.
@@ -336,6 +336,33 @@ Obj GAP_ElmList(Obj list, UInt pos);
 Obj GAP_NewPlist(Int capacity);
 
 
+////
+//// matrix obj
+////
+
+// Returns 1 if <obj> is a GAP matrix obj, 0 if not.
+int GAP_IsMatrixObj(Obj obj);
+
+// Returns the number of rows of the given GAP matrix obj.
+// If <mat> is not a GAP matrix obj, an error may be raised.
+UInt GAP_NrRows(Obj mat);
+
+// Returns the number of columns of the given GAP matrix obj.
+// If <mat> is not a GAP matrix obj, an error may be raised.
+UInt GAP_NrCols(Obj mat);
+
+// Assign <val> at position <pos> into the GAP matrix obj <mat>.
+// If <val> is zero, then this unbinds the list entry.
+// If <mat> is not a GAP matrix obj, an error may be raised.
+void GAP_AssMat(Obj mat, UInt row, UInt col, Obj val);
+
+// Returns the element at the <row>, <col> in the GAP matrix obj <mat>.
+// Returns 0 if <row> or <col> are out of bounds, i.e., if either
+// is zero, or larger than the number of rows respectively columns of the list.
+// If <mat> is not a GAP matrix obj, an error may be raised.
+Obj GAP_ElmMat(Obj mat, UInt row, UInt col);
+
+
 ////
 //// records
 ////

src/libgap_intern.h

@@ -0,0 +1,29 @@
+/****************************************************************************
+**
+**  This file is part of GAP, a system for computational discrete algebra.
+**
+**  Copyright of GAP belongs to its developers, whose names are too numerous
+**  to list here. Please refer to the COPYRIGHT file for details.
+**
+**  SPDX-License-Identifier: GPL-2.0-or-later
+*/
+
+#ifndef LIBGAP_INTERN_H
+#define LIBGAP_INTERN_H
+
+#include "system.h"
+
+
+/****************************************************************************
+**
+*F * * * * * * * * * * * * * initialize module * * * * * * * * * * * * * * *
+*/
+
+/****************************************************************************
+**
+*F  InitInfoLibGapApi() . . . . . . . . . . . . . . . table of init functions
+*/
+StructInitInfo * InitInfoLibGapApi(void);
+
+
+#endif // LIBGAP_INTERN_H

src/modules_builtin.c

@@ -16,6 +16,7 @@
 #include "info.h"
 #include "intfuncs.h"
 #include "iostream.h"
+#include "libgap_intern.h"
 #include "objccoll.h"
 #include "objset.h"
 #include "profile.h"
@@ -139,6 +140,9 @@ const InitInfoFunc InitFuncsBuiltinModules[] = {
     InitInfoThreadAPI,
     InitInfoAObjects,
     InitInfoSerialize,
+#else
+    // libgap API
+    InitInfoLibGapApi,
 #endif
 
     0

tst/testlibgap/api.c

@@ -9,8 +9,11 @@
  * TODO: Test error handling
  *
  */
+
 #include "common.h"
 
+#include <intobj.h>
+
 
 void records(void)
 {
@@ -56,6 +59,29 @@ void lists(void)
     assert(ret == 0);
 }
 
+void matrices(void)
+{
+    Obj mat, val, row, ret;
+
+    mat = GAP_NewPlist(1);
+    val = INTOBJ_INT(42);
+
+    assert(!GAP_IsMatrixObj(mat));   // empty list, not yet a matrix
+    assert(!GAP_IsMatrixObj(0));
+    assert(!GAP_IsMatrixObj(val));
+
+    row = GAP_NewPlist(2);
+    GAP_AssList(row, 1, INTOBJ_INT(1));
+    GAP_AssList(row, 2, INTOBJ_INT(2));
+    GAP_AssList(mat, 1, row);
+    assert(!GAP_IsMatrixObj(row));
+    assert(GAP_IsMatrixObj(mat));
+
+    GAP_AssMat(mat, 1, 1, val);
+    ret = GAP_ElmMat(mat, 1, 1);
+    assert(ret == val);
+}
+
 void strings(void)
 {
     const char *ts = "Hello, world!";
@@ -175,6 +201,10 @@ int main(int argc, char ** argv)
     lists();
     printf("success\n");
 
+    printf("# Testing matrices... ");
+    matrices();
+    printf("success\n");
+
     printf("# Testing integers... ");
     integers();
     printf("success\n");

tst/testlibgap/api.expect

@@ -2,6 +2,7 @@
 # Testing strings... success
 # Testing records... success
 # Testing lists... success
+# Testing matrices... success
 # Testing integers... success
 # Testing operations... success
 # Testing global variables... success