CarpetIO*: Output complex numbers as pair of real numbers

Don't use the (real,imag) notation that C++ would use by default.

1 files changed, 11 insertions, 0 deletions

diff --git a/Carpet/CarpetIOScalar/src/ioscalar.cc b/Carpet/CarpetIOScalar/src/ioscalar.cc
index 1cd0ee102..cfdf1843b 100644
--- a/Carpet/CarpetIOScalar/src/ioscalar.cc
+++ b/Carpet/CarpetIOScalar/src/ioscalar.cc
@@ -471,12 +471,23 @@ namespace CarpetIOScalar {
             file << " ";
             
             switch (specific_cactus_type(vartype)) {
+#define CARPET_NO_COMPLEX
 #define TYPECASE(N,T)                           \
               case N:                           \
                 file << *(T const*)result;      \
               break;
 #include "typecase.hh"
 #undef TYPECASE
+#undef CARPET_NO_COMPLEX
+#define CARPET_COMPLEX
+#define TYPECASE(N,T)                                   \
+              case N:                                   \
+                file << real(*(T const*)result) << " "  \
+                     << imag(*(T const*)result);        \
+              break;
+#include "typecase.hh"
+#undef TYPECASE
+#undef CARPET_COMPLEX
             default:
               UnsupportedVarType (n);
             }