Looks like PA-RISC 1.1 doesn't have a fused multiply-add instruction, just a confusing FMPYADD which does a multiply and an add at the same time but on unrelated inputs, so we'll have to use the portable C versions. (PA-RISC 2.0 has FMPYFADD.)