From 8934d034756bddf474928c8857a280439c4165b1 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Thu, 26 Feb 2026 10:55:17 -0500
Subject: [PATCH 01/23] warnings should only be printed on MPI rank 0

---
 src/MC/fix_gcmc.cpp | 24 ++++++++++--------------
 1 file changed, 10 insertions(+), 14 deletions(-)

diff --git a/src/MC/fix_gcmc.cpp b/src/MC/fix_gcmc.cpp
index e697663a7a3..dc59e9d1e9c 100644
--- a/src/MC/fix_gcmc.cpp
+++ b/src/MC/fix_gcmc.cpp
@@ -280,7 +280,7 @@ void FixGCMC::options(int narg, char **arg)
       if (imol == -1)
         error->all(FLERR, iarg + 1, "Molecule template ID {} for fix gcmc does not exist",
                    arg[iarg + 1]);
-      if (atom->molecules[imol]->nset > 1 && comm->me == 0)
+      if ((atom->molecules[imol]->nset > 1) && (comm->me == 0))
         error->warning(FLERR, "Molecule template for fix gcmc has multiple molecules");
       exchmode = EXCHMOL;
       onemols = atom->molecules;
@@ -712,8 +712,7 @@ void FixGCMC::init()
   // warning if group id is "all"
 
   if ((comm->me == 0) && (groupbit & 1))
-    error->warning(FLERR, "Fix gcmc is being applied "
-                   "to the default group all");
+    error->warning(FLERR, "Fix gcmc is being applied ""to the default group all");
 
   // construct group bitmask for all new atoms
   // aggregated over all group keywords
@@ -790,9 +789,8 @@ void FixGCMC::pre_exchange()
 
   if (full_flag) {
     energy_stored = energy_full();
-    if (overlap_flag && energy_stored > MAXENERGYTEST)
-        error->warning(FLERR,"Energy of old configuration in "
-                       "fix gcmc is > MAXENERGYTEST.");
+    if (overlap_flag && (energy_stored > MAXENERGYTEST) && (comm->me == 0))
+      error->warning(FLERR,"Energy of old configuration in fix gcmc is > MAXENERGYTEST.");
 
     for (int i = 0; i < ncycles; i++) {
       int ixm = static_cast<int>(random_equal->uniform()*ncycles) + 1;
@@ -861,8 +859,8 @@ void FixGCMC::attempt_atomic_translation()
   if (i >= 0) {
     double **x = atom->x;
     double energy_before = energy(i,ngcmc_type,-1,x[i]);
-    if (overlap_flag && energy_before > MAXENERGYTEST)
-        error->warning(FLERR,"Energy of old configuration in fix gcmc is > MAXENERGYTEST.");
+    if (overlap_flag && (energy_before > MAXENERGYTEST) && (comm->me == 0))
+      error->warning(FLERR,"Energy of old configuration in fix gcmc is > MAXENERGYTEST.");
     double rsq = 1.1;
     double rx,ry,rz;
     rx = ry = rz = 0.0;
@@ -1085,9 +1083,8 @@ void FixGCMC::attempt_molecule_translation()
   if (translation_molecule == -1) return;
 
   double energy_before_sum = molecule_energy(translation_molecule);
-  if (overlap_flag && energy_before_sum > MAXENERGYTEST)
-    error->warning(FLERR,"Energy of old configuration in "
-                   "fix gcmc is > MAXENERGYTEST.");
+  if (overlap_flag && (energy_before_sum > MAXENERGYTEST) && (comm->me == 0))
+    error->warning(FLERR,"Energy of old configuration in fix gcmc is > MAXENERGYTEST.");
 
   double **x = atom->x;
   double rx,ry,rz;
@@ -1184,9 +1181,8 @@ void FixGCMC::attempt_molecule_rotation()
   if (rotation_molecule == -1) return;
 
   double energy_before_sum = molecule_energy(rotation_molecule);
-  if (overlap_flag && energy_before_sum > MAXENERGYTEST)
-    error->warning(FLERR,"Energy of old configuration in "
-                   "fix gcmc is > MAXENERGYTEST.");
+  if (overlap_flag && (energy_before_sum > MAXENERGYTEST) && (comm->me == 0))
+    error->warning(FLERR,"Energy of old configuration in fix gcmc is > MAXENERGYTEST.");
 
   int *mask = atom->mask;
   int nmolcoords = 0;

From 3911b91b5511b8270437b6c5c3a1223837439233 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Thu, 26 Feb 2026 11:27:26 -0500
Subject: [PATCH 02/23] avoid getting stuck in region restricted MC moves and
 improve errors

---
 src/MC/fix_gcmc.cpp | 37 +++++++++++++++++++++++++++----------
 1 file changed, 27 insertions(+), 10 deletions(-)

diff --git a/src/MC/fix_gcmc.cpp b/src/MC/fix_gcmc.cpp
index dc59e9d1e9c..a66a461ca1b 100644
--- a/src/MC/fix_gcmc.cpp
+++ b/src/MC/fix_gcmc.cpp
@@ -875,6 +875,7 @@ void FixGCMC::attempt_atomic_translation()
     coord[1] = x[i][1] + displace*ry;
     coord[2] = x[i][2] + displace*rz;
     if (region) {
+      int region_attempt = 0;
       while (region->match(coord[0],coord[1],coord[2]) == 0) {
         rsq = 1.1;
         while (rsq > 1.0) {
@@ -886,10 +887,14 @@ void FixGCMC::attempt_atomic_translation()
         coord[0] = x[i][0] + displace*rx;
         coord[1] = x[i][1] + displace*ry;
         coord[2] = x[i][2] + displace*rz;
+        ++region_attempt;
+        if (region_attempt >= max_region_attempts) break;
       }
+      if (!region->match(coord[0],coord[1],coord[2]))
+        error->one(FLERR,"Fix gcmc translation put atom outside region");
     }
     if (!domain->inside_nonperiodic(coord))
-      error->one(FLERR,"Fix gcmc put atom outside box");
+      error->one(FLERR,"Fix gcmc translation put atom outside box");
 
     double energy_after = energy(i,ngcmc_type,-1,coord);
 
@@ -1115,6 +1120,7 @@ void FixGCMC::attempt_molecule_translation()
     coord[0] = com[0] + displace*rx;
     coord[1] = com[1] + displace*ry;
     coord[2] = com[2] + displace*rz;
+    int region_attempt = 0;
     while (region->match(coord[0],coord[1],coord[2]) == 0) {
       rsq = 1.1;
       while (rsq > 1.0) {
@@ -1126,10 +1132,14 @@ void FixGCMC::attempt_molecule_translation()
       coord[0] = com[0] + displace*rx;
       coord[1] = com[1] + displace*ry;
       coord[2] = com[2] + displace*rz;
+      ++region_attempt;
+      if (region_attempt >= max_region_attempts) break;
     }
     com_displace[0] = displace*rx;
     com_displace[1] = displace*ry;
     com_displace[2] = displace*rz;
+    if (!region->match(coord[0],coord[1],coord[2]))
+      error->one(FLERR,"Fix gcmc translation put molecule COM outside region");
   }
 
   double energy_after = 0.0;
@@ -1139,7 +1149,7 @@ void FixGCMC::attempt_molecule_translation()
       coord[1] = x[i][1] + com_displace[1];
       coord[2] = x[i][2] + com_displace[2];
       if (!domain->inside_nonperiodic(coord))
-        error->one(FLERR,"Fix gcmc put atom outside box");
+        error->one(FLERR,"Fix gcmc translation put molecule atom outside box");
       energy_after += energy(i,atom->type[i],translation_molecule,coord);
     }
   }
@@ -1330,8 +1340,7 @@ void FixGCMC::attempt_molecule_insertion()
       (region_yhi-region_ylo);
     com_coord[2] = region_zlo + random_equal->uniform() *
       (region_zhi-region_zlo);
-    while (region->match(com_coord[0],com_coord[1],
-                                           com_coord[2]) == 0) {
+    while (region->match(com_coord[0],com_coord[1],com_coord[2]) == 0) {
       com_coord[0] = region_xlo + random_equal->uniform() *
         (region_xhi-region_xlo);
       com_coord[1] = region_ylo + random_equal->uniform() *
@@ -1542,6 +1551,7 @@ void FixGCMC::attempt_atomic_translation_full()
     coord[1] = x[i][1] + displace*ry;
     coord[2] = x[i][2] + displace*rz;
     if (region) {
+      int region_attempt = 0;
       while (region->match(coord[0],coord[1],coord[2]) == 0) {
         rsq = 1.1;
         while (rsq > 1.0) {
@@ -1553,10 +1563,14 @@ void FixGCMC::attempt_atomic_translation_full()
         coord[0] = x[i][0] + displace*rx;
         coord[1] = x[i][1] + displace*ry;
         coord[2] = x[i][2] + displace*rz;
+        ++region_attempt;
+        if (region_attempt >= max_region_attempts) break;
       }
+      if (!region->match(coord[0],coord[1],coord[2]))
+        error->one(FLERR,"Fix gcmc translation put atom outside region");
     }
     if (!domain->inside_nonperiodic(coord))
-      error->one(FLERR,"Fix gcmc put atom outside box");
+      error->one(FLERR,"Fix gcmc translation put atom outside box");
     xtmp[0] = x[i][0];
     xtmp[1] = x[i][1];
     xtmp[2] = x[i][2];
@@ -1791,6 +1805,7 @@ void FixGCMC::attempt_molecule_translation_full()
         mask[i] &= molecule_group_inversebit;
       }
     }
+    int region_attempt = 0;
     double com[3];
     com[0] = com[1] = com[2] = 0.0;
     group->xcm(molecule_group,gas_mass,com);
@@ -1808,7 +1823,11 @@ void FixGCMC::attempt_molecule_translation_full()
       coord[0] = com[0] + displace*rx;
       coord[1] = com[1] + displace*ry;
       coord[2] = com[2] + displace*rz;
+      ++region_attempt;
+      if (region_attempt >= max_region_attempts) break;
     }
+    if (!region->match(coord[0],coord[1],coord[2]))
+      error->one(FLERR,"Fix gcmc translation put molecule COM outside region");
     com_displace[0] = displace*rx;
     com_displace[1] = displace*ry;
     com_displace[2] = displace*rz;
@@ -1820,7 +1839,7 @@ void FixGCMC::attempt_molecule_translation_full()
       x[i][1] += com_displace[1];
       x[i][2] += com_displace[2];
       if (!domain->inside_nonperiodic(x[i]))
-        error->one(FLERR,"Fix gcmc put atom outside box");
+        error->one(FLERR,"Fix gcmc put molecule atom outside box");
     }
   }
 
@@ -2062,8 +2081,7 @@ void FixGCMC::attempt_molecule_insertion_full()
       (region_yhi-region_ylo);
     com_coord[2] = region_zlo + random_equal->uniform() *
       (region_zhi-region_zlo);
-    while (region->match(com_coord[0],com_coord[1],
-                                           com_coord[2]) == 0) {
+    while (region->match(com_coord[0],com_coord[1],com_coord[2]) == 0) {
       com_coord[0] = region_xlo + random_equal->uniform() *
         (region_xhi-region_xlo);
       com_coord[1] = region_ylo + random_equal->uniform() *
@@ -2491,8 +2509,7 @@ void FixGCMC::update_gas_atoms_list()
 
       for (int i = 0; i < nlocal; i++) {
         if (mask[i] & groupbit) {
-          if (region->match(comx[molecule[i]],
-             comy[molecule[i]],comz[molecule[i]]) == 1) {
+          if (region->match(comx[molecule[i]],comy[molecule[i]],comz[molecule[i]]) == 1) {
             local_gas_list[ngas_local] = i;
             ngas_local++;
           }

From fe09513dd8d0ee6b3472b3f859fb957e709fdf5f Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Thu, 26 Feb 2026 12:05:59 -0500
Subject: [PATCH 03/23] use PPM format images instead of PNG for better
 compatibility

---
 examples/GRAPHICS/in.breakable          |  2 +-
 examples/GRAPHICS/in.cubes-and-pyramids |  2 +-
 examples/GRAPHICS/in.obstacle-lines     | 84 ++++++++++++-------------
 examples/GRAPHICS/in.peptide-hbonds     |  2 +-
 examples/GRAPHICS/in.water-arrows       |  2 +-
 5 files changed, 46 insertions(+), 46 deletions(-)

diff --git a/examples/GRAPHICS/in.breakable b/examples/GRAPHICS/in.breakable
index d483a000313..f92eadee201 100644
--- a/examples/GRAPHICS/in.breakable
+++ b/examples/GRAPHICS/in.breakable
@@ -86,7 +86,7 @@ fix labels all graphics/labels  500 text "LAMMPS Tutorial 2 - Stretching CNT usi
                                     text "Atoms with bonds  1: $(c_count[1])  2: $(c_count[2])  3: $(c_count[3])  4: $(c_count[4])" &
                                     32.5 7.5 5.0 size 24 fontcolor white backcolor silver
 
-dump viz all image 500 breakable.*.png c_nbond type size 1600 300 zoom 12.5 shiny 0.2 autobond 1.8 0.4 &
+dump viz all image 500 breakable.*.ppm c_nbond type size 1600 300 zoom 12.5 shiny 0.2 autobond 1.8 0.4 &
         adiam 0.8 box no 0.01 view 0 90 fsaa yes ssao yes 1325123 0.6 axes yes 0.5 0.1 &
         fix arrows const 0 0  fix labels const 0 0 fix vel const 0 0
 
diff --git a/examples/GRAPHICS/in.cubes-and-pyramids b/examples/GRAPHICS/in.cubes-and-pyramids
index bd924e73954..bed142c5151 100644
--- a/examples/GRAPHICS/in.cubes-and-pyramids
+++ b/examples/GRAPHICS/in.cubes-and-pyramids
@@ -45,7 +45,7 @@ thermo       1000
 variable progress equal step/v_steps
 fix pbar all graphics/objects ${out} progbar 3 4 y 24.0 12.0 -1.4 30.0 0.5 v_progress 10
 
-dump         2 all image ${out} cubes-and-pyramids.*.png type type size 600 600 shiny 0.4 &
+dump         2 all image ${out} cubes-and-pyramids.*.ppm type type size 600 600 shiny 0.4 &
                          box yes 0.025 axes no 0.5 0.1 zoom 1.5 view 80 0 fsaa yes ssao yes 315123 0.4 &
                          body type 0.2 3 fix pbar type 0 0
 #            set cylinder diameter--^  ^-- set viz style (1 only faces, 2 cylinders, 3 both)
diff --git a/examples/GRAPHICS/in.obstacle-lines b/examples/GRAPHICS/in.obstacle-lines
index c9d06bfa788..6c9bd5b54cc 100644
--- a/examples/GRAPHICS/in.obstacle-lines
+++ b/examples/GRAPHICS/in.obstacle-lines
@@ -1,70 +1,70 @@
 # 2d LJ obstacle flow
 
-dimension	2
-boundary	p s p
+dimension       2
+boundary        p s p
 
-atom_style	atomic
-neighbor	0.3 bin
-neigh_modify	delay 5
+atom_style      atomic
+neighbor        0.3 bin
+neigh_modify    delay 5
 
 # create geometry
 
-lattice		hex 0.7
-region		box block 0 40 0 10 -0.25 0.25
-create_box	4 box
-create_atoms	1 box
+lattice         hex 0.7
+region          box block 0 40 0 10 -0.25 0.25
+create_box      4 box
+create_atoms    1 box
 
-mass		* 1.0
+mass            * 1.0
 
 # LJ potentials
 
-pair_style	lj/cut 1.12246
-pair_coeff	* * 1.0 1.0 1.12246
+pair_style      lj/cut 1.12246
+pair_coeff      * * 1.0 1.0 1.12246
 
 # define groups
 
-region	        1 block INF INF INF 1.25 INF INF
-group		lower region 1
-region		2 block INF INF 8.75 INF INF INF
-group		upper region 2
-group		boundary union lower upper
-group		flow subtract all boundary
+region          1 block INF INF INF 1.25 INF INF
+group           lower region 1
+region          2 block INF INF 8.75 INF INF INF
+group           upper region 2
+group           boundary union lower upper
+group           flow subtract all boundary
 
-set		group lower type 2
-set		group upper type 3
+set             group lower type 2
+set             group upper type 3
 
 # initial velocities
 
-compute	        mobile flow temp
-velocity	flow create 1.0 482748 temp mobile
-fix		1 all nve
-fix		2 flow temp/rescale 200 1.0 1.0 0.02 1.0
-fix_modify	2 temp mobile
+compute         mobile flow temp
+velocity        flow create 1.0 482748 temp mobile
+fix             1 all nve
+fix             2 flow temp/rescale 200 1.0 1.0 0.02 1.0
+fix_modify      2 temp mobile
 
 # Poiselle flow
 
-velocity	boundary set 0.0 0.0 0.0
-fix		3 lower setforce 0.0 0.0 0.0
-fix		4 upper setforce 0.0 NULL 0.0
-fix		5 upper aveforce 0.0 -0.5 0.0
-fix		6 flow addforce 1.0 0.0 0.0
+velocity        boundary set 0.0 0.0 0.0
+fix             3 lower setforce 0.0 0.0 0.0
+fix             4 upper setforce 0.0 NULL 0.0
+fix             5 upper aveforce 0.0 -0.5 0.0
+fix             6 flow addforce 1.0 0.0 0.0
 
 # 2 obstacles
 
-region		void1 sphere 10 4 0 3
-delete_atoms	region void1
-region		void2 sphere 20 7 0 3
-delete_atoms	region void2
+region          void1 sphere 10 4 0 3
+delete_atoms    region void1
+region          void2 sphere 20 7 0 3
+delete_atoms    region void2
 
-fix		7 flow indent 100 sphere 10 4 0 4
-fix		8 flow indent 100 sphere 20 7 0 4
-fix		9 all enforce2d
+fix             7 flow indent 100 sphere 10 4 0 4
+fix             8 flow indent 100 sphere 20 7 0 4
+fix             9 all enforce2d
 
 # Run
 
-timestep	0.003
-thermo		1000
-thermo_modify	temp mobile
+timestep        0.003
+thermo          1000
+thermo_modify   temp mobile
 
 # select atoms to trace and assign atom type 4 so they have a different color (yellow)                
 region   trace block 2 3 1.25 8.75 INF INF
@@ -74,8 +74,8 @@ set      group trace type 4
 # create trajectory lines by averaging positions over 10 MD steps each 10 steps apart
 fix      lines trace graphics/lines 10 10 100 25
 # visualize the trace in silver with a diameter of 0.6 which is half of that of the atoms
-dump     viz   all image 100 obstacle-*.png type type size 800 600 zoom 2.5 center s 0.50 0.7 0 &
+dump     viz   all image 100 obstacle-*.ppm type type size 800 600 zoom 2.5 center s 0.50 0.7 0 &
                                            shiny 0.5 fsaa yes box no 0.025 fix lines const 1 0.6
 dump_modify  viz  pad 6 backcolor darkgray adiam * 1.2 fcolor lines silver
 
-run		30000
+run             30000
diff --git a/examples/GRAPHICS/in.peptide-hbonds b/examples/GRAPHICS/in.peptide-hbonds
index 228f54e671c..df1341fb0bc 100644
--- a/examples/GRAPHICS/in.peptide-hbonds
+++ b/examples/GRAPHICS/in.peptide-hbonds
@@ -53,7 +53,7 @@ fix obj   all graphics/objects 100 arrow 5 80.0 61.0 39 80.0 67.0 39 0.3 0.2 &
                                    arrow 6 80.0 61.0 37.2 80.0 67.0 37.2 0.3 0.2
 
 # combine the graphics into visualization
-dump viz viz image 100 hbonds-*.png element element size 600 600 zoom 2.1 view 80 0 bond atom 0.3 &
+dump viz viz image 100 hbonds-*.ppm element element size 600 600 zoom 2.1 view 80 0 bond atom 0.3 &
                        fsaa yes ssao yes 12384 0.6 shiny 0.1 box no 0.1 center s 0.5 0.52 0.52 &
                        compute hb1 const -0.4 0.3 compute hb2 const -0.4 0.3 fix label const 1 0 fix obj type 0.0 0.0
 dump_modify viz pad 5 boxcolor white backcolor darkgray backcolor2 silver &
diff --git a/examples/GRAPHICS/in.water-arrows b/examples/GRAPHICS/in.water-arrows
index 8b9acfb99b4..370419004bf 100644
--- a/examples/GRAPHICS/in.water-arrows
+++ b/examples/GRAPHICS/in.water-arrows
@@ -58,7 +58,7 @@ variable dip2y equal v_scale*c_dip[2]
 variable dip2z equal v_scale*c_dip[3]
 fix dipole all graphics/objects 1 arrow 1  v_dip1x v_dip1y v_dip1z v_dip2x v_dip2y v_dip2z 0.3 0.2
 
-dump viz all image 10 water-arrows-*.png element type size 600 600 zoom 1.3 shiny 0.2 bond atom 0.2 &
+dump viz all image 10 water-arrows-*.ppm element type size 600 600 zoom 1.3 shiny 0.2 bond atom 0.2 &
                                      view 70 20 box yes 0.025 fsaa yes ssao yes 315465 0.8 &
                                      fix dipole const 0 0  fix vec const 0 0  fix vel const 0 0 fix bin const 0 0
 #                                    all arrows from fixes use the "const" style so we set the color with dump_modify below

From f647ba8c280307877e82a876b6806fa3cc7f66c9 Mon Sep 17 00:00:00 2001
From: jtclemm <jtclemm@sandia.gov>
Date: Thu, 26 Feb 2026 14:28:23 -0700
Subject: [PATCH 04/23] Fix set structure mismatch for box

---
 src/EXTRA-FIX/fix_deform_pressure.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/EXTRA-FIX/fix_deform_pressure.cpp b/src/EXTRA-FIX/fix_deform_pressure.cpp
index 69fa1760d6d..f9169833fac 100644
--- a/src/EXTRA-FIX/fix_deform_pressure.cpp
+++ b/src/EXTRA-FIX/fix_deform_pressure.cpp
@@ -406,7 +406,7 @@ void FixDeformPressure::init()
 
   // reset cumulative counters to match resetting "start" variables in parent
 
-  set[6].cumulative_shift = 0.0;
+  set_box.cumulative_shift = 0.0;
   for (int i = 0; i < 7; i++) {
     set_extra[i].cumulative_vshift[0] = 0.0;
     set_extra[i].cumulative_vshift[1] = 0.0;

From 2ea33836d13138f41ca3417cea36877140263fe2 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Thu, 26 Feb 2026 17:52:15 -0500
Subject: [PATCH 05/23] suggest to check for overlapping atoms also for NaN and
 Inf values

---
 doc/src/Errors_details.rst | 23 ++++++++++++++++++-----
 1 file changed, 18 insertions(+), 5 deletions(-)

diff --git a/doc/src/Errors_details.rst b/doc/src/Errors_details.rst
index 1ab91838b07..2191b857ea8 100644
--- a/doc/src/Errors_details.rst
+++ b/doc/src/Errors_details.rst
@@ -142,11 +142,24 @@ contacts or bad geometries (for the given force styles in use) leading
 to forces that can no longer be represented as numbers.  Those will show
 as "NaN" or "Inf".  On most machines, the program will continue, but
 there is no way to recover from it and those NaN or Inf values will
-propagate.  So-called :doc:`"soft-core" potentials <pair_fep_soft>` or
-the :doc:`"soft" repulsive-only pair style <pair_soft>` are less prone
-for this behavior (depending on the settings in use) and can be used at
-the beginning of a simulation.  Also, single precision numbers can
-overflow much faster, so for the GPU, KOKKOS, or INTEL package it may be
+propagate.
+
+If the "NaN" or "Inf" appears in the first simulation step, the most
+common cause is overlapping atoms.  Note that when atoms are *very*
+close, this cannot be seen when visualizing the geometry, since the
+atoms are effectively sitting on top of each other.  A good test is to
+insert a command like :doc:`delete_atoms 0.1 all all <delete_atoms>` and
+then monitor the output to see how many atoms are deleted, if any.  A
+non-zero number would be an indication of overlapping atoms. Note that
+atoms can also overlap through periodic boundaries when the box
+dimensions are too small (e.g. determined by min/max position of atoms
+without padding).
+
+So-called :doc:`"soft-core" potentials <pair_fep_soft>` or the
+:doc:`"soft" repulsive-only pair style <pair_soft>` are less prone for
+this behavior (depending on the settings in use) and can be used at the
+beginning of a simulation.  Also, single precision numbers can overflow
+much faster, so for the GPU, KOKKOS, or INTEL package it may be
 beneficial to run with double precision initially before switching to
 mixed or single precision for faster execution when the system has
 relaxed.

From fa415a3ef6bfb83f8a5ec3e9049cefd3dca5eef5 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Fri, 27 Feb 2026 00:08:19 -0500
Subject: [PATCH 06/23] don't print last line for group bitmask lookup errors

---
 src/group.cpp | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/src/group.cpp b/src/group.cpp
index 5eb9bf9de95..ff515878cd2 100644
--- a/src/group.cpp
+++ b/src/group.cpp
@@ -658,7 +658,8 @@ int Group::get_bitmask_by_id(const std::string &file, int line, const std::strin
 {
   int igroup = find(name);
   if (igroup < 0)
-    error->all(file, line, "Group ID {} requested by {} does not exist", name, caller);
+    error->all(file, line, Error::NOLASTLINE, "Group ID {} requested by {} does not exist", name,
+               caller);
   return bitmask[igroup];
 }
 

From 79f401e49285fc8ac1972458d5df0d7292d78685 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Fri, 27 Feb 2026 00:08:30 -0500
Subject: [PATCH 07/23] re-apply clang-format

---
 src/group.cpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/group.cpp b/src/group.cpp
index ff515878cd2..4ee2dd79c3d 100644
--- a/src/group.cpp
+++ b/src/group.cpp
@@ -668,7 +668,7 @@ int Group::get_bitmask_by_id(const std::string &file, int line, const std::strin
 ------------------------------------------------------------------------- */
 
 int Group::get_inversemask_by_id(const std::string &file, int line, const std::string &name,
-                             const std::string &caller)
+                                 const std::string &caller)
 {
   int igroup = find(name);
   if (igroup < 0)
@@ -704,7 +704,7 @@ void Group::add_molecules(int /*igroup*/, int bit)
   memory->create(list, n, "group:list");
 
   n = 0;
-  for(const auto pos : hash) list[n++] = pos;
+  for (const auto pos : hash) list[n++] = pos;
 
   molbit = bit;
   comm->ring(n, sizeof(tagint), list, 1, molring, nullptr, (void *) this);

From d3c25783be34e8573801c29f5aed5bae21b0b584 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Fri, 27 Feb 2026 01:21:37 -0500
Subject: [PATCH 08/23] improve error message

---
 src/neighbor.cpp | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/src/neighbor.cpp b/src/neighbor.cpp
index 6496a8c39af..a532f09b1ab 100644
--- a/src/neighbor.cpp
+++ b/src/neighbor.cpp
@@ -856,7 +856,8 @@ int Neighbor::init_pair()
     if (requests[i]->cut) {
       if (comm_cutoff < (requests[i]->cutoff + (requests[i]->occasional ? 0.0 : skin)))
         error->all(FLERR,  Error::NOLASTLINE,
-                   "Custom neighbor list cutoff too large for communication cutoff");
+                   "Custom neighbor list cutoff {} is too large for communication cutoff {}",
+                   (requests[i]->cutoff + (requests[i]->occasional ? 0.0 : skin)), comm_cutoff);
     }
   }
 

From 06fc47e2d6f46b31c8270c13bb9cbc3da9cfab4b Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Fri, 27 Feb 2026 10:03:31 -0500
Subject: [PATCH 09/23] must clamp hi/lo values of regions to box boundaries in
 case they were set to INF

---
 src/GRAPHICS/dump_image.cpp | 93 ++++++++++++++++++++++++-------------
 1 file changed, 62 insertions(+), 31 deletions(-)

diff --git a/src/GRAPHICS/dump_image.cpp b/src/GRAPHICS/dump_image.cpp
index e943e11b478..f631ffa2028 100644
--- a/src/GRAPHICS/dump_image.cpp
+++ b/src/GRAPHICS/dump_image.cpp
@@ -1689,7 +1689,7 @@ void DumpImage::create_image()
                vec3{objarray[i][4], objarray[i][5], objarray[i][6]}, color, opacity);
       } else if (objvec[i] == Graphics::PIXMAP) {
         // get pointer to pixmap buffer and get background transparency color
-        const auto *pixmap = (const unsigned char *) ubuf(objarray[i][6]).i; // NOLINT
+        const auto *pixmap = (const unsigned char *) ubuf(objarray[i][6]).i;    // NOLINT
         double transcolor[3] = {objarray[i][7], objarray[i][8], objarray[i][9]};
         if (iobj.flag1 == 0.0)    // coordinates are in box coordinates
           image->draw_pixmap(&objarray[i][1], (int) objarray[i][4], (int) objarray[i][5], pixmap,
@@ -1820,11 +1820,17 @@ void DumpImage::create_image()
           // inconsistent style. should not happen.
           if (!myreg) continue;
 
-          corners = cornerdata{
-              vec3{myreg->xlo, myreg->ylo, myreg->zlo}, vec3{myreg->xlo, myreg->ylo, myreg->zhi},
-              vec3{myreg->xlo, myreg->yhi, myreg->zhi}, vec3{myreg->xlo, myreg->yhi, myreg->zlo},
-              vec3{myreg->xhi, myreg->ylo, myreg->zlo}, vec3{myreg->xhi, myreg->ylo, myreg->zhi},
-              vec3{myreg->xhi, myreg->yhi, myreg->zhi}, vec3{myreg->xhi, myreg->yhi, myreg->zlo}};
+          // clamp region boundaries to box boundaries
+          double xlo = MAX(myreg->xlo, domain->boxlo[0]);
+          double ylo = MAX(myreg->ylo, domain->boxlo[1]);
+          double zlo = MAX(myreg->zlo, domain->boxlo[2]);
+          double xhi = MIN(myreg->xhi, domain->boxhi[0]);
+          double yhi = MIN(myreg->yhi, domain->boxhi[1]);
+          double zhi = MIN(myreg->zhi, domain->boxhi[2]);
+
+          corners = cornerdata{vec3{xlo, ylo, zlo}, vec3{xlo, ylo, zhi}, vec3{xlo, yhi, zhi},
+                               vec3{xlo, yhi, zlo}, vec3{xhi, ylo, zlo}, vec3{xhi, ylo, zhi},
+                               vec3{xhi, yhi, zhi}, vec3{xhi, yhi, zlo}};
         }
 
         if (regstyle == "prism") {
@@ -1832,24 +1838,32 @@ void DumpImage::create_image()
           // inconsistent style. should not happen.
           if (!myreg) continue;
 
+          // clamp region boundaries to box boundaries
+          double xlo = MAX(myreg->xlo, domain->boxlo[0]);
+          double ylo = MAX(myreg->ylo, domain->boxlo[1]);
+          double zlo = MAX(myreg->zlo, domain->boxlo[2]);
+          double xhi = MIN(myreg->xhi, domain->boxhi[0]);
+          double yhi = MIN(myreg->yhi, domain->boxhi[1]);
+          double zhi = MIN(myreg->zhi, domain->boxhi[2]);
+
           corners = cornerdata{
-              vec3{myreg->xlo, myreg->ylo, myreg->zlo},
-              vec3{myreg->xlo + myreg->xz, myreg->ylo + myreg->yz, myreg->zhi},
-              vec3{myreg->xlo + myreg->xy + myreg->xz, myreg->yhi + myreg->yz, myreg->zhi},
-              vec3{myreg->xlo + myreg->xy, myreg->yhi, myreg->zlo},
-              vec3{myreg->xhi, myreg->ylo, myreg->zlo},
-              vec3{myreg->xhi + myreg->xz, myreg->ylo + myreg->yz, myreg->zhi},
-              vec3{myreg->xhi + myreg->xy + myreg->xz, myreg->yhi + myreg->yz, myreg->zhi},
-              vec3{myreg->xhi + myreg->xy, myreg->yhi, myreg->zlo}};
+              vec3{xlo, ylo, zlo},
+              vec3{xlo + myreg->xz, ylo + myreg->yz, zhi},
+              vec3{xlo + myreg->xy + myreg->xz, yhi + myreg->yz, zhi},
+              vec3{xlo + myreg->xy, yhi, zlo},
+              vec3{xhi, ylo, zlo},
+              vec3{xhi + myreg->xz, ylo + myreg->yz, zhi},
+              vec3{xhi + myreg->xy + myreg->xz, yhi + myreg->yz, zhi},
+              vec3{xhi + myreg->xy, yhi, zlo}};
         }
 
         for (int i = 0; i < 8; ++i)
           reg.ptr->forward_transform(corners[i][0], corners[i][1], corners[i][2]);
 
-#define DRAW_CYLINDER(i, j) \
-  image->draw_cylinder(corners[i].data(), corners[j].data(), reg.color, reg.diameter, 3, 1.0)
-#define DRAW_TRIANGLE(i, j, k) \
-  image->draw_triangle(corners[i].data(), corners[j].data(), corners[k].data(), reg.color, opacity)
+#define DRAW_CYLINDER(j, k) \
+  image->draw_cylinder(corners[j].data(), corners[k].data(), reg.color, reg.diameter, 3, 1.0)
+#define DRAW_TRIANGLE(j, k, l) \
+  image->draw_triangle(corners[j].data(), corners[k].data(), corners[l].data(), reg.color, opacity)
 
         if (reg.style == FRAME) {
           DRAW_CYLINDER(0, 1);
@@ -1907,21 +1921,29 @@ void DumpImage::create_image()
           // inconsistent style. should not happen.
           if (!myreg) continue;
 
-          length = myreg->hi - myreg->lo;
           radiuslo = myreg->radiuslo;
           radiushi = myreg->radiushi;
           if (myreg->axis == 'x') {
             xdir = 1.0;
-            lo = {myreg->lo, myreg->c1, myreg->c2};
-            hi = {myreg->hi, myreg->c1, myreg->c2};
+            double xlo = MAX(myreg->lo, domain->boxlo[0]);
+            double xhi = MIN(myreg->hi, domain->boxhi[0]);
+            length = xhi - xlo;
+            lo = {xlo, myreg->c1, myreg->c2};
+            hi = {xhi, myreg->c1, myreg->c2};
           } else if (myreg->axis == 'y') {
             ydir = 1.0;
-            lo = {myreg->c1, myreg->lo, myreg->c2};
-            hi = {myreg->c1, myreg->hi, myreg->c2};
+            double ylo = MAX(myreg->lo, domain->boxlo[1]);
+            double yhi = MIN(myreg->hi, domain->boxhi[1]);
+            length = yhi - ylo;
+            lo = {myreg->c1, ylo, myreg->c2};
+            hi = {myreg->c1, yhi, myreg->c2};
           } else {    // myreg->axis == 'z'
             zdir = 1.0;
-            lo = {myreg->c1, myreg->c2, myreg->lo};
-            hi = {myreg->c1, myreg->c2, myreg->hi};
+            double zlo = MAX(myreg->lo, domain->boxlo[2]);
+            double zhi = MIN(myreg->hi, domain->boxhi[2]);
+            length = zhi - zlo;
+            lo = {myreg->c1, myreg->c2, zlo};
+            hi = {myreg->c1, myreg->c2, zhi};
           }
         }
 
@@ -1936,16 +1958,25 @@ void DumpImage::create_image()
           radiushi = myreg->radius;
           if (myreg->axis == 'x') {
             xdir = 1.0;
-            lo = {myreg->lo, myreg->c1, myreg->c2};
-            hi = {myreg->hi, myreg->c1, myreg->c2};
+            double xlo = MAX(myreg->lo, domain->boxlo[0]);
+            double xhi = MIN(myreg->hi, domain->boxhi[0]);
+            length = xhi - xlo;
+            lo = {xlo, myreg->c1, myreg->c2};
+            hi = {xhi, myreg->c1, myreg->c2};
           } else if (myreg->axis == 'y') {
             ydir = 1.0;
-            lo = {myreg->c1, myreg->lo, myreg->c2};
-            hi = {myreg->c1, myreg->hi, myreg->c2};
+            double ylo = MAX(myreg->lo, domain->boxlo[1]);
+            double yhi = MIN(myreg->hi, domain->boxhi[1]);
+            length = yhi - ylo;
+            lo = {myreg->c1, ylo, myreg->c2};
+            hi = {myreg->c1, yhi, myreg->c2};
           } else {    // myreg->axis == 'z'
             zdir = 1.0;
-            lo = {myreg->c1, myreg->c2, myreg->lo};
-            hi = {myreg->c1, myreg->c2, myreg->hi};
+            double zlo = MAX(myreg->lo, domain->boxlo[2]);
+            double zhi = MIN(myreg->hi, domain->boxhi[2]);
+            length = zhi - zlo;
+            lo = {myreg->c1, myreg->c2, zlo};
+            hi = {myreg->c1, myreg->c2, zhi};
           }
         }
 

From 6e15c1d172520373c68f0de2e3e324577a08b0dc Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Fri, 27 Feb 2026 12:22:18 -0500
Subject: [PATCH 10/23] add color style "atom" to the body particle
 visualization

---
 doc/src/dump_image.rst      | 49 ++++++++++++++++++++++---------------
 src/GRAPHICS/dump_image.cpp | 43 ++++++++++++++++++++++----------
 2 files changed, 59 insertions(+), 33 deletions(-)

diff --git a/doc/src/dump_image.rst b/doc/src/dump_image.rst
index 689cd33c36a..6d98cde8e68 100644
--- a/doc/src/dump_image.rst
+++ b/doc/src/dump_image.rst
@@ -56,7 +56,7 @@ Syntax
          level = mesh refinement level, value between 1 (low resolution) and 6 (ultra high resolution)
          width = diameter of wireframe edges (distance units) (ignored for triangles)
        *body* = color bflag1 bflag2
-         color = *type* or *index*
+         color = *type* or *index* or *atom*
          bflag1,bflag2 = 2 numeric flags to affect how bodies are drawn
        *compute* = computeID color cflag1 cflag2
          computeID = ID of computes that generates objects to draw
@@ -585,25 +585,34 @@ is used to define body particles with internal state
 body style.  If this keyword is not used, such particles will be drawn
 as spheres, the same as if they were regular atoms.
 
-The :doc:`Howto body <Howto_body>` page describes the body styles
-LAMMPS currently supports, and provides more details as to the kind of
-body particles they represent and how they are drawn by this dump
-image command.  For all the body styles, individual atoms can be
-either a body particle or a usual point (non-body) particle.  Non-body
-particles will be drawn the same way they would be as a regular atom.
-The *bflag1* and *bflag2* settings are numerical values which are
-passed to the body style to affect how the drawing of a body particle
-is done.  See the :doc:`Howto body <Howto_body>` page for a
-description of what these parameters mean for each body style.
+The :doc:`Howto body <Howto_body>` page describes the body styles LAMMPS
+currently supports, and provides more details as to the kind of body
+particles they represent and how they are drawn by this dump image
+command.  For all the body styles, individual atoms can be either a body
+particle or a usual point (non-body) particle.  Non-body particles will
+be drawn the same way they would be as a regular atom.  The *bflag1* and
+*bflag2* settings are numerical values which are passed to the body
+style to affect how the drawing of a body particle is done.  See the
+:doc:`Howto body <Howto_body>` page for a description of what these
+parameters mean for each body style.
 
 .. versionchanged:: 11Feb2026
 
-The there are currently two supported settings for the *color* value:
-*type*, or *index*.  With the *type* setting the body particles will be
-colored according to the atom type of the particle.  With the *index*
-setting the coloring follows the body index instead.  For both settings,
-the value (type or index) is mapped to the colors of atom types.  The
-list of colors is by default as follows:
+   added *index* color style
+
+.. versionchanged:: TBD
+
+   added *atom* color style
+
+The there are currently three supported settings for the *color* value:
+*type*, *index*, or *atom*.  With the *atom* setting, the color follows
+the coloring selected for coloring atoms (including using color maps).
+With the *type* setting the body particles will be colored according to
+the atom type of the particle.  With the *index* setting the coloring
+follows the body index instead.  For both settings, the value (type or
+index) is mapped to the colors of atom types thus the coloring style for
+atoms **must** be set to *type*.  The list of colors is by default as
+follows:
 
 * type 1 = red
 * type 2 = green
@@ -614,9 +623,9 @@ list of colors is by default as follows:
 
 and repeats itself for types > 6.  This list can by changed with the
 :doc:`dump_modify acolor <dump_image>` command.  If more different
-colors than atom types are desired, the number of atom types must be
-increased when using either the :doc:`create_box <create_box>` or the
-:doc:`read_data <read_data>` command.
+colors than atom types are desired, the *number of atom types* must be
+*increased* correspondingly when using either the :doc:`create_box
+<create_box>` or the :doc:`read_data <read_data>` command.
 
 ----------
 
diff --git a/src/GRAPHICS/dump_image.cpp b/src/GRAPHICS/dump_image.cpp
index f631ffa2028..9c5219728a2 100644
--- a/src/GRAPHICS/dump_image.cpp
+++ b/src/GRAPHICS/dump_image.cpp
@@ -294,11 +294,13 @@ DumpImage::DumpImage(LAMMPS *lmp, int narg, char **arg) :
       if (iarg+4 > narg) utils::missing_cmd_args(FLERR,"dump image body", error);
       bodyflag = YES;
       if (strcmp(arg[iarg+1],"type") == 0) bodycolor = TYPE;
+      else if (strcmp(arg[iarg+1],"atom") == 0) bodycolor = ATOM;
       else if (strcmp(arg[iarg+1],"index") == 0) bodycolor = INDEX;
       else
-        error->all(FLERR, iarg+1, "Dump image body only supports color by type or index");
-      if (acolor != TYPE)
-        error->all(FLERR, iarg+1, "Must color atoms by type with body particles");
+        error->all(FLERR, iarg+1, "Dump image body only supports color by type, atom, or index");
+      if ((bodycolor != ATOM) && (acolor != TYPE))
+        error->all(FLERR, iarg+1,
+                   "Must color atoms by type with body particles colored by type or index");
       bodyflag1 = utils::numeric(FLERR,arg[iarg+2],false,lmp);
       bodyflag2 = utils::numeric(FLERR,arg[iarg+3],false,lmp);
       iarg += 4;
@@ -1332,6 +1334,18 @@ void DumpImage::create_image()
       } else if (bodycolor == INDEX) {
         itype = (body[j] % atom->ntypes) + 1;
         color = colortype[itype];
+      } else if (bodycolor == ATOM) {
+        if (acolor == TYPE) {
+          itype = static_cast<int>(buf[m]);
+          color = colortype[itype];
+        } else if (acolor == ELEMENT) {
+          itype = static_cast<int>(buf[m]);
+          color = colorelement[itype];
+        } else if (acolor == ATTRIBUTE) {
+          color = image->map_value2color(0,buf[m]);
+        } else {
+          color = image->color2rgb("white");
+        }
       } else {
         color = image->color2rgb("white");
       }
@@ -1344,8 +1358,12 @@ void DumpImage::create_image()
           image->draw_sphere(bodyarray[k],color,bodyarray[k][3],opacity);
         else if (bodyvec[k] == Graphics::LINE)
           image->draw_cylinder(&bodyarray[k][0],&bodyarray[k][3],color,bodyarray[k][6],3,opacity);
-        else if (bodyvec[k] == Graphics::TRI)
+        else if (bodyvec[k] == Graphics::TRI) {
+          // brighten flat surfaces a little bit
+          image->ambientColor[0] = image->ambientColor[1] = image->ambientColor[2] = 0.2;
           image->draw_triangle(&bodyarray[k][0],&bodyarray[k][3],&bodyarray[k][6],color,opacity);
+          image->ambientColor[0] = image->ambientColor[1] = image->ambientColor[2] = 0.0;
+        }
       }
 
       m += size_one;
@@ -1846,15 +1864,14 @@ void DumpImage::create_image()
           double yhi = MIN(myreg->yhi, domain->boxhi[1]);
           double zhi = MIN(myreg->zhi, domain->boxhi[2]);
 
-          corners = cornerdata{
-              vec3{xlo, ylo, zlo},
-              vec3{xlo + myreg->xz, ylo + myreg->yz, zhi},
-              vec3{xlo + myreg->xy + myreg->xz, yhi + myreg->yz, zhi},
-              vec3{xlo + myreg->xy, yhi, zlo},
-              vec3{xhi, ylo, zlo},
-              vec3{xhi + myreg->xz, ylo + myreg->yz, zhi},
-              vec3{xhi + myreg->xy + myreg->xz, yhi + myreg->yz, zhi},
-              vec3{xhi + myreg->xy, yhi, zlo}};
+          corners = cornerdata{vec3{xlo, ylo, zlo},
+                               vec3{xlo + myreg->xz, ylo + myreg->yz, zhi},
+                               vec3{xlo + myreg->xy + myreg->xz, yhi + myreg->yz, zhi},
+                               vec3{xlo + myreg->xy, yhi, zlo},
+                               vec3{xhi, ylo, zlo},
+                               vec3{xhi + myreg->xz, ylo + myreg->yz, zhi},
+                               vec3{xhi + myreg->xy + myreg->xz, yhi + myreg->yz, zhi},
+                               vec3{xhi + myreg->xy, yhi, zlo}};
         }
 
         for (int i = 0; i < 8; ++i)

From 7858b25eeb813787749f12227b202cb8545a89fb Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Fri, 27 Feb 2026 12:33:53 -0500
Subject: [PATCH 11/23] tweak lighting for triangles of body particles

---
 src/GRAPHICS/dump_image.cpp | 8 +++++++-
 1 file changed, 7 insertions(+), 1 deletion(-)

diff --git a/src/GRAPHICS/dump_image.cpp b/src/GRAPHICS/dump_image.cpp
index 9c5219728a2..133723d5ff3 100644
--- a/src/GRAPHICS/dump_image.cpp
+++ b/src/GRAPHICS/dump_image.cpp
@@ -1360,9 +1360,15 @@ void DumpImage::create_image()
           image->draw_cylinder(&bodyarray[k][0],&bodyarray[k][3],color,bodyarray[k][6],3,opacity);
         else if (bodyvec[k] == Graphics::TRI) {
           // brighten flat surfaces a little bit
-          image->ambientColor[0] = image->ambientColor[1] = image->ambientColor[2] = 0.2;
+          image->ambientColor[0] = image->ambientColor[1] = image->ambientColor[2] = 0.3;
+          image->keyLightColor[0] = image->keyLightColor[1] = image->keyLightColor[2] = 0.8;
+          image->fillLightColor[0] = image->fillLightColor[1] = image->fillLightColor[2] = 0.45;
+          image->backLightColor[0] = image->backLightColor[1] = image->backLightColor[2] = 0.8;
           image->draw_triangle(&bodyarray[k][0],&bodyarray[k][3],&bodyarray[k][6],color,opacity);
           image->ambientColor[0] = image->ambientColor[1] = image->ambientColor[2] = 0.0;
+          image->keyLightColor[0] = image->keyLightColor[1] = image->keyLightColor[2] = 0.9;
+          image->fillLightColor[0] = image->fillLightColor[1] = image->fillLightColor[2] = 0.45;
+          image->backLightColor[0] = image->backLightColor[1] = image->backLightColor[2] = 0.9;
         }
       }
 

From f29e03e0e5c7a94025f9e7dba14c8ce1950fe181 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Fri, 27 Feb 2026 19:51:46 -0500
Subject: [PATCH 12/23] add support for entering colors as hex numbers

---
 doc/src/dump_image.rst      | 18 +++++++++++++++---
 src/GRAPHICS/dump_image.cpp | 24 ++++++++++++++++++------
 2 files changed, 33 insertions(+), 9 deletions(-)

diff --git a/doc/src/dump_image.rst b/doc/src/dump_image.rst
index 6d98cde8e68..27d8279a2db 100644
--- a/doc/src/dump_image.rst
+++ b/doc/src/dump_image.rst
@@ -181,9 +181,10 @@ Syntax
          color = name of color for simulation box lines and processor subdomain lines
        *subboxtrans* arg = transparency
          transparency = transparency for simulation subbox lines (value between 0 (invisible) and 1 (fully opaque))
-       *color* args = name R G B
+       *color* args = name R G B *or* name hex
          name = name of color
          R,G,B = red/green/blue numeric values from 0.0 to 1.0
+         hex = 24-bit RGB color in hexadecimal
        *ccolor* args = computeID color
          computeID = ID of the compute
          color = name of color for image objects provided by this compute when using "const" color style
@@ -218,6 +219,7 @@ Examples
 
    labelmap atom 1 C 2 H 3 O 4 N
    dump_modify 1 acolor C gray acolor H white acolor O red acolor N blue
+   dump_modify 1 color gray80 0.8 0.8 0.8 color gray20 0x333333
 
 Description
 """""""""""
@@ -1134,12 +1136,22 @@ dump_modify color option.
 
 ----------
 
+.. versionchanged:: TBD
+
+   add support for entering colors in hexadecimal
+
 The *color* keyword allows definition of a new color name, in addition
 to the 140-predefined colors (see below), and associates three
 red/green/blue RGB values with that color name.  The color name can
 then be used with any other dump_modify keyword that takes a color
-name as a value.  The RGB values should each be floating point values
-between 0.0 and 1.0 inclusive.
+name as a value.  The RGB values should be either specified as three
+floating point values between 0.0 and 1.0 inclusive or as a single
+24-bit hexadecimal number. The following two commands are equivalent.
+
+.. code-block:: LAMMPS
+
+   dump_modify 1 color mygray 0.431 0.498 0.502
+   dump_modify 1 color mygray 0x6e7f80
 
 When a color name is converted to RGB values, the user-defined color
 names are searched first, then the 140 pre-defined color names.  This
diff --git a/src/GRAPHICS/dump_image.cpp b/src/GRAPHICS/dump_image.cpp
index 133723d5ff3..57115fd4648 100644
--- a/src/GRAPHICS/dump_image.cpp
+++ b/src/GRAPHICS/dump_image.cpp
@@ -2582,12 +2582,24 @@ int DumpImage::modify_param(int narg, char **arg)
   }
 
   if (strcmp(arg[0],"color") == 0) {
-    if (narg < 5) utils::missing_cmd_args(FLERR, "dump_modify color", error);
-    int flag = image->addcolor(arg[1],utils::numeric(FLERR,arg[2],false,lmp),
-                               utils::numeric(FLERR,arg[3],false,lmp),
-                               utils::numeric(FLERR,arg[4],false,lmp));
-    if (flag) error->all(FLERR, argoff + 1 + flag, "Incorrect dump_modify color command");
-    return 5;
+    if (narg < 3) utils::missing_cmd_args(FLERR, "dump_modify color", error);
+    if (utils::strmatch(arg[2], "^0x[0-9a-fA-F]+$")) {
+      char *ptr = nullptr;
+      auto val = strtol(arg[2], &ptr, 16);
+      double rval =  ((val & 0xff0000) >> 16) / 255.0;
+      double gval =  ((val & 0x00ff00) >> 8) / 255.0;
+      double bval =  (val & 0x0000ff) / 255.0;
+      int flag = image->addcolor(arg[1], rval, gval, bval);
+      if (flag) error->all(FLERR, argoff + 1 + flag, "Incorrect dump_modify color command");
+      return 3;
+    } else {
+      if (narg < 5) utils::missing_cmd_args(FLERR, "dump_modify color", error);
+      int flag = image->addcolor(arg[1],utils::numeric(FLERR,arg[2],false,lmp),
+                                 utils::numeric(FLERR,arg[3],false,lmp),
+                                 utils::numeric(FLERR,arg[4],false,lmp));
+      if (flag) error->all(FLERR, argoff + 1 + flag, "Incorrect dump_modify color command");
+      return 5;
+    }
   }
 
   if (strcmp(arg[0],"ccolor") == 0) {

From 6e4e4f12cb9e3cd97272d82e3417910e2f494e15 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Sat, 28 Feb 2026 09:46:44 -0500
Subject: [PATCH 13/23] add missing references

---
 doc/src/Packages_details.rst | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/doc/src/Packages_details.rst b/doc/src/Packages_details.rst
index c21e14cb6c8..8e2636b2f8a 100644
--- a/doc/src/Packages_details.rst
+++ b/doc/src/Packages_details.rst
@@ -263,6 +263,8 @@ particle models including ellipsoids, 2d lines, and 3d triangles.
 * :doc:`pair_style gayberne <pair_gayberne>`
 * :doc:`pair_style resquared <pair_resquared>`
 * :doc:`pair_style ylz <pair_ylz>`
+* :doc:`pair_style line/lj <pair_line_lj>`
+* :doc:`pair_style tri/lj <pair_tri_lj>`
 * `doc/PDF/pair_gayberne_extra.pdf <PDF/pair_gayberne_extra.pdf>`_
 * `doc/PDF/pair_resquared_extra.pdf <PDF/pair_resquared_extra.pdf>`_
 * ``examples/ASPHERE``

From 7cb64691568808b4ddb57a163bf9396021866f6b Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Sat, 28 Feb 2026 09:46:57 -0500
Subject: [PATCH 14/23] improve error message

---
 src/random_mars.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/random_mars.cpp b/src/random_mars.cpp
index fb6fddc1eea..da646401234 100644
--- a/src/random_mars.cpp
+++ b/src/random_mars.cpp
@@ -36,7 +36,7 @@ RanMars::RanMars(LAMMPS *lmp, int seed) : Pointers(lmp),
   double s,t;
 
   if (seed <= 0 || seed > 900000000)
-    error->one(FLERR, Error::NOLASTLINE, "Invalid seed for Marsaglia random # generator");
+    error->one(FLERR, Error::NOLASTLINE, "Invalid seed {} for the Marsaglia random # generator", seed);
 
   save = 0;
   u = new double[97+1];

From dc41a502d239caaccd4607d3b7303693e3b3720d Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Mon, 2 Mar 2026 17:46:52 -0500
Subject: [PATCH 15/23] improve visualization now that we have a smarter
 selection

---
 examples/GRAPHICS/in.peptide-hbonds | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/examples/GRAPHICS/in.peptide-hbonds b/examples/GRAPHICS/in.peptide-hbonds
index 37b1752f297..6c7c22953c6 100644
--- a/examples/GRAPHICS/in.peptide-hbonds
+++ b/examples/GRAPHICS/in.peptide-hbonds
@@ -62,14 +62,14 @@ fix label all graphics/labels ${vizsteps} text " Hydrogen bonds donated:  $(c_hb
                                           text "Hydrogen bonds accepted: $(c_hb2:%02.0f)" 210 30 0.0 &
                                                size 24 backcolor darkgray
 # create colored arrows to go with text labels
-fix obj   all graphics/objects ${vizsteps} arrow 5 80.0 61.0 39 80.0 67.0 39 0.3 0.2 &
-                                   arrow 6 80.0 61.0 37.2 80.0 67.0 37.2 0.3 0.2
+fix obj   all graphics/objects ${vizsteps} arrow 5 80.0 57.0 35.0 80.0 66.0 35.0 0.3 0.2 &
+                                           arrow 6 80.0 57.0 33.5 80.0 66.0 33.5 0.3 0.2
 
 # combine the graphics into visualization
 dump viz viz image ${vizsteps} hbonds-*.ppm element element size 800 800 zoom 2.1 view 80 0 bond atom 0.3&
-                       fsaa yes ssao yes 12384 0.6 shiny 0.1 box no 0.1 center s 0.5 0.52 0.52 &
+                       fsaa yes ssao yes 12384 0.6 shiny 0.1 box no 0.1 center s 0.65 0.49 0.42 &
                        compute hb1 const -0.4 0.3 compute hb2 const -0.4 0.3 &
-                       fix label const 1 0 fix obj type 0.0 0.0 fix graphics element 1 0.1
+                       fix label const 1 0 fix obj type 0.0 0.0
 dump_modify viz pad 5 boxcolor white backcolor darkgray backcolor2 silver &
     element C C O H N C C C O H H S O H ccolor hb1 cyan ccolor hb2 magenta
 

From 8f09a286a2d8b58bdb504937594c23b68f4caf22 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Mon, 2 Mar 2026 18:18:28 -0500
Subject: [PATCH 16/23] improve peptide example with added imported image
 showing hbond definition

---
 examples/GRAPHICS/README            |   3 ++-
 examples/GRAPHICS/hbond-diagram.tga | Bin 0 -> 134647 bytes
 examples/GRAPHICS/in.peptide-hbonds |  16 +++++++++-------
 3 files changed, 11 insertions(+), 8 deletions(-)
 create mode 100644 examples/GRAPHICS/hbond-diagram.tga

diff --git a/examples/GRAPHICS/README b/examples/GRAPHICS/README
index 1779b4af2e6..7f4024f3cf1 100644
--- a/examples/GRAPHICS/README
+++ b/examples/GRAPHICS/README
@@ -21,7 +21,8 @@ Here is a list of the inputs and which graphics features they demonstrate
 
 - in.peptide-hbonds:        peptide with surrounding water molecules
                             visualize donated and accepted hydrogen bonds
-                            https://youtube.com/shorts/1QEjIITapwQ
+                            overlay graphics from image without background
+                            https://youtube.com/shorts/Q6KBE08AG8c
 
 - in.obstacle-lines         trajectory lines with obstacle flow example
                             https://youtu.be/9HEsGaOsdik
diff --git a/examples/GRAPHICS/hbond-diagram.tga b/examples/GRAPHICS/hbond-diagram.tga
new file mode 100644
index 0000000000000000000000000000000000000000..aa91f6a2c6d16302bebe920f4f129dd56e4321ca
GIT binary patch
literal 134647
zcmc${NwY3ld7fRT3O+}CBM&_I!V`xhH1O0O_z(C8c;dg{nFk(u;Mt*Q4t?JUD2}S2
z=&C4$=q1n*ecyLNNR$wr=!^g%2_aP9b=~<q`L1{E%)N6LrSOW#oo}wSp6-5{wKCsk
z-*wmjdsnxIp8UW)|KDAY(SFSF{6G4mKl1fndn*Gy=O;DkFF4=z=}&+9DbU{gzy8<%
zddh$KFaPDPfB*0Q{YyXo@sHp2umAPG&h>x)@BjTNjNSF-H^13=d;INx;p61@&0_i$
z&Nu$zFa82CHxYRQb<pB%X!B&wT%im7XoAv3vw6CCT}EFmHlO4CB*g#j@BS`={@Z{1
zZ~y5({ii$12NmSgck9+I#{bX%`9GV2>HOR4u~4Gz;~Y=MESMt&Ll@&ejqoQOwEX+m
zWAO>a(mu%XBqT33BMY0quGq=uOy2nOjcd?<sObNg^F7D}O+d!E-$CHr3&J3<Atha3
zx|-}x(f*-Ri6yK8OeFjW3HNYD9#Ss}pidp5Y#rZqHXF!(aa4|2oHX2{d^fXff534M
zVMp*|53}El9Fx@~O@d;T33?YsKbOg#b8SB%m1M@X(uo6Hz=lb3-)dU;FDLQsCGwGz
z6qH222hs1zW5i<<-g=06nsk%Pp*XxwYe3;jYm`^kJYjOIsWx@H8^<h+AJ70)TA~^B
z(HO@WkaKBqYH1JoB}L`oCscy^+TO@<H`t;#0r493=DcFR0RJ&2O+^{woP75ZMQARw
zt@Z{wpD+P$DhUtDlt=(s!sS(GXcN&3+FxdpXfIw_3b|A7Mt97;B>nGy;@H2wG=4yq
zco*%^<dvv#7P5~rd={F{xEtHB@RFE#)o;!9wQ#Wg9*6O}6v|hqAFf5!r9il@s2^8L
z>G@q|+*pPs-)b+X^S5BB6KwZ!=xbbt<6eH!3~2D-PZ>Tdq0XqV#QAdAYVK2hJFc%O
zwy%^E&@8bS>EN;`_xFsIa`XY{%T#WOze~#Qvwu65MuqJa95<GLi`=Ce6)0y*ue&{?
zjmeip+tYswe(GX-DaYNYu*BI*PH>c3&u|>Kf*<Otmr-^}eQ}@nxEa0-V%v*3L<-MD
zB|LSi5U!+EY<iqC3hqGlPrnT_HL)3=gGM*Nt3Ysu2rUS{SOgb!chvcZ$IVF{Y|rN?
zT_O(0oB}05hH7}}YF*qIKR!<CV0#XS*|t=J7oFm#(kaz~aTFw;Grpj+lWvxe--e&T
zVS5%w@j>1T%PaV(mKp2psotoUF5OY;N52g-i-+wN$GJvBCb^KFvPN|;QODx{_&7Zi
zN}H9VtOuPJh|XEHaR+Ez5JXq#^6kgX$;x7TDu<gH;Ey6_;o%lhQ^{^}+;Q(fch-?d
zA2+9`L1}vm$3<N#k>BIsQk1Sx<QtEV&ApJ>?&cT(!SGXhzmGAMz+G1;@ik{&pN;kV
zpNkatN64Z9osWIM5#YxOSBU@8<KXcdX{aQG!1HzWN$y7`xlAgrFmTpB&((E!DD&6S
zU={KrcPXf{^@hCc%?O$Sz9Vz{3?0|u;i~j2X|RLPZ@CrP;kG!hN2+kma--ojc>E&_
zvILP*So6H|l=^}i4RZH!{-`sr!9$sUCk+-7ypc<MmA(C#2_!OtCwBH2AG#(Fw|rUV
z7L3CERO!FY=&2Asx4-}Z`Y*IQ^pv}fviT(>zmSZm0qqFD>Jp$78G>>c?xWplf0TzZ
ze<pXLhG>plFfFvd<p??qBRgon_nJI@&U}~fvMO@X8N=v&?7NO22~8#>939&`xwsCG
zzhS<Q|KxrVBNxnVZVa3?*PH9`&{iMyPG}a+?i;o}z`$91D_7Uy@e}6D{jQ6l6By3J
zjSxGNRu?|OILtt4Kj09@jvLMu$%(Zkk~`LR>eZXB$>WF2|7&T$6wf}V)`f~!^Q$U~
za_^!b(Q~1_fu6$*w9<$Xniz9R%b{>a5mI>v+t*!(fzpy2ph;~EYp^W@vpOM$s@2x-
zV0%A&jH07FzQ_FUNV8<kv&9yShZU}{Rqh_Px39xME8Q}RZHzf#OYgG(q6;OW%s#%a
zzAh7`cHIM?u%`|UNIs1jm4h3j46e-$K;gVGbDPg!LGzc=4OpL>XXG|6dtmzv=Qs!L
zeRPMbw<VPlS3v!eYjIFs*0m)L!yCK@oV`V>DzbqdgKP7Kcf|P!pKn6*=h6*0lQ~10
zzJSrIaDM)^IB2hXe!0e^^CWt)%$D+%QBr$QhrzXZH+?N8N_7vjkJ3vWMI1s~wOXAR
zIBU=3>biXX44OZcZWrxWFnSf*&$t!`?e*s7F6va9G?c7Y$?SyoT{xef3(YJ&&OsZ$
zI#Qg6;(ZJaAn(ca25*&lB?oQ%iWFE7JI}+|C&X*4?1$71rX6OWI3Lq1@>v&GgV!=7
zwA9UvV!}osd-r%P21@**tYJn;Y}P}`c<|zAZ9EuQTQ8K?;qfJ?bJK%t%!%xCoZrX{
z%W>l|sMV5x#S4``%s^{%qk1A5EB-_BTQcgI<%JEY-q{RNuV$c~pOXk?R45HH8tF-U
z3x8{<r4FMt>9Op%uERh(KYJ<bY<{;{wp8zD_I2A}pq;si-i7gPvfReMrT)w**I}TY
zpXLzOjd;o3#UPSEYf*jv>9we}@KYlN5jZ5qI~f|1YO5^$<XY6)_DLzaKAcy8_l!f^
zCCTmlacX6KVx%}{>if4l%os`Y)?uA;CAG3X#<2t1`(@UyJ%^f+wD0owIQ7Swod+q~
zL2-xSRUR2TxR#MLx3yyw4<U7jp@~}Cl-G!iHRA*HkHLQKS_aeH#vZ3u)<;H)p`d(-
z{xPMu<+){E<y}dwtPe}UMr@;oIXsBGi+G^ttUb!racUWAhcLIip;C7!VVa}b!Q`{%
z*6}#CwtaA)!tkwN*pruS>bABU)Y@iT+{L25VV|XEl(%vB?Pl#bwYL3<OqbQ**2Tpy
zHGRpJm^g$4bZPFJuEjtL@0^IjckrT1iB^(;dm`O8je&JJwHE$yF=Z<?GtG8n7Xawe
z+_<RKVFq^~uz-yY8_)~CE2T!aZTPntXlEW?4n!m;y=emH#~8)HM-LknuHo+_x|~`I
zm1nL6*dgcorXj+%%Z&^8K~kms(a!WleP`LWq-+y5l-wXyN|w5+YInhruC(NHv)QO;
z*&Q69Q^tEa)>?F~wRhB8xI_p3WRc!BcqY}7Y!tSqd;fkHJNVgoEy<d?qvTAgJ@0n<
zcJ+Fcp`B!4B+dGMoLZ~i#Zkp#qk4U-MgceWp%^`9Ez9CkZH_U}LZ!Nzh3QayzD@t6
z4Xn@Pagi0}FtxJWSoSKZGn(Ufc`=mUtpzsANeATRq)K{+lFEOhr0-9XHkO_{$4RxP
zH|JC#f!92G(|*9;u2OzDySmp?m(>HL+M^_)i5Nn?`TWR!42j>V>l?00X*))$ls}Lo
zhHP;Bu|o)-m7qJ#z&uE*B=4qHKMfFhe_;e-P9$ykDP^Z!d^xH1{C;t4um_I0d4j{7
z22iA#DnmMqyWg2qdz7^0r=QusjVdE)oscuS8>C8kE5{DUpF6bS@rdD1r{ac-XHq3S
zBD)Q+e?k8aG>oJrY>$y@&s#WlIR28MEhQp`KZ1XdRC|=P1<~Kzf4QIbT)aW5l!x~@
z{=<$Ebf;b7;+a%Q4{>Zk^eg&zpotj%1da_qZIEiuoA){XnxQQv=uT6~9c5-xrz&Zy
z<ow3|%l&i>Qm4L2j@y;=pB=g+pRPfwB%coY(7dbBP`?;CLkE}GHMTc!ahzJ!@-(u`
z_E7rddXqzUrRY#_?HilYR{XD^RaWLVzhRw{=e_bOpLw=#Y!7m=>T(6GHboXXKjY)e
z=YFQO@2RrGYZ(UXjr)YxIpta$l=*t8+K$OXPh+gSqg<h*=d8V!t80)eQDiRR`izky
zcSc;Ewfnic2D$dd7bw32-YuqhOR)pJ((_)!wKym<9>Ijf;}j#$RA|T6G^WJcjvu5|
zmX5F#RRD~)NQ^x8Y!2^VqRK|_+;=^4ExXM1imSqy$RTp+IcrMXi<}K|B|abt%EyS|
zzE-}^b*w6$xY;+hS8;KaR?8x5NVq8<pC`M^W9F4kxfTay8s4klryBKq$0)E{@x;($
z-i92dRhChGNzA-Zc3m2#y3WvJFTWnSmL-px<y37HcSX!7an<N*awXoDWCybsJH{<o
z@EmT%=uUene~;2?*=^_Vd;BGr4>S7M6;n2bZ$tGOjy*drvjMVv)8HuWZRU5OU=7;J
z?8`i)xGS7IlWW;op_syO4#AB&b6FH`kSp=k9BZwZYV%o+O(6(bzMF8AR?8v_8CiS)
zb*eLC5iztBF}9S1<K#-T5{`K`70a_ixnZD%_8R2c7nY^uY(#r^E0`c%>iR)iWeM1*
z-eiH|E*WMtOXVd~HaKWs9JQn8ZN@EIdp2fq-Wi{1m1PWGqG3pSncP_G-FjSuT+0&m
ztVXebWsk14WyNAO&e*gnJ1IMzgS4xx-o{K-IkGM>EnB&RwA%Cv4y#{uhfaOFWq)&t
zUD>D}q+Mn0=oVR**!6Uj>b<@CwvL`>T5U=l-O(+(;Y;j#8uWXZ%(Sbl9lJ9B61$$N
ztb?@L^s-BI^WLIYU$c$Nvvx06dv?vVtE}4M&6sCCrNy}MZJViTEB7giJ5H`;S*!QC
z=GM2i7h4A2RN`^+RpJgy?|?61$LM1>IF3_mUUo41I=XGY%!ilgGn1~;b`*R_QTv+o
zW0F*cDVjD&wJY1neU5o{)&Vi??0<8)GBfEat!_!@Xvh<=o%D~WgLKMyk*akHPjQ}y
zshz&lKH}HcpcA>s>1ASCM)K^i?Ll{@uPt&8(yek@XY5)T@gh${du+1x9~oau)fIHw
zl6@!J)vaamZSrcJr*i%)#}3l1a(eSH`5<zzBRWjWcH|tSTjfydo=(a0hj>tPZExq5
zPlqCB!X^)Q57LR*3*?0zTZtDr9S2iV{{1!RRyn;4PPY~BkXP$WJaD<T$T>)-ob1qd
zeDF@jeJT@MsT`{iUi{LrgLGmRIgQO!ULvO|lkzUeIY_t4=`uF?`d#V_URyB|dk)em
zCr?=5A8e+af46cvcBlOLg<}Wl#4K`lcty?*ugEz_x5_DOHnt4!QSF9jDm_1U>^P|w
z<#E_f6?w1hc0B$w#}3l1ayo*i>JvFT5ga)O=~g*A8EM+A)o04&zj5p!otQnx`8(Ao
za(2dkdJVc&&Q8UA-z7Rn&Otim<koV>2k%$Tjt_p~*g-lmi=1B3rs@+pJ24eG2kBNh
zJ8AOA>b##OKX&*ysTMs;lAZj&<5+7xMgEVjLAT23jGXFO<g^Y`?T(y-bgLYAY}Mxv
zs55)zjDaGTYm1zNbgP{5#Qd&N<TP&BwF=?I4>{U{SJCO4JX&aMrl^aYj=BlcA1LP_
z-72Sbo+``-<<%8t@=xR(q+8{X-CI8RkaBj$B4=l8rd#Fg=$!ej^W>kuzQop<PT%Bl
z%#IH-e`oCb%Gnv4=~g*AI)C^Qog?QU-704X=Z`3-$H11$wM9;kO}RDGt#Wq!llh(V
z6^mlx@9I5RXR&k7Zj9fg-?Pem)4qdz_1!a^zvGAKG3Ai^^PBW(L;WNW);^UxhWd>u
zRumu7`paFwPPIXqr`0gILm9Qm`|}k=XNC1iid;D@-$6U5(Q3hthy8H=Ty5S)-#s9P
zvrx+|GaWuGMLR!LRbLOjDJfb$g)Sc=Jgu(p<yx`z!w~t_8kS&}khs_;wWZ`*K@cE|
z7hEr2Jv*rFYTZuF{@5t6;_(ss{52WE2d!FpA1>C17X3={-AfVMu`VIO41bj)n+L_T
zL&`1@k-1LC8U!*a-!j`<U}k08xUHMO%$*Y8A^7_<7AS_P23&ik4D;tDZtj?rhpUYF
zOFwauCUGsB{6Zzd@-z>#63d7{P(|eDn0-V|fr})#Wd5{gh*|K=7^pElALH>WDg6iw
zoD*s6VmhhN4SbD99lOjm3PrjdJVs|;BgU-rCzjCW-Xew-L;(*zBavW6_j0tqK~(C0
z8e)f!n6yFwxS3y(0gJSc^LI|j*VA|U0VK2M^mBEsVpPT!`Rlmk7&`U1c#N-a5SLgY
z#7+-fKsruwOtB0zc@SGFf1+rQn6a7-YEcbYkgqMgP(NdpseatS(A6qE1*su#FB&i@
zN;8_ef|mNTW4?`wK?uuS@NnrjE^<u_<mr0>;I_?BqZmaXtPxLVz7MwM?G>uGseuK~
z(szrQ8c~|se?Nym4O7U}&m^K+%wnK@qk3-M<LV99p3TFIHhDDB`=Qfs-H3L<N`Jnk
zbwq2^Ah_<te5}MpBQM`mz8Fkc$GjK@=&f+BF<{Hk7=u?b7YAi7dgL0kM`L}2G9k+h
zkfne9N>h}qM7G~H8?YJ-dFEFJEVJwDka#B1iy!TGjrV)f9le~n@;*MQl2Q7?RByA#
z+>04VxfdFDE-{me*mFlS-Q`R^U4+s9*c}dY;fUI+nQh7x>7N>S@p5^e85h>{Q&I9+
zkBQh$iKXb&+K-mLl)1G)Wf<R%#$|YrfvhEHR3J|Y4k+7(jR7gWEk@!0jG?oAFXyPQ
zK@KUB-6X=JCKR!%ub~{Rda)cZ#fUW5aq~E=fYa|qW*wrQ)rzMIM>LgTphEVjJ>6-s
ztkbY~E7)*QsZ*u#J&rDzZe(Iv`){&L-`4Mx;?Y_dd|Wiz_q$9o+d+pJL7>v&4;6J#
z;nl!$52S1R46iJPNjVI1ppDZ(H;{`ZzQNa)nR{zgPF9&;^;^0WB7Bjr4`867(+pn4
z>-y<~#K36M0TowsP&aC$A{^HilaM}Ekv1tI!Wd#kXX%WnGMVGjp{7vuuVbhny={x{
zTA%h5Svf~@4yw4CG^q{qw3syU(xh6}we(aqe5at*i?8tCmT|oUDC17QS~+fJbo?p9
z6jVDs^-2&PRCzTKc;cA~Hu>0M;^6&+J{6m3_;(PTqDXH{A~&83@m@F@b#8J9@LQ;U
z4`*VLZ0e<;_!tYjD_pAM0j=A@)R3wZ@J)+#Q@SP%(4|<P?r^+sZ3k*uECiY@Cnm43
z$;k6pw=syBm^{K{L`*9iNFZ9>?r<DY*%ImYetp=5W?ZiAv%Cr)<lA{4N5!Yzh7Ez3
zePi1KNU^QG2}eO}+w&SvS@&+BFvG3~J*<V5oPz*5^R?$AYOW^1IjSQ}`fi%E-bA3J
z{kJ)$4cLMun5QswwviVxj?pd0#W-Bxn8raqRwQHKR9A)sXC6b#_Y4oJy_yN$nkjv^
z&+%g2k>eu!AK<84#R#U1GXgdVw!UHL6e<T9s27~qj>X|_K-zN#!BkX8f{MqFIPrj{
zt5pa|L4gINYl~&m*y|WQ-P@<AEP`d`{dznFk?WY(S#mG4jIrEbw58KrOX*Zt7AWao
z?uI4}?YGU^uhv$zm>+5P{Pyqb;c*^6sgd-$F^O56qdMoTTa)TzfVwZq+=fURd6n4-
zeCRt)abLyhC!BRa=heO{_Hu@FZLtptNz>x}Mn<auw4dTf#3ieURP@X;$C;~LGZ~KO
z$-^vu71N3L$;1>Ra=WTro>?5!+3M*L8;ZUSjD3jB7<Srk89kNln`Ev&jEN{5MKB7i
zb%jyJRWgBDqjzFPgMx{O2?b2LIN{eVW0hIIzz&f*%~jq$+pNEK_Jgj`rcCC~=Q+}k
z`u07{Q5#2eY6G;*`iNmb?n5`15#Wxu%A~pS#1J_Mm<B@TxxF8?uszS?ir2mH@_lT_
ze3VKzR~2VxvB#kLb{p+jZN|6S)^j|6fn(lhztuiuL}Bd84`iY!U9qoG7$Xp^^0?GK
z$(iHMwA1CQ_8kzvstTWk3<g8Gwx~Pp?EF|kMVY_H<bH7hXiHh0T9fmGm&i%rZs~ha
zhgD~TRh#kIsxU92^Tf-~4fqW&Nc6}J%p|7t@|n`0Y6Mu~YHc^MBme_KCE2tSGSDn<
zB^*?1H3NY$1C`~`EoPulN08y>9u}{f09u_h>ZE9j1gLsD;D3a{JZ7OwCk%8@XF+2?
zvz7D>j<qhA2EU>2>571jb8_QbnzYtq&N{~yjiJdHYT5TOj)vIW?v=tZWsiI!S2IeF
z$Q}6UE}u%TA0jCI=t}z{f5Q(`?YK<#pD{l3ndNIh^IK)NK2-;mV^j9}HdyI@rTmH4
zmi&?Dssp}NuL!O`Us8TY5?;6aO8EmHKbty6p|+3wR@vB=d1SpZx<!7Ep@;pW{A5n0
zmrrQNIEXKrw)k%9({tdz>eDt8=B<n8Gkr*6`!47MKVx%KpNa3K`oRB!KJZp%x9C&;
zw(nA1nCD{UykOJ1RNtm9#e3C*3OD*wUbdrpEVC=$CSFJ6CvVWH`K+?Zn{n);40L@+
zjFeBB&U_!|Z^|BcI)A45Q~KQP-b5eIxse~PmGv%;QBL#%cju3>>s;F`BWpN%_?EfP
zg1qgVU;F&LeUe3xdh8tQkoSCf4c!Z+pY(IQ_1m)X%Z%4QqQYBlaeR(IA)pq%<B~Ck
z@{1<T+v{%PHugWnb?rmR(m@6`XxpV@*hskn(WWC|sRz%2t})TB(=E08TK5(N1hmg%
z%zZ=&wlukv?U~F_KK0w`%lZZfsHTgvQ{KSgTeSDU4t)@XDMGC%@au<tXr2carM}X0
z8w7#MtI5?AOv=NMC{Rn3eYUT0=<-(XpS(tDph2@dE2qrl+Ift-pFxr6$QjKAbG!Gn
zo9Um>M0HJ);XS<OC2TSI79I78gqQ{PFt=BWFIwp&@yM>4^FMOk7B=drW)~y*t%dw2
zE{s7e+@DXc`x<ETZ#nQVhPiJX*^Y2Wkt->bcAs2e7}ypo?M@5irI{=5`k)BTMzh=J
z&;d%|CL~bWcl}+w1Wz>PE>J{Z{`CZBuCzoxcSFP_^RrPWH4{?w`fNNBF4z>SCvZo5
zb^!Jtw3QEi^yJewbjg1mLOfSW+e24<uclzqf-c%!g|0Pw{TLIF1zqe9ok8w6{rzEy
z_<p)pBXrL{$T`^{hZ`V>i}J4uU>afq$T*krMN=|k26ZrF2UPzXppNgX08#%Lm^kCN
z*f201<F$QzIR(02K=j@8`Q{K};XQ0|9(4fdV96Mf>P+|R=2`Ftzu*f|@J7vC15MU5
z=BmE9NB8Nvo|fmIcrAE4D;&7Ytta9WM{p2qxfa!c%N%Kz2O}?>P`$zCWxDVrGg2n}
z_mG`whL34-6VkkuB#UCzm1l<h2#vB*Bw8gR=<W3R4-K&osaWgNj~gh4o1Jg_OAgJ*
zN0K<BK5mjyjoXqq0;<$Q&MsnWS;2J;8N6>H!_^L6TuD%Y<wikRFPjsGA~$kilK+zn
zr?6QRcS@H9Yzua+X;FmB2F{nd5WHp=A{b$Dq5q21s7%ULblyNALdF|OA1-qhwdsli
zP;WW3W3){=#z8&e-ZmBP?$qiD6Urk#;_#l6ZQtS$aXjxH<n+&!!7baoftmBi(iRWH
zHW5Gojh?Yk!PXsm6{SvsYz7&PahK21cLoP~gJ*(+F3f6PN;E@GkL%(e?8CJ|TJ==B
z&44~GzS2K)!X=lR3o2X^2nu%d=esCsr~Y6Ox#$&>2t&HmU-Y?(Nld&lF0Rc185iQ>
zVfx%oVi2-#p`;ryXWG8akw~Wh^j~e!^C1TY5S+Z`IEp1A2&EVU)G?;7E9lGK1fo$V
zivv@TDv5=neTcyBLN^&{h%*R)h(=_hpj9a*%1f4BiREt@1=+kIoG4VaK&QQQB?{03
z%?%7u!TSLY?Sl#-6qGoR<vx~On*%Z~M8O;B^DQZ0K}L3JAmU8hcRAG9{<=3IVx*Zv
z;0Aq|qVmrrvDUpjp%rBzi($af$(L*TEf8x2ZH!wvbytAEow`k|hOTyeZYi8zUR8~V
zZ?9ux#9xVBs0Q@+a-S)&0|Cg*B4D^>(*zC$Nv&;Rib_opmts+0U+MB|(5g+Af~G;0
z*U@*j(&x$hQktA^`yPiYy*v85IYUl(NI7m}oT<xk*Q=R0;v8eZH`%589S5Gcne(ON
zh{`zbeHrYNxn_>y2AAs?x@6nw5x2Ij^v~FCI<)75O&9YWhkC2o`_0wbjxsdMy-i0X
zJg#)*HE0pNxUyC(e<$f|MbC`{&^gieeU5_58;G<2`kkMTnvB$w@Qx?3<qRr3*p;WD
zjwf}2p4}S<+8%#7?BhvP3@<#Ymp5H_2`*=PAwM};m|Wq$d?KO5s&ZJl)M^4Eb5~#E
zSGw*X9af*Y&1^g^xc*OS1IbyoA8;68_Sb?-ZNwX|c3^Ut*1EOGErW=0oxCCNeVjAC
z_U-)koSw)Sk(!9@vW(!EJtx$TsYg5kX{&+MwKB#Nl!pWaXT|osc%gF=(^L67&^UJ>
zI$srz(kVZ+u!1F@u_73;7OvhV7z56-{fNWLl7IV7tnd;#Y39$w#0?C`4c2kyL?>Ra
zZm8F~u%m5t1ATErEk~|`TW;91bB@$GA?`w-;EevfTz!x|S|=L%tv=v&e8Bw3<b22n
z2k4X^AK=&4NsP|cBcD0cTIqb-k2%Cnj8E&4pQ>pt&Wb1SXgtv+Y_Q_Rav83bEZ0})
zdy-*A0B?#wBfjg64U0FuQP%o#4MlQ{YQM3r<B$$~@-?7iPCMZ`PKg@VphboJ`d#bg
z)04B+<vuZ5^mp!KKjEm=#TDahb@5StGjeX03UVkoSD_mHu;&*8;iyltLBGV{*)P8?
z9gp~=*%B)MH+^E!JS`o+r$Y!-0k-{=Ls0EMTRJ4W;EfJ5O!6J`_3<hZ9VCY=+4D}k
zNqzsn>>;0|Yu>8HMp8bJ<?A&AoGwM)H-RQVzE+R?$N%^rUjTT%yXD^X<xlWyt<~$w
zs(UdcxO`l2aW5w%Tz28LidSWgmwO5<;^pD`tTtlY3hxtG8xQ8-1rL`BZk()*b7~*i
zAe_-qngH;u{<*PmQ6iYYkG0`FIv1Sc)F=3L9#w<!>vgIbnY%4LG&(gwjH+1<2(Pn>
z$8w>QGX^~CzN22dvq`>=p5%4*9bMkAF0eq;TSf1unrD`ug8Xdjhfde<<hqLm-I&z7
zOVo*o-l^7tyl$BopFQhVtX18kZMd)`TaSj+&4Mc&yDcQb=zUMHSI0=v!?jwz5iw`0
zrK@8r#mnVm=GI~L>3WP$a*ei0X26NzsGdxnuLmE-$30H}n#uHKtYr}gTmHWc<c@;|
z*}bQ>T(PRAis8?Y5q60KHwK!FPtH_k-dVPy_Jue&KpWneD#UI<RrBs=KJ%F%LID4X
zlK+L$JD+TT_rppg(C`^<yx?+ES}n8%{Z`o9UT}!Z(YEdH*gteHj!gN7wq=2b*tRYi
zD=g#MA&l~9<r19FLIJ+NO(=9kB0ee&*BuIh%B5GRK-BA~FwDJ-_6Cjs1*r2S%TI~?
zzmWg)F$PR|SIJOj%oxzjbr)jHiR#eAOe$r}#Pjhu{qgu!<Lzba3VlOUi5{;bD`UO%
z8)fKs;h|d7B6Egy^nAhkD_F?>@zjCR#C-}c4ZEtG?een_K;LOf)H%!~S^{jv!mAMX
zFCGi{{X#4>zMNxU`<AQ?-a~986Y**AV18oZlJVX+^nww8cQqjm1ngGD=2@1r*l)(d
zmEd+n4B-C5H6j9c#5Q~~=1F}~NajXF&oV3Ma-R@+L^Px6s=AQ78>dvlUCl^1>YJ?s
z4A?%Efkl2=Rx$TH^2j4};U*BgQ|OIL{uA{&s|TKVNTK8eHldQ4Q9{)!)-o6uV`X5|
zwa%eGn<K6}g=P!^toLG9O+k337=TFbO?V0h*SC;p>JpFbf+&i@xdyN(g|`xiSGpuZ
zoDJiPF8Mhk<8{<go4T}=^KyOG)wd8)fiX#BhB26K0@Y`{8q`cIsW$))2N(s;syJKM
z>Fcn|=cr(yZsXG{HK^MS{EK4#oXx+~dX(D@{PQM30Do-bO8>X@s@LA@Oq{DsLaKtH
z-ga*&L8OzjWAVtuNUWeWgkY06Mn9FP8&_xf$(tt=NzC91v1DSq#hj;eT)+}&N8Hr0
z)^0i&U3qYa7J&-ITfXo4-vkNIx$RLe^vttn5J|DcjJ<g#inU6X#R3jzfupE;0s++j
zk~gWQkYyJKJ%}==hYwnDA#iX=RL05#4p!!xK+BtGSanP5;07Q8s4K=LLvW2GruZ$8
z`Z`p)3=$N{O)j;jqqV>0)P?X~q$SQaV2?W-SFv;0uM-FaT6;^4YN;AN9hUGYITq)9
zu8gx`diiXBW~GtjCG*dM!P}<05ziFVf1u33XHLvpF1W+bb$O3v-#0m|V@~y%A%sB8
zpqxgZMrBkCj1UA~2y)F?=JjDkKRjph#ddrORd)qo@-x?P@j6g72#{qN!Uj`NM_mNu
z%uVNjT*pnSSgOjg<m`N0w17`tAky5CT~6R*946t0K(yg3A*<&Pr~-z!)dqCvUz`O(
zaTaLU_ZJ+lw$3}-2v~kDIIGO_BhU=vhc?#P>*r3<FU66#$lJ$tUAE<Or*y%iUFFBQ
z=*j0!F)H=Byfp;-fMVD)Qy4CSavg^uHZ~BUohcx+^Xcgl1)TW;9%t7p?b@fgXYmr%
zgLw2qaMo(5&T`D{*x5>*_shZ0i79Om$88{G-=A~%i4@;JIjy!YqW6N&%1d`I+{%N}
zruG6etPMfy9+o5uw|g07&M;TRM9<G-Gf24kX`Z|*tbXOfPK0xea{o~-+{9?li+Bp=
z>n1D}2uZOUB`h-aVMvv^fl`FGy#n8P({Hc3Mz{rFiXnWA(C-BhH}`65IWv#T=x+4g
zhMk___j`|kPXF?|wV06~fGY;yQ2K(Qvrh}5JN|GjyHltzDxnXg4A+&<b_^uh7@B20
zbVks{-jr8HA<H)zICCXZ`n0~u053GWB`;tihElZ4vx$|3LrVk5?BOP?C&+$=kp7ZD
z-Ob;itCO*U{Ec*JdRz>{Bq}C*eR|6(KtPGfWomY5%=he?QeURQzTt~7I7_gO@^J-m
zpN^p9A?)cn`-d5?lpV+*m02$nGR^?##DqI4#^)eupFVe?bn4kcnGXa&8-ojs_?_43
zLrRbN7E8U+jZScmjo7;2;kNve;f8?84bMty)I41nF<TflA8qCDcX+u<WkRRUDl5t!
zpmCye%P|{(CMqhGO;AIE*=;FTYt~0Hv6jq06YX&SaT9XSx4Aa|xWv)Ha~d`5Dj$=)
zCD(_zXAki6$$sFWe{muZC-aN=>-4*TI`3=&%_kBj_2Wboib(*SGuM92Evmy*WyjZL
z1|z;nXe=*A3HEqDlK0vMBmFd18F>w0lf1Yoz@i~c(z3u8M)%WGs}*qo9xmxdi8Es&
ziI}CX`uz6W@U4^6f+Z(Q=;-3C?F4_8Sk5v{-w<f_J<Q<}qoy}`eS@P}mA|@xNG@==
z36pLO+A^^Ii-C5^gU$wSX2UO^W5mXCcgi2^V;&mG;A<R~SCn=3sdme0HLrAvDX@-t
zoMn){E%WGWt=h!yLN!U>$dRuG)U|yb86l<rP4t}#Q))e`qOg;n#*vSyh)Th&L9XIX
zJu6H~%QoD}B{9mwECvBGQ?51$+Di=bqc$!#g}EGOoAmx!o+e-m<7^-y766IA%3*vu
z=PXAPLldH@)Wo&MHz4V|%t2S%)A^x3c~&xr{eAk*5^78y3+0MA@h*ceb5ue(^^{O@
zS|^u-qCD2epo+Nb7z9Zo35O`#81RO;_5*g=Qh||}so3@o`c7HonmKs1@ANXcRQeJG
zl|{}zWs#iL2ZE(U^e5<IQw3ey*((fWOOY6~>Y%H-;S}x3*Lj#W1e|;6HFvS+JzhiR
zBLZHo9}7#|5ia&=!aX2O7*YXI@7;8CjdDV^p^V<bz|J50?{JL5cGUFSo-VyGm&*1M
zj*-&7_uGe>QEjN<$)UmWXalckQ>9w^JJMeb+6$SI2PUwrVj{I0R=1)Z2=+yzF1^Ft
zaaD1it53e$GR(O{-^j1A!Rs*)>C9N$mw7ez_czq@G-v=Wq7U#;oX5zsg>A|roFl8^
z#JK{075%m4P^%*cFSI__mNVcX594?I>a=}=>sCFPR@rC&D(ea!ZP)r?umKnOMCK5b
z@fR58PIAb{_zRoIu#(<3^02fha|~eTbjw_H)*tpo9(*HjNP@^4lA!55k33d+<Yjal
z{ISl(=8P47jcJz50guQ$$Y+(?yld@lJ)<L~X4ENi!FJqQIX`j_@>%7ARbzN2bHXe#
zvB{lL&x^bR^sB6j4`Z`=c>W?s@~SZWL$2NZ=ZV@-7jo}$fX6Bi9$WePh2lSsMJBfN
zQ#_`mMdm?1%Jpy`wiVxz+xpjP7rEdXxttfN<;#P7R=Mr7+CMd~$R%Fd&y8=MEe{LJ
zS*~2d!>vs0oZ|9K-V}$8Z|*UMWg++;<g?0c-a{QzcI<2n3d_bfa)%lhxd-{Ia@%)<
z{>uG4$GEiF7rE#cxfR#;pT+PMd{(*S%NVQju_-tGRqx)*wQoh|Y3)$IBJ%(b<#~)e
zxlH>Uf7dmoax(JpUE~eEFpiG!XnA{UN@N!1D$$tcErOwj#J-(1W#r-ggfiDPWxz=O
znLEWqW18m&g|l-bvt#WFKC9g3J=OTgopPYDjodkY4)R&$Zq(?=Z2gg~O_7NoBXiI(
zwjSWI%IiF-8YgetXF2BDnEg1ubWfZ~{E57SJXV=IHTpBkoW>#(43+7&zU4k=V{VVl
ze3a`k+jBO*ZKzm{T<UG)QjdFX>;?C&H6wBdKC9gJ-4Gk{%lyGtjjs<UeMosG^8gQ@
zE5yzzKATVEP4U?HM&_`-#Lk0!R=Lf4s4vQn+}5A7eT{G4LLXu;em%%%m0NpDV@KZI
z+_=oOKDUSXIPWhHYfR)G<g?0!Pip7DN7<iJE;fyIHFAlg$Q|M!au4!Ru1EL|Jm@<6
zmBZe%MJBZ&G7ER->f0kcTHacNhxm|B?CcmRpGPLP`#i}xw!~8V{|X*H>{G5!AEM(B
z8}j<3WAI4+RgN^a`MOx;uXBU#K|ZV8_QkM%$m<gvt52&QG`^8L%s=z#(7(k;xp~0Z
zHKppib8V4ZHKP6cF8b9a`;~wDZ~txXNe8;%eURTOA1s6E<e$jzvC+rqbjv#R@BZDt
zldj&LNY}g^HTe~+eT*F_c3fsWttGt7jSGsuAIYZrR&dmtX>_2qkpy*G&^ah#RcR+G
zJ_g5ir&XSgcc?_7Qt#3-E3r(S@IJZTnG4+L7p?enS=h{iRLo<>qj_ciWMmXWNgrSO
zV6e35hh?6<Q}BTxh~^EZp~zeoRI{{IyB*JG{xsH}^o<^a2`zHpQZ+d>0d?}2mf<!u
zYE7(qY!>4pu4-#-8ES<zA9ZZXoA#z}A{J;_Vo{~j#hCStIw<p2;}GU&v37;DRl8m}
zs>-OD?IRqUv4-e8>!rd%8lV~wx-__+TSvR0o+jW9N>e+Jx%SeP7VS0y&i|bUFpe;9
zT}gmXZo|Xp{#>ry0$b(rxDj$}7PKmcpm+pn&g+)>jYql0So4!W*e-c;6AgVhtG>8;
z4*Y}6eI6c^@JYB%KL2n_gC8%Ke>>G4Apd+Adjn0-&0$-xIKXZd-Rxkqa9I9+h|`+t
zwPE?d*hj)8Acl$**$4S4-(xR8A|F@-3$mEHm|uNf+wS8y%43ze<DSR_%fv;kN*|Oz
zGSfHU9JvSi7$~uI$1lo=ofSrr*?N@gnX=e=fX6DYeNlYnVLS6<`8G1Um<{t|=RrQo
z^;mEE2fS9f@V?+v{F}ZZmV77uD4$jC#9NuRKjBzoWb`ZYB6mo;$UVr%^ozZ<=!|b5
zD)$b@T7nDH;9Nd*j_o#w^(?wDW?#2Eycl_`G8?mD_asf`xBjDV=^NwPxwkpi8|W+e
zDA&XK8a;^#F4y+QyjmB-s;OuBG{l0RdpOF+8WXu4_l08%e!sjLyBdq!UJr)(k$aGj
za#P<rCx&HO*^x_bj`5@1*@ySjx8iaIk5y*#9r8h$c~C$3J+>0Rkvq(f+=G0q+p)Ft
zVW>H)+{UfeiO4P7onw1|V_h?@;GxW{7wxBM9geNU%ouBt*>OC~f30#4@>%7=XZ96r
zmIq~B?}&B#T!xL;c@2o%gM55&8ha<+^0(yy-NY+$JN_73W<K^<<R0XsTn}q$?5(x5
zdG&)8ZF9p~!-2X*3?=Dx*JHIxpQ5LwXBXemTM-khdp;`k+=%Y}X}2=-c|Jt3ddimO
zmHC|=RV5=Al;h<gx*mjRuS(8*l<Topw_#;bc5V;aqg4&E<LuSGBF?$5=Qzse1<Z}y
zsS?V^_D+tuN-8%IF+^5m9^|1+kG<*^nbT?!jJwJV)#?q`;Iqn|{0ED)xAS-QS1=}K
z<_w73gM3~H-`G3F>mv2B$0B!%&&ch{F|4p}RL6sSl<UzIZphLW{Ra+D2_Dg16f3oy
z|0d^F{{dF3^hOU0Wz+3#yjUM}RrE&mM)8Q~gM9pOPpob}!~Q|pkvrAP$nEGIGB0uu
z@=>mbx1NcXAxf3|Rt|4NJBh1M<XON*@p$NZtX|CYNbfir>H$NMzJ6|~oEp(xKL(T{
z`XHZGE|%S5CD5bF>-H92l+`0Bik?eI#2#m(WDh^g7*<25(*2l{`-6dot$ne$vuE(n
z!_Ga(=OxUITrAzPIJ}S-ohs%dx8oNR)y0;%2l*)1V^7aBua&(YtNHW({gG?%@#_we
zyQAlu<u%zGx%5Tu;Frif$Va&zd-f(bJ4#1;!8h20H{@yN-cdd;V{YUUC-Lv#FPWtR
zZ`r$1KO*-aALV+i^<!9Fm75i3viIZ8uUh6DTk2capeuO1ocWPiE9*vCe?3!Hl?u{`
z-6$LpeVmojJy!OUwhuaXk)6`#+VKb1W3x)W$Q?#$*)&ct#_ld3xVz^@`f*mTfX-vj
z?(EVpvKyT}z61I8>#<s;UzE`H<aTVb{Szf6DP8I2{CC)YoRwPukGW#dm6p5Di|k^1
zC-x8~?}X?YtXAn$9+mLNBE5H>B_@V<ZVyh6=!1M-i2#p1Rc|{s)y2rIHHqQ0ck}ls
zk5%T5j<@ACt$Cj;{&6w^z-3zG9^~^XC_MJ`yw9;MTRX9ay50Af<60Vp_%ria<?iU2
z`R&_UlVWeL%gAruk$aHOt0D4O?H$5Xe!+r>9`_27+XZw;rhJ&G*On{zta6zrKGtqv
zRb;obkQ+_zOsYQTNAz)4qNxVs-fr4*Ce7h^!SIb_Qu;MmDLwvd_TWp`+N|%s`3WQm
zozm=6bgv;lEbG#`j7HO|S^1>5pE3t|JT`D+bXim@*apW$0cOwn^?nto+prslxUD*L
zWTC!-yJZJ6jYaxKmPPtBcaP6)20VmubSzzI?u>OLMg)1z_M2H-?}p~#?{QW>hrP`)
zhRRT;7+JRFlP@EEBi-&?kCoE%B$$d3+lSSFq4lmaW0hJFy^(Hz!h8Wa!pBu9E7N|2
z<naO)CLD?S3dcy_NVg9t{aV~s9d_cRpCQE_%tFBOaV%VsKo}Av(hu@+E0XO<yL711
z%J9kP7;hM5_h&gQ$aA9Q>Ev-%?!=USNxbFJY;(jYW-o?QeqwlmVZYL^!Aj{K-o*x2
z&RnkT*}U>uP{vTfff$a+h#s)Y^Rt6|-1?5&-;KM1XDRw(yhS=?FR<4G4z<et@5Sj3
z**3VXI!yL01gCdJM_fL~VC){2wZ7vYG1}wd{YLIN+MNS_(maD>!+nZfACV3&%1!Wb
zR^DhSeUwnqOmCVBtm{-M7wH&0Cz;Z(!Aj{K8WOCAu%&zXsNCcWMVE*;e9*Qx%|}?h
z0X80EJWa^oTas#383Hoq1Nxk5y~DEK>4zAun$-4m_>Mmcq&<bdK3k7rRsTO<ebHgq
zf<=eptlkJ4kJZCN7}LFc^9i>^G|6ONMECMC;Ge?Zv3CWZH^E1_x)OlHeh6d5SY`6z
zzp^f38}(tXVh^%Wsz(Rh<fPveCv1vBRo=cJMeVD`d*IJpi_@E#?a@e2e*Pzi0Br<1
zg0Vh6MggXzCH@NNE7&{)frkrkSPy>rwfz8Te?|AV=n_kU?5POw>;7DQ&xu{%2jt_C
zV~ksm+8{@jfzz*)i@!;l3?cU~?cz+(MM3O`;5({<$UuT>T&AhYao!Iz+hZpV{$6G<
zZ__~H(dt4Qzw9RJs4HCa+X-|n|8cyE+7x3OynSfwu@f18=g<%tk8)kw`CYBzFI=mS
zsz%VYe6?^pR%XVJFx_J-R(@&!5GzFGqjrIZpJ}Y0JVdU)KUhE6fiv@YN6tOHcSm$V
zo%R*}9+cP$?l}4va2kUMZ61|5YDW)zhEt#JEY2WmugHGWS-3MnH6)0t9hdc1W_#=e
z)z4&h9#rueaor!k1f^%nfUc;<^<Q|6`=Z+5)~{Ym(C;(bV<*D?#-WZdP%^9=n`b$O
z9=S%5D74|mKRwRgC$=72UjM26=W!O4@Xo!&nh*V+?G&`l2R_mE&%JK%p&^%PU+3>Z
z5v|dVvws4oDX@&S4y-7ELh%ubL^^*|;5bL0=jNF%UcN<On%9f}!|z8jI}e}K47yA?
z-E{h47ptu+)lEh`Kh<kq$BCR@<u#gS93goYfxPUfa4y4X-{9{U&X?=sGhvU8GQkfU
z{;2*KV)t`>6L!zwtlq$f-2C<gRL0rz@0R=}Hx!HFgoV4s)4w6N*U!Z>M`IPixsCLB
zjI0@jWq#F(2x1vyj&LcyHv@4U6X{NiOS#q&mzm*%=sb)@Lgr#eu~_~~R&AZinYS6c
z2mDEJ81v$3bm&z?`F;5M$+TDkWq(TDo#^PkWB?>uswIA?7LYXJ25{(DdV>SAl%Y%3
z+1L#bkp~qc79G`s9?f`>o$*FQC;5+=+Q^IvFo5%YmBOLf;1Tinajj$N!h)oB!6W=G
zx&F=yUCd*`5oA(I5vifBPrw547Ans4m)e+=U*|z0su-H0t};g15;?g?W#WAE;Tm$#
z)^{@1!_12PxPZ&F@A7xWh5lD_?Jo^s2SW4x08~moGa4Ag%hJd><MbzwAmJsmFAE=^
znV6T^)GQ!vA(UyT>z8tB$LTsMJdp(@OwiN3R{tavh;f<rD1VRG`!42r7_h+xUR<Vq
zo4;cj#!VJ7J%{~WePb<}sFz>0@FlbAH{05CsMZF=UJMztD!zwrn82Ku-h`c)_?-kJ
za|<kqD-=QxGQ^&Gm8$l-$VUT)lDLbM+ss<;HBg<)v@C+hFnKow9;!D(_$~Se2&*jl
z?wrX7C<Kv_T~)QPb6*)l(OHjb9Ly7akX;6(V-`4Y4+q@YUCjZj#u63i0mJDxP)B@V
zVo(J)JSB(IDlq7aAS1$6AmM-rF4O*szsC`23=6BkxYUm@dz8QB;S3QTt~Sob+3il#
zfZlOCm)#U7YDvaSPorUqzd#Bw8RHz(7Z?M=gjuE+D=}kcsvjDJ18{bYLwNLCg3{bF
z)q>WLuA+6^cHg=3D2iKpZmA`7rWwbzd+CsTnO!J|Zt>t+LZUgnSp?`cFQvsUjkyvw
zNF<$b*o6&1&OFf9_%fR+VjX~Bs+NAQ2Ay^G`qmq&ahWz%sH;JEFJvD5tYe~Z|24t(
z61pzNK02XBDp(^2T+fj+%3lRr(vKJp4$Kc9Ebi+uJ@Xn_86$`y1B`M_1U*-k&)o_~
zCcIc>7)YoAGM8!J<?nIxea|lM8p6p)njU}rnan^SL56k;^&vuMC1<Xl=dAW~<5uZ&
zS$M^ESj1LL3Qm2El^1YoOgW2-b@VCBb0K^_*Qezqv9&1oLMLKEa6&0A)4s>w-WgnB
zwKN_(({rnaum<TxOx0yn%b}KH-wWi?wZG%yC{A!o26r4HU^kxN?bz!lMPLCH4hncE
z6ii^qytKdO??EHq$5fAXV==4^U$egfAue=5?FWaXPXMDtDha5hmLEhTmgCf93zcOT
zi6ML6h!%BgL@v76S{)D^Ca60X?H?R8@BPg4SWSi(muYz#gNz}!5|DUEyX!d*C|JJ5
zGQGG28tfXgXM8=cZ)QqkpmY26g$!yuIKP>Br(1-0f!<ynx<GR#V3?XFBreqd;XnL`
zj`kXPY{ODP={JYx>0e@O3N`bl%j*pw!2YbKM3%03650_=U_=lz<~^A>XzOLzWOLsY
z7V1OzApXV7^+`>3BFLdDuog%vpNJW}ZT>t7QBZ;#dTrpeVYnFG(HblbPSA6i_67c~
zI4#ad;xHv43?fuNcbr0(+l~dKL$DNll3oQ%etU!Yy<ifcAd)?6Bf@nM!8wca_3y}0
zkfHo&Bo+YaF+1<T#A8~q=}G$>e=qKkRL@?#a$bB*u+yENzmWAr22G=544;b3uncu|
z3qSzIyfl{~^`%RXoqpB3F(S-rYJSG(qQ!s0u*VKepD6(#e+vP%BQmIKxx_Vto*6(N
z#1a%YQOo=N%b>bW?P-`?2{(5&XSH}R{X1nKZnaX~Kxy*K5^ycEKLCM8p_pVFqV-eL
z5P|zDyHErmI-ZDpWWSODbDfl67`*^k)eb>LipElMAvMKphF@2|yhdVg*|Wi7PUb>m
zu1=dB!UZ)lkqL=16k=Xl;6rb^lCc<;`)6#Y>i`TNgpJ1t+~A&1TpW5$b%Bi9_$Ur(
zK7Y|;`0ql2LaFQk%g{3rFfV`&1?vJ}MHdW^CVy$%uY&iWMD}6K4oaMvl$N_RW>-eI
zYgEW_n?J;C5BJ}ujfFd$`SUT3iPQabP4^y)8!NneuXq+x=z9UBvfTUham;YKfRe6P
zHO_{SYfkSgy;D@5rL=dr!-T#p&csSEBf^-P=7SV>1YmVbwpgCK=VJ(b7-|oJ3SbBe
zTmF1x?8l!Ea$OI6q8iH3Yp^@fz<@aBSCFjz>V=sjG=Zg|KzPQO6KJ+z1A>ep@+@<W
zEaXo$A2eX<TC{Exn8;<?qx>}qkki{DA4nVlK~QF3=t@iEHIV3;4KMic60B0=7rkOU
z!ebD6^xEWggO>duKxPaMj>WOCO1!|T=`wc$EpxFsGEia!-?&zXv97M(be-P1U&{kw
ztYe~}xsIW_A-m#+ngN18GZr7E&%=~Rpbc|f5Y~QFeu(SyG%&iJ>MZ2PAW<%=3!pVH
zRFMEhR0WIzQC`uSdAq<>1dD@ms%}IyMn@H<bcQYqk}J3AmkdQXmuaS|H4t9;4P9op
zGp<77V~~1yXXD)IKDqby!_Xro3>Z6tqAhx<{SEZO0Hb|l70A+rgOUuW8GBJX@DW4d
zfZsJztqeJs&SfrRoy%w*X^;<H$pA1>$NIr;*K_R*eE65kv~Td&0>_ZMYadgkZXSk2
z06)&rAvS@So9Jf78F=2t?9=j|N8{)kosk|Cu*CpSH_mMtsK<D^3_T?QmkjEiu<*dR
z_y&X<#zq!G-KB3(#rWC;7Y776d;>BvHI)(7e4uZLd)tq*n49esqs&b~;>;Bc(tO-h
z?@;>2k*Si^ypr~6<iy#%NS6dDP&Ul@^Ii^ZqqjrfH;(*21?W#N0;4x8gOMdUIQWc7
z00BSo#;h}#fq^a&uYZXM`AvF!9}y0r!FFj%dhll`obehk#9^F~m9oJ(s&E$7<OE8F
zSyP_tWHk7;BWf}co9?5-TAPq&E|+Ov=dUX%6Z_YRvn7J6Q~7exJ_!Sl6=w7vDBdmq
z0cK1PdABD?Wv+eW%WR8T4IVf$4K3n)W?3iAJruQc)w_-iY$4Q3L4N!gCT~(Cl=F_|
zYcfS!q(uw(=aEcXf2A#-f`P}do13098&r<sr^FeoU63kCVjZ?e;j-OIf0HR9nktQv
zR_TJ#70)+p7m#-V$V!1*nJ0@I*v|aGfRw<nAPnE-e{RSfX8D<4&OTVUmQf~_M-jEH
z^rXF=zthl1=sInAvyc=1sHgI-UU}<P1B(`o)!a7<m}NB(MNCWi07Xn_;6lsjm^XrW
z5*StxGwC4(S}@@*>RhJz9IHM6(Erd55&c*JY!6uJJ-3I(4>1t_{IR}IIAhx<sux}I
zUmLI(^YW5WbZ}?%uZUJA01#tn&kGV9G1;Fp&-wSP1I+FMY7z1EzhW{ZOk);IcmYEu
zKFr|y;X5$mGVRO!b&;W~e=Au{u4UV2A@G>00zGN@m={{4`mABTd?O3}R6QFy)Q!97
zYtk?B!V-ZoFVT~r=EJb#f_#jWM@b;uzCup(e+w_xG6w{p=TJ&(spfsLz}-?lXUA=y
z=co+J59WfQQ8ub}xgZ(j`5)5~b;I<8FA4-Q^hSGTFJEaI^@1dj4W!Wvv;jJLG30G|
zZ}L3L%C;dzg#^8)PfP3@7}7%gF{yT2qEhVxzTSW=F#=9_(MLs=Rb~Zjv}%Us19;>!
z=NHAPU!JozK=8p}32YB?7{0hOEMZ;n{T?G{Ix;lFm&8UZhWt-}PQP7b;dVrp%9&SP
z?rVESPuc_gb>YrRfE+v1SMd3`Z-s9_Bej@0s-S8Zc)wGzJuAY=(3RsL#meWgPA9Ae
z=C>ItsC$K?XZdo=s8>UV4;{CCila6huTn1-shHDxT`B0RFz~3Vxq-7WQa+7#<?h$f
z-44Cu>kc}CjLWnxyeOQ1>&aQNYsqp`(HjfU-#|O}7oeS+h|<nLJkt&w6K(@`ZfS*`
zq1T8V9_jc!CdjuZJ+ud)quM#(FE0_~Y)0W*@VuT8!{J^O@@`K>7>$ktf**4T+Zc+D
z=$blkO^IKJjkf?o(ux#P^yUX01{b>}sI>i25zW|rl>WLu^8e@7GlT4it|LO7avL~o
zuFg{_t8Y6O)(|}QV})<y^<jsIL;$~EpF(M+t@lVbf@kLc=gjty$+l0Q?<mqQ`t&D^
zi0}+ncBoJ=^nglHJY4*)E5hZ^XEi%ce6B2s3%JM5lM~fz2QOFHsvoaaH9{|2hHWQq
zi%_RWi_doAWD7bu@uyJif%Rm-LK7&^13NG&(GTR=;q6AtEXZf5{Ty}Mn|<qe(VHr`
zgW4w<dNNR>BMK6*=!oTc4~Dq;11WdoXItX8p!q9~(Ye@9PukP?TU*7+xqrjR7LM*F
zy!}Q;GW|dQ=l`640fM^kEP}2BS4vT>Z*z=EbYA9eI=W<CiZD{>dF*F)qrygA%W4R4
zc?9!2zw<l!G%NF%{Vu5n4__(jI}msj<awtS4PZs*WB=?3mmQg$fy9xYbF9b_V)og-
z$l+>suV-}Z+urwZeH|I!W&Zc%f%GVa=qkFumKADXTSiW~s!}6nLzUK}{x56tl)tTA
z?s=itq{EU_>-$PtaBdw4X*wVKdq;o|+xvNe-G9rmix-x}mxE?<wV?PM?y#W0X0k^C
zowztTQx$4CY(&eMzhq>6CNY4fGqq*bkA!fqqksSSn&N)I{2wX~y`mqn1CI185x@8~
zBUF_N!{|q<Wo;OOE%$)u`aMP7n9vX73B_znkBU<_*KILJQWPOTU_<9)y%;n%g0d>W
z-${ci^hC=6#m_Nh=z$^=RB*!73XtdTVhZaII^`PjUDA75hdf^w<6Pir(eQ|;J}gaV
zLUk)9RoEOhChOtmq_vNMnS_K}>I_iXJ8*qezDQbP3{~*$P^f@RdXZ4)Bf+UWk&xhZ
z<PTtb7peaz-F?;q{tP5Q+Dm|(&Db@RaAkAZsHE&jEwxhW7Z@^vAP`Ze|3(3<D%MjD
zux^g%l+SbZUGjW>kjLJfb6)@(SRWGZLeh%~It*$t@YuIFE}A0J1}WkKB(kc!Sz?XV
z;9kbimFG$~Q9PCP+<>P`LaB%Wu=ym{GC!(}B7Y&xc|ah9&c}>HR|8ayX-hpxyfPo3
z`<Uhiy-0MvP^`5MBKvHg;qVUdtPqzWzD*k|P&rS!IXnM48tHL^c1%UIDl5FlC@Y5a
zjTMWlBevAdt_8yQp!@>;@@x+?(RQ@rwCu$Ou*$n6_>{Tk`gK&&gsyd<pt6380Bs82
z)MPGEB9@>txzS+}m5Bw_P5MKHB1wblb|_k1w8Wy6bin9-ugHbHqUoqIF46xW8Bmm(
z^$tndU=a}%42un%@q*jPURpeJnFXOsXx+!zZR@8C@QfeN&3O(pj7&o-kX=ugUqSwl
z9JQXua&W{WWrBN=GQlmFn>a=S9xq)Ty=7<8@|`_^jcXXQ&!+zVS7N53<F*ygW{$~}
zH$Jf7t6X1GC9TMg4hckn%APypV@$AoFaC%z#1zUuOX0j!@d10+yY&e=x^bDNu^TIc
z+kEVU2LfUbER#O$Kk|U$+KT-Xl77Qc@s)2=VOlK_A7Z2){o}FN!x`a*gprKSF=;`a
z&g|*dqRUMCY^hmXtM0T=6R+|!He(FJeO_=)mE8Z^>K<US<PIYzOES2G+4Cg9L)fJy
z7JqscGq5Vz4$Ah~ToBfKhI4aXz|h@nvETyh)OybFJ*tZQZTE11EyTb(5xe_8*)c>Q
zpuLbWmxW|Z?e)HAA9}_%qMpk%SLRc(?UFY*=lM>602R>Xn_fL0ofPH@ZBt6+d8{I6
zOQUx?f(3x@$SJ8_!hmrRlWKBq0hidDYjM1T*AhlpJ=!HQ8fDn>@~^J13>2}!vjujr
z<!(X*d9Zfy4;)0s1#^5%Xuj(W%RrU2Z`+NPnqzlmAH9h;HJ*&RiLl(gajqaasDeyY
z-e!p%jX*Riqqu`_;ixYwR}7!Na%xtWRx#z2{@J6UKO}JuR>{m}QJYN3J^Nt7VFl!;
z_O@qo)PBUb8|#aL1??Lc^bUY@K*_uv30{c2PpU;xP0%a!WRGbLh%1)Zx7(KD4&`F(
z#yaVw>nkGFv;0G|rWXgF;@piW6tAedpFwjF6;Rgk8;dEw$G~_`#^qy=R$pv`!akcu
zTz<dgs59^8Ama$^Yp5cjwcPexjxxq~J-wMLMa14;<N#+p9!s#4!Phb{Z(A%qE}aP{
zrg52O1>0C7MQCB>tdyL^ilaJ6=AvLry&#fMRC4PQFRKoj%=H-~{1Eo2$>8f5Fowi_
z$a0@HZ<XQ{EsP(EFEbMpE=<>xEu{;Q&Gn$Fs&DL`zCqiAVznfO<Mjc7H!<o?FWR3n
z;2U;x88=k2$mgRRR->`h*81*nBAST8ifiZ~<%|47dof2@;Xj8pClFUTk$ni4DCD3q
zh65gwgc-O?E3`&VZYnIZ>8mKW|0DJjj+?7S(I%J|JC150Yu7K^OF1fciIH&Nj4oBC
ztjPE98gF0%dN){fAoUT3H#RA5`P>@EprO+0ld00`p3{2IMPmE6P)Qaq^AGJ6993kA
z5Uul*HW-5TsP`L&@4*WEL>SbJTH_gZwVadwYxV(s$}qgj-Snn{Z5Y>O%W!%8x+tJ)
z<*Vd`e&~z_zR!PVm3NnCR^{z>CKE~TVr<>Ez=F%=ec_Hips>%Dr%}*g;?A1#N@hJy
zMHKLA|Il8;;jRR~_#bNT!4YVSCb;G&jCps#Fe)<Qj8m$;m*MgD1qynW4+!x37Lz;4
z?AN5LZpWObW{*=N$u-F4`d5)`ahPrUF2g6=YR0J;zK`LGC+!$&*B2RStNRX=uJz;g
zL&4T_XLo1LQ@_WrmmFO3ZLW)_rN;k%tM~BK5;rc9?$znQJcnzgBJ_Nsyiq!Q6ycku
znrx(UnP%zOSXBiGSx$L?<8kUD`Rn{cbCFc%CA<7ggG(hv@g~~1Z6K~x^2waJLR0z~
z#%@1iP#-xq(&k&9JdlA3c3G^TU{>0P_`7rW=D17WZ?BV_+;-aJ=nk8l`#rRNGVWRG
zK7CjPs#@?QUM`u(JjR`R+WTe8hW)t*Dkppvsx>sB<}&Sr{5_TVoi#bj<8P^y(lpK{
zzq17M&F?2K-o(B54Gjs!cQbFpoz6ZhMxSEd#{S7Fs;_S`e+=0;6<{})@50NQnI}?3
z*yZQ?JinFWY$13M1uY{L%<B9U5O|PzmJwek$Z<BaoPAnYh}*<3k-?az4h6}bqM2Bc
zeS^s6iq14~s6sYE;pWV7m)R#$O?m#le`s&#IOF{=iu-@9@_ABXpnqfFdD!lkjM9(#
z&N~?oIC;LfB(KIB;|j<g8NAX55BqFx-6osucKnIaM{e)%56$?ME4hd&Sod=@^tPK0
zRR21|`Zsxsjgt7$qcvtBvd{Jb4)YWr<H4PE|Bsn;i9Gqd$3L_fjO<mlAC_uM^0`3F
zTddwB;)RV}R^Euj)@h9!v1U}!nsT;vJs0g1msr~!4mUDr32NU~MF3b-{5E47Nw0&~
zJ$ofHUHK;E#yaJ43(E?5P27_w6OFW)IeyFo{wPDZq8;x2Dhi&yGS47^=Udk6-4Syt
z|KzYf^&I@1^$BeT7@3*(3=Rr(x?F6Z;HZzxz!O!`6Kg+#o=?V&T+s)(=!tVI4El?a
zTg>PUnoC4=Z}|NVqnF&ynV!crU~)3ov+q-Ee3a&Wh4HdZEExYCQ#ZCO2HF}5SjlDD
z9sZhQW6|r-XrsTuoP*Nk>BhtMRgV7$2my2f9M>+GzePkr_9mD2tf({9o1Hn;?$~GC
zO9TOzY1x}WgAYgLPYbxSzQu7(9VGJ(%J!EWZh2gF-0)e6Q+*NY<Kw&~f=;Q*&OuM<
zgjcjv&FHt!f+V6HH&p=~(m7*s0IqwTT~mHsudl0x#*ghqT#OG;i@gPAD>PKt8w@bs
zsj<!6w|2&+pt47=4U6y}GO{`oFq$gW8nc;A@6{52)r57<a}L0#s3_4|=!myogCgx7
zhvE(-4#h<)FWiSX*EUOh34QZ(Ll=}T!y+o3E$hKWS!xEp%bex!jHHh~47>zV&3%`B
zypEv*hVQo@3=%r$(iY~+-+vx2c>OYtR(g2g=DIyT;;kr_*np`t1LrHmnbKoZ+~ATM
zlmTFHhWnsi9Y*wG7j`DgCxPN1ZMRw4%NbBw$3@i-w{vUA4-h8gnLbJDU0_R@A@mWB
z$j3|NcXxN~rRi;nI}8;BZ4*6bO;E@7yIjnL@e%q)+Xg5sIC==PQESfLN{*4gQ^w{~
z6N4fa!NM_b1J>K7frr8@Rmr4niN!%K-EQ`aI=_l@a~H1)3U80vQ$rn*D0Qpj{3Q2$
z7vsckvyqAvz-f&)pn5wpB!pbfipC<+?y)^%yx{e#IX2egH`+Y}CU!KTHPYK&0_iNb
zw&gv&$?AL`=QUDlq4$Qm3xB^?<MdQx!lt8e6z2_ybH+BydJTg*Cs9%bA~$NtwK@&S
zEt?p2F}u!J?91fcX0YuTJo8$+JN<@GGeD+1i-p6=%Z#+yYq#kw(ilQPsDQk`MmlH;
zV!33s0ZvpnkaEwKWN5#~>#^Y0Mp$`m^XNN;)nmd_?2pzC>rYF0R==ooc;mORZrx_g
z_=(wYF0-Xjq(f7p?=G3r%pxC<2be`<1-&6VrpU0=q_$iXshvZ)Q;=<jG=V{-YFjEc
zI`46+BjZ|=6`W~fCJd3bh*|#n_SV>5$3-t+L*REMjih1ws-&GobUV?XHwxpm^2CAJ
zBX7&dj{kalvYGB1aLy8w+JF<^)dJuGyKaL99w{x%8D-7OD5QBK#dEjU^I|Pu6(4Sb
z#`qW_J8QrorrGq_UE0`E3A<udIAre#kR1-3o5Y6}r*$+zafvV)jm;Xrv81^>9=nHf
z<krp~b;73v<`M^IQcLAgitNVLD@zrqVbOd6)7C|O@Rf1#BskAf{Z?SPU=%Q$na8vX
zc>)nCr<)(ECx<F&WX#f>Kektl7rf4rjK+z);^mu)Yvh$o$;8m*lP;mVN|XUSw%%|=
zRDon&83S;L5T&Wb80S9FS+N{>`+Ao5unf>u=!Q!3OuEMf2yYS>GO!t!5{*)QR)f?H
zhi!=*d1FrT&V4IIYEnj#DgKj4e^J4)vCXb?*SN8)x!`Wa;yay=4H$D4_(&e5_}ZP-
z*dF0x4T!1{PqPnpW8g@mn6uUMDq*e}WAnN7%D6BOk6um(g_G4|Oq6kVGY7A4cx;xI
zn+ax5esK@I6PJ5AEkL6kzMg@68+{uq*%&hC;pn}hjF04-+nli}OwNR^)+yKFu|NV>
zQz+?t)x0XdZ#QY3|A6y|8K`>GxEG6amYs5Y!)>=qMxSy=+tJ7Q*9z+Ep(JjlvpQ~w
znBfcby`bTOUnycE1~U(suZnZm__3KNz5TR2u~_{*T_b`+HEj1Vv0{A0Z{+K1cQ9j1
z?xQ0N6X;7_`Xy5aal_V1S|Pr4#oUPZ*Gf3FtQhFqAxKF4EyJVG#h|89(K+q}c4xty
z?Tov7HN$s5mvl*ian(_pQ)fx`kFDQH>^@U}-C`r=Z@8}9oy(x;U1P0uYcaN76Na+n
zdN2wD%ce*+;Er+8I<pN3dTyWHuZ>NK#eN78ACT?{P5h$7D0F|NoDIw!D<R*J!P2NK
zi&UT-tlP@rt@UyCmeqK@c_kQr=rw%9f4E<0R~^VJm)x;sliX#8^ItRAB8KguBv!|3
zt&477#`b<L<UNGM_o_>aCUME$CG!|ORgL|Xah8`SGF>oWKn#6nU*2BgZmo(lib^ub
zYF}*(d|OEc>{OQo^5~N25;(RyTnxoWIKM^LD(So!Dd~DkqPvN)_2cQT4v#XpBP=<L
zK4Usvx+KoBJ|OR5z50uNSzD<qEE}uhjNk|5$^(>+LQ?n}c0rveyojKwdez$KIx#k-
zTr5Ukw|9x@06M7Kt~v#}zHKZ=*1p#Mv3-b(F)Xm}PfuN2vaY)@x_poA!(6PJ37qy-
zu6In6^sJbJz$HuKY}-ea)Q^Q&vmzN5&)b){wX0_KS+jnR?W0^w)v29GTIY@ucRh<*
zG(Hs4ZI`T)R&3vzbD7Uubm;wJ1rc>O>F>Ak=t`R(iK)~WMrcAXnkas#7+s#jOY<>$
z5isGxpVSzHbpVO?nGay&@^TYJTu5E_Acw!-sjCn1G3Ql``4A8M5lC?Q9Ylwjj4y=6
z?*$c|ZhA^_<+L`XkL_bzXcqi#!OsnS3Jsi8*n8?NqYR`mD`b0Guyv@vF-YvggBMSA
zI_G;8!8vRZkuRztvhSlYr86HSy>ZgR{dKmd)wNu_Gp(<-<;Nc|8xVAwYUR7FQgW>#
z5TFY$xL%v=$M$h9y7nkGixx~3AI2NbHN0mWEK=_7>T{83)&^V!!!4*PbNxvG*H?OW
z5O{it&O6-9lvwqj0cBA08j}@nm=F*sMrYbT_lr7zf^)AXp52Yq0|6Q6Ann`=?wtn#
z27w=<go@zG_e0U<DhT3NtboLCQTidWMsIX~zs`>mB0(~ZR!rXd#{``(07)bGF%5O-
zO#7#PQRh!`mP0Bu3S?6WU~C|ax1hTNz&!*BQxE~2nJ`9#J5Jx9C(7d)rX0b*jZ5c0
zwxEdNYDg`-l1cE~0w?@2)dJzo4j3dqA%+!1aRz)6LfG7CnQIhOQL7qL5SsRn{i4pF
z;ym8{!@KV@<Q{))z;*Rj(#oQlgeyU7BA-6T9ItB)y<26Z?TLLz@2E&QfUcuP`hZln
z-#`*($Hxr^_%}PJJEmF;pQ}Y*L6RSE=+~ab$l40QEm+j9I=3MBLGMID!$e<G@9u!n
zQV1;)+9y~11Z8x;1R71RIv@0)PEKQwiiqX~bSm^6OZ`JqXbmmcnJr#h%<CiWQ3%$k
zM`zkE`$e4{s2X=m!OTt(g+w>$7WG7SYRv#j7I;FEQB(VD4TW0FEJbLUs97^Ux?cf}
zB`>xiH;WVJM9~H=XY<e}*V8;#<N5?SdZ<tNiT<=-^ou%whO-|;Nckc5O@oRZsm<R@
zV5yryl->cyt@p+7#c1Yd<8l=~ET3&>@Ilord{X^zM&{z$%x&Jrnr2qqyg-}I9mdE|
zqJ1PXw$F00mS~uo3DzK^&WP9nDAVBVex`CT`vTO&H*}(d4~qRnJ}P$%bCowoT6@m1
z22Lru2-ccdN`_}!%DGV*Ee$FleJ?!RXNC_>-kSR&cutXnQV`Ogww9Z^{v7A=E(){s
zwm{v&IG2qArngqiDB9&eqQZZ|#stN~k4OtT#8pxX^A-=iPyv?!bGE`k5^xZ+W$;TB
zviKsyyd1H6>7)a?EMZH3+K>80oyDUJ*8cLi4(|NAo7aoL4MehC$4BwmIWpcQTj<U+
zTK%uY0+nMN)WM>!!EiKI=Zgi9CA$E8T!dcEaNQoDLzZ{W1QPveKkOHE{ygXDZJqcO
z`yh<<mjY2wEYucGy{#D`ElumoIPiU$6R!#lsM^<oW0*g}@jQ6d=dSUDVBk=_ID+D$
zX=y`AW&}1;lg_jmkuT7zda1Lh@i5n8-Dp@igUrICKJ-Rc?%^gh)F;tZM%_RQNMBKc
zI^^ftG44OKbHr@Kz}3s!KBQ8EOFcc&n;QjqZ)o{KbmUzEz!#N8$OHhGKIAp*^4-#9
z0ih=?YtRxOqv76VLTENvWD<s`1`OBn{qrGjY({&HWle!&|EeFHv%M*RMfS8q%<kn9
zd%vWb*q2(cc=b(oVZnyUZbhmXwaM96fC4<ui{Oi8Gequ{9)`)5#Rjt>SQekYlGqzG
ziG=o)`YKX8d-~JnK>sqm-bBQ1;sOD0unQ$N?6&He>`MS-F+F#@z?{AkH^@J<6qM(@
zGrla}jf0_d--9-wj#|`}>{5zPmJGUadt8EB8DSFiX+Q24bruX6Gf_Om<@MDCG=34e
zYFg7}6##*abrBv34a9kuKki`QMP?yqPUy=*!w?E6g!xu@BsqY=KtKCR<O`7Adnvsz
z?HaHVh}xi^_KP}yg>waZ5JA9$bh(`p4l|}HOaQMZEKZMm+?vRXIlXz}Dh|6UJ=$hq
z3kDr1(8cb$PO$`VtGk|J*cfQ*ttC}61_kI$`&qxJ^H({Kk2+rKHG*cl$M=dVJU+<O
zX&b%3jL(y;Z7=LLeGV@1F}WUlq2BU&?8SSIm1y@tZLiRtD&SC+P=legU-XMQe~q*B
zsr#sLzt_M&TkUT7rlK7`_xgwFsz(`}c1*x6S&G_pOqxxPE6JOf%KVI%I}QH}t}eH$
zpM|gOb}gKWAn`gx`QI^NY=6qd_>=;fx9O@ak*1$zRFh^oM|Lj(BNo<<k22A0!+NYr
zF<kVBAsCch2oqG;z4>^<AgBSNV}$F>zYGY*#??U3#ywCLO>h_r;^LXM=pXt;odss7
zv=w9boqGL_$LV{2P5rPF2%fHi6(^@R0zb}#^T<@ygtKRu3P)g}Sf0TG+l+ILshqR6
zLw6xQlq)ZMl(`KlG$Sk>mz^O?ac*PgM_OM$8Fx$OFcpxTx%NApbRl}kvNxS*&VmI+
zKJ5_7V{f~+afn?glU)s`*rM+g$Jy)pl$j<J;=5=dD;$W-aIDOV^_5ADgqqa%%3+-L
zn|@K}KjR!=#Q^2pTUst3;QI7k#wL1z7uEzY+1NF>**q;vL1+kvBn}Ef8xJEaNO2Al
z*cE*No-U`j*jC2Pbt4=aJCFc`P{4$%;Dh(JFj1}gsSb$Ru=)HAI@bD$!LF=srrsc;
za2|r<H4I&NA~U6p+tS+Bquu~Y$AY<0T3Lyh;l7)>2?G?3DdM66fI^)4v|L?Ct!9%m
z#DS-RA+PO&(sJha;yibNrQBg!9;6qqdx6#{(8<1|Lic$l(7*l)ny+nt&SCo8Diz`B
zEo|+?9;El;8xyysfwf+#vJ{!nVkDi6Wd#TD&L@pD)ez;6j>b(TIzK{+ZOq<+B5?-C
zyomw~&Cy#A*}!FO=GNG>?qUtidIeaoQzsNcttQGw<DzLs1}--83~JT9d!a&sR>UQ<
zqTlDDCgKutTUHqwHQD!(G^Bmr{8}f^vtep7JvHdvM~VQxhvSed*-LPpz&SVYTkvYD
zJm+x~XAsI(4Li7UFQi29N-q+KZ2f*_9;}T+^{ri)r2z{qX`t9uV=Hi>S=(&zHyIh<
z`OBRhF1cUU2;^nc>ci=&$hpCPDKAWF#l__G-;^*ARY4E)P}ps|A!LTrY3<LMAtSMQ
z7*qa$21LwScqG@MgE~8H2As~!wTbH39_3=(_(nRS0MUO*aIpqHC;=Wy;THU2kZWYm
zwu(~g>{h|RyGp6>z@%GZjV9=m?L2lN`9AgHjPNDg;d=UiK9g!UnKj4T<54L*;N2=A
zIX84@ncrpMjF!=aKM!$K8EBloN7onwk5>uVnhJoR8*;}$4CL2umJ}I4c_TunTX;qX
zh`##-X!TbbyM#=%0*Y+r>2k+dKTgJ^H0@r$bBHGBV@{7~k@Vkj){rqHS}9vb9+`2!
z^A1GRx19$~KjCtRs7Qj#U${EGok)M$Tll-N<XB=*8cvrX7Y$3RqZ-CyXriHctpZq>
z2Q!=zxoEch7r05cA3FOoyhd%D$91#wJeP%7M^)0H_2J<ZQpkp11RTo!BQP=7#&S?g
z^$!?2qMLejdq=;hv-(tsM5QHgvNBV1UZe{bx8R>9-MnK%#bhM2p9sCo>M8a`pt@cl
zdZYt@w^BYt8e_d6J1rm4iqtW>U4Ih!^*V7S1=r2c^NiN_>07t!xa=m^B{V-neW6bH
zIA=qv`SVZA$`31e?)KNb=x4|NxFwfZoGljdvvg#xN{)90OD1)B%#`7e=`gmy0^@j}
zto<MV<A1nqxwt-Q%EgWL<-)8EV^dhc(Z6pA>OI3Jd{g>mH@e6`sbOx@%o5C3C16Sa
zm@&rT?j|r`HiSqsVA`{1w0<Nt+;oGr4UHh|+3$69du!HH`mA@y_G2#Ydg@c3>V9*G
zY-3(|&#Fjw{X2zpMo~duS4|js`r3#cn+Et)hZDH;^Qph*(k)>Xu-?+;13w^{Yo_Oa
z<1(k<^w#UL9Sg%}ZPrEH#s@|q`8RaA`NAmDi#Vs8yaDUxpQ9<z`v!+u9q(ddy*!i5
z#SD0afAgE)3<OoQtS$)8ZjVZU<?lfDlX!{o=tF>%=`Vugl-IMspG+MJG>`xH|Nh@F
zz6E3N)(kLv;f@j?p!dQ<G8e2%_ZpL+{6CbjPf!i{n97#sQLeq|K_;esLYc^3>?>ds
zSp>uoSsR;)QHZY%uCe(r*iA$PBhaTEC&6Dp@cA<i%jrgLp7sVKb6=Y^WBWN5cm3D@
z`d@Fsx>nb(&{c115>6ohr;5dH^g^Y7mg4dhnjdKkr$3{`3A8{i(l~SV1v*a;P5)9g
zmSU}+y=%tR?zneWbVkmaqQ>TP=5bYo_h;!ky;tNdquDzC`pz?%VO7cpadH3pi#>pG
zzhH1C5{c|@$)2hPG|!U=WyiE?oXmNK&)Tf%-!kgnB(8_wC%M*-Kn=BbKpQFj0F&nI
z&&|@Zu403r*KguuyGl-3XyBd-Bo}9G*7BE(_8O~oALH7k6yEOp&wg`|ZXj~+Q<ZWk
z6S+wMc35EJ*O}wvCUAub-NQZB8J>$h%s#a_*#3^8>8Z|#=rYDY6$$z2RUlo$8L57l
z)Fi{xYJ>)^CurYqO8Btr*-SdVXAq{bxm*|ow`7iQka~wM{fksO{OQ`Wn(|?XQ3@?l
z%Hh)n!>4nuak-fV665u{8x6y4Y+2A27pSDG%=;NXJsx%Uf2N@;VrhtF-n$4J>X_2f
z_q}$~>AYfGxl5VV^Bn?BPIe&X54o;KGW|XYy$C^mKfJ_)poebG@IeW1n*3&IKgoTT
z*B5o$+1$)-{Z&y%E4aM-q*d-sa4Z_2l8<MF5@(}N@8r-&5d8tL=)>IH+AuN8$+=$4
zJZlOcn{ODbYX?wemkTR!$z91(w2*1hs&qiIV+GD9ibI~Xf|r~8Q)xJ$()>QouVo0M
zTgpH_@%1*kEbH_}D+&i@GT{Bo=_l#)e*&1p*hi(n4(Ov8I_K^KAkS(;p`~}BIBQ#a
z>0($#HCS==2sEXd;$3O+7m~9DWC9XUDXKxql}T!?@u{FXZIsrnoz?6)KkPj&QgD3u
zq>Hz1^;7__TJ9>6(7~}Ai!L8v!1d!6U1rLLKFm04F6CqU4HxUp8m7OAYfVC*2*P;4
zO{GA<(96{wX+;e2Pg-06Knt+*Mu-7spEX5~ExQczp*KQrV6eVW;ay4vtnpP<APPc&
zy4G`ZTG|^+{uwC4<NDE0_;(cSYyB~|X-Rn<<HiaGMjyhMD~U><^b)c_ky?dbT&Yrh
zZvq7;t@d|qGdht2$HKcTyI<*xyCNM2m2p|)c?}pS;f=Qd(vcfTdh=msx^BW^Du>FN
zn+~1Z-E_ba{<q-o`>4N1&kf9Mktcng!HCul`tXvDcQWA4%l8EoazGbVo?FboABi3~
z-^T<K3=3oYNfVd#lO}mqKsw<xtEZ(a?@*xGsLj@PP8GIPX}YKa7_VfWPC~2P5(5q|
zF<Ay~dro_TXtVgqs%diDb=Uv=pZ_xsg&UXrm2_@fx}<_*`z^@oVfHJab6bNN5r^AK
z7vE(V6m>p|{yx2~SAh{l8kr)Rwb1&xJ(Z5G6xJ?Syp*nr+-N~>Wz)+YSv;2!A|6qB
zBuasD1v$IfEUH=LU2_N5s-aBm5)U<ciCqZeWrSfb6QqsTmQXez1Ow0>I<at5Iq%Om
zFqtdAyi=%=+KS+y+$;~`mLPC)A{nJBd|dBc?W|7zX8AqL#D(;3gvGq4l3(^4%<b^(
zjW_&aD4xJjXWq|I;ACX+gmu}~=Cz`7a+?I*w!~Th#=VLUeMjae<Wp#IHR$-1t&o4c
z^w@fTbuppE_1XJ?xjAXvkKm!~>ldiQAEuZiVU`@AE+|x2Wn3U5_N1e;ySX~fO|ga=
zklj*Y&Rlgb{%Xt0X^AUBnRyvmppFBe1Lgc9_>-oiysZdx11`Elgy{ggU^a-5k-1g5
z*lJ<1wOj<exzx`U6cTiEF`iqt*UargoQ)gKN>zoO^Rpv7c_32gt}S2Q(n9f?JPO!a
z6?Nx&e~REPd4nmkcsobMTbmEEua+2h_uFV22|RC^WcHjQ)3;AT)g^x5*PdPc=Tz6X
z5!_|O>SA@z|6ri&8dXf$Luk2Y<q;i!HA^TNUw%_?K{!uK7ad(Tta`cOrK^HZg0XFt
z@T4o~K$Uwd&ZV<kpaOOfWfgZRk0_J^n^8o+)>(6V9vv2#zGo`Z_tRZKwj9{^I)Xm=
zUofP|h7T2CAkNL&r<QK;8@dnBS9_T3kr;NBa5#6oq{5E3%&@t=fDQvJtLF){579k<
zX8gJWxDm>GZp_>$a56Wu=k`K6^lL9vOE``{)E`SwBQwD$KID4GEhh2YBo+UE^Kbr5
z@JC)D{jnM$_*+hZ1*Z04X*BaRmgg1~QLTlPil1bjK&BfD=r$(Zne^@1Mun2O@iez~
zOSqwT32_xzg<v4vU0%G~sOK>?dFH9ja|;fowTvDA9EU+S7?^SmTiibF<7bAVGq)rR
z;VHuw@&zWf5Z6B?IF>lXvYili3plsHT>8q=aUZzUCNhN{&PKm!n|&09)J$iO&4is>
z%IPHGD@-y&;M5YbbHSE``F>4#+035ncoCarRT~OteT`!jLe6EcN+e~GK-f+-aTh%|
zv!=M+DmfMRUv~um0(^|YR*S4%#BY?^FDtGUbDP$8O)62}V5)o@H$-`3VEJmLb}EZ&
z;oN%1Qd+kAy>%HZ^^>TG2P$6pWW>jU=rD8h+U@MzvNbt}$0?6;3{o@_hxDhL%U&mr
z?e=(4uMM0L{sQ1zbi21p%3<M7{csmHw}c`_c4#{OZ4Lo4JTI>l56tYi&OT}4_}${z
zRca>|-J{MeVT>|S5Y!z1K8H0ReZ<2|<Mo3Zoi5U&AC|)C#J}cPIqtx>rCkR<oD!5L
zjz5$m1-WB&$K=;4eQGjQnThmAbeHoTah;q9bGgMai35nuu|n=A6UQAWx4e8eq;MgU
zqf0^dQyE_Dj$Xr`tw6lNc%j%En4i&Q9Jbaw5W2sgTUUt&=@(3BNV{*Qnfc^79)?Zg
zWi`5C2ejK0dw2YGUoax}OS`W|>@^$Qv?Hua18?6Bz%|z`pBT<?c_(%iTr(f<q52`N
z6|TG`s{?S&H9OoWWigN#fXnjV(-*8#n@>U<9owtM3ts;t$7&;)<**KExZGISH@BHj
zjKMG@M}t$kCs-<Hw;tomd0zj5G9{vo!L(MO5X%TC>?}b+rAavZUV*~LsBk#x&!%DS
zcdXmY%V<&yq^rf4G?vjdI(CwB-c}?y3Dl@C?Ahd{m_l^V5i`5LQw2naT{)DlW$|HE
zO#e<Df4LKf7G4BEyI}x3t0PX{7#bL?|H`!cOzPDS!zQJ7meFpfZp@oLA3==(p-qd)
zR=6=8f~Utl&%0FFt(h>;3>){jdnX}gcB!lFhryz3Q>%2FWTQ&`+MCg_b&cKtJNxQ#
zkvB>sS-+my{hWa{L<SXd7l7q)!qDG*JimCn;B|NrW-~OoE5hvb@8~u2%Bt3o4IYr!
zm~cC)y=YP`>{_)U6JE$X9X+c&lFz*S#Ka2Okl6AXV{pjdk)6}V7Iqsdh{#o{&AhU+
z80y5Jg6s}x=biG{<hZ4UA->f=($2iRLD*42b~{5ON4-BbIc{lT@!IMiX=h%&&FdOF
zU>9~BjDyEVcHWhbO?E!#Nn{P>$Ru}JNdFLtBkd-yUSfub53eiqGLgsjV=k;Ub27S!
zuX+tNYNXxdRZH_&zrqW>A_F`}9{qae<-XnT=g)e%fm-$2S*t93V+*gwe!wfdC?=O-
zzsai>H}z{d?lmp4;WhO4xM%D8*mAeixPV)a-R8v_D}%m0?s=0}myH3J+sw5(MX#Op
zX_FVWc&|PFM<4G%+-rmvJh{Bkd<rtQA8_HSK4TnS9e2ZOKhkdU>RcHX0eMxmh{`)_
z%FNE&*Bvr%2iGVYB#qj+_>V1j!)?)!G2yj=z(~8vt21dp+Pi6^Z1&7$ymH^a5dq%o
zuSIc)n30zcfW{VH8<`Ydok>F^kF=Y-I+})zc|AmKoQ;@~zF*ILzUTa|6hj4lgZ-`R
zP`#UvXZMz43$OOofJu0@jKQ}f?Iy3v$e#BmUT9nVx08*Vd^TvpYvSK>JkNZ-tLzOM
z!l!NC!O1QC*urarCcHLejI^7)COaO2V8sdBNBVxf$!7zn@GAae#cn>HpEh3b`a2w3
zScMlDrJ4=#F_Lfc>bfxC`!Gb|gKi^zzux3Se9yW676`1$>EFR>lULV~Az#Dm3aos#
zIJR$dv92&fB@eGFu-fF+SPgae5oF9O$f(06pAD?us;-UI&N8vdt7!(T!mH}bs7t@z
z;u9O7`7)FI+B#69CCu|lX#$Ge@dm(WUJO5taX)m%)z|3tK~Xpi$h0)iJ{B0;UvaT|
z23n<B@r{)FK;4h%s*ghfhmVF&-$(okV-dF#=$qnNMO1Eax8SPSAHdI#=dlCK`%qYp
zxnEK#x)Q#Oa(SU`N-F5pQto-Iji#1xS9L(mJA<$Wt}i;PorvHIBIwwZ0_SK3?U^Ld
zQj$d~?{Y3;zF^lta#&2OMtmve9(aHnz3`0>NQ&<ibFXh89{_FwcYQaeAw;3p!ik49
zfML5Jh616F)D1sa=7YE>7ae*kI9gZ=?NFoU=T)=A04`e^$a!o=YKPjN@miKhiBw*8
zAV&oKE^SU}#W+W5?{bhna_`r()Ne81dzO(5%{ROzgd!G8R<U>^4{04|UerLmGDC|e
zczL01%7ec&Zk{6t-G9kzQy|yyk1h1KKQ1t)W0QAB9uUdj;0#Z2d8yafYQ(s*y?wl<
zKlhF`u2W5x*2i08>xf^Ed7rNdgy_Z`E-$oA-3XSZ`Ew}bmus-|v4uql$l>TG<4;A%
zRkA)MUIF14%)rqnF@9D9btLj*R$M84)c)DU9YJN29H3$Ddz?F<#^r^!X$_p(JFZPD
z4|qL~Yd!28t^3c+L$w!hU2mG-V%T?0(SXj9?|x?08~&oAbtu~(W`djCed_~;;ldf@
z=q%fe&u4OZ`iyJdK13K|92iC9y2J}+KM3IEg|-Dlf*<-zAwa{-!{I9LZB^dF$Fv>G
zjkhW)d2<#H5bH-36bwK1gxQ15CSc*xA94UC3&o`dgI1JdE>kbGn^VtUa=q>gq1`jw
zrv^s2E;0rZ3mC3y>SIjI{uZ3FU|6`~B4YjjFYioyY%QuN%;kdrM*jsSBt{db_;-K+
zHzq<<5<eJo9T)`_P@KhCG#YV2#0dcrWe6yu;?zECoxRT4U0rp$Y47EOzDegjHLcov
z?V;)&&S||@aTQu8CY=J28pG<{lzIhV#Owyv9Y5=e+<ZK8Omz)E=(1N=EO`zwbZu(g
zDeHx-Q6BM7i94$%W=uLq7rp-E6uaxFcYaH~Lx~9*NA-@|N9<G4+qrT`g|o~dx~vXr
zd4B=nx{FYcBv%?LKdSpx5>X0qwkP5NdWkFwY8qPvyXq^4?nnC>rh^u3AXm;SW5$sL
zL-ELQ8f_h9(o$lAD{XTl<Ta>)CvMUWeQ3sFQgAf0IV;zZGDNL5SIT~rtTI)$CXS`D
zUZtk0W1BfZ(oHcVFsj6iY=A(?!WV@sG&DZ5uH3_fVOJAIKf9!SS_D!r>D!~~DO0ta
z)@#BZ0Rm2u$xQ)ZreKr$Ou;7pJynQ8WTl9DmTa#BK2^ZV5GV2~=A-M8W9n5Q3WMw=
z3V-o?stB%xt#oc36CkQb`jI7*jJ>Tt7~+@V$|O*31V-8gXieEDxZELrXm<)_W-l2)
zw5nstX^I|UhKkdY<1j;1iKZ-{ZwH#Bau2U3A`Qbxr720-Zy={VF$7b#ZqQ_Uz%X_q
zsDes6M7KkWcZ{@~F3v}i?ZxB*`Rudk!!)J?Gg;)AOG#FAAo7xqP=`9EIAc-Abc8N2
zjExDbe(0dx>SGc;V%{4!r=cVCGj<g>Nhp#_i9~_UM_s4b17a~OZ51mZrrtrD$3jc;
zmS9Z<jO(*74pC5+km|=$5IUh5m?OtDI!;Of`dM~#5RR65yM^Rd8D1i{c_`D*cyK=`
zHwTd!7ddq0fLzAG^qw~^L47zpmJi0rHIEQM#ec;Gx>*KvFc*{)BG$KbvVP^3OCI3W
z{ThE$9-?>SLEnLr9gfYKuat)*^*v$XaIGalv0xtJBg&-miPC)^54Le)2suy7g}&)<
zp?S!{wMq@<ASntv*luuQ-%#AX${g4qu-&mdAY*i!#l|H=rLMZte=iq_RdlxuDHGe?
z9K4!wLA>na;ijh3pstIj;lR7S#00V6`H2Dh17&m2;@rNEoXUrX@&mjn4rY1Mn8`o%
zGe*{x<ObU93$GRj%yR>J%16vW#l5c7EiZB<XPbj?gg-?7&VlvLmLa+^hjT;L!{og-
z*P)ygF2X<KJs0DWyDzqX=(~I7qv_zB5E&8|DI?=-OJzRrAIpzz5AoB+hHJNzgOi_J
z4j<aLAXnek=rZ6C7ht-*xV;duA};LX$@R*c__(gs55-6H!$&MX_Lg5IJ|z1Y_=Gx(
zi0D=iQOB}XqMzHL)^H4UIdr6*vfG1)n4!e|DMuWYFDNCb=s8f?lpupCL1u?bWTI?$
zlB|sq(Nov17O|qMU4rUD4Vy+0gNA3o%q;{a$6JWxZiF)`lUpbV^mxj~E1}vH7nGa0
zx9lc&N$TKiF`|~z0_~GaeNk0EAO(@ShpzpKjawWnF8sxSrmTpmdtkAMV72Lu#Y{%X
ztc6Y+ppa!NAtk2pkQkj=E)-#1T1{<di5(HiqLiI4rye?>5fyQ6s&xw`VU`=ErU|Vi
zy+~vT3TUHKx5nc$9$3MpAEg3sgUw)6&CR5AlPif(S}j!(rYi8m**WzBlvdeIQ_9__
zRob|qw8b1yktNLT9w$JG@phx!26+<&VOOf^9)eY_=#Zlvy7P!nSq=L6(KFLCgR4;1
zhGeK_RC2x$np#YzxdBl$kHv*+*Q17moHmyeU9fRPYvjVXU2DZYinKe*2BVkV2*pgL
zCRRKlgzhQ(DVE9rx*Z|(xItGmy7yr+7{kXS$8vWDud5xZ6^6_Z;pKw4Ja<T^-7;5Z
zNI0mEevuDBP}CzNHO11sJu&bqO)o!W0XOQE`Zd504d^8c0231QT0nSXtPGzUEILy>
za?H5czPkg~UOuu4z?=D?Z0(TDV!_9rhec<R7J2X%c=Q!86&fccff}A3Rju+&{U#6H
zJssT0Wz<0~N{9LmqaCJtl>z7=k>p0cMJ8S;AK`_H(_Z8|>!Foh5Sb3`EHHyr<P1ya
zm|Bk#A&?0xFuFx$S;?XM<@EP&v{7Ys8ANdS=^|(BeWN;Z2GDbZ=L_q`liaIr<cr8N
z42<#8t7g;^pd_TE3a0OtTr9E>;X0<(l=gp=bs7JWrh0;mKOl(oub<R}SS!~_<MXwJ
zBO@r|Lz;fKNjwU=UbmpuLCG0M0A9~e>!EAWT)#f<tMcGtG(8;IJ<*oe?w)ApPDNzU
z*xgT8x-Z~>p`YJEKP{igyo$@aN)e!u)=NP#BzM>QAXMw=m@cpKI)T2WRZV&mTE8~8
z>DWrh{C*89wAEKqra_Z27-1#%eGQwB{yCTJT#QrQc>G_@MGR)KeMKC?oyTE=kw2BE
z9F|}Hug#-1aP=H~=6^||C*kceI=MgGopg>{8PnfKIO#n7LAmE5&m~_fw6DkJ)~D))
z>$%FVx(qt;T=5@rPR97{Bm7bNrmmlr_xd4jX0dFX<OVWhACu2}u~-q(&sdR_#Ul0n
zJjDtpq{H;Rmk-<b^cQh6#lpYh10BnC*dlfPJmNe~<JkWm>bCnh<n810d4|*Xap`9)
z7QB`1v~Ko8@@VnQjy-M0iS3#l(ysFqD`T3yjfH=Wa~z8Wr(F-@MAn7frZU9`uwUl9
zn-4wCQs2*0Ed0p!Yj$!Toz(BcdSWrZu_p2!ixpltr*HCm=P+$L>zbz@mI+>ae`5r@
z*|d+F;*<~S+BSXt&UtuR*l&t88Dp%>!8vpE^N96Z8n-dicrw<YcgN8#s~^T0`Wpj&
z_WAe}XQIQMOtxJ&7CcR{COjD{b2M9YeMYP&HUGrs#v1hQIJ&UuhjFUE^hukHHM2w7
zb)I4+E}FgX5Y|NReOSMv(J|b(HBRN5)wT2}lW=A{87uj2kvnzi+7Dw*xVoKA+qG)3
z?oF}6Ra=*C6Bc^|>(2V#=g;R6>nR%Czr~AjB5TX3!O(ezqkB5WYWtI0>ATf$!PO|J
zWYUeVQIA_xAHCathU)Dx#ftxMwVaN|0>g7GbR4mM&)~fO(|zb8)a&L2>N!2MX)>M8
zx{uHg!zFZ^s|s$beyi_eN31_+4D)*aMzgUdW;qrs=5L|+A<Yr@S}Ya1M(>B^%C@SR
z%viMrm9BkEK95*W)1b6sleKj7!9{&=be>{O=08v0?N)6Oj=D|I56d;-<vg9VZMM9P
z$3Aw%dPZX=_R1HWNBedx-P-7fu@Xkjk3K*rJ=U;Xt-g<)VKo=8JM|f%U&2wM+!stE
z;<Fl)0BlnDCPjjL5lQdL#C0}SjTuDKt&y#0i<~&ZJ*P2DZFQP+`db?FX}UmDxY(g-
zWVqP6G0?GeIoJ=au`|(5^P>pHn#g6WSf<5nm*d7>;PdkfxpdP?Ipdt};k9L<51(R9
z<TBPVa_MUx{dj@Gj4!{zeT-aFEbL^rXv?Lsj*;u#T_`?Fu8Jqa;Z^COuE$X}!u?TW
zCZZW@QlGhW)~(Zi2x}plELCsck5YGp5%EtnM!GWuK+E0o?`^(ooOU?z+YiBY0bA%O
z^t3~;yERVx@Dc7sjR9kHQKe>0lSHtPP7i7NA-FDG8z#i<Obh8oxIb$QdCStZMMN8b
zt(Po<Q~}eEmndWb9qgd3PSWI%K|8_}F?=^WXm_Y3#!@e`uAWHr<1ZAnAb3pMkzGrD
zx5HxRp{&&p!7c0{#N+HR!u?faNJ^HDbEdjR^+@abg&oFU0{ffB;1`Oafo$CdS~O)w
zou`PAqS@gp!A+#@SX))9@@0)lJTx!uf#00EL(P6IuGygJ`fp8wB}ZfLL;JhN#HmVx
zf2l6QHpM%Rb`<^i2L*BLxoeE;3cHEfhw1Yaie<!4`->~_S6q397PK91{Pn<p(wIPR
z64+Nl)WW~Cq3ef9VPoOxizxeV5mj%z^db2?gR0=`bUBkWx>~pfNoVa^`(dlox6{?~
zrDtfd)rX$%h59_9%S%mcyP!+g3*&Kyih;4Lsgi0BV>7f%LEEexoiy)}Dp|_H3(|CH
zG#-xQ5LXOom)c%<T3tnbc*LUYEJIqh9ZC#cEp)eqjmgLZT~LmPb6~O86w|eze5bs3
zSaCYY3bc7%!`1D9!tDxlKz+V>te~trhRW6&*KF#WGbmuV#d4(^2IJv{1ak~<zJG}G
z8E)Z(x-T*w*HNKl=5e*mJ)yc_Eo^ez2t?T{aK@Dkj|uCjxZ~&<oc9L<+@mLhyM4r`
zEVypM%EG@AB5B2Z5+QA)*g@{-h5D@hP)DpF`Gu-6IAD(Bb=ZA$0;;<P<%70Kfj)9#
z&^srf-dAo!GHB`kut1pr30=KdHtw1r>e9ND^w2W*LsZM$dQoJR9X>XJ%O=|63hT$h
zA~N=2eSCuD{e8*M0`3Qb1G-N!j{d}0L0K0dQK`XwU)9IJ-93TR6<&v{RQO%0Wp>t)
zeFEMyLG%Sz^xZkzcM6f9He%U)1h{vsp!^;UWz{IZXnF7fwXpdJDAP{H`=76j_3zct
zjSqH<esf+uj{ekGL0Q)=LX0q_;n!7-D>v!ema^mA*9-Momm3w-5>v0oRZ92@ElQQH
z-tjzQe0r>)tm}>r%IFg5OKRB>ss*b&Vubn3SV8%H8jX1Pp}%EY417V=8wM;a1TE6q
zXD6Vts65+miS5)^SE=vnHR9(ch`u3~p*!70f9Y;e=3I#R`3b1@#R??9npcW{Qi#+s
zhxq_;GV<VylS=K4<xCi0P+Uo&k*sJ~xV?-y9@^b%H0!8wcj_w-mixK|P8X)Xe*~p`
zo<C(rzbuLR#hjLu`pG4Ot~vcOJ;y^E%TB^hRV~m;2?!ZA79@Oi1f{I~WszWc@v2f)
z+pr?nEo&0KHo?_skudx(H>*K4;&<;D0lqO-P}UBv7_dCJL#>K7BvK=xTL53Imr7~l
zaSIj2fQ5wtqF%ieZfe1XZUMhq)Y^}4^+J7?UBrOjVBNWa6V4_O1gf~53_kewq*8k)
zMI{ni9vG-xl?;BgWrWe%GamNLO;QR2g{u`-Gvcdmf!T$;S|!Ir`<>#?Bp_XRrguHe
zL5*7o!3DG58$l^2JTn>|3io==*l>3~GQwz+Fdj}P8`D^-9qu)BEnL6Hyff#B@!(iN
zS(}ZDB#qW9rp$t4=?Yj$5yvtt3N2{(;RIIug9?1ZZj1x0ZjA_*c?`Rb2W{lWL+g3H
zcido+>m?fWRxLgavR_(X^YD`qma<lr-EScYKyUwUsCa=kEsS`A-=D|h7F7Inf~*yI
zha3Rlo>2p-0XSGl!vgtd6J#x#hmq?a$^iPL+Q5PZ2=Q;@ti}HM1YEm;0xtmh2K)m0
zXthBPZXxiRcp>~RC-_>O$`SO1481hJVOzq8ZNzoM)^#Qz2fv<FYR6T`H-UWUf0O@+
zw))@XOsM9P9N(Eh9sFiesU@ZG^<9~IOB{)nXLjIj+fcRue|!S2)kPQ^r_h&L#0*xX
zPFvKo*cKz;itwUly*&<b6Q~3Ez<6jikXad`fdHJ{;sY(0x(q8a+jbUxMOR8vu~6ix
zaTdx-=DYuFuK!@DuX1hXVsa(`OxK8)@e%S+5rJq|SGgJ$3wL^YGM#eq?tho-Tb%lK
z7jO`{VGr%j(@9GPGy2<1G&7Mp%M0j<924w==O&d(rsbwL1E4qft_DwAVE;xFQ14`(
zYna^iOr(C7Zp^hH^7+X`O0PSMJt_?5*#zRGyZ@cSSpSePkOpBG`!E=XqGIf>vNGy;
zJVM0|=E0M|ygXu1RysHf0~@R36%Pmrg5YkG115v~FJ;Dhh392kBviu+_X~`V8?c27
z2}rjp$3r1=Ie>75BpbRA@E$Nv@Lzm`rrTrV@hlbEA~xScSh>VW{!y{ebhl$XUZmn`
zAhns3+u79+Si5aKO;*`1s6Gp)6Q(w4ujEjeTvEqd-~6^K-hSm<u6)xQU;B>NU-^bB
IUU$iV0A<DP+W-In

literal 0
HcmV?d00001

diff --git a/examples/GRAPHICS/in.peptide-hbonds b/examples/GRAPHICS/in.peptide-hbonds
index 6c7c22953c6..34d0f831079 100644
--- a/examples/GRAPHICS/in.peptide-hbonds
+++ b/examples/GRAPHICS/in.peptide-hbonds
@@ -30,7 +30,8 @@ fix             2 all shake 0.0001 10 1000 b 4 6 8 10 12 14 18 a 31
 
 group           peptide type <= 12
 
-# compute center of mass of peptide and center of box and displace atoms so they will match
+# recenter system so the peptide is in the center of the box for that we compute the center
+# of mass of the peptide and the center of the box and displace atoms so tjat they will match
 variable comx equal xcm(peptide,x)
 variable comy equal xcm(peptide,y)
 variable comz equal xcm(peptide,z)
@@ -57,17 +58,18 @@ compute hb1 all hbond/local 3.5 30.0 pdonor woxygen hydrogen
 compute hb2 all hbond/local 3.7 30.0 woxygen pacceptor hydrogen
 
 # create donor/acceptor hydrogen bonds info text
-fix label all graphics/labels ${vizsteps} text " Hydrogen bonds donated:  $(c_hb1:%02.0f)" 207 72 0.0 &
+fix label all graphics/labels ${vizsteps} text " Hydrogen bonds donated:  $(c_hb1:%02.0f)" 206 62 0.0 &
                                                size 24 backcolor darkgray &
-                                          text "Hydrogen bonds accepted: $(c_hb2:%02.0f)" 210 30 0.0 &
-                                               size 24 backcolor darkgray
+                                          text "Hydrogen bonds accepted: $(c_hb2:%02.0f)" 209 20 0.0 &
+                                               size 24 backcolor darkgray &
+                                          image hbond-diagram.tga 190 210 0 scale 0.25 transcolor auto
 # create colored arrows to go with text labels
-fix obj   all graphics/objects ${vizsteps} arrow 5 80.0 57.0 35.0 80.0 66.0 35.0 0.3 0.2 &
-                                           arrow 6 80.0 57.0 33.5 80.0 66.0 33.5 0.3 0.2
+fix obj   all graphics/objects ${vizsteps} arrow 5 80.0 56.0 34.6 80.0 64.0 34.6 0.3 0.2 &
+                                           arrow 6 80.0 56.0 33.2 80.0 64.0 33.2 0.3 0.2
 
 # combine the graphics into visualization
 dump viz viz image ${vizsteps} hbonds-*.ppm element element size 800 800 zoom 2.1 view 80 0 bond atom 0.3&
-                       fsaa yes ssao yes 12384 0.6 shiny 0.1 box no 0.1 center s 0.65 0.49 0.42 &
+                       fsaa yes ssao yes 12384 0.6 shiny 0.1 box no 0.1 center s 0.65 0.45 0.42 &
                        compute hb1 const -0.4 0.3 compute hb2 const -0.4 0.3 &
                        fix label const 1 0 fix obj type 0.0 0.0
 dump_modify viz pad 5 boxcolor white backcolor darkgray backcolor2 silver &

From b1deefdd1fe2c847e0e484c20cfec4c729f3577c Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Mon, 2 Mar 2026 20:43:09 -0500
Subject: [PATCH 17/23] fix incorrect bond opacity lookup for bonds constrained
 with fix shake

---
 src/GRAPHICS/dump_image.cpp | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/src/GRAPHICS/dump_image.cpp b/src/GRAPHICS/dump_image.cpp
index 57115fd4648..0e4bbf1b11d 100644
--- a/src/GRAPHICS/dump_image.cpp
+++ b/src/GRAPHICS/dump_image.cpp
@@ -1502,15 +1502,15 @@ void DumpImage::create_image()
           xmid[2] = x[atom1][2] + 0.5*delz;
           if (bcolor == ATOM)
             image->draw_cylinder(x[atom1],xmid,color1,diameter,3,aopacity[type[atom1]]);
-          else image->draw_cylinder(x[atom1],xmid,color,diameter,3,bopacity[btype]);
+          else image->draw_cylinder(x[atom1],xmid,color,diameter,3,bopacity[itype]);
           xmid[0] = x[atom2][0] - 0.5*delx;
           xmid[1] = x[atom2][1] - 0.5*dely;
           xmid[2] = x[atom2][2] - 0.5*delz;
           if (bcolor == ATOM)
             image->draw_cylinder(xmid,x[atom2],color2,diameter,3,aopacity[type[atom1]]);
-          else image->draw_cylinder(xmid,x[atom2],color,diameter,3,bopacity[btype]);
+          else image->draw_cylinder(xmid,x[atom2],color,diameter,3,bopacity[itype]);
 
-        } else image->draw_cylinder(x[atom1],x[atom2],color,diameter,3,bopacity[btype]);
+        } else image->draw_cylinder(x[atom1],x[atom2],color,diameter,3,bopacity[itype]);
       }
     }
   }

From a97698af6845ad191d2a43800305ecde08f1c191 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Mon, 2 Mar 2026 20:50:49 -0500
Subject: [PATCH 18/23] cleaner version of bond opacity bugfix

---
 src/GRAPHICS/dump_image.cpp | 15 ++++++---------
 1 file changed, 6 insertions(+), 9 deletions(-)

diff --git a/src/GRAPHICS/dump_image.cpp b/src/GRAPHICS/dump_image.cpp
index 0e4bbf1b11d..e1c0b560399 100644
--- a/src/GRAPHICS/dump_image.cpp
+++ b/src/GRAPHICS/dump_image.cpp
@@ -1448,6 +1448,7 @@ void DumpImage::create_image()
         if (atom2 < 0 || !chooseghost[atom2]) continue;
         if (newton_bond == 0 && tag[atom1] > tag[atom2]) continue;
         if (btype == 0) continue;
+        if (btype < 0) btype = -btype;
 
         if (bcolor == ATOM) {
           if (acolor == TYPE) {
@@ -1464,9 +1465,7 @@ void DumpImage::create_image()
             color2 = image->color2rgb("white");
           }
         } else if (bcolor == TYPE) {
-          itype = btype;
-          if (itype < 0) itype = -itype;
-          color = bcolortype[itype];
+          color = bcolortype[btype];
         }
 
         if (bdiam == NUMERIC) {
@@ -1482,9 +1481,7 @@ void DumpImage::create_image()
             diameter = MIN(bufcopy[atom1][1],bufcopy[atom2][1]);
           }
         } else if (bdiam == TYPE) {
-          itype = btype;
-          if (itype < 0) itype = -itype;
-          diameter = bdiamtype[itype];
+          diameter = bdiamtype[btype];
         }
 
         // draw cylinder in 2 pieces if bcolor = ATOM
@@ -1502,15 +1499,15 @@ void DumpImage::create_image()
           xmid[2] = x[atom1][2] + 0.5*delz;
           if (bcolor == ATOM)
             image->draw_cylinder(x[atom1],xmid,color1,diameter,3,aopacity[type[atom1]]);
-          else image->draw_cylinder(x[atom1],xmid,color,diameter,3,bopacity[itype]);
+          else image->draw_cylinder(x[atom1],xmid,color,diameter,3,bopacity[btype]);
           xmid[0] = x[atom2][0] - 0.5*delx;
           xmid[1] = x[atom2][1] - 0.5*dely;
           xmid[2] = x[atom2][2] - 0.5*delz;
           if (bcolor == ATOM)
             image->draw_cylinder(xmid,x[atom2],color2,diameter,3,aopacity[type[atom1]]);
-          else image->draw_cylinder(xmid,x[atom2],color,diameter,3,bopacity[itype]);
+          else image->draw_cylinder(xmid,x[atom2],color,diameter,3,bopacity[btype]);
 
-        } else image->draw_cylinder(x[atom1],x[atom2],color,diameter,3,bopacity[itype]);
+        } else image->draw_cylinder(x[atom1],x[atom2],color,diameter,3,bopacity[btype]);
       }
     }
   }

From b602118fc80d1d66f262baa2c2eb37a7f7b2f092 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Tue, 3 Mar 2026 16:24:33 -0500
Subject: [PATCH 19/23] add a preset for compiling a library for use with
 LAMMPS-GUI

---
 cmake/presets/lammps-gui.cmake | 70 ++++++++++++++++++++++++++++++++++
 1 file changed, 70 insertions(+)
 create mode 100644 cmake/presets/lammps-gui.cmake

diff --git a/cmake/presets/lammps-gui.cmake b/cmake/presets/lammps-gui.cmake
new file mode 100644
index 00000000000..5a1792100ff
--- /dev/null
+++ b/cmake/presets/lammps-gui.cmake
@@ -0,0 +1,70 @@
+set(WIN_PACKAGES
+  AMOEBA
+  ASPHERE
+  BODY
+  BPM
+  BROWNIAN
+  CG-DNA
+  CG-SPICA
+  CLASS2
+  COLLOID
+  CORESHELL
+  DIELECTRIC
+  DIFFRACTION
+  DIPOLE
+  DPD-BASIC
+  DPD-MESO
+  DPD-REACT
+  DPD-SMOOTH
+  DRUDE
+  EFF
+  EXTRA-COMMAND
+  EXTRA-COMPUTE
+  EXTRA-DUMP
+  EXTRA-FIX
+  EXTRA-MOLECULE
+  EXTRA-PAIR
+  FEP
+  GRAPHICS
+  GRANULAR
+  INTERLAYER
+  KSPACE
+  LEPTON
+  MACHDYN
+  MANYBODY
+  MC
+  MEAM
+  MISC
+  ML-IAP
+  ML-POD
+  ML-SNAP
+  ML-UF3
+  MOFFF
+  MOLECULE
+  OPENMP
+  OPT
+  ORIENT
+  PERI
+  PHONON
+  QEQ
+  QTB
+  REACTION
+  REAXFF
+  RHEO
+  RIGID
+  SHOCK
+  SPH
+  SPIN
+  SRD
+  TALLY
+  UEF
+  YAFF)
+
+foreach(PKG ${WIN_PACKAGES})
+  set(PKG_${PKG} ON CACHE BOOL "" FORCE)
+endforeach()
+
+set(FFT KISS CACHE STRING "" FORCE)
+set(BUILD_MPI OFF CACHE BOOL "" FORCE)
+set(BUILD_TOOLS OFF CACHE BOOL "" FORCE)
+set(CMAKE_INSTALL_PREFIX "${CMAKE_CURRENT_BINARY_DIR}/lammps-installer")

From 8b45ceeb97c9b1df95584acc22fe0bbc696a98a4 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Wed, 4 Mar 2026 17:43:40 -0500
Subject: [PATCH 20/23] this is no longer needed

---
 src/.gitignore | 1 -
 1 file changed, 1 deletion(-)

diff --git a/src/.gitignore b/src/.gitignore
index b6da4da02e2..ef50bf2b835 100644
--- a/src/.gitignore
+++ b/src/.gitignore
@@ -845,7 +845,6 @@
 /ewald.h
 /ewald_cg.cpp
 /ewald_cg.h
-/ewald_const.h
 /ewald_dipole.cpp
 /ewald_dipole.h
 /ewald_dipole_spin.cpp

From 1279e336fef73ea953543637f9387592c52f56d1 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Wed, 4 Mar 2026 17:45:04 -0500
Subject: [PATCH 21/23] some more updates for the LAMMPS-GUI custom LAMMPS
 library build

---
 cmake/presets/lammps-gui.cmake | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/cmake/presets/lammps-gui.cmake b/cmake/presets/lammps-gui.cmake
index 5a1792100ff..af18bdac9d0 100644
--- a/cmake/presets/lammps-gui.cmake
+++ b/cmake/presets/lammps-gui.cmake
@@ -8,6 +8,7 @@ set(WIN_PACKAGES
   CG-SPICA
   CLASS2
   COLLOID
+  COLVARS
   CORESHELL
   DIELECTRIC
   DIFFRACTION
@@ -64,6 +65,9 @@ foreach(PKG ${WIN_PACKAGES})
   set(PKG_${PKG} ON CACHE BOOL "" FORCE)
 endforeach()
 
+set(BUILD_SHARED_LIBS ON CACHE BOOL "" FORCE)
+set(CMAKE_BUILD_TYPE Release CACHE STRING "" FORCE)
+set(USE_INTERNAL_LINALG ON CACHE BOOL "" FORCE)
 set(FFT KISS CACHE STRING "" FORCE)
 set(BUILD_MPI OFF CACHE BOOL "" FORCE)
 set(BUILD_TOOLS OFF CACHE BOOL "" FORCE)

From f562372072f501225ee7952a243904bcbf0b6cd5 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Wed, 4 Mar 2026 17:45:38 -0500
Subject: [PATCH 22/23] use new feature, so we can color bodies by atom
 properties also from LAMMPS-GUI

---
 examples/GRAPHICS/in.cubes-and-pyramids | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/examples/GRAPHICS/in.cubes-and-pyramids b/examples/GRAPHICS/in.cubes-and-pyramids
index bed142c5151..01e55496758 100644
--- a/examples/GRAPHICS/in.cubes-and-pyramids
+++ b/examples/GRAPHICS/in.cubes-and-pyramids
@@ -47,7 +47,7 @@ fix pbar all graphics/objects ${out} progbar 3 4 y 24.0 12.0 -1.4 30.0 0.5 v_pro
 
 dump         2 all image ${out} cubes-and-pyramids.*.ppm type type size 600 600 shiny 0.4 &
                          box yes 0.025 axes no 0.5 0.1 zoom 1.5 view 80 0 fsaa yes ssao yes 315123 0.4 &
-                         body type 0.2 3 fix pbar type 0 0
+                         body atom 0.2 3 fix pbar type 0 0
 #            set cylinder diameter--^  ^-- set viz style (1 only faces, 2 cylinders, 3 both)
 dump_modify  2 pad 6 backcolor darkgray backcolor2 gray boxcolor white &
                      acolor 1 steelblue acolor 2 orange  acolor 3 slategray  acolor 4 red

From 320f26bc96acc3bfc224425cb8c4c244920131c5 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Thu, 5 Mar 2026 08:42:05 -0500
Subject: [PATCH 23/23] flush small values for stress or global scalars or
 vectors to zero to avoid false positives

---
 unittest/force-styles/test_fix_timestep.cpp | 14 ++++++++++++--
 1 file changed, 12 insertions(+), 2 deletions(-)

diff --git a/unittest/force-styles/test_fix_timestep.cpp b/unittest/force-styles/test_fix_timestep.cpp
index eb1f1c80362..83d2acace79 100644
--- a/unittest/force-styles/test_fix_timestep.cpp
+++ b/unittest/force-styles/test_fix_timestep.cpp
@@ -197,6 +197,9 @@ void generate_yaml_file(const char *outfile, const TestConfig &config)
         // run_stress, if enabled
         if (ifix->thermo_virial) {
             auto *stress = ifix->virial;
+            // avoid false positives on tiny stresses. force to zero instead.
+            for (int i = 0; i < 6; ++i)
+                if (fabs(stress[i]) < 1.0e-13) stress[i] = 0.0;
             block = fmt::format("{:23.16e} {:23.16e} {:23.16e} {:23.16e} {:23.16e} {:23.16e}",
                                 stress[0], stress[1], stress[2], stress[3], stress[4], stress[5]);
             writer.emit_block("run_stress", block);
@@ -205,6 +208,8 @@ void generate_yaml_file(const char *outfile, const TestConfig &config)
         // global scalar
         if (ifix->scalar_flag) {
             double value = ifix->compute_scalar();
+            // avoid false positives on tiny values. force to zero instead.
+            if (fabs(value) < 1.0e-13) value = 0.0;
             writer.emit("global_scalar", value);
         }
 
@@ -212,8 +217,13 @@ void generate_yaml_file(const char *outfile, const TestConfig &config)
         if (ifix->vector_flag) {
             int num = ifix->size_vector;
             block   = std::to_string(num);
-            for (int i = 0; i < num; ++i)
-                block += fmt::format(" {}", ifix->compute_vector(i));
+            double value;
+            for (int i = 0; i < num; ++i) {
+                // avoid false positives on tiny values. force to zero instead.
+                value = ifix->compute_vector(i);
+                if (fabs(value) < 1.0e-13) value = 0.0;
+                block += fmt::format(" {:23.16e}", value);
+            }
             writer.emit_block("global_vector", block);
         }
     }