langevin 2025.11.8a2__tar.gz → 2025.11.10a2__tar.gz

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/.github/workflows/doxygen-gh-pages.yml

@@ -3,10 +3,10 @@ name: Doxygen gh-pages deploy action
 on:
   release:
     types: [published]
-  push:
-    branches:
-      # - main
-      - none
+  # push:
+  #   branches:
+  #     - main
+  #     - none
   # workflow_dispatch:
 
 jobs:
@@ -19,6 +19,6 @@ jobs:
         with:
           github_token: ${{ secrets.GITHUB_TOKEN }}
           branch: gh-pages
-          folder: docs/html
+          # folder: docs/html
           config_file: Doxyfile
           doxygen_version: 1.14.0

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/.github/workflows/testpypi-publish.yml

@@ -4,8 +4,8 @@ name: Upload Python Package to TestPyPI
 #   https://cibuildwheel.pypa.io/en/stable/ci-services/#
 
 on:
-  release:
-    types: [published]
+  # release:
+  #   types: [published]
   push:
     branches:
       - main

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: langevin
-Version: 2025.11.8a2
+Version: 2025.11.10a2
 Summary: Tools to integrate Langevin equations of absorbing phase transition type
 Keywords: Langevin equation,directed percolation,operator splitting,stochastic differential equations,phase transition
 Author-Email: "Colin P. Stark" <cstarkjp@gmail.com>
@@ -31,8 +31,11 @@ Description-Content-Type: text/markdown
 
 Tools to integrate Langevin equations of absorbing phase transition (APT) type — with a focus on simulation of directed percolation (DP).
 
-![](https://raw.githubusercontent.com/cstarkjp/Langevin/main/images/density_grid.png
- "Density grid")
+![](https://raw.githubusercontent.com/cstarkjp/Langevin/main/images/ρ_a1p18950_b1_D0p04_η1_x100_y50_Δx1_Δt0p1_rs1.gif
+ "Density field evolution over time")
+
+ <!-- ![](https://raw.githubusercontent.com/cstarkjp/Langevin/main/images/density_grid.png
+ "Density grid") -->
 
 The `langevin` package implements the operator-splitting method originally developed by Dornic et al (2005), Pechenik & Levine (1999) and others, and improved upon by Weissmann et al (2018).
 It provides a Python wrapper around core C++ heavily adapted from a code base written by [Paula Villa Martín](https://github.com/pvillamartin), extended by [Victor Buendía](https://github.com/VictorSeven) ("VMB"), and arising from earlier efforts by Ivan Dornic and Juan Bonachela. The wrapper provides easy access to the Langevin integrator, and broad opportunity to experiment, adapt, and extend it further. 

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/README.md

@@ -2,8 +2,11 @@
 
 Tools to integrate Langevin equations of absorbing phase transition (APT) type — with a focus on simulation of directed percolation (DP).
 
-![](https://raw.githubusercontent.com/cstarkjp/Langevin/main/images/density_grid.png
- "Density grid")
+![](https://raw.githubusercontent.com/cstarkjp/Langevin/main/images/ρ_a1p18950_b1_D0p04_η1_x100_y50_Δx1_Δt0p1_rs1.gif
+ "Density field evolution over time")
+
+ <!-- ![](https://raw.githubusercontent.com/cstarkjp/Langevin/main/images/density_grid.png
+ "Density grid") -->
 
 The `langevin` package implements the operator-splitting method originally developed by Dornic et al (2005), Pechenik & Levine (1999) and others, and improved upon by Weissmann et al (2018).
 It provides a Python wrapper around core C++ heavily adapted from a code base written by [Paula Villa Martín](https://github.com/pvillamartin), extended by [Victor Buendía](https://github.com/VictorSeven) ("VMB"), and arising from earlier efforts by Ivan Dornic and Juan Bonachela. The wrapper provides easy access to the Langevin integrator, and broad opportunity to experiment, adapt, and extend it further. 

langevin-2025.11.10a2/TODO.md

@@ -0,0 +1,7 @@
+[ ] Implement video output properly
+
+[ ] Complete bcs for 1d
+
+[ ] Develop DP scaling-validation tests
+
+[ ] Implement PCP

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/cplusplus/dp/wrapper_dplvn.cpp

@@ -25,7 +25,7 @@
  */
 PYBIND11_MODULE(dplvn, module)
 {
-    module.attr("__version__") = "2025.11.8a2";
+    module.attr("__version__") = "2025.11.10a2";
     module.doc() = 
         "Operator-splitting method of integrating DP-type Langevin equations"; 
   

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/docs/doxygen/wrapper__dplvn_8cpp_source.html

@@ -101,7 +101,7 @@ $(function(){initNavTree('wrapper__dplvn_8cpp_source.html','',''); });
 <div class="foldopen" id="foldopen00026" data-start="{" data-end="}">
 <div class="line"><a id="l00026" name="l00026"></a><span class="lineno"><a class="line" href="wrapper__dplvn_8cpp.html#a260583c98d3e428eb1fc79bf477e2528">   26</a></span><a class="code hl_function" href="wrapper__dplvn_8cpp.html#a260583c98d3e428eb1fc79bf477e2528">PYBIND11_MODULE</a>(dplvn, module)</div>
 <div class="line"><a id="l00027" name="l00027"></a><span class="lineno">   27</span>{</div>
-<div class="line"><a id="l00028" name="l00028"></a><span class="lineno">   28</span>    <span class="keyword">module</span>.attr(&quot;__version__&quot;) = &quot;2025.11.8a2&quot;;</div>
+<div class="line"><a id="l00028" name="l00028"></a><span class="lineno">   28</span>    <span class="keyword">module</span>.attr(&quot;__version__&quot;) = &quot;2025.11.10a2&quot;;</div>
 <div class="line"><a id="l00029" name="l00029"></a><span class="lineno">   29</span>    <span class="keyword">module</span>.doc() = </div>
 <div class="line"><a id="l00030" name="l00030"></a><span class="lineno">   30</span>        &quot;Operator-splitting method of integrating DP-type Langevin equations&quot;; </div>
 <div class="line"><a id="l00031" name="l00031"></a><span class="lineno">   31</span>  </div>

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/docs/how-to-run.md

@@ -1,5 +1,8 @@
 # How to run
 
+![](https://raw.githubusercontent.com/cstarkjp/Langevin/main/images/ρ_a1p18950_b1_D0p04_η1_x100_y50_Δx1_Δt0p1_rs1.gif
+ "Density field evolution over time")
+
 1. Navigate to [`Langevin/simulation/dp/`](https://github.com/cstarkjp/Langevin/tree/main/simulation/dp). There you'll find Jupyter notebooks and Python scripts to run DP-type Langevin simulations.
 
 2. Don't run the noteboooks _in-situ_: their output will be written to [`Langevin/experiments/`](https://github.com/cstarkjp/Langevin/tree/main/experiments) which will generate `git` conflicts. 

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/docs/index.md

@@ -1,6 +1,6 @@
 # [**langevin**](https://pypi.org/project/langevin/)
 
-###  _Tools for integrating APT-type Langevin equations._
+###  _Tools for integrating APT-type Langevin equations_
 
 The  [`langevin` package ](https://pypi.org/project/langevin/) provides software tools to help integrate a time-dependent density field described by Langevin equations of absorbing phase transition (APT) type.
 [Directed percolation (DP)](references.md) is the _type-example_ of such an absorbing phase transition — its Langevin equation is:
@@ -19,6 +19,6 @@ where $\rho(\mathbf{x},t)$ is a fluctuating meso-scale field  evolving nonlinear
 
 ![Plot of grid-averaged density $\overline{\rho}(t)$ versus time, for an ensemble of simulations with $a$ taking values ranging symmetrically about criticality $a_c \approx 1.8857$ by up to $\Delta{a}=\pm 0.01$:](images/ρ_t_loglog_reduced.png)
 
-The `langevin` Langevin integrator employs the operator-splitting method originated largely by [Dornic et al (2005)](references.md). The software tools are implemented as a [`pip`-installable Python package](https://pypi.org/project/langevin/) with a C++ core, a set of [Jupyter notebooks](https://github.com/cstarkjp/Langevin/tree/main/simulation/dp), and related [Python scripts](https://github.com/cstarkjp/Langevin/tree/main/python).
+The `langevin` integrator employs the operator-splitting method originated largely by [Dornic et al (2005)](references.md). The software tools are implemented as a [`pip`-installable Python package](https://pypi.org/project/langevin/) with a C++ core, a set of [Jupyter notebooks](https://github.com/cstarkjp/Langevin/tree/main/simulation/dp), and related [Python scripts](https://github.com/cstarkjp/Langevin/tree/main/python).
 
 

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/docs/navigation.md

@@ -10,4 +10,5 @@ Look under under "Simulation tools" in the sidebar for more complete examples an
 
 The key driver of a simulation is the [`Info.json`](info-reference.md) file: care must be taken to match the "job name" implied by this file (a string constructed from the model coefficients and parameters specified by it) with its parent folder name, such that output files are placed correctly.
 
-Refer to the links under "Python modules" to see documentation of the `langevin` Python package. The underlying `C++` core is documented under "C++ source" using `Doxygen`.
+Refer to the links under "Python modules" to see documentation of the 
+[`langevin` Python package](https://pypi.org/project/langevin/). The underlying `C++` core is documented under ["C++ source"](cplusplus.md) using `Doxygen`.

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/experiments/a1p18855_b1_D0p04_/316/2671_x31_y31_/316/224x1_/316/224t0p1/rs1/Outfo.json

@@ -75,7 +75,7 @@
             "a1p18855_b1_D0p04_η1_x31_y31_Δx1_Δt0p1"
         ],
         "name": "a1p18855_b1_D0p04_η1_x31_y31_Δx1_Δt0p1",
-        "dplvn_version": "2025.11.8a2",
-        "date_time": "2025-11-08 09:46:26"
+        "dplvn_version": "2025.11.10a2",
+        "date_time": "2025-11-10 18:24:11"
     }
 }

langevin-2025.11.10a2/experiments/a1p18950_b1_D0p04_/316/2671_x100_y50_/316/224x1_/316/224t0p1/Info.json

@@ -0,0 +1,41 @@
+{
+    "Parameters": {
+        "linear": 1.18950,
+        "quadratic": 1,
+        "diffusion": 0.04,
+        "noise": 1,
+        "dx": 1,
+        "dt": 0.1,
+        "t_final": 300.0,
+        "random_seed": 1,
+        "grid_dimension": "D2",
+        "grid_size": [100, 50],
+        "grid_topologies": ["BOUNDED", "BOUNDED"],
+        "boundary_conditions": ["FLOATING","FLOATING","FLOATING","FLOATING"],
+        "initial_condition": "RANDOM_UNIFORM",
+        "ic_values": [0, 1.8],
+        "integration_method": "RUNGE_KUTTA"
+    },
+    "Analysis": {
+        "a_c": 1.18857
+    },
+    "Misc": {
+        "n_sims": 15,
+        "Δa_range": 0.01,
+        "n_round_Δt_summation": 5,
+        "n_segments": 300,
+        "n_digits": 6,
+        "do_export_combo_graphs": false,
+        "do_export_combo_data": false,
+        "do_export_graphs": true,
+        "do_export_data": true,
+        "do_export_images": true,
+        "do_make_video": true,
+        "ylimits_log": [5e-4, 2e0],
+        "xlimits_log": [5e-1, 1.2e4],
+        "ylimits_rescaled": [3e-2, 1e1],
+        "xlimits_rescaled": [5e-5, 2e3],
+        "video_frame_rate": 50,
+        "video_format": "mp4"
+    }
+}

{langevin-2025.11.8a2/experiments/a1p18855_b1_D0p04_/316/2671_x31_y31_ → langevin-2025.11.10a2/experiments/a1p18950_b1_D0p04_/316/2671_x100_y50_}/316/224x1_/316/224t0p1/Outfo.json

@@ -1,17 +1,17 @@
 {
     "Parameters": {
-        "linear": 1.18855,
+        "linear": 1.1895,
         "quadratic": 1,
         "diffusion": 0.04,
         "noise": 1,
         "dx": 1,
         "dt": 0.1,
-        "t_final": 10000.0,
+        "t_final": 300.0,
         "random_seed": 1,
         "grid_dimension": "D2",
         "grid_size": [
-            31,
-            31
+            100,
+            50
         ],
         "grid_topologies": [
             "BOUNDED",
@@ -26,7 +26,7 @@
         "initial_condition": "RANDOM_UNIFORM",
         "ic_values": [
             0,
-            10
+            1.8
         ],
         "integration_method": "RUNGE_KUTTA"
     },
@@ -36,20 +36,20 @@
         "dp_ν_pp": 0.734,
         "dp_ν_ll": 1.295,
         "dp_δ": 0.451,
-        "n_epochs": 100001
+        "n_epochs": 3001
     },
     "Misc": {
         "n_sims": 15,
         "Δa_range": 0.01,
         "n_round_Δt_summation": 5,
-        "n_segments": 500,
+        "n_segments": 300,
         "n_digits": 6,
         "do_export_combo_graphs": false,
         "do_export_combo_data": false,
         "do_export_graphs": true,
         "do_export_data": true,
-        "do_export_images": false,
-        "do_make_video": false,
+        "do_export_images": true,
+        "do_make_video": true,
         "ylimits_log": [
             0.0005,
             2.0
@@ -68,8 +68,8 @@
         ],
         "video_frame_rate": 50,
         "video_format": "mp4",
-        "dplvn_version": "2025.11.8a2",
-        "date_time": "2025-11-08 09:46:22",
-        "computation_time": "0:00:04"
+        "dplvn_version": "2025.11.10a2",
+        "date_time": "2025-11-10 18:22:31",
+        "computation_time": "0:00:01"
     }
 }

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/experiments/ac1p18857/b1_D0p04_/316/2671_x31_y31_/316/224x1_/316/224t0p1/Outfo.json

@@ -151,7 +151,7 @@
         ],
         "n_cores": 16,
         "computation_time": "0:00:05",
-        "dplvn_version": "2025.11.8a2",
+        "dplvn_version": "2025.11.10a2",
         "date_time": "2025-11-08 09:46:28"
     }
 }

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/meson.build

@@ -5,7 +5,7 @@
 project(
     'Langevin',
     'cpp',
-    version: '2025.11.8a2',
+    version: '2025.11.10a2',
     default_options: ['cpp_std=c++11'],
 )
 add_project_arguments(

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/mkdocs.yml

@@ -37,13 +37,11 @@ nav:
           - dp/batch_ensemble.py: dp-batch_ensemble-py-reference.md
       - JSON info files: 
           - Example: info-reference.md
-    - Python modules:
+    - Python modules: 
         - base/file.py: file-reference.md 
-        # - initialize.py: initialize-reference.md 
         - base/serialize.py: serialize-reference.md 
         - base/utils.py: utils-reference.md 
         - base/viz.py: viz-reference.md
-        # - dp:
         - dp/ensemble.py: ensemble-reference.md
         - dp/simulation.py: simulation-reference.md
         - dp/vizdp.py: vizdp-reference.md

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/pyproject.toml

@@ -9,7 +9,7 @@ requires = ["meson-python", "wheel", "pybind11", "ninja", "numpy"]
 
 [project]
 name = 'langevin'
-version = '2025.11.8a2'
+version = '2025.11.10a2'
 # dynamic = ["version"]
 description = 'Tools to integrate Langevin equations of absorbing phase transition type'
 readme = 'README.md'

langevin-2025.11.10a2/python/__init__.py

@@ -0,0 +1,3 @@
+__version__ = "2025.11.10a2"
+
+__all__ = []

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/python/base/__init__.py

@@ -1,4 +1,4 @@
-__version__ = "2025.11.8a2"
+__version__ = "2025.11.10a2"
 
 __all__ = [
     "file",

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/python/dp/simulation.py

@@ -10,8 +10,14 @@ from datetime import datetime, timedelta
 from numpy.typing import NDArray
 import numpy as np
 from numpy.lib.npyio import NpzFile
+try:
+    import ffmpeg
+except:
+    # Quietly fail
+    pass
 import sys, os
-sys.path.insert(0, os.path.join(os.path.pardir, "Packages"))
+from os.path import join, pardir
+sys.path.insert(0, join(pardir, "Packages"))
 from langevin.base.file import (
     create_directories, export_info, export_plots,
 )
@@ -199,12 +205,12 @@ class Simulation:
             do_verbose: report how the exporting is going
         """
         if self.do_verbose | do_verbose:
-            print(f"Outfo/graph/data path:  {self.misc["path"]}")
+            print(f"Outfo/graphs/videos/data path:  {self.misc["path"]}")
         seed_dir_name: str = f"rs{self.parameters["random_seed"]}"
     
         outfo_path: str = \
             create_directories(
-                (os.path.pardir, *self.misc["path"]), seed_dir_name,
+                (pardir, *self.misc["path"]), seed_dir_name,
             )
         outfo: dict = {
             "Parameters" : self.parameters,
@@ -214,30 +220,56 @@ class Simulation:
         if not do_dummy:
             _ = export_info(outfo_path, "Outfo", outfo, module,)
 
-        if self.misc["do_export_graphs"]:
-            graphs_path: str = \
-                create_directories(
-                    (os.path.pardir,  *self.misc["path"], seed_dir_name,), ".",
-                )
-            if not do_dummy:
-                _ = export_plots(
-                        self.graphs.fdict, 
-                        graphs_path,
-                        do_verbose=self.do_verbose,
-                    )
-
         if self.misc["do_export_data"]:
             data_path: str = \
                 create_directories(
-                    (os.path.pardir, *self.misc["path"], seed_dir_name,), ".", 
+                    (pardir, *self.misc["path"], seed_dir_name,), ".", 
                 )
             if not do_dummy:
                 np.savez_compressed(
-                    os.path.join(data_path, "ρ_t",), 
+                    join(data_path, "ρ_t",), 
                     t_epochs=self.t_epochs,
                     mean_densities=self.mean_densities,
                 )
                 data_npz: NpzFile = np.load(
-                    os.path.join(data_path, "ρ_t"+".npz",), 
+                    join(data_path, "ρ_t"+".npz",), 
                 )
                 data_npz["t_epochs"][-10:], data_npz["mean_densities"][-10:]
+
+        if self.misc["do_export_graphs"]:
+            graphs_path: str = \
+                create_directories(
+                    (pardir,  *self.misc["path"], seed_dir_name,), ".",
+                )
+            if not do_dummy:
+                _ = export_plots(
+                        self.graphs.fdict, 
+                        graphs_path,
+                        do_verbose=self.do_verbose,
+                    )
+
+        if self.misc["do_make_video"]:
+            print("Video export not working yet")
+            # videos_path: str = create_directories(
+            #     (pardir, pardir, "experiments", sim_dir_name,), "Videos",
+            # )
+
+            # video_frame_rate: int = misc["video_frame_rate"]
+            # video_format: str = misc["video_format"]
+            # video_images_wildcard: str = "ρ_t"+"?"*n_digits+".png"
+            # input = ffmpeg.input( 
+            #     join(images_path, video_images_wildcard), 
+            #     pattern_type="glob", 
+            #     framerate=video_frame_rate, 
+            #     pix_fmt="yuv420p",
+            # )
+            # output = ffmpeg.output(
+            #     input.video,
+            #     join(
+            #         videos_path, 
+            #         f"ρ_{sim_dir_name}_rs{parameters["random_seed"]}.{video_format}"
+            #     ),
+            #     vf="crop=floor(iw/2)*2:floor(ih/2)*2",
+            #     vcodec="libx264",
+            #     format=video_format,
+            # )

{langevin-2025.11.8a2 → langevin-2025.11.10a2}/simulation/dp/Simulation.ipynb

@@ -35,7 +35,7 @@
     {
      "data": {
       "text/plain": [
-       "'2025.11.8a2'"
+       "'2025.11.10a1'"
       ]
      },
      "execution_count": 1,
@@ -181,7 +181,7 @@
      "name": "stderr",
      "output_type": "stream",
      "text": [
-      "100%|\u001b[32m██████████\u001b[0m| 501/501 [00:04<00:00, 107.39it/s]\n"
+      "100%|\u001b[32m██████████\u001b[0m| 501/501 [00:03<00:00, 131.24it/s]\n"
      ]
     }
    ],
@@ -245,19 +245,19 @@
      "name": "stdout",
      "output_type": "stream",
      "text": [
-      "Outfo/graph/data path:  ['..', 'experiments', 'a1p18855_b1_D0p04_η1_x31_y31_Δx1_Δt0p1']\n"
+      "Outfo/graphs/videos/data path:  ['..', 'experiments', 'a1p18855_b1_D0p04_η1_x31_y31_Δx1_Δt0p1']\n"
      ]
     },
     {
      "name": "stderr",
      "output_type": "stream",
      "text": [
-      "100%|\u001b[32m██████████\u001b[0m| 2/2 [00:00<00:00,  4.97it/s]\n"
+      "100%|\u001b[32m██████████\u001b[0m| 2/2 [00:00<00:00,  8.45it/s]\n"
      ]
     }
    ],
    "source": [
-    "sim.save(dplvn, do_dummy=False,)"
+    "sim.save(dplvn, do_verbose=True,)"
    ]
   },
   {