grid-simulate 0.1.0__tar.gz

grid_simulate-0.1.0/.github/workflows/publish.yml

@@ -0,0 +1,48 @@
+name: Publish to PyPI
+
+on:
+  release:
+    types: [published]
+
+permissions:
+  contents: read
+  id-token: write
+
+jobs:
+  build:
+    runs-on: ubuntu-latest
+
+    steps:
+      - uses: actions/checkout@v4
+
+      - name: Build wheels with maturin
+        uses: PyO3/maturin-action@v1
+        with:
+          command: build
+          args: --release --sdist
+
+      - name: Upload artifacts
+        uses: actions/upload-artifact@v4
+        with:
+          name: dist
+          path: target/wheels/
+
+  publish:
+    needs: build
+    runs-on: ubuntu-latest
+
+    permissions:
+      id-token: write
+
+    environment:
+      name: pypi
+
+    steps:
+      - name: Download artifacts
+        uses: actions/download-artifact@v4
+        with:
+          name: dist
+          path: dist/
+
+      - name: Publish to PyPI
+        uses: pypa/gh-action-pypi-publish@release/v1

grid_simulate-0.1.0/.gitignore

@@ -0,0 +1 @@
+/target

grid_simulate-0.1.0/.idea/.gitignore

@@ -0,0 +1,8 @@
+# Default ignored files
+/shelf/
+/workspace.xml
+# Editor-based HTTP Client requests
+/httpRequests/
+# Datasource local storage ignored files
+/dataSources/
+/dataSources.local.xml

grid_simulate-0.1.0/.idea/GRiD.iml

@@ -0,0 +1,8 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<module type="PYTHON_MODULE" version="4">
+  <component name="NewModuleRootManager">
+    <content url="file://$MODULE_DIR$" />
+    <orderEntry type="jdk" jdkName="Python 3.13" jdkType="Python SDK" />
+    <orderEntry type="sourceFolder" forTests="false" />
+  </component>
+</module>

grid_simulate-0.1.0/.idea/inspectionProfiles/profiles_settings.xml

@@ -0,0 +1,6 @@
+<component name="InspectionProjectProfileManager">
+  <settings>
+    <option name="USE_PROJECT_PROFILE" value="false" />
+    <version value="1.0" />
+  </settings>
+</component>

grid_simulate-0.1.0/.idea/misc.xml

@@ -0,0 +1,7 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<project version="4">
+  <component name="Black">
+    <option name="sdkName" value="Python 3.13" />
+  </component>
+  <component name="ProjectRootManager" version="2" project-jdk-name="Python 3.13" project-jdk-type="Python SDK" />
+</project>

grid_simulate-0.1.0/.idea/modules.xml

@@ -0,0 +1,8 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<project version="4">
+  <component name="ProjectModuleManager">
+    <modules>
+      <module fileurl="file://$PROJECT_DIR$/.idea/GRiD.iml" filepath="$PROJECT_DIR$/.idea/GRiD.iml" />
+    </modules>
+  </component>
+</project>

grid_simulate-0.1.0/.idea/vcs.xml

@@ -0,0 +1,6 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<project version="4">
+  <component name="VcsDirectoryMappings">
+    <mapping directory="" vcs="Git" />
+  </component>
+</project>

grid_simulate-0.1.0/Cargo.lock

@@ -0,0 +1,280 @@
+# This file is automatically @generated by Cargo.
+# It is not intended for manual editing.
+version = 4
+
+[[package]]
+name = "approx"
+version = "0.5.1"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "cab112f0a86d568ea0e627cc1d6be74a1e9cd55214684db5561995f6dad897c6"
+dependencies = [
+ "num-traits",
+]
+
+[[package]]
+name = "autocfg"
+version = "1.5.1"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "f2032f911046de80f0a198e0901378627c33f59ea0ac00e363d481118bd70a53"
+
+[[package]]
+name = "cfg-if"
+version = "1.0.4"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "9330f8b2ff13f34540b44e946ef35111825727b38d33286ef986142615121801"
+
+[[package]]
+name = "grid"
+version = "0.1.0"
+dependencies = [
+ "approx",
+ "ndarray",
+ "numpy",
+ "pyo3",
+]
+
+[[package]]
+name = "heck"
+version = "0.5.0"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "2304e00983f87ffb38b55b444b5e3b60a884b5d30c0fca7d82fe33449bbe55ea"
+
+[[package]]
+name = "indoc"
+version = "2.0.7"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "79cf5c93f93228cf8efb3ba362535fb11199ac548a09ce117c9b1adc3030d706"
+dependencies = [
+ "rustversion",
+]
+
+[[package]]
+name = "libc"
+version = "0.2.186"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "68ab91017fe16c622486840e4c83c9a37afeff978bd239b5293d61ece587de66"
+
+[[package]]
+name = "matrixmultiply"
+version = "0.3.10"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "a06de3016e9fae57a36fd14dba131fccf49f74b40b7fbdb472f96e361ec71a08"
+dependencies = [
+ "autocfg",
+ "rawpointer",
+]
+
+[[package]]
+name = "memoffset"
+version = "0.9.1"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "488016bfae457b036d996092f6cb448677611ce4449e970ceaf42695203f218a"
+dependencies = [
+ "autocfg",
+]
+
+[[package]]
+name = "ndarray"
+version = "0.16.1"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "882ed72dce9365842bf196bdeedf5055305f11fc8c03dee7bb0194a6cad34841"
+dependencies = [
+ "matrixmultiply",
+ "num-complex",
+ "num-integer",
+ "num-traits",
+ "portable-atomic",
+ "portable-atomic-util",
+ "rawpointer",
+]
+
+[[package]]
+name = "num-complex"
+version = "0.4.6"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "73f88a1307638156682bada9d7604135552957b7818057dcef22705b4d509495"
+dependencies = [
+ "num-traits",
+]
+
+[[package]]
+name = "num-integer"
+version = "0.1.46"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "7969661fd2958a5cb096e56c8e1ad0444ac2bbcd0061bd28660485a44879858f"
+dependencies = [
+ "num-traits",
+]
+
+[[package]]
+name = "num-traits"
+version = "0.2.19"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "071dfc062690e90b734c0b2273ce72ad0ffa95f0c74596bc250dcfd960262841"
+dependencies = [
+ "autocfg",
+]
+
+[[package]]
+name = "numpy"
+version = "0.22.1"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "edb929bc0da91a4d85ed6c0a84deaa53d411abfb387fc271124f91bf6b89f14e"
+dependencies = [
+ "libc",
+ "ndarray",
+ "num-complex",
+ "num-integer",
+ "num-traits",
+ "pyo3",
+ "rustc-hash",
+]
+
+[[package]]
+name = "once_cell"
+version = "1.21.4"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "9f7c3e4beb33f85d45ae3e3a1792185706c8e16d043238c593331cc7cd313b50"
+
+[[package]]
+name = "portable-atomic"
+version = "1.13.1"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "c33a9471896f1c69cecef8d20cbe2f7accd12527ce60845ff44c153bb2a21b49"
+
+[[package]]
+name = "portable-atomic-util"
+version = "0.2.7"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "c2a106d1259c23fac8e543272398ae0e3c0b8d33c88ed73d0cc71b0f1d902618"
+dependencies = [
+ "portable-atomic",
+]
+
+[[package]]
+name = "proc-macro2"
+version = "1.0.106"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "8fd00f0bb2e90d81d1044c2b32617f68fcb9fa3bb7640c23e9c748e53fb30934"
+dependencies = [
+ "unicode-ident",
+]
+
+[[package]]
+name = "pyo3"
+version = "0.22.6"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "f402062616ab18202ae8319da13fa4279883a2b8a9d9f83f20dbade813ce1884"
+dependencies = [
+ "cfg-if",
+ "indoc",
+ "libc",
+ "memoffset",
+ "once_cell",
+ "portable-atomic",
+ "pyo3-build-config",
+ "pyo3-ffi",
+ "pyo3-macros",
+ "unindent",
+]
+
+[[package]]
+name = "pyo3-build-config"
+version = "0.22.6"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "b14b5775b5ff446dd1056212d778012cbe8a0fbffd368029fd9e25b514479c38"
+dependencies = [
+ "once_cell",
+ "target-lexicon",
+]
+
+[[package]]
+name = "pyo3-ffi"
+version = "0.22.6"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "9ab5bcf04a2cdcbb50c7d6105de943f543f9ed92af55818fd17b660390fc8636"
+dependencies = [
+ "libc",
+ "pyo3-build-config",
+]
+
+[[package]]
+name = "pyo3-macros"
+version = "0.22.6"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "0fd24d897903a9e6d80b968368a34e1525aeb719d568dba8b3d4bfa5dc67d453"
+dependencies = [
+ "proc-macro2",
+ "pyo3-macros-backend",
+ "quote",
+ "syn",
+]
+
+[[package]]
+name = "pyo3-macros-backend"
+version = "0.22.6"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "36c011a03ba1e50152b4b394b479826cad97e7a21eb52df179cd91ac411cbfbe"
+dependencies = [
+ "heck",
+ "proc-macro2",
+ "pyo3-build-config",
+ "quote",
+ "syn",
+]
+
+[[package]]
+name = "quote"
+version = "1.0.45"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "41f2619966050689382d2b44f664f4bc593e129785a36d6ee376ddf37259b924"
+dependencies = [
+ "proc-macro2",
+]
+
+[[package]]
+name = "rawpointer"
+version = "0.2.1"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "60a357793950651c4ed0f3f52338f53b2f809f32d83a07f72909fa13e4c6c1e3"
+
+[[package]]
+name = "rustc-hash"
+version = "1.1.0"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "08d43f7aa6b08d49f382cde6a7982047c3426db949b1424bc4b7ec9ae12c6ce2"
+
+[[package]]
+name = "rustversion"
+version = "1.0.22"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "b39cdef0fa800fc44525c84ccb54a029961a8215f9619753635a9c0d2538d46d"
+
+[[package]]
+name = "syn"
+version = "2.0.117"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "e665b8803e7b1d2a727f4023456bbbbe74da67099c585258af0ad9c5013b9b99"
+dependencies = [
+ "proc-macro2",
+ "quote",
+ "unicode-ident",
+]
+
+[[package]]
+name = "target-lexicon"
+version = "0.12.16"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "61c41af27dd6d1e27b1b16b489db798443478cef1f06a660c96db617ba5de3b1"
+
+[[package]]
+name = "unicode-ident"
+version = "1.0.24"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "e6e4313cd5fcd3dad5cafa179702e2b244f760991f45397d14d4ebf38247da75"
+
+[[package]]
+name = "unindent"
+version = "0.2.4"
+source = "registry+https://github.com/rust-lang/crates.io-index"
+checksum = "7264e107f553ccae879d21fbea1d6724ac785e8c3bfc762137959b5802826ef3"

grid_simulate-0.1.0/Cargo.toml

@@ -0,0 +1,26 @@
+[package]
+name = "grid"
+version = "0.1.0"
+edition = "2021"
+readme = "README.md"
+
+[lib]
+name = "grid"
+crate-type = ["cdylib"]
+
+[dependencies]
+pyo3 = { version = "0.22", features = ["extension-module"] }
+numpy = "0.22"
+ndarray = "0.16"
+
+[profile.release]
+opt-level = 3
+lto = true
+codegen-units = 1
+
+[dev-dependencies]
+approx = "0.5"
+
+[[bin]]
+name = "grid_sim_demo"
+path = "src/main.rs"

grid_simulate-0.1.0/PKG-INFO

@@ -0,0 +1,5 @@
+Metadata-Version: 2.4
+Name: grid-simulate
+Version: 0.1.0
+Summary: Rust + PyO3 simulation library
+Requires-Python: >=3.8

grid_simulate-0.1.0/README.md

@@ -0,0 +1,42 @@
+# GRiD — Rigidity-Informed Deep Graph Embeddings
+
+Dimension reduction via rigidity theory, with a learned GNN-based 
+independence scorer replacing classical Gaussian elimination.
+
+## Requirements
+- Rust
+- Python 3.13
+
+## Setup
+
+```bash
+# 1. Create and activate a Python 3.13 virtual environment
+python3.13 -m venv .venv
+.venv\Scripts\activate        # Windows
+source .venv/bin/activate     # macOS / Linux
+
+# 2. Install build tool
+pip install maturin
+
+# 3. Compile the Rust simulation and install the Python bridge
+maturin develop
+
+# 4. Verify everything works
+python python/test_bridge.py
+```
+
+## Running the standalone Rust demo (no Python needed)
+```bash
+cargo run
+cargo test
+```
+
+## If you are having problems compiling the code, it might be because of:
+-  Your current python version (pyo3 mismatch)
+
+### To solve it:
+Just before compiling the code:
+```bash
+$env:PYO3_PYTHON = "Your Python 3.13 path"
+```
+And you should be fine

grid_simulate-0.1.0/pyproject.toml

@@ -0,0 +1,12 @@
+[build-system]
+requires = ["maturin>=1.5"]
+build-backend = "maturin"
+
+[project]
+name = "grid-simulate"
+version = "0.1.0"
+description = "Rust + PyO3 simulation library"
+requires-python = ">=3.8"
+
+[tool.maturin]
+features = ["pyo3/extension-module"]

grid_simulate-0.1.0/python/test_bridge.py

@@ -0,0 +1,158 @@
+import grid
+import numpy as np
+
+"""
+This file is both a bridge and a documentary type of file, so that it is easy for
+everyone to understand which methods are available for our custom lib (grid) and
+what each callable does.
+
+author: Sezer
+date: 17/06/2026
+"""
+
+
+# ---------------------------------------------------------------------------
+# Small helpers so every section prints a readable, labeled block.
+# ---------------------------------------------------------------------------
+def banner(title):
+    line = "=" * 78
+    print(f"\n{line}\n  {title}\n{line}")
+
+
+def kv(label, value):
+    """Print a single labeled value."""
+    print(f"  {label:<20}: {value}")
+
+
+################################################################################
+### grid.simulate ###
+"""
+Runs the rigidity simulation on one set of node coordinates
+
+Args:
+    coords    : np.ndarray shape (n, d)    [node coordinates, float64]
+    d         : int                        [embedding dimension, must be 2 or 3]
+    tol       : float                      [rank tolerance, optional(def=1e-9)]
+
+Returns:
+    "labels"             np.ndarray[uint8]       shape (n*(n-1)/2)
+                         1 = independent (keep), 0 = redundant (skip)
+
+    "accepted_edges"     list[(i, j)]            list of kept edges
+    "distance_matrix"    np.ndarray[float64]     shape (n, n)
+    "final_rank"         int
+    "max_rank"           int                     nd - d(d+1)/2
+    "is_rigid"           bool
+    "n_nodes"            int
+    "d"                  int
+"""
+
+banner("grid.simulate  -  rigidity simulation on a single triangle")
+
+coords = np.array([[0, 0], [1, 0], [0.5, 0.866]], dtype=np.float64)
+result = grid.simulate(coords, d=2)
+result = grid.simulate(coords, d=2, tol=1e-6)  # custom tolerance (overrides above)
+
+kv("n_nodes", result["n_nodes"])
+kv("d", result["d"])
+kv("labels", result["labels"])          # 1 = independent (keep), 0 = redundant
+kv("accepted_edges", result["accepted_edges"])
+kv("final_rank", result["final_rank"])
+kv("max_rank", result["max_rank"])      # nd - d(d+1)/2
+kv("is_rigid", result["is_rigid"])
+print("  distance_matrix:")
+print(result["distance_matrix"])
+
+
+################################################################################
+### grid.generate_training_batch ###
+"""
+Generate many random graphs at once for GNN training data
+All heavy work runs with GIL released so torch can train in parallel
+
+Args:
+    n_graphs : int    (how many graphs to generate)
+    n_nodes  : int    (nodes per graph)
+    d        : int    (embedding dimension)
+    seed     : int    (RNG seed (graph k gets seed+k, so the results are reproducible))
+
+Returns:
+    list[dict]    (same dict as simulate(), one per graph)
+"""
+
+banner("grid.generate_training_batch  -  batched random graphs for GNN data")
+
+batch = grid.generate_training_batch(n_graphs=100, n_nodes=20, d=2, seed=42)
+
+kv("graphs generated", len(batch))
+
+# Inspect the first graph so we see the shape of a single sample.
+first = batch[0]
+kv("graph[0] labels shape", first["labels"].shape)
+kv("graph[0] dist shape", first["distance_matrix"].shape)
+kv("graph[0] is_rigid", first["is_rigid"])
+print("  graph[0] labels:")
+print(first["labels"])
+
+
+################################################################################
+### grid.gram_from_dist ###
+"""
+Compute the Gram matrix from a Euclidean distance matrix
+We will use this to verify isometry such that G should be positive semi definite with rank < d
+
+Args:
+    dist : np.ndarray    shape(n, n)    Euclidean distances
+
+Returns:
+    np.ndarray    shape(n, n)    (Gram matrix G = -1/2 H D^2 H)
+"""
+
+banner("grid.gram_from_dist  -  Gram matrix + isometry (PSD) check")
+
+gram = grid.gram_from_dist(result["distance_matrix"])
+eigenvalues = np.linalg.eigvalsh(gram)
+is_isometric = np.all(eigenvalues >= -1e-9)
+
+print("  gram matrix:")
+print(gram)
+kv("eigenvalues", eigenvalues)
+kv("is_isometric (PSD)", is_isometric)  # True => valid Euclidean distance matrix
+
+
+################################################################################
+### grid.max_independent ###
+"""
+Returns the maximum number of independent edges for n nodes in d dimension
+(nd - d(d+1)/2)
+We are going to use this as the stopping condition or to pre-allocate label tensors.
+"""
+
+banner("grid.max_independent  -  max independent edges = nd - d(d+1)/2")
+
+kv("n=10, d=2", grid.max_independent(n=10, d=2))
+kv("n=5,  d=3", grid.max_independent(n=5, d=3))
+
+
+################################################################################
+### grid.reconstruct ###
+"""
+Reconstructs the original data points from the given Euclidean distance matrix
+"""
+
+banner("grid.reconstruct  -  recover coordinates from a distance matrix")
+
+coords = np.array([[0, 0], [1, 0], [0.5, 0.866]], dtype=np.float64)
+sim = grid.simulate(coords, d=2)
+rec = grid.reconstruct(sim["distance_matrix"], d=2)
+
+kv("stress", rec["stress"])              # should be ~0.0
+kv("is_isometric", rec["is_isometric"])  # True
+kv("eigenvalues", rec["eigenvalues"])    # top-2 eigenvalues of Gram matrix
+print("  recovered coords:")
+print(rec["coords"])
+
+
+
+
+