void pendientes(double *s, double *f, int n, int dim) {
  int i;
  double x1, x2, y1, y2;
  for (i = 0; i < n; i++) {
    x1 = s[i * dim + 0];
    y1 = s[i * dim + 1];
    x2 = s[i * dim + 2];
    y2 = s[i * dim + 3];
    f[i * 2] = (y2 - y1) / (x2 - x1);
    f[i * 2 + 1] = -1 * ((y2 - y1) / (x2 - x1) * x1) + y1;
  }
}


x86 assembly (clang, -O0):
<pendientes>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1204 <pendientes+0x104>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
subsd  -0x38(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm1
subsd  -0x28(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm1
subsd  -0x38(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
subsd  -0x28(%rbp),%xmm0
divsd  %xmm0,%xmm1
mulsd  -0x28(%rbp),%xmm1
movsd  -0x38(%rbp),%xmm0
subsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1119 <pendientes+0x19>
pop    %rbp
ret

x86 assembly (clang, -O2):
<pendientes>:
test   %edx,%edx
jle    120d <pendientes+0x10d>
movslq %ecx,%r8
mov    %edx,%eax
xor    %r9d,%r9d
cmp    $0x1,%edx
je     11ac <pendientes+0xac>
cmp    $0x1,%ecx
jne    11ac <pendientes+0xac>
lea    (%rdi,%rax,8),%rcx
add    $0x18,%rcx
cmp    %rsi,%rcx
jbe    1143 <pendientes+0x43>
mov    %rax,%rcx
shl    $0x4,%rcx
add    %rsi,%rcx
cmp    %rdi,%rcx
jbe    1143 <pendientes+0x43>
xor    %r9d,%r9d
jmp    11ac <pendientes+0xac>
mov    %eax,%r9d
and    $0xfffffffe,%r9d
lea    0x18(%rdi),%rcx
mov    %r8,%r10
shl    $0x4,%r10
mov    %r9,%r11
mov    %rsi,%rdx
nopl   0x0(%rax,%rax,1)
movupd -0x18(%rcx),%xmm0
movupd -0x10(%rcx),%xmm1
movupd -0x8(%rcx),%xmm2
movupd (%rcx),%xmm3
subpd  %xmm1,%xmm3
subpd  %xmm0,%xmm2
divpd  %xmm2,%xmm3
mulpd  %xmm3,%xmm0
subpd  %xmm0,%xmm1
movapd %xmm3,%xmm0
unpcklpd %xmm1,%xmm0
movhlps %xmm3,%xmm1
movups %xmm1,0x10(%rdx)
movupd %xmm0,(%rdx)
add    %r10,%rcx
add    $0x20,%rdx
add    $0xfffffffffffffffe,%r11
jne    1160 <pendientes+0x60>
cmp    %rax,%r9
je     120d <pendientes+0x10d>
mov    %r9,%rcx
imul   %r8,%rcx
lea    (%rdi,%rcx,8),%rcx
add    $0x18,%rcx
shl    $0x3,%r8
sub    %r9,%rax
shl    $0x4,%r9
lea    (%r9,%rsi,1),%rdx
add    $0x8,%rdx
xchg   %ax,%ax
movsd  -0x18(%rcx),%xmm0
movsd  -0x10(%rcx),%xmm1
movsd  -0x8(%rcx),%xmm2
movsd  (%rcx),%xmm3
subsd  %xmm1,%xmm3
subsd  %xmm0,%xmm2
divsd  %xmm2,%xmm3
movsd  %xmm3,-0x8(%rdx)
mulsd  %xmm0,%xmm3
subsd  %xmm3,%xmm1
movsd  %xmm1,(%rdx)
add    %r8,%rcx
add    $0x10,%rdx
add    $0xffffffffffffffff,%rax
jne    11d0 <pendientes+0xd0>
ret

LLVM IR (-O0):
define void @pendientes(double* %0, double* %1, i32 %2, i32 %3) {
%5 = alloca double*
%6 = alloca double*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca double
store double* %0, double** %5
store double* %1, double** %6
store i32 %2, i32* %7
store i32 %3, i32* %8
store i32 0, i32* %9
br label %14
14: ; preds = %80, %4
%15 = load i32, i32* %9
%16 = load i32, i32* %7
%17 = icmp slt i32 %15, %16
br i1 %17, label %18, label %83
18: ; preds = %14
%19 = load double*, double** %5
%20 = load i32, i32* %9
%21 = load i32, i32* %8
%22 = mul nsw i32 %20, %21
%23 = add nsw i32 %22, 0
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds double, double* %19, i64 %24
%26 = load double, double* %25
store double %26, double* %10
%27 = load double*, double** %5
%28 = load i32, i32* %9
%29 = load i32, i32* %8
%30 = mul nsw i32 %28, %29
%31 = add nsw i32 %30, 1
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %27, i64 %32
%34 = load double, double* %33
store double %34, double* %12
%35 = load double*, double** %5
%36 = load i32, i32* %9
%37 = load i32, i32* %8
%38 = mul nsw i32 %36, %37
%39 = add nsw i32 %38, 2
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %35, i64 %40
%42 = load double, double* %41
store double %42, double* %11
%43 = load double*, double** %5
%44 = load i32, i32* %9
%45 = load i32, i32* %8
%46 = mul nsw i32 %44, %45
%47 = add nsw i32 %46, 3
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds double, double* %43, i64 %48
%50 = load double, double* %49
store double %50, double* %13
%51 = load double, double* %13
%52 = load double, double* %12
%53 = fsub double %51, %52
%54 = load double, double* %11
%55 = load double, double* %10
%56 = fsub double %54, %55
%57 = fdiv double %53, %56
%58 = load double*, double** %6
%59 = load i32, i32* %9
%60 = mul nsw i32 %59, 2
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds double, double* %58, i64 %61
store double %57, double* %62
%63 = load double, double* %13
%64 = load double, double* %12
%65 = fsub double %63, %64
%66 = load double, double* %11
%67 = load double, double* %10
%68 = fsub double %66, %67
%69 = fdiv double %65, %68
%70 = load double, double* %10
%71 = fmul double %69, %70
%72 = load double, double* %12
%73 = call double @llvm.fmuladd.f64(double -1.000000e+00, double %71, double %72)
%74 = load double*, double** %6
%75 = load i32, i32* %9
%76 = mul nsw i32 %75, 2
%77 = add nsw i32 %76, 1
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds double, double* %74, i64 %78
store double %73, double* %79
br label %80
80: ; preds = %18
%81 = load i32, i32* %9
%82 = add nsw i32 %81, 1
store i32 %82, i32* %9
br label %14
83: ; preds = %14
ret void
}

LLVM IR (-O2):
define void @pendientes(double* nocapture %0, double* nocapture writeonly %1, i32 %2, i32 %3) {
%5 = icmp sgt i32 %2, 0
br i1 %5, label %6, label %82
6: ; preds = %4
%7 = sext i32 %3 to i64
%8 = zext i32 %2 to i64
%9 = icmp ne i32 %2, 1
%10 = icmp eq i32 %3, 1
%11 = and i1 %9, %10
br i1 %11, label %12, label %55
12: ; preds = %6
%13 = shl nuw nsw i64 %8, 1
%14 = getelementptr double, double* %1, i64 %13
%15 = add nuw nsw i64 %8, 3
%16 = getelementptr double, double* %0, i64 %15
%17 = icmp ugt double* %16, %1
%18 = icmp ugt double* %14, %0
%19 = and i1 %17, %18
br i1 %19, label %55, label %20
20: ; preds = %12
%21 = and i64 %8, 4294967294
%22 = getelementptr inbounds double, double* %1, i64 -1
br label %23
23: ; preds = %23, %20
%24 = phi i64 [ 0, %20 ], [ %51, %23 ]
%25 = mul nsw i64 %24, %7
%26 = getelementptr inbounds double, double* %0, i64 %25
%27 = bitcast double* %26 to <2 x double>*
%28 = load <2 x double>, <2 x double>* %27
%29 = or i64 %25, 1
%30 = getelementptr inbounds double, double* %0, i64 %29
%31 = bitcast double* %30 to <2 x double>*
%32 = load <2 x double>, <2 x double>* %31
%33 = add nsw i64 %25, 2
%34 = getelementptr inbounds double, double* %0, i64 %33
%35 = bitcast double* %34 to <2 x double>*
%36 = load <2 x double>, <2 x double>* %35
%37 = add nsw i64 %25, 3
%38 = getelementptr inbounds double, double* %0, i64 %37
%39 = bitcast double* %38 to <2 x double>*
%40 = load <2 x double>, <2 x double>* %39
%41 = fsub <2 x double> %40, %32
%42 = fsub <2 x double> %36, %28
%43 = fdiv <2 x double> %41, %42
%44 = shl nuw nsw i64 %24, 1
%45 = fmul <2 x double> %28, %43
%46 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %45, <2 x double> <double -1.000000e+00, double -1.000000e+00>, <2 x double> %32)
%47 = or i64 %44, 1
%48 = getelementptr inbounds double, double* %22, i64 %47
%49 = bitcast double* %48 to <4 x double>*
%50 = shufflevector <2 x double> %43, <2 x double> %46, <4 x i32> <i32 0, i32 2, i32 1, i32 3>
store <4 x double> %50, <4 x double>* %49
%51 = add nuw i64 %24, 2
%52 = icmp eq i64 %51, %21
br i1 %52, label %53, label %23
53: ; preds = %23
%54 = icmp eq i64 %21, %8
br i1 %54, label %82, label %55
55: ; preds = %12, %6, %53
%56 = phi i64 [ 0, %12 ], [ 0, %6 ], [ %21, %53 ]
br label %57
57: ; preds = %55, %57
%58 = phi i64 [ %80, %57 ], [ %56, %55 ]
%59 = mul nsw i64 %58, %7
%60 = getelementptr inbounds double, double* %0, i64 %59
%61 = load double, double* %60
%62 = add nsw i64 %59, 1
%63 = getelementptr inbounds double, double* %0, i64 %62
%64 = load double, double* %63
%65 = add nsw i64 %59, 2
%66 = getelementptr inbounds double, double* %0, i64 %65
%67 = load double, double* %66
%68 = add nsw i64 %59, 3
%69 = getelementptr inbounds double, double* %0, i64 %68
%70 = load double, double* %69
%71 = fsub double %70, %64
%72 = fsub double %67, %61
%73 = fdiv double %71, %72
%74 = shl nuw nsw i64 %58, 1
%75 = getelementptr inbounds double, double* %1, i64 %74
store double %73, double* %75
%76 = fmul double %61, %73
%77 = tail call double @llvm.fmuladd.f64(double %76, double -1.000000e+00, double %64)
%78 = or i64 %74, 1
%79 = getelementptr inbounds double, double* %1, i64 %78
store double %77, double* %79
%80 = add nuw nsw i64 %58, 1
%81 = icmp eq i64 %80, %8
br i1 %81, label %82, label %57
82: ; preds = %57, %53, %4
ret void
}