solweig-gpu 2.0.0__tar.gz → 2.0.1__tar.gz

{solweig_gpu-2.0.0/solweig_gpu.egg-info → solweig_gpu-2.0.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: solweig-gpu
-Version: 2.0.0
+Version: 2.0.1
 Summary: GPU-accelerated SOLWEIG model for urban thermal comfort simulation
 Home-page: https://github.com/nvnsudharsan/SOLWEIG-GPU
 Author: Harsh Kamath, Naveen Sudharsan
@@ -56,7 +56,7 @@ Dynamic: summary
   <a href="https://www.gnu.org/licenses/gpl-3.0"><img src="https://img.shields.io/badge/License-GPLv3-%230ab5b3.svg" alt="License: GPL v3"></a>
   <a href="https://pepy.tech/projects/solweig-gpu"><img src="https://static.pepy.tech/personalized-badge/solweig-gpu?period=total&units=INTERNATIONAL_SYSTEM&left_color=GRAY&right_color=GREEN&left_text=downloads" alt="PyPI Downloads"></a>
   <a href="https://joss.theoj.org/papers/27faa2bf5f6058d981df8b565f8e9a34"><img src="https://joss.theoj.org/papers/27faa2bf5f6058d981df8b565f8e9a34/status.svg"></a>
-  <a href="https://doi.org/10.5194/gmd-19-7389-2026"><img src="https://doi.org/10.5194/gmd-19-7389-2026.svg" alt="EMD"></a>
+  <a href="https://doi.org/10.5194/gmd-19-7389-2026"><img src="https://img.shields.io/badge/DOI-10.5194%2Fgmd--19--7389--2026-blue" alt="GMD"></a>
   <a href="https://github.com/nvnsudharsan/solweig-gpu/actions/workflows/tests.yml"><img src="https://img.shields.io/badge/Tests-Passing-%23ffb703.svg" alt="Tests"></a>
 </p>
 
@@ -77,7 +77,7 @@ Dynamic: summary
 
 1. Kamath, H. G., Sudharsan, N., Singh, M., Wallenberg, N., Lindberg, F., & Niyogi, D. (2026). SOLWEIG-GPU: GPU-Accelerated Thermal Comfort Modeling Framework for Urban Digital Twins. Journal of Open Source Software, 11(118), 9535. https://doi.org/10.21105/joss.09535
 
-2. Zonato, A., Kamath, H.G., Sudharsan, N., Monaco, L., Kittner, J., Wolf, L., Demuzere, M.A., Middel, A., Bechtel, B. and Milelli, M., 2026. GLIDE-SOL: A GPU-accelerated Global Lightweight Infrastructure for Diagnostic Environmental Modeling with SOLWEIG. EGUsphere, 2026, pp.1-30. https://doi.org/10.5194/egusphere-2026-776
+2. Zonato, A., Kamath, H. G., Sudharsan, N., Monaco, L., Kittner, J., Wolf, L., Demuzere, M., Middel, A., Bechtel, B., and Milelli, M.: GLIDE-SOL: a GPU-accelerated global lightweight infrastructure for diagnostic environmental modeling with SOLWEIG, Geosci. Model Dev., 19, 7389–7413, https://doi.org/10.5194/gmd-19-7389-2026, 2026.
 
 **SOLWEIG** was originally developed by Dr. Fredrik Lindberg's group. Journal reference: Lindberg, F., Holmer, B. & Thorsson, S. SOLWEIG 1.0 – Modelling spatial variations of 3D radiant fluxes and mean radiant temperature in complex urban settings. *Int J Biometeorol* 52, 697–713 (2008). https://doi.org/10.1007/s00484-008-0162-7
 

{solweig_gpu-2.0.0 → solweig_gpu-2.0.1}/README.md

@@ -13,7 +13,7 @@
   <a href="https://www.gnu.org/licenses/gpl-3.0"><img src="https://img.shields.io/badge/License-GPLv3-%230ab5b3.svg" alt="License: GPL v3"></a>
   <a href="https://pepy.tech/projects/solweig-gpu"><img src="https://static.pepy.tech/personalized-badge/solweig-gpu?period=total&units=INTERNATIONAL_SYSTEM&left_color=GRAY&right_color=GREEN&left_text=downloads" alt="PyPI Downloads"></a>
   <a href="https://joss.theoj.org/papers/27faa2bf5f6058d981df8b565f8e9a34"><img src="https://joss.theoj.org/papers/27faa2bf5f6058d981df8b565f8e9a34/status.svg"></a>
-  <a href="https://doi.org/10.5194/gmd-19-7389-2026"><img src="https://doi.org/10.5194/gmd-19-7389-2026.svg" alt="EMD"></a>
+  <a href="https://doi.org/10.5194/gmd-19-7389-2026"><img src="https://img.shields.io/badge/DOI-10.5194%2Fgmd--19--7389--2026-blue" alt="GMD"></a>
   <a href="https://github.com/nvnsudharsan/solweig-gpu/actions/workflows/tests.yml"><img src="https://img.shields.io/badge/Tests-Passing-%23ffb703.svg" alt="Tests"></a>
 </p>
 
@@ -34,7 +34,7 @@
 
 1. Kamath, H. G., Sudharsan, N., Singh, M., Wallenberg, N., Lindberg, F., & Niyogi, D. (2026). SOLWEIG-GPU: GPU-Accelerated Thermal Comfort Modeling Framework for Urban Digital Twins. Journal of Open Source Software, 11(118), 9535. https://doi.org/10.21105/joss.09535
 
-2. Zonato, A., Kamath, H.G., Sudharsan, N., Monaco, L., Kittner, J., Wolf, L., Demuzere, M.A., Middel, A., Bechtel, B. and Milelli, M., 2026. GLIDE-SOL: A GPU-accelerated Global Lightweight Infrastructure for Diagnostic Environmental Modeling with SOLWEIG. EGUsphere, 2026, pp.1-30. https://doi.org/10.5194/egusphere-2026-776
+2. Zonato, A., Kamath, H. G., Sudharsan, N., Monaco, L., Kittner, J., Wolf, L., Demuzere, M., Middel, A., Bechtel, B., and Milelli, M.: GLIDE-SOL: a GPU-accelerated global lightweight infrastructure for diagnostic environmental modeling with SOLWEIG, Geosci. Model Dev., 19, 7389–7413, https://doi.org/10.5194/gmd-19-7389-2026, 2026.
 
 **SOLWEIG** was originally developed by Dr. Fredrik Lindberg's group. Journal reference: Lindberg, F., Holmer, B. & Thorsson, S. SOLWEIG 1.0 – Modelling spatial variations of 3D radiant fluxes and mean radiant temperature in complex urban settings. *Int J Biometeorol* 52, 697–713 (2008). https://doi.org/10.1007/s00484-008-0162-7
 

{solweig_gpu-2.0.0 → solweig_gpu-2.0.1}/setup.py

@@ -18,7 +18,7 @@ long_description = (this_directory / "README.md").read_text(encoding="utf-8")
 
 setup(
     name="solweig-gpu",
-    version="2.0.0",
+    version="2.0.1",
     description="GPU-accelerated SOLWEIG model for urban thermal comfort simulation",
     long_description=long_description,
     long_description_content_type='text/markdown',

{solweig_gpu-2.0.0 → solweig_gpu-2.0.1}/solweig_gpu/__init__.py

@@ -11,7 +11,7 @@
 #MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
 #GNU General Public License for more details.
 
-__version__ = "2.0.0"
+__version__ = "2.0.1"
 from .solweig_gpu import (
     thermal_comfort,
     preprocess,

{solweig_gpu-2.0.0 → solweig_gpu-2.0.1}/solweig_gpu/landcoverclasses_2016a.txt

@@ -5,4 +5,4 @@ Cobble_stone_2014a 0   0.20 0.95 0.37   -3.41  15.0
 Water              7   0.05 0.98 0.00    0.00  12.0
 Grass_unmanaged    5   0.16 0.94 0.21   -3.38  14.0
 bare_soil          6   0.25 0.94 0.33   -3.01  14.0
-Walls             99   0.20 0.90 0.58   -3.41  15.0
+Walls             99   0.20 0.90 0.37   -3.41  15.0

{solweig_gpu-2.0.0 → solweig_gpu-2.0.1}/solweig_gpu/shadow.py

@@ -32,6 +32,25 @@ gdal.UseExceptions()
 
 device = torch.device("cuda" if torch.cuda.is_available() else "cpu")
 
+# Increment when cached shadow/SVF outputs become incompatible with the
+# current shadow algorithm. Version 2 adds UMEP's pergola-aware vegetation
+# ray tracking.
+SVF_CACHE_VERSION = 2
+
+
+def is_current_svf_cache(npz_path):
+    """Return whether a shadow-matrix cache uses the current algorithm."""
+    if not os.path.isfile(npz_path):
+        return False
+
+    try:
+        with np.load(npz_path) as npz:
+            if "cache_version" not in npz.files:
+                return False
+            return int(np.asarray(npz["cache_version"]).item()) == SVF_CACHE_VERSION
+    except (OSError, ValueError, TypeError):
+        return False
+
 def load_raster_to_tensor(dem_path, device=device):
     dataset = gdal.Open(dem_path)
     if dataset is None:
@@ -158,6 +177,7 @@ def save_svf_zip_npz_outputs(output_dir, gdal_dsm, svf, svfE, svfS, svfW, svfN,
 
     np.savez_compressed(
         npz_path,
+        cache_version=np.asarray(SVF_CACHE_VERSION, dtype=np.int32),
         shadowmat=tensor_to_numpy(shmat).astype(np.float32),
         vegshadowmat=tensor_to_numpy(vegshmat).astype(np.float32),
         vbshmat=tensor_to_numpy(vbshvegshmat).astype(np.float32),
@@ -191,14 +211,14 @@ def shadow(amaxvalue, a, vegdem, vegdem2, bush, azimuth, altitude, scale):
         - Implements anisotropic shadow casting
         - Accounts for vegetation transmittance
     """
-    device = torch.device("cuda" if torch.cuda.is_available() else "cpu")
+    device = a.device
     degrees = torch.pi / 180.
-    if azimuth == 0.0:
-        azimuth = 1e-12
-    azimuth = ensure_tensor(azimuth)
-    altitude = ensure_tensor(altitude)
-    azimuth = azimuth * degrees #torch.tensor(azimuth * degrees, device=a.device)
-    altitude = altitude * degrees #torch.tensor(altitude * degrees, device=a.device)
+    # Use float64 for ray angles. In float32, tan(90 degrees) is a large
+    # negative value and makes the zenith patch cast horizontal shadows.
+    azimuth = ensure_tensor(azimuth, device=device).to(device=device, dtype=torch.float64)
+    altitude = ensure_tensor(altitude, device=device).to(device=device, dtype=torch.float64)
+    azimuth = azimuth * degrees
+    altitude = altitude * degrees
 
     dx = 0.
     dy = 0.
@@ -206,8 +226,6 @@ def shadow(amaxvalue, a, vegdem, vegdem2, bush, azimuth, altitude, scale):
     sizex = a.shape[0]
     sizey = a.shape[1]
 
-    device = a.device
-
     dx = torch.tensor(dx, device=device)
     dy = torch.tensor(dy, device=device)
     dz = torch.tensor(dz, device=device)
@@ -217,10 +235,10 @@ def shadow(amaxvalue, a, vegdem, vegdem2, bush, azimuth, altitude, scale):
     tempvegdem2 = torch.zeros((sizex, sizey), device=device)
     sh = torch.zeros((sizex, sizey), device=device)
     vbshvegsh = torch.zeros((sizex, sizey), device=device)
-    tempbush = torch.zeros((sizex, sizey), device=device)
+    templastfabovea = torch.zeros((sizex, sizey), device=device)
+    templastgabovea = torch.zeros((sizex, sizey), device=device)
 
     f = a.clone()
-    g = torch.zeros((sizex, sizey), device=device)
     bushplant = (bush > 1.).float()
     vegsh = torch.zeros((sizex, sizey), device=device) + bushplant
 
@@ -237,9 +255,14 @@ def shadow(amaxvalue, a, vegdem, vegdem2, bush, azimuth, altitude, scale):
     dscos = torch.abs((1. / cosazimuth))
     tanaltitudebyscale = torch.tan(altitude) / scale
 
-    index = 1
+    # The previous ray height is required to detect vegetation layers that
+    # fall between two ray steps (for example, thin pergola canopies).
+    index = 0
+    dzprev = torch.tensor(0.0, device=device)
     fabovea = None
     gabovea = None
+    lastfabovea = None
+    lastgabovea = None
     vegsh2 = None
 
     while (amaxvalue >= dz and torch.abs(dx) < sizex and torch.abs(dy) < sizey):
@@ -257,6 +280,8 @@ def shadow(amaxvalue, a, vegdem, vegdem2, bush, azimuth, altitude, scale):
         tempvegdem.zero_()
         tempvegdem2.zero_()
         temp.zero_()
+        templastfabovea.zero_()
+        templastgabovea.zero_()
         absdx = torch.abs(dx)
         absdy = torch.abs(dy)
         xc1 = int((dx + absdx) / 2.)
@@ -278,45 +303,41 @@ def shadow(amaxvalue, a, vegdem, vegdem2, bush, azimuth, altitude, scale):
 
         fabovea = tempvegdem > a
         gabovea = tempvegdem2 > a
-        vegsh2 = fabovea.float() - gabovea.float()
+
+        templastfabovea[xp1:xp2, yp1:yp2] = vegdem[xc1:xc2, yc1:yc2] - dzprev
+        templastgabovea[xp1:xp2, yp1:yp2] = vegdem2[xc1:xc2, yc1:yc2] - dzprev
+        lastfabovea = templastfabovea > a
+        lastgabovea = templastgabovea > a
+        dzprev = dz
+
+        # Match UMEP's pergola-aware shadowingfunction_20. A ray intersects
+        # vegetation when any, but not all, of the current/previous canopy
+        # and trunk-base tests are above the target surface.
+        vegetation_states = (
+            fabovea.float()
+            + gabovea.float()
+            + lastfabovea.float()
+            + lastgabovea.float()
+        )
+        vegsh2 = ((vegetation_states > 0) & (vegetation_states < 4)).float()
 
         vegsh = torch.max(vegsh, vegsh2)
         vegsh[(vegsh * sh > 0.)] = 0.
 
         vbshvegsh = vegsh + vbshvegsh
 
-        if index == 1.:
-            firstvegdem = tempvegdem - temp
-            firstvegdem[firstvegdem <= 0.] = 1000.
-            vegsh[firstvegdem < dz] = 1.
-            vegsh = vegsh * (vegdem2 > a).float()
-            vbshvegsh.zero_()
-
-        if bush.max() > 0. and torch.max(fabovea * bush) > 0.:
-            tempbush.zero_()
-            tempbush[int(xp1):int(xp2), int(yp1):int(yp2)] = bush[int(xc1):int(xc2), int(yc1):int(yc2)] - dz
-            g = torch.max(g, tempbush)
-            g *= bushplant
-
         index += 1.
 
     sh = 1. - sh
     vbshvegsh[vbshvegsh > 0.] = 1.
     vbshvegsh = vbshvegsh - vegsh
 
-    if bush.max() > 0.:
-        g = g - bush
-        g[g > 0.] = 1.
-        g[g < 0.] = 0.
-        vegsh = vegsh - bushplant + g
-        vegsh[vegsh < 0.] = 0.
-
-    vegsh[vegsh > 0.] = 1.
     vegsh = 1. - vegsh
     vbshvegsh = 1. - vbshvegsh
 
-    # Changed here
-    del tempvegdem, tempvegdem2, temp, tempbush, fabovea, gabovea, vegsh2
+    del tempvegdem, tempvegdem2, temp
+    del templastfabovea, templastgabovea
+    del fabovea, gabovea, lastfabovea, lastgabovea, vegsh2
     torch.cuda.empty_cache()
     # Changed here
     return sh, vegsh, vbshvegsh

{solweig_gpu-2.0.0 → solweig_gpu-2.0.1}/solweig_gpu/solweig_gpu.py

@@ -293,7 +293,7 @@ def calculate_svf(base_path: str, patch_option: int = 2, overwrite: bool = False
     import torch
     from tqdm import tqdm
 
-    from .shadow import svf_calculator
+    from .shadow import is_current_svf_cache, svf_calculator
 
     building_dsm_dir = os.path.join(base_path, "Building_DSM")
     dem_dir = os.path.join(base_path, "DEM")
@@ -369,7 +369,7 @@ def calculate_svf(base_path: str, patch_option: int = 2, overwrite: bool = False
             not overwrite
             and os.path.exists(expected_tif)
             and os.path.exists(expected_zip)
-            and os.path.exists(expected_npz)
+            and is_current_svf_cache(expected_npz)
         ):
             continue
 

{solweig_gpu-2.0.0 → solweig_gpu-2.0.1}/solweig_gpu/utci_process.py

@@ -30,7 +30,7 @@ import pytz
 import datetime
 from .Tgmaps_v1 import Tgmaps_v1
 from .sun_position import Solweig_2015a_metdata_noload
-from .shadow import svf_calculator, create_patches
+from .shadow import create_patches, is_current_svf_cache, svf_calculator
 from .solweig import Solweig_2022a_calc, clearnessindex_2013b
 from .calculate_utci import utci_calculator
 from .calculate_wbgt import isobaric_wet_bulb_temperature_from_rh, black_globe_temperature
@@ -245,7 +245,7 @@ def _svf_cache_exists(building_dsm_path: str, number: str) -> bool:
     return (
         os.path.isfile(svftotal_path)
         and os.path.isfile(zip_path)
-        and os.path.isfile(npz_path)
+        and is_current_svf_cache(npz_path)
     )
 
 def _load_raster_to_tensor_from_path(raster_path: str):
@@ -411,7 +411,7 @@ def compute_utci(building_dsm_path, tree_path, dem_path, walls_path, aspect_path
                 f"Building DSM shape {tuple(a.shape)}"
             )
  
-    # Added
+    # Landcover warnings
     landcover = 0
 
     if landcover_path is not None:
@@ -439,7 +439,7 @@ def compute_utci(building_dsm_path, tree_path, dem_path, walls_path, aspect_path
         lc_class = np.empty((len(lines), 6), dtype=float)
         for i, ln in enumerate(lines):
             lc_class[i, :] = [float(x) for x in ln.split()[1:]]  # cols 1-6
-    # Added
+    # Landcover warnings
     
     base_date = datetime.datetime.strptime(selected_date_str, "%Y-%m-%d")
     rows, cols = a.shape
@@ -518,7 +518,7 @@ def compute_utci(building_dsm_path, tree_path, dem_path, walls_path, aspect_path
     Tgmap1N = torch.zeros((rows, cols), device=device)
     TgOut1 = torch.zeros((rows, cols), device=device)
     
-    # Added
+    # Landcover classification based calculations 
     if landcover == 1:                                     
         (TgK_np, Tstart_np, alb_np, emis_np, TgK_wall_np, Tstart_wall_np, TmaxLST_np,
          TmaxLST_wall_np) = Tgmaps_v1(lcgrid_np, lc_class)
@@ -540,7 +540,6 @@ def compute_utci(building_dsm_path, tree_path, dem_path, walls_path, aspect_path
         Tstart_wall = -3.41
         TmaxLST = 15.
         TmaxLST_wall = 15.
-    # Added
     
     transVeg = 3. / 100.
     # landcover = 1 # Modified
@@ -929,4 +928,3 @@ def compute_utci(building_dsm_path, tree_path, dem_path, walls_path, aspect_path
     time_taken = end_time - start_time
     print(f"Time taken to execute tile {number}: {time_taken:.2f} seconds")
 
-

{solweig_gpu-2.0.0 → solweig_gpu-2.0.1/solweig_gpu.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: solweig-gpu
-Version: 2.0.0
+Version: 2.0.1
 Summary: GPU-accelerated SOLWEIG model for urban thermal comfort simulation
 Home-page: https://github.com/nvnsudharsan/SOLWEIG-GPU
 Author: Harsh Kamath, Naveen Sudharsan
@@ -56,7 +56,7 @@ Dynamic: summary
   <a href="https://www.gnu.org/licenses/gpl-3.0"><img src="https://img.shields.io/badge/License-GPLv3-%230ab5b3.svg" alt="License: GPL v3"></a>
   <a href="https://pepy.tech/projects/solweig-gpu"><img src="https://static.pepy.tech/personalized-badge/solweig-gpu?period=total&units=INTERNATIONAL_SYSTEM&left_color=GRAY&right_color=GREEN&left_text=downloads" alt="PyPI Downloads"></a>
   <a href="https://joss.theoj.org/papers/27faa2bf5f6058d981df8b565f8e9a34"><img src="https://joss.theoj.org/papers/27faa2bf5f6058d981df8b565f8e9a34/status.svg"></a>
-  <a href="https://doi.org/10.5194/gmd-19-7389-2026"><img src="https://doi.org/10.5194/gmd-19-7389-2026.svg" alt="EMD"></a>
+  <a href="https://doi.org/10.5194/gmd-19-7389-2026"><img src="https://img.shields.io/badge/DOI-10.5194%2Fgmd--19--7389--2026-blue" alt="GMD"></a>
   <a href="https://github.com/nvnsudharsan/solweig-gpu/actions/workflows/tests.yml"><img src="https://img.shields.io/badge/Tests-Passing-%23ffb703.svg" alt="Tests"></a>
 </p>
 
@@ -77,7 +77,7 @@ Dynamic: summary
 
 1. Kamath, H. G., Sudharsan, N., Singh, M., Wallenberg, N., Lindberg, F., & Niyogi, D. (2026). SOLWEIG-GPU: GPU-Accelerated Thermal Comfort Modeling Framework for Urban Digital Twins. Journal of Open Source Software, 11(118), 9535. https://doi.org/10.21105/joss.09535
 
-2. Zonato, A., Kamath, H.G., Sudharsan, N., Monaco, L., Kittner, J., Wolf, L., Demuzere, M.A., Middel, A., Bechtel, B. and Milelli, M., 2026. GLIDE-SOL: A GPU-accelerated Global Lightweight Infrastructure for Diagnostic Environmental Modeling with SOLWEIG. EGUsphere, 2026, pp.1-30. https://doi.org/10.5194/egusphere-2026-776
+2. Zonato, A., Kamath, H. G., Sudharsan, N., Monaco, L., Kittner, J., Wolf, L., Demuzere, M., Middel, A., Bechtel, B., and Milelli, M.: GLIDE-SOL: a GPU-accelerated global lightweight infrastructure for diagnostic environmental modeling with SOLWEIG, Geosci. Model Dev., 19, 7389–7413, https://doi.org/10.5194/gmd-19-7389-2026, 2026.
 
 **SOLWEIG** was originally developed by Dr. Fredrik Lindberg's group. Journal reference: Lindberg, F., Holmer, B. & Thorsson, S. SOLWEIG 1.0 – Modelling spatial variations of 3D radiant fluxes and mean radiant temperature in complex urban settings. *Int J Biometeorol* 52, 697–713 (2008). https://doi.org/10.1007/s00484-008-0162-7
 

{solweig_gpu-2.0.0 → solweig_gpu-2.0.1}/solweig_gpu.egg-info/SOURCES.txt

@@ -38,6 +38,7 @@ tests/test_gpu_smoke.py
 tests/test_integration.py
 tests/test_integration_end_to_end.py
 tests/test_preprocessor.py
+tests/test_shadow.py
 tests/test_sun_position.py
 tests/test_tiling_and_metfiles.py
 tests/test_walls_aspect.py

solweig_gpu-2.0.1/tests/test_shadow.py

@@ -0,0 +1,74 @@
+import numpy as np
+import torch
+
+from solweig_gpu.shadow import SVF_CACHE_VERSION, is_current_svf_cache, shadow
+
+
+def _isolated_tree_scene(size=31):
+    dsm = torch.zeros((size, size), dtype=torch.float32)
+    canopy = torch.zeros_like(dsm)
+    trunk = torch.zeros_like(dsm)
+    bush = torch.zeros_like(dsm)
+
+    center = size // 2
+    canopy[center, center] = 10.0
+    trunk[center, center] = 2.5
+    return dsm, canopy, trunk, bush
+
+
+def test_shadow_tracks_previous_ray_step_for_thin_canopies():
+    dsm, canopy, trunk, bush = _isolated_tree_scene()
+
+    _, vegetation_visible, _ = shadow(
+        torch.tensor(10.0),
+        dsm,
+        canopy,
+        trunk,
+        bush,
+        torch.tensor(5.0),
+        torch.tensor(42.0),
+        1.0,
+    )
+
+    shaded = set(map(tuple, torch.nonzero(vegetation_visible == 0).tolist()))
+    assert shaded == {
+        (18, 15),
+        (19, 15),
+        (20, 15),
+        (21, 14),
+        (22, 14),
+        (23, 14),
+        (24, 14),
+        (25, 14),
+        (26, 14),
+        (27, 14),
+    }
+
+
+def test_zenith_patch_does_not_cast_horizontal_building_shadows():
+    dsm = torch.zeros((9, 9), dtype=torch.float32)
+    dsm[4, 4] = 10.0
+    vegetation = torch.zeros_like(dsm)
+
+    building_visible, _, _ = shadow(
+        torch.tensor(10.0),
+        dsm,
+        vegetation,
+        vegetation,
+        vegetation,
+        torch.tensor(0.0),
+        torch.tensor(90.0),
+        1.0,
+    )
+
+    assert torch.all(building_visible == 1)
+
+
+def test_pre_pergola_shadow_cache_is_invalidated(tmp_path):
+    old_cache = tmp_path / "old.npz"
+    current_cache = tmp_path / "current.npz"
+    np.savez_compressed(old_cache, shadowmat=np.ones((1, 1, 1)))
+    np.savez_compressed(current_cache, cache_version=np.asarray(SVF_CACHE_VERSION))
+
+    assert not is_current_svf_cache(old_cache)
+    assert is_current_svf_cache(current_cache)