    .text
    .type fmadd, %function
    .global fmadd
    /*
     * Scalar floating-point multiply-add: d = a * b + c.
     * @param X0 pointer to a (4x float)
     * @param X1 pointer to b (4x float)
     * @param X2 pointer to c (4x float)
     * @param X3 pointer to d (4x float)
     */
fmadd:
    ldr q0, [x0]
    ldr q1, [x1]
    ldr q2, [x2]
    ldr q3, [x3]
    fmadd s3, s0, s1, s2
    str q3, [x3]
    ret
    .size fmadd, (. - fmadd)
