void f(int *neq, float *t, float *y, float *ydot) {
  ydot[0] = -0.04f * y[0] + 1e4f * y[1] * y[2];
  ydot[2] = 3e7f * y[1] * y[1];
  ydot[1] = -ydot[0] - ydot[2];
}


x86 assembly (clang, -O0):
<f>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  0xee0(%rip),%xmm1
mulss  0x4(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulss  0x8(%rax),%xmm1
movss  0xec6(%rip),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
movss  0xea6(%rip),%xmm0
mulss  0x4(%rax),%xmm0
mov    -0x18(%rbp),%rax
mulss  0x4(%rax),%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,0x8(%rax)
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mov    -0x20(%rbp),%rax
subss  0x8(%rax),%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,0x4(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<f>:
movss  (%rdx),%xmm0
movss  0x4(%rdx),%xmm1
mulss  0xeef(%rip),%xmm1
mulss  0x8(%rdx),%xmm1
mulss  0xee6(%rip),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rcx)
movss  0x4(%rdx),%xmm1
movss  0xed5(%rip),%xmm2
mulss  %xmm1,%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,0x8(%rcx)
xorps  0xec9(%rip),%xmm0
subss  %xmm2,%xmm0
movss  %xmm0,0x4(%rcx)
ret

LLVM IR (-O0):
define void @f(i32* %0, float* %1, float* %2, float* %3) {
%5 = alloca i32*
%6 = alloca float*
%7 = alloca float*
%8 = alloca float*
store i32* %0, i32** %5
store float* %1, float** %6
store float* %2, float** %7
store float* %3, float** %8
%9 = load float*, float** %7
%10 = getelementptr inbounds float, float* %9, i64 0
%11 = load float, float* %10
%12 = load float*, float** %7
%13 = getelementptr inbounds float, float* %12, i64 1
%14 = load float, float* %13
%15 = fmul float 1.000000e+04, %14
%16 = load float*, float** %7
%17 = getelementptr inbounds float, float* %16, i64 2
%18 = load float, float* %17
%19 = fmul float %15, %18
%20 = call float @llvm.fmuladd.f32(float 0xBFA47AE140000000, float %11, float %19)
%21 = load float*, float** %8
%22 = getelementptr inbounds float, float* %21, i64 0
store float %20, float* %22
%23 = load float*, float** %7
%24 = getelementptr inbounds float, float* %23, i64 1
%25 = load float, float* %24
%26 = fmul float 3.000000e+07, %25
%27 = load float*, float** %7
%28 = getelementptr inbounds float, float* %27, i64 1
%29 = load float, float* %28
%30 = fmul float %26, %29
%31 = load float*, float** %8
%32 = getelementptr inbounds float, float* %31, i64 2
store float %30, float* %32
%33 = load float*, float** %8
%34 = getelementptr inbounds float, float* %33, i64 0
%35 = load float, float* %34
%36 = fneg float %35
%37 = load float*, float** %8
%38 = getelementptr inbounds float, float* %37, i64 2
%39 = load float, float* %38
%40 = fsub float %36, %39
%41 = load float*, float** %8
%42 = getelementptr inbounds float, float* %41, i64 1
store float %40, float* %42
ret void
}

LLVM IR (-O2):
define void @f(i32* nocapture readnone %0, float* nocapture readnone %1, float* nocapture %2, float* nocapture writeonly %3) {
%5 = load float, float* %2
%6 = getelementptr inbounds float, float* %2, i64 1
%7 = load float, float* %6
%8 = fmul float %7, 1.000000e+04
%9 = getelementptr inbounds float, float* %2, i64 2
%10 = load float, float* %9
%11 = fmul float %8, %10
%12 = tail call float @llvm.fmuladd.f32(float %5, float 0xBFA47AE140000000, float %11)
store float %12, float* %3
%13 = load float, float* %6
%14 = fmul float %13, 3.000000e+07
%15 = fmul float %13, %14
%16 = getelementptr inbounds float, float* %3, i64 2
store float %15, float* %16
%17 = fneg float %12
%18 = fsub float %17, %15
%19 = getelementptr inbounds float, float* %3, i64 1
store float %18, float* %19
ret void
}