diff --git a/include/dlaf/auxiliary/norm/api.h b/include/dlaf/auxiliary/norm/api.h
index 6b0be3a2e6..31772ebda5 100644
--- a/include/dlaf/auxiliary/norm/api.h
+++ b/include/dlaf/auxiliary/norm/api.h
@@ -9,6 +9,8 @@
 //
 #pragma once
 
+#include <complex>
+
 #include <pika/execution.hpp>
 
 #include <dlaf/communication/communicator_grid.h>
diff --git a/include/dlaf/eigensolver/band_to_tridiag/api.h b/include/dlaf/eigensolver/band_to_tridiag/api.h
index a61671bfb7..8bae30b155 100644
--- a/include/dlaf/eigensolver/band_to_tridiag/api.h
+++ b/include/dlaf/eigensolver/band_to_tridiag/api.h
@@ -10,6 +10,8 @@
 
 #pragma once
 
+#include <complex>
+
 #include <dlaf/matrix/matrix.h>
 #include <dlaf/types.h>
 
diff --git a/include/dlaf/eigensolver/bt_band_to_tridiag/api.h b/include/dlaf/eigensolver/bt_band_to_tridiag/api.h
index fce15eb4c3..df07fc8662 100644
--- a/include/dlaf/eigensolver/bt_band_to_tridiag/api.h
+++ b/include/dlaf/eigensolver/bt_band_to_tridiag/api.h
@@ -10,6 +10,8 @@
 
 #pragma once
 
+#include <complex>
+
 #include <dlaf/communication/communicator_grid.h>
 #include <dlaf/matrix/matrix.h>
 #include <dlaf/matrix/matrix_ref.h>
diff --git a/include/dlaf/eigensolver/bt_reduction_to_band/api.h b/include/dlaf/eigensolver/bt_reduction_to_band/api.h
index 3f2dddfd52..b261f7c278 100644
--- a/include/dlaf/eigensolver/bt_reduction_to_band/api.h
+++ b/include/dlaf/eigensolver/bt_reduction_to_band/api.h
@@ -9,6 +9,8 @@
 //
 #pragma once
 
+#include <complex>
+
 #include <dlaf/matrix/matrix.h>
 #include <dlaf/matrix/matrix_ref.h>
 #include <dlaf/types.h>
diff --git a/include/dlaf/eigensolver/eigensolver/api.h b/include/dlaf/eigensolver/eigensolver/api.h
index 9b7653e96f..dfe39fe499 100644
--- a/include/dlaf/eigensolver/eigensolver/api.h
+++ b/include/dlaf/eigensolver/eigensolver/api.h
@@ -10,6 +10,8 @@
 
 #pragma once
 
+#include <complex>
+
 #include <dlaf/blas/tile.h>
 #include <dlaf/common/vector.h>
 #include <dlaf/communication/communicator_grid.h>
diff --git a/include/dlaf/eigensolver/gen_eigensolver/api.h b/include/dlaf/eigensolver/gen_eigensolver/api.h
index 431f514953..7f9869796e 100644
--- a/include/dlaf/eigensolver/gen_eigensolver/api.h
+++ b/include/dlaf/eigensolver/gen_eigensolver/api.h
@@ -9,6 +9,8 @@
 //
 #pragma once
 
+#include <complex>
+
 #include <dlaf/communication/communicator_grid.h>
 #include <dlaf/eigensolver/eigensolver/api.h>
 #include <dlaf/matrix/matrix.h>
diff --git a/include/dlaf/eigensolver/gen_to_std/api.h b/include/dlaf/eigensolver/gen_to_std/api.h
index 3965ebd2da..b2d88d2a8d 100644
--- a/include/dlaf/eigensolver/gen_to_std/api.h
+++ b/include/dlaf/eigensolver/gen_to_std/api.h
@@ -9,6 +9,8 @@
 //
 #pragma once
 
+#include <complex>
+
 #include <dlaf/matrix/matrix.h>
 #include <dlaf/types.h>
 
diff --git a/include/dlaf/eigensolver/reduction_to_band/api.h b/include/dlaf/eigensolver/reduction_to_band/api.h
index 5ead6f0325..89233d5c48 100644
--- a/include/dlaf/eigensolver/reduction_to_band/api.h
+++ b/include/dlaf/eigensolver/reduction_to_band/api.h
@@ -10,6 +10,8 @@
 
 #pragma once
 
+#include <complex>
+
 #include <pika/execution.hpp>
 
 #include <dlaf/matrix/matrix.h>
diff --git a/include/dlaf/eigensolver/tridiag_solver/api.h b/include/dlaf/eigensolver/tridiag_solver/api.h
index d2ec6d00db..889a2befba 100644
--- a/include/dlaf/eigensolver/tridiag_solver/api.h
+++ b/include/dlaf/eigensolver/tridiag_solver/api.h
@@ -9,6 +9,8 @@
 //
 #pragma once
 
+#include <complex>
+
 #include <dlaf/matrix/matrix.h>
 #include <dlaf/types.h>
 
diff --git a/include/dlaf/eigensolver/tridiag_solver/gpu/kernels.h b/include/dlaf/eigensolver/tridiag_solver/gpu/kernels.h
index f56a69c8dc..f8be53a83f 100644
--- a/include/dlaf/eigensolver/tridiag_solver/gpu/kernels.h
+++ b/include/dlaf/eigensolver/tridiag_solver/gpu/kernels.h
@@ -12,6 +12,8 @@
 
 #ifdef DLAF_WITH_GPU
 
+#include <complex>
+
 #include <whip.hpp>
 
 #include <dlaf/eigensolver/tridiag_solver/coltype.h>
diff --git a/include/dlaf/eigensolver/tridiag_solver/impl.h b/include/dlaf/eigensolver/tridiag_solver/impl.h
index 38667f1d12..68ae560e89 100644
--- a/include/dlaf/eigensolver/tridiag_solver/impl.h
+++ b/include/dlaf/eigensolver/tridiag_solver/impl.h
@@ -11,6 +11,7 @@
 
 #include <algorithm>
 #include <atomic>
+#include <complex>
 #include <sstream>
 #include <tuple>
 #include <utility>
diff --git a/include/dlaf/eigensolver/tridiag_solver/kernels.h b/include/dlaf/eigensolver/tridiag_solver/kernels.h
index f51edf5fd9..fb6d58389c 100644
--- a/include/dlaf/eigensolver/tridiag_solver/kernels.h
+++ b/include/dlaf/eigensolver/tridiag_solver/kernels.h
@@ -10,6 +10,8 @@
 
 #pragma once
 
+#include <complex>
+
 #include <dlaf/common/callable_object.h>
 #include <dlaf/eigensolver/tridiag_solver/coltype.h>
 #include <dlaf/lapack/tile.h>
diff --git a/include/dlaf/factorization/cholesky/api.h b/include/dlaf/factorization/cholesky/api.h
index 1d5e57472e..4d1f518213 100644
--- a/include/dlaf/factorization/cholesky/api.h
+++ b/include/dlaf/factorization/cholesky/api.h
@@ -9,6 +9,8 @@
 //
 #pragma once
 
+#include <complex>
+
 #include <dlaf/matrix/matrix.h>
 #include <dlaf/types.h>
 
diff --git a/include/dlaf/factorization/qr/api.h b/include/dlaf/factorization/qr/api.h
index 140b0406bc..e4fef98c37 100644
--- a/include/dlaf/factorization/qr/api.h
+++ b/include/dlaf/factorization/qr/api.h
@@ -10,6 +10,8 @@
 
 #pragma once
 
+#include <complex>
+
 #include <pika/execution.hpp>
 
 #include <dlaf/communication/communicator_pipeline.h>
diff --git a/include/dlaf/gpu/blas/gpublas.h b/include/dlaf/gpu/blas/gpublas.h
index 18bb0dfa31..f0730274b3 100644
--- a/include/dlaf/gpu/blas/gpublas.h
+++ b/include/dlaf/gpu/blas/gpublas.h
@@ -13,6 +13,7 @@
 /// Provides gpublas wrappers for BLAS operations.
 
 #ifdef DLAF_WITH_GPU
+#include <complex>
 #include <cstddef>
 #include <utility>
 
diff --git a/include/dlaf/lapack/gpu/add.h b/include/dlaf/lapack/gpu/add.h
index b5a02f58b5..b357c87ef1 100644
--- a/include/dlaf/lapack/gpu/add.h
+++ b/include/dlaf/lapack/gpu/add.h
@@ -12,6 +12,8 @@
 
 #ifdef DLAF_WITH_GPU
 
+#include <complex>
+
 #include <blas.hh>
 #include <whip.hpp>
 
diff --git a/include/dlaf/lapack/gpu/lacpy.h b/include/dlaf/lapack/gpu/lacpy.h
index 92642a43d4..3798ee83e9 100644
--- a/include/dlaf/lapack/gpu/lacpy.h
+++ b/include/dlaf/lapack/gpu/lacpy.h
@@ -12,6 +12,8 @@
 
 #ifdef DLAF_WITH_GPU
 
+#include <complex>
+
 #include <blas.hh>
 #include <whip.hpp>
 
diff --git a/include/dlaf/lapack/gpu/larft.h b/include/dlaf/lapack/gpu/larft.h
index 5e08f59683..34d2035e40 100644
--- a/include/dlaf/lapack/gpu/larft.h
+++ b/include/dlaf/lapack/gpu/larft.h
@@ -12,6 +12,8 @@
 
 #ifdef DLAF_WITH_GPU
 
+#include <complex>
+
 #include <blas.hh>
 #include <whip.hpp>
 
diff --git a/include/dlaf/lapack/gpu/laset.h b/include/dlaf/lapack/gpu/laset.h
index f4c6ea490d..7246cfce90 100644
--- a/include/dlaf/lapack/gpu/laset.h
+++ b/include/dlaf/lapack/gpu/laset.h
@@ -12,6 +12,8 @@
 
 #ifdef DLAF_WITH_GPU
 
+#include <complex>
+
 #include <blas.hh>
 #include <whip.hpp>
 
diff --git a/include/dlaf/lapack/tile.h b/include/dlaf/lapack/tile.h
index 8f36976731..7ef832763a 100644
--- a/include/dlaf/lapack/tile.h
+++ b/include/dlaf/lapack/tile.h
@@ -12,6 +12,7 @@
 
 /// @file
 
+#include <complex>
 #include <cstddef>
 #include <utility>
 
diff --git a/include/dlaf/matrix/matrix.h b/include/dlaf/matrix/matrix.h
index f0c260c240..33e0d07024 100644
--- a/include/dlaf/matrix/matrix.h
+++ b/include/dlaf/matrix/matrix.h
@@ -12,6 +12,7 @@
 
 /// @file
 
+#include <complex>
 #include <cstddef>
 #include <exception>
 #include <utility>
diff --git a/include/dlaf/matrix/matrix_mirror.h b/include/dlaf/matrix/matrix_mirror.h
index 6bd12a8d4f..e46e043d22 100644
--- a/include/dlaf/matrix/matrix_mirror.h
+++ b/include/dlaf/matrix/matrix_mirror.h
@@ -12,6 +12,7 @@
 
 /// @file
 
+#include <complex>
 #include <exception>
 #include <vector>
 
diff --git a/include/dlaf/matrix/matrix_ref.h b/include/dlaf/matrix/matrix_ref.h
index 9d4aae2d9c..af452bf67c 100644
--- a/include/dlaf/matrix/matrix_ref.h
+++ b/include/dlaf/matrix/matrix_ref.h
@@ -12,6 +12,7 @@
 
 /// @file
 
+#include <complex>
 #include <utility>
 
 #include <dlaf/matrix/distribution.h>
diff --git a/include/dlaf/matrix/print_numpy.h b/include/dlaf/matrix/print_numpy.h
index c7a3a98db1..1d5caa6ab5 100644
--- a/include/dlaf/matrix/print_numpy.h
+++ b/include/dlaf/matrix/print_numpy.h
@@ -12,6 +12,7 @@
 
 /// @file
 
+#include <complex>
 #include <iostream>
 #include <sstream>
 #include <string>
diff --git a/include/dlaf/matrix/tile.h b/include/dlaf/matrix/tile.h
index 47e7412524..aec09292b8 100644
--- a/include/dlaf/matrix/tile.h
+++ b/include/dlaf/matrix/tile.h
@@ -12,6 +12,7 @@
 
 /// @file
 
+#include <complex>
 #include <exception>
 #include <memory>
 #include <ostream>
diff --git a/include/dlaf/memory/memory_view.h b/include/dlaf/memory/memory_view.h
index 942cfe404c..8217f39fc1 100644
--- a/include/dlaf/memory/memory_view.h
+++ b/include/dlaf/memory/memory_view.h
@@ -12,6 +12,7 @@
 
 /// @file
 
+#include <complex>
 #include <cstdlib>
 #include <memory>
 #include <utility>
diff --git a/include/dlaf/multiplication/general/api.h b/include/dlaf/multiplication/general/api.h
index 25ee0f0d98..a2c2701212 100644
--- a/include/dlaf/multiplication/general/api.h
+++ b/include/dlaf/multiplication/general/api.h
@@ -10,6 +10,8 @@
 
 #pragma once
 
+#include <complex>
+
 #include <blas.hh>
 
 #include <dlaf/common/index2d.h>
diff --git a/include/dlaf/multiplication/hermitian/api.h b/include/dlaf/multiplication/hermitian/api.h
index b5ed1a04d9..4cfc4bf242 100644
--- a/include/dlaf/multiplication/hermitian/api.h
+++ b/include/dlaf/multiplication/hermitian/api.h
@@ -9,6 +9,8 @@
 //
 #pragma once
 
+#include <complex>
+
 #include <dlaf/matrix/matrix.h>
 #include <dlaf/matrix/matrix_ref.h>
 #include <dlaf/types.h>
diff --git a/include/dlaf/multiplication/triangular/api.h b/include/dlaf/multiplication/triangular/api.h
index f360ff66cf..ca1ce48db3 100644
--- a/include/dlaf/multiplication/triangular/api.h
+++ b/include/dlaf/multiplication/triangular/api.h
@@ -9,6 +9,8 @@
 //
 #pragma once
 
+#include <complex>
+
 #include <blas.hh>
 
 #include <dlaf/matrix/matrix.h>
diff --git a/include/dlaf/permutations/general/perms.h b/include/dlaf/permutations/general/perms.h
index 149b736fc0..9fc2f7a30c 100644
--- a/include/dlaf/permutations/general/perms.h
+++ b/include/dlaf/permutations/general/perms.h
@@ -12,6 +12,7 @@
 
 #ifdef DLAF_WITH_GPU
 
+#include <complex>
 #include <vector>
 
 #include <whip.hpp>
diff --git a/include/dlaf/solver/triangular/api.h b/include/dlaf/solver/triangular/api.h
index 4d64eb2496..ca3212ccca 100644
--- a/include/dlaf/solver/triangular/api.h
+++ b/include/dlaf/solver/triangular/api.h
@@ -9,6 +9,8 @@
 //
 #pragma once
 
+#include <complex>
+
 #include <blas.hh>
 
 #include <dlaf/matrix/matrix.h>
diff --git a/include/dlaf/util_matrix.h b/include/dlaf/util_matrix.h
index 939249ebef..582c0fe088 100644
--- a/include/dlaf/util_matrix.h
+++ b/include/dlaf/util_matrix.h
@@ -10,6 +10,7 @@
 #pragma once
 
 #include <cmath>
+#include <complex>
 #include <cstddef>
 #include <optional>
 #include <random>
diff --git a/include/dlaf/util_rocblas.h b/include/dlaf/util_rocblas.h
index 932b3d5758..7668d82bce 100644
--- a/include/dlaf/util_rocblas.h
+++ b/include/dlaf/util_rocblas.h
@@ -12,6 +12,8 @@
 
 #ifdef DLAF_WITH_HIP
 
+#include <complex>
+
 #include <blas.hh>
 
 #include <dlaf/gpu/blas/api.h>
diff --git a/src/auxiliary/norm/mc.cpp b/src/auxiliary/norm/mc.cpp
index e10b4045e9..57b616f09d 100644
--- a/src/auxiliary/norm/mc.cpp
+++ b/src/auxiliary/norm/mc.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/auxiliary/norm/mc.h>
 
 namespace dlaf::auxiliary::internal {
diff --git a/src/blas/scal.cpp b/src/blas/scal.cpp
index 6ace168dbb..bbd0a869b9 100644
--- a/src/blas/scal.cpp
+++ b/src/blas/scal.cpp
@@ -8,6 +8,7 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
 #include <cstdint>
 
 #include <blas/config.h>
diff --git a/src/c_api/factorization/cholesky.cpp b/src/c_api/factorization/cholesky.cpp
index f930284496..69ae1901b3 100644
--- a/src/c_api/factorization/cholesky.cpp
+++ b/src/c_api/factorization/cholesky.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf_c/factorization/cholesky.h>
 #include <dlaf_c/utils.h>
 
diff --git a/src/communication/datatypes.cpp b/src/communication/datatypes.cpp
index 1cf4f9e042..b5be924d3f 100644
--- a/src/communication/datatypes.cpp
+++ b/src/communication/datatypes.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/communication/datatypes.h>
 
 namespace dlaf {
diff --git a/src/eigensolver/band_to_tridiag/mc.cpp b/src/eigensolver/band_to_tridiag/mc.cpp
index 1261a01fe0..ac92b0118f 100644
--- a/src/eigensolver/band_to_tridiag/mc.cpp
+++ b/src/eigensolver/band_to_tridiag/mc.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/eigensolver/band_to_tridiag/mc.h>
 
 namespace dlaf::eigensolver::internal {
diff --git a/src/eigensolver/bt_band_to_tridiag/mc.cpp b/src/eigensolver/bt_band_to_tridiag/mc.cpp
index 3f288b514b..d411deabda 100644
--- a/src/eigensolver/bt_band_to_tridiag/mc.cpp
+++ b/src/eigensolver/bt_band_to_tridiag/mc.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/eigensolver/bt_band_to_tridiag/impl.h>
 
 namespace dlaf::eigensolver::internal {
diff --git a/src/eigensolver/bt_reduction_to_band/gpu.cpp b/src/eigensolver/bt_reduction_to_band/gpu.cpp
index aea50bf041..0bf72663b2 100644
--- a/src/eigensolver/bt_reduction_to_band/gpu.cpp
+++ b/src/eigensolver/bt_reduction_to_band/gpu.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/eigensolver/bt_reduction_to_band/impl.h>
 
 namespace dlaf::eigensolver::internal {
diff --git a/src/eigensolver/bt_reduction_to_band/mc.cpp b/src/eigensolver/bt_reduction_to_band/mc.cpp
index 1b01e54934..e337ae2338 100644
--- a/src/eigensolver/bt_reduction_to_band/mc.cpp
+++ b/src/eigensolver/bt_reduction_to_band/mc.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/eigensolver/bt_reduction_to_band/impl.h>
 
 namespace dlaf::eigensolver::internal {
diff --git a/src/eigensolver/eigensolver/gpu.cpp b/src/eigensolver/eigensolver/gpu.cpp
index 52da7518c2..1cad3840dc 100644
--- a/src/eigensolver/eigensolver/gpu.cpp
+++ b/src/eigensolver/eigensolver/gpu.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/eigensolver/eigensolver/impl.h>
 
 namespace dlaf::eigensolver::internal {
diff --git a/src/eigensolver/eigensolver/mc.cpp b/src/eigensolver/eigensolver/mc.cpp
index d65f41790d..932411d461 100644
--- a/src/eigensolver/eigensolver/mc.cpp
+++ b/src/eigensolver/eigensolver/mc.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/eigensolver/eigensolver/impl.h>
 
 namespace dlaf::eigensolver::internal {
diff --git a/src/eigensolver/gen_eigensolver/gpu.cpp b/src/eigensolver/gen_eigensolver/gpu.cpp
index eccbbda5dd..3a8061fe20 100644
--- a/src/eigensolver/gen_eigensolver/gpu.cpp
+++ b/src/eigensolver/gen_eigensolver/gpu.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/eigensolver/gen_eigensolver/impl.h>
 
 namespace dlaf::eigensolver::internal {
diff --git a/src/eigensolver/gen_eigensolver/mc.cpp b/src/eigensolver/gen_eigensolver/mc.cpp
index 0d93a6362d..7f480ad0e7 100644
--- a/src/eigensolver/gen_eigensolver/mc.cpp
+++ b/src/eigensolver/gen_eigensolver/mc.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/eigensolver/gen_eigensolver/impl.h>
 
 namespace dlaf::eigensolver::internal {
diff --git a/src/eigensolver/gen_to_std/gpu.cpp b/src/eigensolver/gen_to_std/gpu.cpp
index f6fe03f6c0..d0e9eec3b7 100644
--- a/src/eigensolver/gen_to_std/gpu.cpp
+++ b/src/eigensolver/gen_to_std/gpu.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/eigensolver/gen_to_std/impl.h>
 
 namespace dlaf {
diff --git a/src/eigensolver/gen_to_std/mc.cpp b/src/eigensolver/gen_to_std/mc.cpp
index a6a3c5aba1..854ee72620 100644
--- a/src/eigensolver/gen_to_std/mc.cpp
+++ b/src/eigensolver/gen_to_std/mc.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/eigensolver/gen_to_std/impl.h>
 
 namespace dlaf {
diff --git a/src/eigensolver/reduction_to_band/gpu.cpp b/src/eigensolver/reduction_to_band/gpu.cpp
index 7e0cce7e77..3ab273b07c 100644
--- a/src/eigensolver/reduction_to_band/gpu.cpp
+++ b/src/eigensolver/reduction_to_band/gpu.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/eigensolver/reduction_to_band/impl.h>
 
 namespace dlaf::eigensolver::internal {
diff --git a/src/eigensolver/reduction_to_band/mc.cpp b/src/eigensolver/reduction_to_band/mc.cpp
index 93c8be8806..0cda66bab4 100644
--- a/src/eigensolver/reduction_to_band/mc.cpp
+++ b/src/eigensolver/reduction_to_band/mc.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/eigensolver/reduction_to_band/impl.h>
 
 namespace dlaf::eigensolver::internal {
diff --git a/src/eigensolver/tridiag_solver/gpu/kernels.cu b/src/eigensolver/tridiag_solver/gpu/kernels.cu
index 3ef6cf918d..e574e44396 100644
--- a/src/eigensolver/tridiag_solver/gpu/kernels.cu
+++ b/src/eigensolver/tridiag_solver/gpu/kernels.cu
@@ -8,6 +8,7 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
 #include <type_traits>
 
 #include <whip.hpp>
diff --git a/src/factorization/cholesky/gpu.cpp b/src/factorization/cholesky/gpu.cpp
index d923da217d..b6bbafda65 100644
--- a/src/factorization/cholesky/gpu.cpp
+++ b/src/factorization/cholesky/gpu.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/factorization/cholesky/impl.h>
 
 namespace dlaf {
diff --git a/src/factorization/cholesky/mc.cpp b/src/factorization/cholesky/mc.cpp
index da181dbdba..12665bb83e 100644
--- a/src/factorization/cholesky/mc.cpp
+++ b/src/factorization/cholesky/mc.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/factorization/cholesky/impl.h>
 
 namespace dlaf {
diff --git a/src/factorization/qr/gpu.cpp b/src/factorization/qr/gpu.cpp
index 603a0977a3..74a7a61ba3 100644
--- a/src/factorization/qr/gpu.cpp
+++ b/src/factorization/qr/gpu.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/factorization/qr/t_factor_impl.h>
 
 namespace dlaf::factorization::internal {
diff --git a/src/factorization/qr/mc.cpp b/src/factorization/qr/mc.cpp
index 614d2d2153..ceb81f3bee 100644
--- a/src/factorization/qr/mc.cpp
+++ b/src/factorization/qr/mc.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/factorization/qr/t_factor_impl.h>
 
 namespace dlaf::factorization::internal {
diff --git a/src/lapack/gpu/add.cu b/src/lapack/gpu/add.cu
index 4ef3329b64..3494ff4aac 100644
--- a/src/lapack/gpu/add.cu
+++ b/src/lapack/gpu/add.cu
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <whip.hpp>
 
 #include <dlaf/gpu/assert.cu.h>
diff --git a/src/lapack/gpu/lacpy.cu b/src/lapack/gpu/lacpy.cu
index b5abd62211..597b7ff6cb 100644
--- a/src/lapack/gpu/lacpy.cu
+++ b/src/lapack/gpu/lacpy.cu
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <whip.hpp>
 
 #include <dlaf/common/assert.h>
diff --git a/src/lapack/gpu/larft.cu b/src/lapack/gpu/larft.cu
index aa34d7e9c6..932774244c 100644
--- a/src/lapack/gpu/larft.cu
+++ b/src/lapack/gpu/larft.cu
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <whip.hpp>
 
 #include <dlaf/common/assert.h>
diff --git a/src/lapack/gpu/laset.cu b/src/lapack/gpu/laset.cu
index c0e68ca1d3..f66a051520 100644
--- a/src/lapack/gpu/laset.cu
+++ b/src/lapack/gpu/laset.cu
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <whip.hpp>
 
 #include <dlaf/gpu/assert.cu.h>
diff --git a/test/unit/communication/test_all_reduce.cpp b/test/unit/communication/test_all_reduce.cpp
index cab786c2fc..2273ed3d48 100644
--- a/test/unit/communication/test_all_reduce.cpp
+++ b/test/unit/communication/test_all_reduce.cpp
@@ -8,6 +8,7 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
 #include <utility>
 
 #include <dlaf/common/data_descriptor.h>
diff --git a/test/unit/communication/test_broadcast_matrix.cpp b/test/unit/communication/test_broadcast_matrix.cpp
index 4071402417..06e694182e 100644
--- a/test/unit/communication/test_broadcast_matrix.cpp
+++ b/test/unit/communication/test_broadcast_matrix.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/communication/functions_sync.h>
 #include <dlaf/matrix/matrix.h>
 
diff --git a/test/unit/communication/test_broadcast_tile.cpp b/test/unit/communication/test_broadcast_tile.cpp
index e273db918a..c0c3fd9042 100644
--- a/test/unit/communication/test_broadcast_tile.cpp
+++ b/test/unit/communication/test_broadcast_tile.cpp
@@ -8,6 +8,8 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
+
 #include <dlaf/communication/functions_sync.h>
 #include <dlaf/matrix/tile.h>
 
diff --git a/test/unit/communication/test_reduce.cpp b/test/unit/communication/test_reduce.cpp
index 6e10b47c19..874d0a9abc 100644
--- a/test/unit/communication/test_reduce.cpp
+++ b/test/unit/communication/test_reduce.cpp
@@ -8,6 +8,7 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
 #include <utility>
 
 #include <dlaf/common/data_descriptor.h>
diff --git a/test/unit/matrix/test_matrix_local.cpp b/test/unit/matrix/test_matrix_local.cpp
index 29eef594b6..c1594bb09a 100644
--- a/test/unit/matrix/test_matrix_local.cpp
+++ b/test/unit/matrix/test_matrix_local.cpp
@@ -8,6 +8,7 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
 #include <sstream>
 #include <string>
 #include <utility>
diff --git a/test/unit/matrix/test_matrix_output.cpp b/test/unit/matrix/test_matrix_output.cpp
index 34cf1466b0..4d95ff8751 100644
--- a/test/unit/matrix/test_matrix_output.cpp
+++ b/test/unit/matrix/test_matrix_output.cpp
@@ -8,6 +8,7 @@
 // SPDX-License-Identifier: BSD-3-Clause
 //
 
+#include <complex>
 #include <sstream>
 #include <string>
 #include <tuple>