Merge pull request #7 from wegank/raglib-install-1

Add installation script

Author: mohabsafey

README.md

@@ -45,14 +45,14 @@ Source code and installation instructions for
 ## Basic usage
 
 At the moment, `RAGlib` provides two main functions:
-- `HasRealSolutions(eqs, pos, ineqs)` where eqs, pos, ineqs are lists
+- `RAG[HasRealSolutions](eqs, pos, ineqs)` where eqs, pos, ineqs are lists
     of polynomials encoding equality, positivity and non-vanishing
     constraints respectively.   
     It decides whether the set of real solutions to the input
     constraints is empty or not. In case of emptiness, it returns an
     empty list, otherwise, it returns a list of witness points.  
     All such points are encoded by isolating boxes. 
-- `PointsPerComponents(eqs, pos, ineqs)` where eqs, pos, ineqs are lists
+- `RAG[PointsPerComponents](eqs, pos, ineqs)` where eqs, pos, ineqs are lists
     of polynomials encoding equality, positivity and non-vanishing
     constraints respectively.  
     It returns a list of points (encoded with isolating boxes) meeting
@@ -61,7 +61,7 @@ At the moment, `RAGlib` provides two main functions:
 
 For instance, the call 
 
-    PointsPerComponents([x*y-1], [x^2+y^2-4], []);
+    RAG[PointsPerComponents]([x*y-1], [x^2+y^2-4], []);
 
 returns 
 
@@ -80,7 +80,7 @@ connected components.
 
 For instance, the call 
 
-    PointsPerComponents([], [], [(x-1)^2+y^2-1, (x+1)^2+y^2-1-1/2^32]);
+    RAG[PointsPerComponents]([], [], [(x-1)^2+y^2-1, (x+1)^2+y^2-1-1/2^32]);
 
 returns 
 

install.mm

@@ -0,0 +1,37 @@
+# Installation script for the RAG package.
+# Run from the repository root with:
+# maple -q install.mm
+
+homedir:=kernelopts(homedir):
+savelibname:=cat(homedir, "/libs/"):
+mladirname:=cat(savelibname, "rag.mla"):
+ssystem(cat("rm -f ", mladirname)):
+libname:=savelibname,libname:
+
+try
+  with(MSolve):
+catch:
+  printf("Error: failed to load MSolve. Please install the file interface first.\n"):
+  quit:
+end try:
+
+raglib_root:=currentdir():
+raglib_src_dir:=cat(raglib_root, "/src"):
+kernelopts(includepath=raglib_src_dir):
+
+raglib_src:=cat(raglib_src_dir, "/rag.mm"):
+try
+  read raglib_src:
+catch:
+  printf("Error: failed to load RAGlib sources. Please run the installation script from the repository root.\n"):
+  quit:
+end try:
+
+if not(evalb(whattype(eval(RAG)) = package)) then
+  printf("Error: failed to load RAG package.\n"):
+  quit:
+end if:
+
+march(`create`, mladirname):
+savelib(`RAG`):
+printf("Installed RAG package into %s\n", mladirname):

src/rag-main.mm

@@ -0,0 +1,54 @@
+
+(** This file is part of RAGlib (Real Algebraic Geometry Library).
+ *
+ * RAGlib 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.
+ *
+ * RAGlib 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 RAGlib.  If not, see <https://www.gnu.org/licenses/>
+ *
+ * Authors:
+ * Mohab Safey El Din **)
+
+#Exported functions of RAGlib
+
+HasRealSolutions:=proc(eqs, pos, ineqs, opts:={})
+local isempty, _l, newopts;
+
+  if type(subs(opts, "isempty"), integer) then 
+    isempty:=subs(opts, "isempty");
+  else 
+    isempty:=0:
+  end if;
+
+  if isempty = 0 then 
+    return SemiAlgebraicSolve(eqs, pos, ineqs, opts union {"isempty"=1});
+  else 
+    newopts:=map(_l->if lhs(_l)="isempty" then lhs(_l)=1 else _l fi, opts);
+    return SemiAlgebraicSolve(eqs, pos, ineqs, newopts);
+  end if;
+end proc:
+
+PointsPerComponents:=proc(eqs, pos, ineqs, opts:={})
+local isempty, _l, newopts;
+
+  if type(subs(opts, "isempty"), integer) then 
+    isempty:=subs(opts, "isempty");
+  else 
+    isempty:=0:
+  end if;
+
+  if isempty = 0 then 
+    return SemiAlgebraicSolve(eqs, pos, ineqs, opts);
+  else 
+    newopts:=map(_l->if lhs(_l)="isempty" then lhs(_l)=0 else _l fi, opts);
+    return SemiAlgebraicSolve(eqs, pos, ineqs, newopts);
+  end if;
+end proc: