#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
scalar_t__ fabs(double);
int round(double);
int plot_line(int (*P)[2], int w, int h, double a, double b, double c) {
  int r = 0;
  if (fabs(a) < fabs(b)) {
    double p = -a / b;
    double q = -c / b;
    for (int x = 0; x < w; x++) {
      int y = round(p * x + q);
      P[r][0] = x;
      if (y >= 0 && y < h)
        P[r++][1] = y;
    }
  } else {
    double p = -b / a;
    double q = -c / a;
    for (int y = 0; y < h; y++) {
      int x = round(p * y + q);
      P[r][1] = y;
      if (x >= 0 && x < w)
        P[r++][0] = x;
    }
  }
  return r;
}


x86 assembly (clang, -O0):
<plot_line>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
movsd  %xmm2,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
movsd  -0x18(%rbp),%xmm0
call   1030 <fabs@plt>
mov    %rax,-0x68(%rbp)
movsd  -0x20(%rbp),%xmm0
call   1030 <fabs@plt>
mov    %rax,%rcx
mov    -0x68(%rbp),%rax
cmp    %rcx,%rax
jge    1243 <plot_line+0x123>
movsd  -0x18(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
divsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x28(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
divsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    123e <plot_line+0x11e>
movsd  -0x38(%rbp),%xmm0
cvtsi2sdl -0x44(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1040 <round@plt>
mov    %eax,-0x48(%rbp)
mov    -0x44(%rbp),%edx
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %edx,(%rax,%rcx,8)
cmpl   $0x0,-0x48(%rbp)
jl     122b <plot_line+0x10b>
mov    -0x48(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    122b <plot_line+0x10b>
mov    -0x48(%rbp),%edx
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    %ecx,%esi
add    $0x1,%esi
mov    %esi,-0x2c(%rbp)
movslq %ecx,%rcx
mov    %edx,0x4(%rax,%rcx,8)
jmp    1230 <plot_line+0x110>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    11c3 <plot_line+0xa3>
jmp    1316 <plot_line+0x1f6>
movsd  -0x20(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x28(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1311 <plot_line+0x1f1>
movsd  -0x50(%rbp),%xmm0
cvtsi2sdl -0x5c(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1040 <round@plt>
mov    %eax,-0x60(%rbp)
mov    -0x5c(%rbp),%edx
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %edx,0x4(%rax,%rcx,8)
cmpl   $0x0,-0x60(%rbp)
jl     12fe <plot_line+0x1de>
mov    -0x60(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12fe <plot_line+0x1de>
mov    -0x60(%rbp),%edx
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    %ecx,%esi
add    $0x1,%esi
mov    %esi,-0x2c(%rbp)
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,8)
jmp    1303 <plot_line+0x1e3>
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    1296 <plot_line+0x176>
jmp    1316 <plot_line+0x1f6>
mov    -0x2c(%rbp),%eax
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<plot_line>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x40,%rsp
movapd %xmm2,0x10(%rsp)
movapd %xmm1,0x20(%rsp)
mov    %edx,%r14d
mov    %esi,%r15d
mov    %rdi,%r12
movaps %xmm0,0x30(%rsp)
call   1030 <fabs@plt>
mov    %rax,%rbp
movapd 0x20(%rsp),%xmm0
call   1030 <fabs@plt>
cmp    %rax,%rbp
jge    1205 <plot_line+0xe5>
test   %r15d,%r15d
jle    12a1 <plot_line+0x181>
movapd 0xe8d(%rip),%xmm0
movapd 0x30(%rsp),%xmm1
xorpd  %xmm0,%xmm1
movapd 0x20(%rsp),%xmm2
divsd  %xmm2,%xmm1
movapd %xmm1,0x30(%rsp)
movapd 0x10(%rsp),%xmm1
xorpd  %xmm0,%xmm1
divsd  %xmm2,%xmm1
movapd %xmm1,0x10(%rsp)
xorpd  %xmm1,%xmm1
xor    %ebp,%ebp
xor    %ebx,%ebx
jmp    11ca <plot_line+0xaa>
nopl   0x0(%rax,%rax,1)
movsd  0x8(%rsp),%xmm1
addsd  0xe52(%rip),%xmm1
add    $0x1,%ebx
cmp    %ebx,%r15d
je     12a3 <plot_line+0x183>
movapd 0x30(%rsp),%xmm0
movsd  %xmm1,0x8(%rsp)
mulsd  %xmm1,%xmm0
addsd  0x10(%rsp),%xmm0
call   1040 <round@plt>
movslq %ebp,%rcx
mov    %ebx,(%r12,%rcx,8)
test   %eax,%eax
js     11b0 <plot_line+0x90>
cmp    %r14d,%eax
movsd  0x8(%rsp),%xmm1
jge    11b6 <plot_line+0x96>
lea    0x1(%rcx),%ebp
mov    %eax,0x4(%r12,%rcx,8)
jmp    11b6 <plot_line+0x96>
test   %r14d,%r14d
jle    12a1 <plot_line+0x181>
movapd 0xdea(%rip),%xmm0
movapd 0x20(%rsp),%xmm1
xorpd  %xmm0,%xmm1
movapd 0x30(%rsp),%xmm2
divsd  %xmm2,%xmm1
movapd %xmm1,0x20(%rsp)
movapd 0x10(%rsp),%xmm1
xorpd  %xmm0,%xmm1
divsd  %xmm2,%xmm1
movapd %xmm1,0x10(%rsp)
xorpd  %xmm1,%xmm1
xor    %ebp,%ebp
xor    %ebx,%ebx
jmp    1266 <plot_line+0x146>
xchg   %ax,%ax
movsd  0x8(%rsp),%xmm1
addsd  0xdb2(%rip),%xmm1
add    $0x1,%ebx
cmp    %ebx,%r14d
je     12a3 <plot_line+0x183>
movapd 0x20(%rsp),%xmm0
movsd  %xmm1,0x8(%rsp)
mulsd  %xmm1,%xmm0
addsd  0x10(%rsp),%xmm0
call   1040 <round@plt>
movslq %ebp,%rcx
mov    %ebx,0x4(%r12,%rcx,8)
test   %eax,%eax
js     1250 <plot_line+0x130>
cmp    %r15d,%eax
movsd  0x8(%rsp),%xmm1
jge    1256 <plot_line+0x136>
lea    0x1(%rcx),%ebp
mov    %eax,(%r12,%rcx,8)
jmp    1256 <plot_line+0x136>
xor    %ebp,%ebp
mov    %ebp,%eax
add    $0x40,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @plot_line([2 x i32]* %0, i32 %1, i32 %2, double %3, double %4, double %5) {
%7 = alloca [2 x i32]*
%8 = alloca i32
%9 = alloca i32
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca i32
%14 = alloca double
%15 = alloca double
%16 = alloca i32
%17 = alloca i32
%18 = alloca double
%19 = alloca double
%20 = alloca i32
%21 = alloca i32
store [2 x i32]* %0, [2 x i32]** %7
store i32 %1, i32* %8
store i32 %2, i32* %9
store double %3, double* %10
store double %4, double* %11
store double %5, double* %12
store i32 0, i32* %13
%22 = load double, double* %10
%23 = call i64 @fabs(double %22)
%24 = load double, double* %11
%25 = call i64 @fabs(double %24)
%26 = icmp slt i64 %23, %25
br i1 %26, label %27, label %72
27: ; preds = %6
%28 = load double, double* %10
%29 = fneg double %28
%30 = load double, double* %11
%31 = fdiv double %29, %30
store double %31, double* %14
%32 = load double, double* %12
%33 = fneg double %32
%34 = load double, double* %11
%35 = fdiv double %33, %34
store double %35, double* %15
store i32 0, i32* %16
br label %36
36: ; preds = %68, %27
%37 = load i32, i32* %16
%38 = load i32, i32* %8
%39 = icmp slt i32 %37, %38
br i1 %39, label %40, label %71
40: ; preds = %36
%41 = load double, double* %14
%42 = load i32, i32* %16
%43 = sitofp i32 %42 to double
%44 = load double, double* %15
%45 = call double @llvm.fmuladd.f64(double %41, double %43, double %44)
%46 = call i32 @round(double %45)
store i32 %46, i32* %17
%47 = load i32, i32* %16
%48 = load [2 x i32]*, [2 x i32]** %7
%49 = load i32, i32* %13
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds [2 x i32], [2 x i32]* %48, i64 %50
%52 = getelementptr inbounds [2 x i32], [2 x i32]* %51, i64 0, i64 0
store i32 %47, i32* %52
%53 = load i32, i32* %17
%54 = icmp sge i32 %53, 0
br i1 %54, label %55, label %67
55: ; preds = %40
%56 = load i32, i32* %17
%57 = load i32, i32* %9
%58 = icmp slt i32 %56, %57
br i1 %58, label %59, label %67
59: ; preds = %55
%60 = load i32, i32* %17
%61 = load [2 x i32]*, [2 x i32]** %7
%62 = load i32, i32* %13
%63 = add nsw i32 %62, 1
store i32 %63, i32* %13
%64 = sext i32 %62 to i64
%65 = getelementptr inbounds [2 x i32], [2 x i32]* %61, i64 %64
%66 = getelementptr inbounds [2 x i32], [2 x i32]* %65, i64 0, i64 1
store i32 %60, i32* %66
br label %67
67: ; preds = %59, %55, %40
br label %68
68: ; preds = %67
%69 = load i32, i32* %16
%70 = add nsw i32 %69, 1
store i32 %70, i32* %16
br label %36
71: ; preds = %36
br label %117
72: ; preds = %6
%73 = load double, double* %11
%74 = fneg double %73
%75 = load double, double* %10
%76 = fdiv double %74, %75
store double %76, double* %18
%77 = load double, double* %12
%78 = fneg double %77
%79 = load double, double* %10
%80 = fdiv double %78, %79
store double %80, double* %19
store i32 0, i32* %20
br label %81
81: ; preds = %113, %72
%82 = load i32, i32* %20
%83 = load i32, i32* %9
%84 = icmp slt i32 %82, %83
br i1 %84, label %85, label %116
85: ; preds = %81
%86 = load double, double* %18
%87 = load i32, i32* %20
%88 = sitofp i32 %87 to double
%89 = load double, double* %19
%90 = call double @llvm.fmuladd.f64(double %86, double %88, double %89)
%91 = call i32 @round(double %90)
store i32 %91, i32* %21
%92 = load i32, i32* %20
%93 = load [2 x i32]*, [2 x i32]** %7
%94 = load i32, i32* %13
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds [2 x i32], [2 x i32]* %93, i64 %95
%97 = getelementptr inbounds [2 x i32], [2 x i32]* %96, i64 0, i64 1
store i32 %92, i32* %97
%98 = load i32, i32* %21
%99 = icmp sge i32 %98, 0
br i1 %99, label %100, label %112
100: ; preds = %85
%101 = load i32, i32* %21
%102 = load i32, i32* %8
%103 = icmp slt i32 %101, %102
br i1 %103, label %104, label %112
104: ; preds = %100
%105 = load i32, i32* %21
%106 = load [2 x i32]*, [2 x i32]** %7
%107 = load i32, i32* %13
%108 = add nsw i32 %107, 1
store i32 %108, i32* %13
%109 = sext i32 %107 to i64
%110 = getelementptr inbounds [2 x i32], [2 x i32]* %106, i64 %109
%111 = getelementptr inbounds [2 x i32], [2 x i32]* %110, i64 0, i64 0
store i32 %105, i32* %111
br label %112
112: ; preds = %104, %100, %85
br label %113
113: ; preds = %112
%114 = load i32, i32* %20
%115 = add nsw i32 %114, 1
store i32 %115, i32* %20
br label %81
116: ; preds = %81
br label %117
117: ; preds = %116, %71
%118 = load i32, i32* %13
ret i32 %118
}

LLVM IR (-O2):
define i32 @plot_line([2 x i32]* nocapture writeonly %0, i32 %1, i32 %2, double %3, double %4, double %5) {
%7 = tail call i64 @fabs(double %3)
%8 = tail call i64 @fabs(double %4)
%9 = icmp slt i64 %7, %8
br i1 %9, label %10, label %34
10: ; preds = %6
%11 = fneg double %3
%12 = fdiv double %11, %4
%13 = fneg double %5
%14 = fdiv double %13, %4
%15 = icmp sgt i32 %1, 0
br i1 %15, label %16, label %58
16: ; preds = %10, %30
%17 = phi i32 [ %31, %30 ], [ 0, %10 ]
%18 = phi i32 [ %32, %30 ], [ 0, %10 ]
%19 = sitofp i32 %18 to double
%20 = tail call double @llvm.fmuladd.f64(double %12, double %19, double %14)
%21 = tail call i32 @round(double %20)
%22 = sext i32 %17 to i64
%23 = getelementptr inbounds [2 x i32], [2 x i32]* %0, i64 %22, i64 0
store i32 %18, i32* %23
%24 = icmp sgt i32 %21, -1
%25 = icmp slt i32 %21, %2
%26 = and i1 %24, %25
br i1 %26, label %27, label %30
27: ; preds = %16
%28 = add nsw i32 %17, 1
%29 = getelementptr inbounds [2 x i32], [2 x i32]* %0, i64 %22, i64 1
store i32 %21, i32* %29
br label %30
30: ; preds = %27, %16
%31 = phi i32 [ %28, %27 ], [ %17, %16 ]
%32 = add nuw nsw i32 %18, 1
%33 = icmp eq i32 %32, %1
br i1 %33, label %58, label %16
34: ; preds = %6
%35 = fneg double %4
%36 = fdiv double %35, %3
%37 = fneg double %5
%38 = fdiv double %37, %3
%39 = icmp sgt i32 %2, 0
br i1 %39, label %40, label %58
40: ; preds = %34, %54
%41 = phi i32 [ %55, %54 ], [ 0, %34 ]
%42 = phi i32 [ %56, %54 ], [ 0, %34 ]
%43 = sitofp i32 %42 to double
%44 = tail call double @llvm.fmuladd.f64(double %36, double %43, double %38)
%45 = tail call i32 @round(double %44)
%46 = sext i32 %41 to i64
%47 = getelementptr inbounds [2 x i32], [2 x i32]* %0, i64 %46, i64 1
store i32 %42, i32* %47
%48 = icmp sgt i32 %45, -1
%49 = icmp slt i32 %45, %1
%50 = and i1 %48, %49
br i1 %50, label %51, label %54
51: ; preds = %40
%52 = add nsw i32 %41, 1
%53 = getelementptr inbounds [2 x i32], [2 x i32]* %0, i64 %46, i64 0
store i32 %45, i32* %53
br label %54
54: ; preds = %51, %40
%55 = phi i32 [ %52, %51 ], [ %41, %40 ]
%56 = add nuw nsw i32 %42, 1
%57 = icmp eq i32 %56, %2
br i1 %57, label %58, label %40
58: ; preds = %54, %30, %34, %10
%59 = phi i32 [ 0, %10 ], [ 0, %34 ], [ %31, %30 ], [ %55, %54 ]
ret i32 %59
}