Squashed commit of the following:

commit 0e1f80c
Author: xhuang-jpl <118782850+xhuang-jpl@users.noreply.github.com>
Date:   Wed Jun 24 16:38:22 2026 -0700

    Fix the Ampcor issues from Lijun's branch and increase the limit of the search radius  (#282)

    * sm r3.4

    * change the SM commit id for R4.0.2

    * remove the mask description for the wrapped ifgram

    * remove the tmp folder

    * fix minors

    * comment out the checkPixelInImageRange

    * pycuampcor: fix the pixels out of range issue

    * pycuampcor: fix an error in determining the corrwindow size

    * update the maximum search radius to be 256

    ---------

    Co-authored-by: Xiaodong Huang <xhuang@nisar-adt-dev-3.jpl.nasa.gov>
    Co-authored-by: Lijun Zhu <ljzhu@gps.caltech.edu>

commit 813d729
Author: Brian Hawkins <1729052+bhawkins@users.noreply.github.com>
Date:   Wed Jun 24 13:56:09 2026 -0500

    Handle case where no samples in a rangline are valid (#327)

    * Handle case where all subswaths are invalid.

    * Log some details about missing data.

    * Fix typo

    * Use signed integer for negative values...

commit 075795f
Author: Gustavo H. X. Shiroma <52007211+gshiroma@users.noreply.github.com>
Date:   Thu Jun 18 15:59:58 2026 -0700

    Fix STATIC filenames with fractional posting by rounding posting after converting to decimeters. (#314)

    * disable polarimetric symmetrization by default

    * revert changes to `symmetrize_cross_pol_channels`

    * Update GCOV and GSLC specification XMLs

    * Revert changes to the GCOV and GSLC specification XMLs

    * round posting after conversion to decimeters

commit 3f7027d
Author: SamNemo <11642807+nemo794@users.noreply.github.com>
Date:   Thu Jun 18 10:29:43 2026 -0700

    Fix text encoding: Replace en dash (U+2013) with ASCII hyphen in descriptions (#292)

    Replaced Unicode en dash character (U+2013 "–") with standard ASCII hyphen
    (U+002d "-") in 19 locations across 4 Python files. The en dash was causing
    display issues in some programs, appearing as "â€"" instead of "–".

    Files modified:
    - python/packages/isce3/signal/rfi_freq_null.py (1 occurrence in citation)
    - python/packages/nisar/products/insar/GUNW_writer.py (5 occurrences in mask descriptions)
    - python/packages/nisar/products/insar/GOFF_writer.py (4 occurrences in mask descriptions)
    - python/packages/nisar/products/insar/InSAR_L1_writer.py (9 occurrences in mask descriptions)

    This fix ensures HDF5 attribute descriptions display correctly across all
    programs and avoids confusion with special Unicode characters in source code.

    Co-authored-by: Samantha C. Niemoeller <samantha.c.niemoeller@jpl.nasa.gov>

commit db4f65b
Author: Gustavo H. X. Shiroma <52007211+gshiroma@users.noreply.github.com>
Date:   Wed Jun 17 14:59:11 2026 -0700

    Add margins to the STATIC product Doppler LUT grid boundaries to ensure full radargrid coverage. (#298)

    * disable polarimetric symmetrization by default

    * revert changes to `symmetrize_cross_pol_channels`

    * Update GCOV and GSLC specification XMLs

    * Revert changes to the GCOV and GSLC specification XMLs

    * Add margins to the Doppler LUT grid boundaries to ensure full radargrid coverage.

    * Add margins to the Doppler LUT grid boundaries to ensure full radargrid coverage.

    * Add margins to the Doppler LUT grid boundaries to ensure full radargrid coverage.

commit 67f0878
Author: Gustavo H. X. Shiroma <52007211+gshiroma@users.noreply.github.com>
Date:   Mon Jun 8 11:43:17 2026 -0700

    Improve DEM bounding box selection when product coordinates are in a different projection than the DEM (#275)

    * disable polarimetric symmetrization by default

    * revert changes to `symmetrize_cross_pol_channels`

    * Update GCOV and GSLC specification XMLs

    * Revert changes to the GCOV and GSLC specification XMLs

    * Update DEM bounding box selection

    * Improve variable names, docstrings, and handle incorrect starting/ending ordering

    * Improve variable names

    * Ensure latitude values remain within [-90, 90]

    * Ensure latitude values remain within [-90, 90]

    * Simplify code

    * Use parenthesis to make it clear evaluation order

    * improve comments

    * revert changes to constructor of variables `all_x` and `all_y`

    * change margin types from `float` to `double`

Author: Tyler-g-hudson

.github/workflows/build-and-run.yml

@@ -92,11 +92,6 @@ jobs:
         run: |
           micromamba install -c nvidia cuda
 
-      - name: Install CUDA Toolkit
-        if: matrix.deps.label == 'Documentation'
-        run: |
-          micromamba install -c nvidia cuda
-
       # Fix missing conda command on macOS
       - name: macOS conda command fixup
         if: runner.os == 'macOS'

cxx/isce3/cuda/matchtemplate/pycuampcor/cuAmpcorParameter.cu

@@ -88,8 +88,8 @@ void cuAmpcorParameter::setupParameters()
     windowSizeWidth = windowSizeWidthRaw*rawDataOversamplingFactor;  //
     windowSizeHeight = windowSizeHeightRaw*rawDataOversamplingFactor;
 
-    searchWindowSizeWidthRaw =  windowSizeWidthRaw + 2*halfSearchRangeDownRaw;
-    searchWindowSizeHeightRaw = windowSizeHeightRaw + 2*halfSearchRangeAcrossRaw;
+    searchWindowSizeWidthRaw =  windowSizeWidthRaw + 2*halfSearchRangeAcrossRaw;
+    searchWindowSizeHeightRaw = windowSizeHeightRaw + 2*halfSearchRangeDownRaw;
 
     searchWindowSizeWidthRawZoomIn = windowSizeWidthRaw + 2*halfZoomWindowSizeRaw;
     searchWindowSizeHeightRawZoomIn = windowSizeHeightRaw + 2*halfZoomWindowSizeRaw;
@@ -237,14 +237,28 @@ void cuAmpcorParameter::setChunkStartPixels()
                     if(sChunkEA < vpixel) sChunkEA = vpixel;
                 }
             }
+            // Clamp chunk start pixels to valid image range
+            if (mChunkSD < 0) mChunkSD = 0;
+            if (mChunkSA < 0) mChunkSA = 0;
+            if (sChunkSD < 0) sChunkSD = 0;
+            if (sChunkSA < 0) sChunkSA = 0;
+            // Compute end pixels and clamp to image range
+            int mEndD = mChunkED + windowSizeHeightRaw;
+            if (mEndD > referenceImageHeight) mEndD = referenceImageHeight;
+            int mEndA = mChunkEA + windowSizeWidthRaw;
+            if (mEndA > referenceImageWidth) mEndA = referenceImageWidth;
+            int sEndD = sChunkED + searchWindowSizeHeightRaw;
+            if (sEndD > secondaryImageHeight) sEndD = secondaryImageHeight;
+            int sEndA = sChunkEA + searchWindowSizeWidthRaw;
+            if (sEndA > secondaryImageWidth) sEndA = secondaryImageWidth;
             referenceChunkStartPixelDown[idxChunk]   = mChunkSD;
             referenceChunkStartPixelAcross[idxChunk] = mChunkSA;
             secondaryChunkStartPixelDown[idxChunk]    = sChunkSD;
             secondaryChunkStartPixelAcross[idxChunk]  = sChunkSA;
-            referenceChunkHeight[idxChunk] = mChunkED - mChunkSD + windowSizeHeightRaw;
-            referenceChunkWidth[idxChunk]  = mChunkEA - mChunkSA + windowSizeWidthRaw;
-            secondaryChunkHeight[idxChunk]  = sChunkED - sChunkSD + searchWindowSizeHeightRaw;
-            secondaryChunkWidth[idxChunk]   = sChunkEA - sChunkSA + searchWindowSizeWidthRaw;
+            referenceChunkHeight[idxChunk] = mEndD - mChunkSD;
+            referenceChunkWidth[idxChunk]  = mEndA - mChunkSA;
+            secondaryChunkHeight[idxChunk]  = sEndD - sChunkSD;
+            secondaryChunkWidth[idxChunk]   = sEndA - sChunkSA;
             if(maxReferenceChunkHeight < referenceChunkHeight[idxChunk]) maxReferenceChunkHeight = referenceChunkHeight[idxChunk];
             if(maxReferenceChunkWidth  < referenceChunkWidth[idxChunk] ) maxReferenceChunkWidth  = referenceChunkWidth[idxChunk];
             if(maxSecondaryChunkHeight  < secondaryChunkHeight[idxChunk]) maxSecondaryChunkHeight = secondaryChunkHeight[idxChunk];
@@ -256,60 +270,37 @@ void cuAmpcorParameter::setChunkStartPixels()
 /// check whether reference and secondary windows are within the image range
 void cuAmpcorParameter::checkPixelInImageRange()
 {
+#ifdef CUAMPCOR_DEBUG
     int endPixel;
     for(int row=0; row<numberWindowDown; row++)
     {
         for(int col = 0; col < numberWindowAcross; col++)
         {
             int i = row*numberWindowAcross + col;
-            if(referenceStartPixelDown[i] <0)
-            {
-                fprintf(stderr, "Reference Window start pixel out of range in Down, window (%d,%d), pixel %d\n", row, col, referenceStartPixelDown[i]);
-                exit(EXIT_FAILURE); //or raise range error
-            }
-            if(referenceStartPixelAcross[i] <0)
-            {
-                fprintf(stderr, "Reference Window start pixel out of range in Across, window (%d,%d), pixel %d\n", row, col, referenceStartPixelAcross[i]);
-                exit(EXIT_FAILURE);
-            }
+            if(referenceStartPixelDown[i] < 0)
+                fprintf(stderr, "Warning: Reference Window start pixel out of range in Down, window (%d,%d), pixel %d\n", row, col, referenceStartPixelDown[i]);
+            if(referenceStartPixelAcross[i] < 0)
+                fprintf(stderr, "Warning: Reference Window start pixel out of range in Across, window (%d,%d), pixel %d\n", row, col, referenceStartPixelAcross[i]);
             endPixel = referenceStartPixelDown[i] + windowSizeHeightRaw;
             if(endPixel >= referenceImageHeight)
-            {
-                fprintf(stderr, "Reference Window end pixel out of range in Down, window (%d,%d), pixel %d\n", row, col, endPixel);
-                exit(EXIT_FAILURE);
-            }
+                fprintf(stderr, "Warning: Reference Window end pixel out of range in Down, window (%d,%d), pixel %d\n", row, col, endPixel);
             endPixel = referenceStartPixelAcross[i] + windowSizeWidthRaw;
             if(endPixel >= referenceImageWidth)
-            {
-                fprintf(stderr, "Reference Window end pixel out of range in Across, window (%d,%d), pixel %d\n", row, col, endPixel);
-                exit(EXIT_FAILURE);
-            }
+                fprintf(stderr, "Warning: Reference Window end pixel out of range in Across, window (%d,%d), pixel %d\n", row, col, endPixel);
             //secondary
-            if(secondaryStartPixelDown[i] <0)
-            {
-                fprintf(stderr, "Secondary Window start pixel out of range in Down, window (%d,%d), pixel %d\n", row, col, secondaryStartPixelDown[i]);
-                exit(EXIT_FAILURE);
-            }
-            if(secondaryStartPixelAcross[i] <0)
-            {
-                fprintf(stderr, "Secondary Window start pixel out of range in Across, window (%d,%d), pixel %d\n", row, col, secondaryStartPixelAcross[i]);
-                exit(EXIT_FAILURE);
-            }
+            if(secondaryStartPixelDown[i] < 0)
+                fprintf(stderr, "Warning: Secondary Window start pixel out of range in Down, window (%d,%d), pixel %d\n", row, col, secondaryStartPixelDown[i]);
+            if(secondaryStartPixelAcross[i] < 0)
+                fprintf(stderr, "Warning: Secondary Window start pixel out of range in Across, window (%d,%d), pixel %d\n", row, col, secondaryStartPixelAcross[i]);
             endPixel = secondaryStartPixelDown[i] + searchWindowSizeHeightRaw;
             if(endPixel >= secondaryImageHeight)
-            {
-                fprintf(stderr, "Secondary Window end pixel out of range in Down, window (%d,%d), pixel %d\n", row, col, endPixel);
-                exit(EXIT_FAILURE);
-            }
+                fprintf(stderr, "Warning: Secondary Window end pixel out of range in Down, window (%d,%d), pixel %d\n", row, col, endPixel);
             endPixel = secondaryStartPixelAcross[i] + searchWindowSizeWidthRaw;
             if(endPixel >= secondaryImageWidth)
-            {
-                fprintf(stderr, "Secondary Window end pixel out of range in Across, window (%d,%d), pixel %d\n", row, col, endPixel);
-                exit(EXIT_FAILURE);
-            }
-
+                fprintf(stderr, "Warning: Secondary Window end pixel out of range in Across, window (%d,%d), pixel %d\n", row, col, endPixel);
         }
     }
+#endif
 }
 
 

cxx/isce3/cuda/matchtemplate/pycuampcor/cuArraysCopy.cu

@@ -68,7 +68,7 @@ void cuArraysCopyToBatch(cuArrays<float2> *image1, cuArrays<float2> *image2,
 
 // kernel for cuArraysCopyToBatchWithOffset
 template<typename T_in, typename T_out>
-__global__ void cuArraysCopyToBatchWithOffset_kernel(const T_in *imageIn, const int inNY,
+__global__ void cuArraysCopyToBatchWithOffset_kernel(const T_in *imageIn, const int inNX, const int inNY,
     T_out *imageOut, const int outNX, const int outNY, const int nImages,
     const int *offsetX, const int *offsetY)
 {
@@ -77,8 +77,12 @@ __global__ void cuArraysCopyToBatchWithOffset_kernel(const T_in *imageIn, const
     int outy = threadIdx.y + blockDim.y*blockIdx.y;
     if(idxImage>=nImages || outx >= outNX || outy >= outNY) return;
     int idxOut = idxImage*outNX*outNY + outx*outNY + outy;
-    int idxIn = (offsetX[idxImage]+outx)*inNY + offsetY[idxImage] + outy;
-    imageOut[idxOut] = T_out{imageIn[idxIn]};
+    int inx = offsetX[idxImage] + outx;
+    int iny = offsetY[idxImage] + outy;
+    if (inx >= 0 && inx < inNX && iny >= 0 && iny < inNY)
+        imageOut[idxOut] = T_out{imageIn[inx*inNY + iny]};
+    else
+        imageOut[idxOut] = T_out{0.0f};
 }
 
 /**
@@ -98,14 +102,14 @@ void cuArraysCopyToBatchWithOffset(cuArrays<float2> *image1, const int lda1, cuA
     dim3 blockSize(nthreads, nthreads, 1);
     dim3 gridSize(IDIVUP(image2->height,nthreads), IDIVUP(image2->width,nthreads), image2->count);
     cuArraysCopyToBatchWithOffset_kernel<<<gridSize,blockSize, 0 , stream>>> (
-        image1->devData, lda1,
+        image1->devData, image1->height, lda1,
         image2->devData, image2->height, image2->width, image2->count,
         offsetH, offsetW);
     getLastCudaError("cuArraysCopyToBatchAbsWithOffset_kernel");
 }
 
 // same as above, but from complex to real(take amplitudes)
-__global__ void cuArraysCopyToBatchAbsWithOffset_kernel(const float2 *imageIn, const int inNY,
+__global__ void cuArraysCopyToBatchAbsWithOffset_kernel(const float2 *imageIn, const int inNX, const int inNY,
     float2 *imageOut, const int outNX, const int outNY, const int nImages,
     const int *offsetX, const int *offsetY)
 {
@@ -114,8 +118,12 @@ __global__ void cuArraysCopyToBatchAbsWithOffset_kernel(const float2 *imageIn, c
     int outy = threadIdx.y + blockDim.y*blockIdx.y;
     if(idxImage>=nImages || outx >= outNX || outy >= outNY) return;
     int idxOut = idxImage*outNX*outNY + outx*outNY + outy;
-    int idxIn = (offsetX[idxImage]+outx)*inNY + offsetY[idxImage] + outy;
-    imageOut[idxOut] = make_float2(complexAbs(imageIn[idxIn]), 0.0);
+    int inx = offsetX[idxImage] + outx;
+    int iny = offsetY[idxImage] + outy;
+    if (inx >= 0 && inx < inNX && iny >= 0 && iny < inNY)
+        imageOut[idxOut] = make_float2(complexAbs(imageIn[inx*inNY + iny]), 0.0f);
+    else
+        imageOut[idxOut] = make_float2(0.0f, 0.0f);
 }
 
 /**
@@ -135,7 +143,7 @@ void cuArraysCopyToBatchAbsWithOffset(cuArrays<float2> *image1, const int lda1,
     dim3 blockSize(nthreads, nthreads, 1);
     dim3 gridSize(IDIVUP(image2->height,nthreads), IDIVUP(image2->width,nthreads), image2->count);
     cuArraysCopyToBatchAbsWithOffset_kernel<<<gridSize,blockSize, 0 , stream>>> (
-        image1->devData, lda1,
+        image1->devData, image1->height, lda1,
         image2->devData, image2->height, image2->width, image2->count,
         offsetH, offsetW);
     getLastCudaError("cuArraysCopyToBatchAbsWithOffset_kernel");
@@ -158,7 +166,7 @@ void cuArraysCopyToBatchWithOffsetR2C(cuArrays<float> *image1, const int lda1, c
     dim3 blockSize(nthreads, nthreads, 1);
     dim3 gridSize(IDIVUP(image2->height,nthreads), IDIVUP(image2->width,nthreads), image2->count);
     cuArraysCopyToBatchWithOffset_kernel<<<gridSize,blockSize, 0 , stream>>> (
-        image1->devData, lda1,
+        image1->devData, image1->height, lda1,
         image2->devData, image2->height, image2->width, image2->count,
         offsetH, offsetW);
     getLastCudaError("cuArraysCopyToBatchWithOffsetR2C_kernel");

cxx/isce3/geometry/RTC.cpp

@@ -1174,7 +1174,7 @@ void _RunBlock(const int jmax, const int block_size,
     the bounding-box corners.
     */
     int dem_margin_x_in_pixels = 100;
-    int dem_margin_y_in_pixels = 200;
+    int dem_margin_y_in_pixels = 100;
     auto error_code = loadDemFromProj(
             dem_raster, minX, maxX, minY, maxY, &dem_interp_block, proj,
             dem_margin_x_in_pixels, dem_margin_y_in_pixels);