#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 */
typedef int /*<<< orphan*/ FILE;

/* Variables and functions */
int EOF;
double EXIT_FAILURE;
scalar_t__ calloc(int, int);
char fgetc(int /*<<< orphan*/ *);
int /*<<< orphan*/ free(double *);
int fscanf(int /*<<< orphan*/ *, char *, double *);
int /*<<< orphan*/ perror(char *);
int /*<<< orphan*/ printf(char *, double);
scalar_t__ realloc(double *, int);
int /*<<< orphan*/ skipWord(int /*<<< orphan*/ *);
double column(FILE *stream, int *tab) {
  double *tabNumber = NULL;
  double number;
  int lengthTab = 1, countMax = 0, count = 0;
  int test, i, j = 0, col = 1;
  if (!(tabNumber = (double *)calloc(1, sizeof(double)))) {
    perror("num-utils-ng");
    return EXIT_FAILURE;
  }
  while ((test = fscanf(stream, "%lf", &number)) != EOF) {
    if (!test)
      skipWord(stream);
    if (count == lengthTab) {
      lengthTab *= 2;
      if (!(tabNumber =
                (double *)realloc(tabNumber, lengthTab * sizeof(double)))) {
        perror("num-utils-ng");
        return EXIT_FAILURE;
      }
    }
    if (fgetc(stream) != '\n') {
      tabNumber[count] += number;
      count++;
    } else {
      if (count >= countMax)
        countMax = count;
      tabNumber[count] += number;
      count = 0;
    }
  }
  for (i = 0; i < countMax + 1; i++) {
    if (!tab || (col == tab[j])) {
      printf("%lf\n", tabNumber[i]);
      j++;
    }
    col++;
  }
  free(tabNumber);
  return 0;
}


x86 assembly (clang, -O0):
<column>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
movq   $0x0,-0x20(%rbp)
movl   $0x1,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
movl   $0x0,-0x40(%rbp)
movl   $0x1,-0x44(%rbp)
mov    $0x1,%edi
mov    $0x8,%esi
call   1070 <calloc@plt>
mov    %rax,-0x20(%rbp)
cmp    $0x0,%rax
jne    11f9 <column+0x79>
lea    0xe21(%rip),%rdi
call   10a0 <perror@plt>
mov    0x2ded(%rip),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1383 <column+0x203>
jmp    11fe <column+0x7e>
mov    -0x10(%rbp),%rdi
lea    0xe04(%rip),%rsi
lea    -0x28(%rbp),%rdx
call   1050 <fscanf@plt>
mov    %eax,-0x38(%rbp)
mov    0x2dcc(%rip),%rcx
cmp    (%rcx),%eax
je     1302 <column+0x182>
cmpl   $0x0,-0x38(%rbp)
jne    1237 <column+0xb7>
mov    -0x10(%rbp),%rdi
call   1090 <skipWord@plt>
mov    -0x34(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jne    1293 <column+0x113>
mov    -0x2c(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x20(%rbp),%rdi
movslq -0x2c(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%esi
call   1080 <realloc@plt>
mov    %rax,-0x20(%rbp)
cmp    $0x0,%rax
jne    128e <column+0x10e>
lea    0xd8c(%rip),%rdi
call   10a0 <perror@plt>
mov    0x2d58(%rip),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1383 <column+0x203>
jmp    1293 <column+0x113>
mov    -0x10(%rbp),%rdi
call   1060 <fgetc@plt>
movsbl %al,%eax
cmp    $0xa,%eax
je     12cd <column+0x14d>
movsd  -0x28(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    12fd <column+0x17d>
mov    -0x34(%rbp),%eax
cmp    -0x30(%rbp),%eax
jl     12df <column+0x15f>
mov    -0x34(%rbp),%eax
mov    %eax,-0x30(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x34(%rbp)
jmp    11fe <column+0x7e>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1372 <column+0x1f2>
cmpq   $0x0,-0x18(%rbp)
je     1339 <column+0x1b9>
mov    -0x44(%rbp),%eax
mov    -0x18(%rbp),%rcx
movslq -0x40(%rbp),%rdx
cmp    (%rcx,%rdx,4),%eax
jne    135b <column+0x1db>
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
lea    0xcc4(%rip),%rdi
call   1040 <printf@plt>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1309 <column+0x189>
mov    -0x20(%rbp),%rdi
call   1030 <free@plt>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<column>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rsi,%r12
mov    %rdi,%r15
mov    $0x1,%edi
mov    $0x8,%esi
call   1070 <calloc@plt>
test   %rax,%rax
je     12f8 <column+0x178>
mov    %rax,%rbx
lea    0xe57(%rip),%rsi
lea    0x8(%rsp),%rdx
mov    %r15,%rdi
call   1050 <fscanf@plt>
mov    0x2e1e(%rip),%rcx
xor    %r14d,%r14d
mov    $0x0,%r13d
cmp    (%rcx),%eax
je     12a1 <column+0x121>
mov    %r12,0x10(%rsp)
mov    $0x1,%ebp
xor    %r12d,%r12d
xor    %r13d,%r13d
jmp    122c <column+0xac>
nopl   (%rax)
cmp    %r13d,%r12d
cmovl  %r13d,%r12d
movsd  0x8(%rsp),%xmm0
addsd  (%rbx,%rcx,8),%xmm0
movsd  %xmm0,(%rbx,%rcx,8)
mov    %r12d,%r13d
xor    %r12d,%r12d
mov    %r15,%rdi
lea    0xdf6(%rip),%rsi
lea    0x8(%rsp),%rdx
call   1050 <fscanf@plt>
mov    0x2dc0(%rip),%rcx
cmp    (%rcx),%eax
je     1297 <column+0x117>
test   %eax,%eax
je     1240 <column+0xc0>
cmp    %ebp,%r12d
je     124d <column+0xcd>
jmp    126b <column+0xeb>
nopw   0x0(%rax,%rax,1)
mov    %r15,%rdi
call   1090 <skipWord@plt>
cmp    %ebp,%r12d
jne    126b <column+0xeb>
mov    %r12d,%esi
shl    $0x4,%esi
mov    %rbx,%rdi
call   1080 <realloc@plt>
test   %rax,%rax
je     12f8 <column+0x178>
mov    %rax,%rbx
lea    (%r12,%r12,1),%ebp
mov    %r15,%rdi
call   1060 <fgetc@plt>
movslq %r12d,%rcx
cmp    $0xa,%al
je     11f0 <column+0x70>
movsd  0x8(%rsp),%xmm0
addsd  (%rbx,%rcx,8),%xmm0
movsd  %xmm0,(%rbx,%rcx,8)
add    $0x1,%r12d
jmp    120d <column+0x8d>
test   %r13d,%r13d
mov    0x10(%rsp),%r12
js     12dd <column+0x15d>
add    $0x1,%r13d
lea    0xd65(%rip),%r15
xor    %eax,%eax
jmp    12c9 <column+0x149>
movsd  (%rbx,%rax,8),%xmm0
mov    %r15,%rdi
call   1040 <printf@plt>
add    $0x1,%r14d
mov    %rbp,%rax
cmp    %rbp,%r13
je     12dd <column+0x15d>
lea    0x1(%rax),%rbp
test   %r12,%r12
je     12b0 <column+0x130>
movslq %r14d,%rcx
cmp    (%r12,%rcx,4),%ebp
je     12b0 <column+0x130>
jmp    12c1 <column+0x141>
mov    %rbx,%rdi
call   1030 <free@plt>
xorpd  %xmm0,%xmm0
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
lea    0xd01(%rip),%rdi
call   10a0 <perror@plt>
mov    0x2ccd(%rip),%rax
movsd  (%rax),%xmm0
jmp    12e9 <column+0x169>

LLVM IR (-O0):
@.str = [13 x i8] c"num-utils-ng\00"
@EXIT_FAILURE = global double 0.000000e+00
@.str.1 = [4 x i8] c"%lf\00"
@EOF = global i32 0
@.str.2 = [5 x i8] c"%lf\0A\00"
define double @column(i32* %0, i32* %1) {
%3 = alloca double
%4 = alloca i32*
%5 = alloca i32*
%6 = alloca double*
%7 = alloca double
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
store i32* %0, i32** %4
store i32* %1, i32** %5
store double* null, double** %6
store i32 1, i32* %8
store i32 0, i32* %9
store i32 0, i32* %10
store i32 0, i32* %13
store i32 1, i32* %14
%15 = call i64 @calloc(i32 1, i32 8)
%16 = inttoptr i64 %15 to double*
store double* %16, double** %6
%17 = icmp ne double* %16, null
br i1 %17, label %21, label %18
18: ; preds = %2
%19 = call i32 @perror(i8* getelementptr inbounds ([13 x i8], [13 x i8]* @.str, i64 0, i64 0))
%20 = load double, double* @EXIT_FAILURE
store double %20, double* %3
br label %117
21: ; preds = %2
br label %22
22: ; preds = %81, %21
%23 = load i32*, i32** %4
%24 = call i32 @fscanf(i32* %23, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.1, i64 0, i64 0), double* %7)
store i32 %24, i32* %11
%25 = load i32, i32* @EOF
%26 = icmp ne i32 %24, %25
br i1 %26, label %27, label %82
27: ; preds = %22
%28 = load i32, i32* %11
%29 = icmp ne i32 %28, 0
br i1 %29, label %33, label %30
30: ; preds = %27
%31 = load i32*, i32** %4
%32 = call i32 @skipWord(i32* %31)
br label %33
33: ; preds = %30, %27
%34 = load i32, i32* %10
%35 = load i32, i32* %8
%36 = icmp eq i32 %34, %35
br i1 %36, label %37, label %52
37: ; preds = %33
%38 = load i32, i32* %8
%39 = mul nsw i32 %38, 2
store i32 %39, i32* %8
%40 = load double*, double** %6
%41 = load i32, i32* %8
%42 = sext i32 %41 to i64
%43 = mul i64 %42, 8
%44 = trunc i64 %43 to i32
%45 = call i64 @realloc(double* %40, i32 %44)
%46 = inttoptr i64 %45 to double*
store double* %46, double** %6
%47 = icmp ne double* %46, null
br i1 %47, label %51, label %48
48: ; preds = %37
%49 = call i32 @perror(i8* getelementptr inbounds ([13 x i8], [13 x i8]* @.str, i64 0, i64 0))
%50 = load double, double* @EXIT_FAILURE
store double %50, double* %3
br label %117
51: ; preds = %37
br label %52
52: ; preds = %51, %33
%53 = load i32*, i32** %4
%54 = call signext i8 @fgetc(i32* %53)
%55 = sext i8 %54 to i32
%56 = icmp ne i32 %55, 10
br i1 %56, label %57, label %67
57: ; preds = %52
%58 = load double, double* %7
%59 = load double*, double** %6
%60 = load i32, i32* %10
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds double, double* %59, i64 %61
%63 = load double, double* %62
%64 = fadd double %63, %58
store double %64, double* %62
%65 = load i32, i32* %10
%66 = add nsw i32 %65, 1
store i32 %66, i32* %10
br label %81
67: ; preds = %52
%68 = load i32, i32* %10
%69 = load i32, i32* %9
%70 = icmp sge i32 %68, %69
br i1 %70, label %71, label %73
71: ; preds = %67
%72 = load i32, i32* %10
store i32 %72, i32* %9
br label %73
73: ; preds = %71, %67
%74 = load double, double* %7
%75 = load double*, double** %6
%76 = load i32, i32* %10
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds double, double* %75, i64 %77
%79 = load double, double* %78
%80 = fadd double %79, %74
store double %80, double* %78
store i32 0, i32* %10
br label %81
81: ; preds = %73, %57
br label %22
82: ; preds = %22
store i32 0, i32* %12
br label %83
83: ; preds = %111, %82
%84 = load i32, i32* %12
%85 = load i32, i32* %9
%86 = add nsw i32 %85, 1
%87 = icmp slt i32 %84, %86
br i1 %87, label %88, label %114
88: ; preds = %83
%89 = load i32*, i32** %5
%90 = icmp ne i32* %89, null
br i1 %90, label %91, label %99
91: ; preds = %88
%92 = load i32, i32* %14
%93 = load i32*, i32** %5
%94 = load i32, i32* %13
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds i32, i32* %93, i64 %95
%97 = load i32, i32* %96
%98 = icmp eq i32 %92, %97
br i1 %98, label %99, label %108
99: ; preds = %91, %88
%100 = load double*, double** %6
%101 = load i32, i32* %12
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds double, double* %100, i64 %102
%104 = load double, double* %103
%105 = call i32 @printf(i8* getelementptr inbounds ([5 x i8], [5 x i8]* @.str.2, i64 0, i64 0), double %104)
%106 = load i32, i32* %13
%107 = add nsw i32 %106, 1
store i32 %107, i32* %13
br label %108
108: ; preds = %99, %91
%109 = load i32, i32* %14
%110 = add nsw i32 %109, 1
store i32 %110, i32* %14
br label %111
111: ; preds = %108
%112 = load i32, i32* %12
%113 = add nsw i32 %112, 1
store i32 %113, i32* %12
br label %83
114: ; preds = %83
%115 = load double*, double** %6
%116 = call i32 @free(double* %115)
store double 0.000000e+00, double* %3
br label %117
117: ; preds = %114, %48, %18
%118 = load double, double* %3
ret double %118
}

LLVM IR (-O2):
@.str = [13 x i8] c"num-utils-ng\00"
@EXIT_FAILURE = global double 0.000000e+00
@.str.1 = [4 x i8] c"%lf\00"
@EOF = global i32 0
@.str.2 = [5 x i8] c"%lf\0A\00"
define double @column(i32* %0, i32* %1) {
%3 = alloca double
%4 = bitcast double* %3 to i8*
%5 = tail call i64 @calloc(i32 1, i32 8)
%6 = inttoptr i64 %5 to double*
%7 = icmp eq i64 %5, 0
br i1 %7, label %12, label %8
8: ; preds = %2
%9 = call i32 @fscanf(i32* %0, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.1, i64 0, i64 0), double* %3)
%10 = load i32, i32* @EOF
%11 = icmp eq i32 %9, %10
br i1 %11, label %17, label %23
12: ; preds = %2
%13 = tail call i32 @perror(i8* getelementptr inbounds ([13 x i8], [13 x i8]* @.str, i64 0, i64 0))
%14 = load double, double* @EXIT_FAILURE
br label %91
15: ; preds = %63
%16 = icmp slt i32 %64, 0
br i1 %16, label %88, label %17
17: ; preds = %8, %15
%18 = phi i32 [ %64, %15 ], [ 0, %8 ]
%19 = phi double* [ %44, %15 ], [ %6, %8 ]
%20 = icmp eq i32* %1, null
%21 = add nuw i32 %18, 1
%22 = zext i32 %21 to i64
br label %69
23: ; preds = %8, %63
%24 = phi i32 [ %66, %63 ], [ %9, %8 ]
%25 = phi i32 [ %65, %63 ], [ 0, %8 ]
%26 = phi i32 [ %64, %63 ], [ 0, %8 ]
%27 = phi i32 [ %45, %63 ], [ 1, %8 ]
%28 = phi double* [ %44, %63 ], [ %6, %8 ]
%29 = icmp eq i32 %24, 0
br i1 %29, label %30, label %32
30: ; preds = %23
%31 = call i32 @skipWord(i32* %0)
br label %32
32: ; preds = %30, %23
%33 = icmp eq i32 %25, %27
br i1 %33, label %34, label %43
34: ; preds = %32
%35 = shl nsw i32 %25, 1
%36 = shl i32 %25, 4
%37 = call i64 @realloc(double* %28, i32 %36)
%38 = inttoptr i64 %37 to double*
%39 = icmp eq i64 %37, 0
br i1 %39, label %40, label %43
40: ; preds = %34
%41 = call i32 @perror(i8* getelementptr inbounds ([13 x i8], [13 x i8]* @.str, i64 0, i64 0))
%42 = load double, double* @EXIT_FAILURE
br label %91
43: ; preds = %34, %32
%44 = phi double* [ %38, %34 ], [ %28, %32 ]
%45 = phi i32 [ %35, %34 ], [ %27, %32 ]
%46 = call signext i8 @fgetc(i32* %0)
%47 = icmp eq i8 %46, 10
br i1 %47, label %55, label %48
48: ; preds = %43
%49 = load double, double* %3
%50 = sext i32 %25 to i64
%51 = getelementptr inbounds double, double* %44, i64 %50
%52 = load double, double* %51
%53 = fadd double %49, %52
store double %53, double* %51
%54 = add nsw i32 %25, 1
br label %63
55: ; preds = %43
%56 = icmp slt i32 %25, %26
%57 = select i1 %56, i32 %26, i32 %25
%58 = load double, double* %3
%59 = sext i32 %25 to i64
%60 = getelementptr inbounds double, double* %44, i64 %59
%61 = load double, double* %60
%62 = fadd double %58, %61
store double %62, double* %60
br label %63
63: ; preds = %55, %48
%64 = phi i32 [ %26, %48 ], [ %57, %55 ]
%65 = phi i32 [ %54, %48 ], [ 0, %55 ]
%66 = call i32 @fscanf(i32* %0, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.1, i64 0, i64 0), double* %3)
%67 = load i32, i32* @EOF
%68 = icmp eq i32 %66, %67
br i1 %68, label %15, label %23
69: ; preds = %17, %83
%70 = phi i64 [ 0, %17 ], [ %86, %83 ]
%71 = phi i32 [ 1, %17 ], [ %85, %83 ]
%72 = phi i32 [ 0, %17 ], [ %84, %83 ]
br i1 %20, label %78, label %73
73: ; preds = %69
%74 = sext i32 %72 to i64
%75 = getelementptr inbounds i32, i32* %1, i64 %74
%76 = load i32, i32* %75
%77 = icmp eq i32 %71, %76
br i1 %77, label %78, label %83
78: ; preds = %73, %69
%79 = getelementptr inbounds double, double* %19, i64 %70
%80 = load double, double* %79
%81 = call i32 @printf(i8* getelementptr inbounds ([5 x i8], [5 x i8]* @.str.2, i64 0, i64 0), double %80)
%82 = add nsw i32 %72, 1
br label %83
83: ; preds = %78, %73
%84 = phi i32 [ %82, %78 ], [ %72, %73 ]
%85 = add nuw nsw i32 %71, 1
%86 = add nuw nsw i64 %70, 1
%87 = icmp eq i64 %86, %22
br i1 %87, label %88, label %69
88: ; preds = %83, %15
%89 = phi double* [ %44, %15 ], [ %19, %83 ]
%90 = call i32 @free(double* %89)
br label %91
91: ; preds = %88, %40, %12
%92 = phi double [ %42, %40 ], [ 0.000000e+00, %88 ], [ %14, %12 ]
ret double %92
}