void fuerza(float *f[3], float v[3], float b[3], float q) {
  *f[0] = q * (v[1] * b[2] - v[2] * b[1]);
  *f[1] = q * (v[2] * b[0] - v[0] * b[2]);
  *f[2] = q * (v[0] * b[1] - v[1] * b[0]);
}


x86 assembly (clang, -O0):
<fuerza>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movss  %xmm0,-0x1c(%rbp)
movss  -0x1c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm1
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm3
mov    -0x10(%rbp),%rax
movss  0x8(%rax),%xmm2
mov    -0x18(%rbp),%rax
mulss  0x4(%rax),%xmm2
movd   %xmm2,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm2
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
mulss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    (%rax),%rax
movss  %xmm0,(%rax)
movss  -0x1c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movss  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm3
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm2
mov    -0x18(%rbp),%rax
mulss  0x8(%rax),%xmm2
movd   %xmm2,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm2
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
mulss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    0x8(%rax),%rax
movss  %xmm0,(%rax)
movss  -0x1c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm1
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm3
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm2
mov    -0x18(%rbp),%rax
mulss  (%rax),%xmm2
movd   %xmm2,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm2
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
mulss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    0x10(%rax),%rax
movss  %xmm0,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<fuerza>:
movss  0x4(%rsi),%xmm1
movss  0x4(%rdx),%xmm2
mulss  0x8(%rsi),%xmm2
mulss  0x8(%rdx),%xmm1
subss  %xmm2,%xmm1
mulss  %xmm0,%xmm1
mov    (%rdi),%rax
movss  %xmm1,(%rax)
movss  0x8(%rsi),%xmm1
movss  0x8(%rdx),%xmm2
mulss  (%rsi),%xmm2
mulss  (%rdx),%xmm1
subss  %xmm2,%xmm1
mulss  %xmm0,%xmm1
mov    0x8(%rdi),%rax
movss  %xmm1,(%rax)
movss  (%rsi),%xmm1
movss  (%rdx),%xmm2
mulss  0x4(%rsi),%xmm2
mulss  0x4(%rdx),%xmm1
subss  %xmm2,%xmm1
mulss  %xmm0,%xmm1
mov    0x10(%rdi),%rax
movss  %xmm1,(%rax)
ret

LLVM IR (-O0):
define void @fuerza(float** %0, float* %1, float* %2, float %3) {
%5 = alloca float**
%6 = alloca float*
%7 = alloca float*
%8 = alloca float
store float** %0, float*** %5
store float* %1, float** %6
store float* %2, float** %7
store float %3, float* %8
%9 = load float, float* %8
%10 = load float*, float** %6
%11 = getelementptr inbounds float, float* %10, i64 1
%12 = load float, float* %11
%13 = load float*, float** %7
%14 = getelementptr inbounds float, float* %13, i64 2
%15 = load float, float* %14
%16 = load float*, float** %6
%17 = getelementptr inbounds float, float* %16, i64 2
%18 = load float, float* %17
%19 = load float*, float** %7
%20 = getelementptr inbounds float, float* %19, i64 1
%21 = load float, float* %20
%22 = fmul float %18, %21
%23 = fneg float %22
%24 = call float @llvm.fmuladd.f32(float %12, float %15, float %23)
%25 = fmul float %9, %24
%26 = load float**, float*** %5
%27 = getelementptr inbounds float*, float** %26, i64 0
%28 = load float*, float** %27
store float %25, float* %28
%29 = load float, float* %8
%30 = load float*, float** %6
%31 = getelementptr inbounds float, float* %30, i64 2
%32 = load float, float* %31
%33 = load float*, float** %7
%34 = getelementptr inbounds float, float* %33, i64 0
%35 = load float, float* %34
%36 = load float*, float** %6
%37 = getelementptr inbounds float, float* %36, i64 0
%38 = load float, float* %37
%39 = load float*, float** %7
%40 = getelementptr inbounds float, float* %39, i64 2
%41 = load float, float* %40
%42 = fmul float %38, %41
%43 = fneg float %42
%44 = call float @llvm.fmuladd.f32(float %32, float %35, float %43)
%45 = fmul float %29, %44
%46 = load float**, float*** %5
%47 = getelementptr inbounds float*, float** %46, i64 1
%48 = load float*, float** %47
store float %45, float* %48
%49 = load float, float* %8
%50 = load float*, float** %6
%51 = getelementptr inbounds float, float* %50, i64 0
%52 = load float, float* %51
%53 = load float*, float** %7
%54 = getelementptr inbounds float, float* %53, i64 1
%55 = load float, float* %54
%56 = load float*, float** %6
%57 = getelementptr inbounds float, float* %56, i64 1
%58 = load float, float* %57
%59 = load float*, float** %7
%60 = getelementptr inbounds float, float* %59, i64 0
%61 = load float, float* %60
%62 = fmul float %58, %61
%63 = fneg float %62
%64 = call float @llvm.fmuladd.f32(float %52, float %55, float %63)
%65 = fmul float %49, %64
%66 = load float**, float*** %5
%67 = getelementptr inbounds float*, float** %66, i64 2
%68 = load float*, float** %67
store float %65, float* %68
ret void
}

LLVM IR (-O2):
define void @fuerza(float** nocapture %0, float* nocapture %1, float* nocapture %2, float %3) {
%5 = getelementptr inbounds float, float* %1, i64 1
%6 = load float, float* %5
%7 = getelementptr inbounds float, float* %2, i64 2
%8 = load float, float* %7
%9 = getelementptr inbounds float, float* %1, i64 2
%10 = load float, float* %9
%11 = getelementptr inbounds float, float* %2, i64 1
%12 = load float, float* %11
%13 = fneg float %10
%14 = fmul float %12, %13
%15 = tail call float @llvm.fmuladd.f32(float %6, float %8, float %14)
%16 = fmul float %15, %3
%17 = load float*, float** %0
store float %16, float* %17
%18 = load float, float* %9
%19 = load float, float* %2
%20 = load float, float* %1
%21 = load float, float* %7
%22 = fneg float %20
%23 = fmul float %21, %22
%24 = tail call float @llvm.fmuladd.f32(float %18, float %19, float %23)
%25 = fmul float %24, %3
%26 = getelementptr inbounds float*, float** %0, i64 1
%27 = load float*, float** %26
store float %25, float* %27
%28 = load float, float* %1
%29 = load float, float* %11
%30 = load float, float* %5
%31 = load float, float* %2
%32 = fneg float %30
%33 = fmul float %31, %32
%34 = tail call float @llvm.fmuladd.f32(float %28, float %29, float %33)
%35 = fmul float %34, %3
%36 = getelementptr inbounds float*, float** %0, i64 2
%37 = load float*, float** %36
store float %35, float* %37
ret void
}