gsimplex 0.0.2__tar.gz

gsimplex-0.0.2/PKG-INFO

@@ -0,0 +1,25 @@
+Metadata-Version: 2.4
+Name: gsimplex
+Version: 0.0.2
+Summary: Implementation of simplex algorithm contrelled by the primal-dual gap
+Author-email: Riccardo Ciucci <riccardo@ciucci.dev>
+License-Expression: MIT
+Project-URL: Homepage, https://github.com/Richie314/GapControlledSimplex
+Project-URL: Issues, https://github.com/Richie314/GapControlledSimplex/issues
+Classifier: Programming Language :: Python :: 3
+Classifier: Operating System :: OS Independent
+Requires-Python: >=3.12
+Description-Content-Type: text/markdown
+Requires-Dist: numpy>=2.2.0
+Requires-Dist: aiohttp>=3.9.0
+Requires-Dist: pulp>=3.1.0
+Provides-Extra: dev
+Requires-Dist: pytest; extra == "dev"
+
+# Gap controlled Simplex
+
+[![Test and publish package](https://github.com/Richie314/GapControlledSimplex/actions/workflows/pypi.yml/badge.svg)](https://github.com/Richie314/GapControlledSimplex/actions/workflows/pypi.yml)
+
+![Algorithm diagram](./assets/AlgorithmDiagram.drawio.svg)
+
+## Download test problems

gsimplex-0.0.2/README.md

@@ -0,0 +1,7 @@
+# Gap controlled Simplex
+
+[![Test and publish package](https://github.com/Richie314/GapControlledSimplex/actions/workflows/pypi.yml/badge.svg)](https://github.com/Richie314/GapControlledSimplex/actions/workflows/pypi.yml)
+
+![Algorithm diagram](./assets/AlgorithmDiagram.drawio.svg)
+
+## Download test problems

gsimplex-0.0.2/pyproject.toml

@@ -0,0 +1,44 @@
+[build-system]
+requires = [ "setuptools >= 80.0.1", "wheel" ]
+build-backend = "setuptools.build_meta"
+
+[project]
+name = "gsimplex"
+version = "0.0.2"
+authors = [
+  { name="Riccardo Ciucci", email="riccardo@ciucci.dev" },
+]
+description = "Implementation of simplex algorithm contrelled by the primal-dual gap"
+readme = "README.md"
+requires-python = ">=3.12"
+classifiers = [
+    "Programming Language :: Python :: 3",
+    "Operating System :: OS Independent",
+]
+license = "MIT"
+license-files = ["LICEN[CS]E*"]
+dependencies = [
+    "numpy >= 2.2.0",
+    "aiohttp >= 3.9.0",
+    "pulp >= 3.1.0",
+]
+
+[project.scripts]
+gsimplex = "gsimplex.main:__main"
+gsimplex-demo = "gsimplex.demo:demo"
+gsimplex-download-netlib = "gsimplex.benchmarks.netlib:main"
+gsimplex-emps = "gsimplex.benchmarks.netlib_emps:__main"
+gsimplex-download-plato = "gsimplex.benchmarks.plato:main"
+
+[project.urls]
+Homepage = "https://github.com/Richie314/GapControlledSimplex"
+Issues = "https://github.com/Richie314/GapControlledSimplex/issues"
+
+[project.optional-dependencies]
+dev = ["pytest"]
+
+[tool.setuptools.packages.find]
+where = ["src"]
+
+[tool.setuptools]
+include-package-data = true

gsimplex-0.0.2/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

gsimplex-0.0.2/src/gsimplex/__init__.py

@@ -0,0 +1 @@
+# Simplex library

gsimplex-0.0.2/src/gsimplex/benchmarks/downloader.py

@@ -0,0 +1,76 @@
+import asyncio
+import aiohttp
+from pathlib import Path
+from typing import Dict, List, Tuple, Optional
+
+class Downloader:
+    def __init__(self, benchmark_dir: Optional[str] = None, quiet: bool = False):
+        self._quiet = quiet
+        if benchmark_dir is None:
+            self._benchmark_dir = Path.cwd() / "benchmark"
+        else:
+            self._benchmark_dir = Path(benchmark_dir)
+
+    async def download_async(self, 
+                             url: str, 
+                             filename: str, 
+                             cached_filename: Optional[str] = None,
+                             ) -> Optional[str]:
+        
+        # Make sure the benchmark directory exists
+        self._benchmark_dir.mkdir(parents=True, exist_ok=True)
+
+        if not cached_filename:
+            cached_filename = filename
+
+        filepath = self._benchmark_dir / filename
+        cached_filepath = self._benchmark_dir / cached_filename
+        
+        # If already downloaded, return it
+        if cached_filepath.exists():
+            if not self._quiet:
+                print(f"Using cached: {cached_filename}")
+            return str(cached_filepath)
+        
+        if not self._quiet:
+            print(f"Downloading: {url}...")
+        filepath.parent.mkdir(parents=True, exist_ok=True)
+        try:
+            async with aiohttp.ClientSession() as session:
+                async with session.get(url) as response:
+                    if not response.ok:
+                        if not self._quiet:
+                            print(f"Failed to download {filename}: HTTP {response.status}")
+                        return None
+                    
+                    content = await response.read()
+                    with open(filepath, 'wb') as f:
+                        f.write(content)
+                    
+                    return str(filepath)
+        except Exception as e:
+            if not self._quiet:
+                print(f"Failed to download {url}: {e}")
+            if filepath.exists():
+                try:
+                    filepath.unlink()
+                except:
+                    pass
+            return None
+
+    async def download_many_async(self, 
+                                  files: List[Tuple[str, str, str, Optional[str]]], 
+                                  post_process=None,
+                                  ) -> Dict[str, str]:
+        tasks = [self.download_async(url, filename, cached_filename) for url, problem_name, filename, cached_filename in files]
+        results = await asyncio.gather(*tasks)
+        
+        problem_files = {}
+        for (url, problem_name, filename, cached_filename), path in zip(files, results):
+            if path:
+                if post_process is not None:
+                    path = post_process(path)
+                problem_files[problem_name] = path
+        
+        return problem_files
+

gsimplex-0.0.2/src/gsimplex/benchmarks/netlib.py

@@ -0,0 +1,78 @@
+#!/usr/bin/env python3
+
+import asyncio
+import sys
+import argparse
+from typing import Dict, List, Optional, Tuple
+from pathlib import Path
+
+from gsimplex.benchmarks.downloader import Downloader
+from gsimplex.benchmarks.netlib_emps import expand_mps
+
+class NetLibDownloader(Downloader):
+    BASE_URL = "https://www.netlib.org/lp/data/"
+    
+    async def download_netlib_benchmarks_async(self, problem_names: List[str]) -> Dict[str, str]:
+        files: List[Tuple[str, str, str, Optional[str]]] = [
+            (f"{self.BASE_URL}{name}", name, f"netlib/{name}.mps.netlib", f"netlib/{name}.mps")
+            for name in problem_names
+            if name.strip()
+        ]
+        return await self.download_many_async(files, post_process=NetLibDownloader.post_process_download)
+    
+    @staticmethod
+    def post_process_download(netlib_mps_file: str|Path) -> str:
+        downloaded_file = Path(netlib_mps_file)
+        if not downloaded_file.exists():
+            raise FileNotFoundError(f"Downloaded file path mismatch: {downloaded_file} not found!")
+        
+        download_dir = downloaded_file.parent
+        target_file = download_dir / downloaded_file.name.removesuffix('.netlib')
+
+        expand_mps(str(downloaded_file), str(target_file))
+        downloaded_file.unlink()
+
+        assert target_file.exists()
+        return str(target_file)
+  
+async def download_netlib_benchmarks(dir: Optional[str] = None, quiet: bool = False) -> bool:
+    downloader = NetLibDownloader(benchmark_dir=dir, quiet=quiet)
+    
+    # All Netlib problems
+    problem_names = [
+        "25fv47", "80bau3b", "adlittle", "afiro", "agg", "agg2", "agg3",
+        "bandm", "beaconfd", "blend", "bnl1", "bnl2", "boeing1", "boeing2",
+        "bore3d", "brandy", "capri", "cycle", "czprob", "d2q06c", "d6cube",
+        "degen2", "degen3", "dfl001", "e226", "etamacro", "fffff800",
+        "finnis", "fit1d", "fit1p", "fit2d", "fit2p", "forplan", "ganges",
+        "gfrd-pnc", "greenbea", "greenbeb", "grow7", "grow15", "grow22",
+        "israel", "kb2", "lotfi", "maros", "maros-r7", "modszk1", "nesm",
+        "perold", "pilot", "pilot.ja", "pilot.we", "pilot4", "pilot87",
+        "pilotnov", "recipe", "sc105", "sc205", "sc50a", "sc50b",
+        "scagr25", "scagr7", "scfxm1", "scfxm2", "scfxm3", "scorpion",
+        "scrs8", "scsd1", "scsd6", "scsd8", "sctap1", "sctap2", "sctap3",
+        "seba", "share1b", "share2b", "shell", "ship04l", "ship04s",
+        "ship08l", "ship08s", "ship12l", "ship12s", "sierra", "stair",
+        "standata", "standgub", "standmps", "stocfor1", "stocfor2",
+        "tuff", "vtp.base", "wood1p", "woodw"
+    ]
+    
+    if not quiet:
+        print(f"Downloading {len(problem_names)} Netlib problems...")
+    results = await downloader.download_netlib_benchmarks_async(problem_names)
+    if not quiet:
+        print(f"Downloaded {len(results)} problems successfully")
+
+    return len(results) == len(problem_names)
+
+def main():
+    parser = argparse.ArgumentParser(description="Download Netlib benchmarks")
+    parser.add_argument('--quiet', action='store_true', help='Run in quiet mode')
+    parser.add_argument('--dir', type=str, default=None, help='Directory to save benchmarks')
+    args = parser.parse_args()
+
+    esit = asyncio.run(download_netlib_benchmarks(quiet=args.quiet, dir=args.dir))
+    return 0 if esit else 1
+
+if __name__ == "__main__":
+    sys.exit(main())