diff --git a/src/tribol/search/InterfacePairFinder.cpp b/src/tribol/search/InterfacePairFinder.cpp
index 40b7fb09..fb971939 100644
--- a/src/tribol/search/InterfacePairFinder.cpp
+++ b/src/tribol/search/InterfacePairFinder.cpp
@@ -126,7 +126,9 @@ class CartesianProduct : public SearchBase {
     auto& contactPairs = m_coupling_scheme->getInterfacePairs();
     contactPairs.resize( countArray_host[0] );
 
-    countArray.fill( 0 );
+    int zero = 0;
+    // Workaround for axom::Array::fill() issue for optimized CUDA code, see Axom #1833.
+    axom::copy( countArray.data(), &zero, sizeof( int ) );
     auto pairs_view = m_coupling_scheme->getInterfacePairs().view();
     // fill proximate pairs array
     forAllExec(