Compile c code with float instead of double

Question

I have a lengthy numeric integration scheme written in C. I\'d like to test my algorithm in floating point precision. Is there a way to tell gcc to demote every occurrence o

Answer 1

typedef double float; 

Before any doubles that you want to replace should work, however be warned it may confuse some external libraries.

In the future, the best approach is to define your own float type:

#ifdef USE_FLOATS
    typedef float MyFloatType; 
#else 
    typedef double MyFloatType;
#endif

Or use templates, which has the added benefit of allowing you to change the code at runtime to use one or the other.

Answer 2

You can't safely do this without modifying your source code, but that shouldn't be terribly difficult to do.

Using the preprocessor to force the keyword double in your program to be treated as float is a bad idea; it will make your program difficult to read, and if you happen to use long double anywhere it would be treated as long float, which is a syntax error.

As stix's answer suggests, you can add a typedef, either at the top of your program (if it's a single source file) or in some header that's #includeed by all the relevant source files:

typedef double real; /* or pick a different name */

Then go through your source code and change each occurrence of double to real. (Be careful about doing a blind global search-and-replace.)

Make sure that the program still compiles, runs, and behaves the same way after this change. Then you can change the typedef to:

 typedef float real;

and recompile to use float rather than double.

It's not quite that simple, though. If you're using functions declared in <math.h>, you'll want to use the right function for whatever floating-point type you're using; for example, sqrt() is for double, sqrtf() is for float, and sqrtl() is for long double.

If your compiler supports it, you might use the <tgmath.h> header, which defines type-generic macros corresponding to the math functions from <math.h>. If you use <tgmath.h>, then sqrt(x) will resolve to call the correct square root function depending on the type of the argument.