fix(tidy3d): FXC-4563-use-2-d-intersections-for-adjoint-monitor-sizes-in-2-d-simulations

[marcorudolphflex]

Previously, we just took the geometry bounds to construct adjoint monitors. This PR uses the plane intersection in 2D simulations.

Greptile Overview

Greptile Summary

This PR improves adjoint monitor sizing in 2D simulations by using planar geometry intersections instead of full bounding boxes. In 2D simulations (detected via self.size.count(0.0) == 1), the code now calls geometry.intersections_plane() to compute tight monitor bounds based on the actual geometry cross-section at the simulation plane.

Key changes:

• simulation.py:4887 detects 2D simulations and passes plane info to structure monitor generation
• structure.py:328-352 implements _box_from_plane_intersection() that computes monitor bounds from plane intersections with multi-level fallback (geometry → bounding box → full bounding box)
• Comprehensive tests validate behavior for spheres, meshes, cylinders, boxes, and polyslabs in both 2D and 3D cases
• Changelog entry added under 'Fixed' section

The implementation handles edge cases through a cascading fallback strategy and maintains backward compatibility for 3D simulations.

Confidence Score: 5/5

• This PR is safe to merge with comprehensive test coverage and proper fallback handling
• The changes are well-tested with multiple geometry types, include proper fallback logic for edge cases, maintain backward compatibility for 3D simulations, and fix an actual bug where adjoint monitors were oversized in 2D simulations. Code follows established patterns and includes appropriate documentation.
• No files require special attention

Important Files Changed

File Analysis

Filename	Score	Overview
CHANGELOG.md	5/5	Added changelog entry under 'Fixed' section documenting adjoint monitor cropping in 2D simulations
tests/test_components/autograd/test_adjoint_monitors.py	5/5	New comprehensive test file covering 2D and 3D adjoint monitor sizing with various geometries
tidy3d/components/simulation.py	5/5	Detects 2D simulations and passes plane info to structure monitor generation
tidy3d/components/structure.py	4/5	Implements plane intersection logic for 2D adjoint monitors, with fallback to bounding box

Sequence Diagram

sequenceDiagram
    participant Sim as Simulation
    participant Struct as Structure
    participant Geom as Geometry
    participant BoxCls as Box

    Note over Sim: User calls _make_adjoint_monitors
    Sim->>Sim: Check if 2D via size.count(0.0) == 1
    alt Is 2D simulation
        Sim->>Sim: sim_plane = self
    else Is 3D simulation
        Sim->>Sim: sim_plane = None
    end
    
    loop For each structure
        Sim->>Struct: _make_adjoint_monitors(freqs, index, field_keys, plane)
        
        alt plane is not None (2D case)
            Struct->>Struct: Call _box_from_plane_intersection
            Struct->>Geom: geometry.intersections_plane
            Geom-->>Struct: Return intersection shapes
            
            alt Intersections found
                Struct->>Struct: Compute union of bounds
            else No intersections
                Struct->>Geom: geom_box.intersections_plane
                alt Box intersections found
                    Struct->>Struct: Compute union of bounds
                else Still no intersections
                    Struct->>Struct: Use geom_box as fallback
                end
            end
            
            Struct->>BoxCls: Box.from_bounds(rmin, rmax)
            BoxCls-->>Struct: Return planar box
        else plane is None (3D case)
            Struct->>Struct: box = geom_box
        end
        
        Struct->>Struct: Create FieldMonitor with box
        Struct->>Struct: Create PermittivityMonitor with box
        Struct-->>Sim: Return monitors
    end
    
    Sim-->>Sim: Collect all monitors

Loading

---

Note

Improves adjoint monitor sizing for 2D simulations by cropping to the simulation plane’s geometry intersection; 3D behavior unchanged.

• Simulation passes a planar context to structures in _make_adjoint_monitors; structures compute monitor Box via geometry.intersections_plane(...) with robust fallbacks
• New tests validate 2D plane-cropped bounds (spheres, meshes incl. disjoint components) and confirm 3D still uses bounding_box
• CHANGELOG updated; minor deterministic test data tweak (random seed)

^{Written by Cursor Bugbot for commit 119c999. This will update automatically on new commits. Configure here.}

[marcorudolphflex]

@greptile

[greptile-apps]

_{4 files reviewed, 1 comment}

_{Edit Code Review Agent Settings | Greptile}

[marcorudolphflex]

@greptile

[greptile-apps]

_{4 files reviewed, 2 comments}

_{Edit Code Review Agent Settings | Greptile}

[greptile-apps]

_{4 files reviewed, no comments}

_{Edit Code Review Agent Settings | Greptile}

[yaugenst-flex]

Sort of an edge case but nice to have, thanks @marcorudolphflex

[github-actions]

Diff Coverage

Diff: origin/develop...HEAD, staged and unstaged changes

• tidy3d/components/simulation.py (100%)
• tidy3d/components/structure.py (87.5%): Missing lines 336-337,339

Summary

• Total: 25 lines
• Missing: 3 lines
• Coverage: 88%

tidy3d/components/structure.py

Lines 332-343

  332 
  333             intersections = geometry.intersections_plane(**{axis_char: plane_position})
  334             bounds = [shape.bounds for shape in intersections if not shape.is_empty]
  335             if len(bounds) == 0:
! 336                 intersections = geom_box.intersections_plane(**{axis_char: plane_position})
! 337                 bounds = [shape.bounds for shape in intersections if not shape.is_empty]
  338             if len(bounds) == 0:  # fallback
! 339                 return geom_box
  340 
  341             min_plane = (min(b[0] for b in bounds), min(b[1] for b in bounds))
  342             max_plane = (max(b[2] for b in bounds), max(b[3] for b in bounds))

---