diff --git a/src/include/types.h b/src/include/types.h
index 742a5ff8821919259cc8fa80acc63d5b03211e5f..453514b4d89762c779e7f6fb2ba05888244e465a 100644
--- a/src/include/types.h
+++ b/src/include/types.h
@@ -62,7 +62,21 @@ typedef __complex__ double dcomplex;
 #endif // lapack_int
 #endif // np_int
 
-#define imag(z) ( __imag__ (z) )
-#define real(z) ( __real__ (z) )
 #define dconjg(z) ( (__real__ (z) - I * (__imag__ (z))) )
+
+/*! \brief Get the imaginary part of double precision complex number.
+ *
+ * \param z: `const dcomplex &` Reference to the complex number from
+ * which to extract the imaginary part.
+ * \return Im(z): `double` Imaginary part of z.
+ */
+double inline imag(const dcomplex &z) { return __imag__ z; }
+
+/*! \brief Get the imaginary part of double precision complex number.
+ *
+ * \param z: `const dcomplex &` Reference to the complex number from
+ * which to extract the real part.
+ * \return Re(z): `double` Real part of z.
+ */
+double inline real(const dcomplex &z) { return __real__ z; }
 #endif