On many CPUs, floating-point values are usually stored in special "floating-point registers", and are added, subtracted, etc with special "floating-point instructions", but other than the name these registers and instructions are exactly analogous to regular integer registers and instructions. For example, the integer PowerPC assembly code to add registers 1 and 2 into register 3 is "add r3,r1,r2"; the floating-point code to add floating-point registers 1 and 2 into floating-point register 3 is "fadd fr3,fr1,fr2".

x86 is not like that.

The problem is that the x86 instruction set wasn't designed with floating-point in mind; they added floating-point instructions to the CPU later (with the 8087, a separate chip that handled all floating-point instructions). Unfortunately, there weren't many unused opcode bytes left, and (being the 1980's, when bytes were expensive) the designers really didn't want to make the instructions longer. So instead of the usual instructions like "add register A to register B", x86 floating-point has just "add", which saves the bits that would be needed to specify the source and destination registers!

But the question is, what the heck are you adding? The answer is the "top two values on the floating-point register stack". That's not "the stack" (the memory area used by function calls), it's a separate set of values totally internal to the CPU's floating-point hardware. There are various load functions that push values onto the floating-point register stack, and most of the arithmetic functions read from the top of the floating-point register stack. So to compute stuff, you load the values you want to manipulate onto the floating-point register stack, and then use some arithmetic instructions.

For example, to add together the three values a, b, and c, you'd "load a; load b; add; load c; add;". Or, you could "load a; load b; load c; add; add;". If you've ever used an HP calculator, or written Postscript or Forth code, you've seen this "Reverse Polish Notation".

fldpi ; Push "pi" onto floating-point stackThere are lots of useful floating-point instructions:

sub esp,8 ; Make room on the stack for an 8-byte double

fstp QWORD [esp]; Push printf's double parameter onto the stack

push my_string ; Push printf's string parameter (below)

extern printf

call printf ; Print string

add esp,12 ; Clean up stack

ret ; Done with function

my_string: db "Yo! Here's our float: %f",0xa,0

Assembly |
Description |

fld1 |
Pushes into the floating-point registers the constant 1.0 |

fldz |
Pushes into the floating-point registers the constant 0.0 |

fldpi |
Pushes the constant pi. (Try this in NetRun now!) |

fld DWORD [eax] |
Pushes
into the floating-point registers the 4-byte "float" loaded from memory
at address eax. This is how most constants get loaded into the
program. (Try this in NetRun now!) |

fild DWORD [eax] |
Pushes into the floating-point registers the 4-byte "int" loaded from memory at address eax. |

fld QWORD [eax] |
Pushes an 8-byte "double" loaded from address eax. (Try this in NetRun now!) |

fld st0 |
Duplicates the top float, so there are now two copes of it. (Try this in NetRun now!) |

fstp DWORD [eax] | Pops the top floating-point value, and stores it as a "float" to address eax. |

fst DWORD [eax] | Reads the top floating-point value and stores it as a "float" to address eax. This doesn't change the value stored on the floating-point stack. |

fstp QWORD [eax] | Pops the top floating-point value, and stores it as a "double" to address eax. |

faddp |
Add the top two values, pushes the result. (Try this in NetRun now!) |

fsubp |
Subtract the two values, pushes the result. Note "fld A; fld B; fsubp;" computes A-B. (Try this in NetRun now!) There's also a "fsubrp" that subtracts in the opposite order (computing B-A). |

fmulp |
Multiply the top two values. |

fdivp |
Divide the top two values. Note "fld A; fld B; fdivp;" computes A/B. (Try this in NetRun now!) There's also a "fdivrp" that divides in the opposite order (computing B/A). |

fabs |
Take the absolute value of the top floating-point value. |

fsqrt |
Take the square root of the top floating-point value. |

fsin |
Take the sin() of the top floating-point value, treated as radians. (Try this in NetRun now!) |

In general, the "p" instructions pop a value from the floating-point stack.

The non-"p" instructions don't. For example, there isn't a "fsinp" instruction, since sin only takes one argument, so the stack stays the same height after doing a sin().

x86 has quite a few really bizarre-sounding floating-point instructions. Intel's Reference Volume 2 has the complete list (Section 3, alphabetized under "f"). The "+1" and "-1" versions are designed to decrease roundoff, by shifting the input to the most sensitive region.

F2XM1 |
2^{x} - 1 |

FYL2X |
y*log_{2}(x), where x is on top of the floating-point stack. |

FYL2XP1 | y*log_{2}(x+1), where x is on top |

FCHS |
-x |

FSINCOS |
Computes *both* sin(x) and cos(x). cos(x) ends up on top. |

FPATAN |
atan2(a/b), where b is on top |

FPREM |
fmod(a,b), where b is on top |

FRNDINT |
Round to the nearest integer |

FXCH |
Swap the top two values on the floating-point stack |