Avoid using complex number literal extension, instead relying on

standard-defined memory layout.

Author: BenBrock

include/binsparse/detail/declamp_values.h

@@ -42,7 +42,12 @@ static inline void bsp_matrix_declamp_values(bsp_matrix_t matrix) {
       real_value = bsp_suitesparse_declamp_impl_(real_value);
       imaginary_value = bsp_suitesparse_declamp_impl_(imaginary_value);
 
-      values[i] = real_value + 1j * imaginary_value;
+      double _Complex complex_value;
+
+      ((double*) &complex_value)[0] = real_value;
+      ((double*) &complex_value)[1] = imaginary_value;
+
+      values[i] = complex_value;
     }
   }
 }

include/binsparse/matrix_market/matrix_market_read.h

@@ -89,7 +89,10 @@ static inline bsp_matrix_t bsp_mmread_explicit_array(const char* file_path,
       double real_value, imaginary_value;
       sscanf(buf, "%lf %lf", &real_value, &imaginary_value);
 
-      double _Complex value = real_value + 1j * imaginary_value;
+      double _Complex value;
+
+      ((double*) &value)[0] = real_value;
+      ((double*) &value)[1] = imaginary_value;
 
       size_t i = count % matrix.ncols;
       size_t j = count / matrix.ncols;
@@ -203,12 +206,15 @@ bsp_mmread_explicit_coordinate(const char* file_path, bsp_type_t value_type,
       bsp_array_write(matrix.indices_1, count, j);
     } else if (mm_type == BSP_MM_COMPLEX) {
       unsigned long long i, j;
-      double real_value, complex_value;
-      sscanf(buf, "%llu %llu %lf %lf", &i, &j, &real_value, &complex_value);
+      double real_value, imaginary_value;
+      sscanf(buf, "%llu %llu %lf %lf", &i, &j, &real_value, &imaginary_value);
       i--;
       j--;
 
-      double _Complex value = real_value + 1j * complex_value;
+      double _Complex value;
+
+      ((double*) &value)[0] = real_value;
+      ((double*) &value)[1] = imaginary_value;
 
       bsp_array_write(matrix.values, count, value);
       bsp_array_write(matrix.indices_0, count, i);