#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 /*<<< orphan*/ fclose(int /*<<< orphan*/ *);
int /*<<< orphan*/ *fopen(char *, char *);
scalar_t__ getc(int /*<<< orphan*/ *);
void mp_mul_f_cmuladd(int nfft, double dinout[], char *f_name_in2,
                      char *f_name_in3) {
  int j, k;
  double xr, xi, yr, yi, din2[3], din3[3];
  char *cin2, *cin3;
  FILE *f_in2, *f_in3;
  cin2 = (char *)din2;
  cin3 = (char *)din3;
  f_in2 = fopen(f_name_in2, "rb");
  f_in3 = fopen(f_name_in3, "rb");
  for (k = 0; k < 3 * (int)sizeof(double); k++) {
    cin2[k] = (char)getc(f_in2);
    cin3[k] = (char)getc(f_in3);
  }
  dinout[0] = din3[0];
  dinout[1] = din3[1] + dinout[1] * din2[1];
  dinout[2] = din3[2] + dinout[2] * din2[2];
  for (j = 3; j < nfft; j += 2) {
    for (k = 0; k < 2 * (int)sizeof(double); k++) {
      cin2[k] = (char)getc(f_in2);
      cin3[k] = (char)getc(f_in3);
    }
    xr = dinout[j];
    xi = dinout[j + 1];
    yr = din2[0];
    yi = din2[1];
    dinout[j] = din3[0] + (xr * yr - xi * yi);
    dinout[j + 1] = din3[1] + (xr * yi + xi * yr);
  }
  for (k = 0; k < (int)sizeof(double); k++) {
    cin2[k] = (char)getc(f_in2);
    cin3[k] = (char)getc(f_in3);
  }
  dinout[nfft + 1] = din3[0] + dinout[nfft + 1] * din2[0];
  fclose(f_in2);
  fclose(f_in3);
}


x86 assembly (clang, -O0):
<mp_mul_f_cmuladd>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
lea    -0x60(%rbp),%rax
mov    %rax,-0x88(%rbp)
lea    -0x80(%rbp),%rax
mov    %rax,-0x90(%rbp)
mov    -0x18(%rbp),%rdi
lea    0xe95(%rip),%rsi
call   1040 <fopen@plt>
mov    %rax,-0x98(%rbp)
mov    -0x20(%rbp),%rdi
lea    0xe7e(%rip),%rsi
call   1040 <fopen@plt>
mov    %rax,-0xa0(%rbp)
movl   $0x0,-0x28(%rbp)
cmpl   $0x18,-0x28(%rbp)
jge    11e5 <mp_mul_f_cmuladd+0xb5>
mov    -0x98(%rbp),%rdi
call   1050 <getc@plt>
mov    %al,%dl
mov    -0x88(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0xa0(%rbp),%rdi
call   1050 <getc@plt>
mov    %al,%dl
mov    -0x90(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1195 <mp_mul_f_cmuladd+0x65>
movsd  -0x80(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x78(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  -0x58(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movsd  -0x70(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  -0x50(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x10(%rax)
movl   $0x3,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    136f <mp_mul_f_cmuladd+0x23f>
movl   $0x0,-0x28(%rbp)
cmpl   $0x10,-0x28(%rbp)
jge    12a4 <mp_mul_f_cmuladd+0x174>
mov    -0x98(%rbp),%rdi
call   1050 <getc@plt>
mov    %al,%dl
mov    -0x88(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0xa0(%rbp),%rdi
call   1050 <getc@plt>
mov    %al,%dl
mov    -0x90(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1254 <mp_mul_f_cmuladd+0x124>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x80(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm3
movsd  -0x38(%rbp),%xmm2
mulsd  -0x48(%rbp),%xmm2
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x78(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm1
movsd  -0x48(%rbp),%xmm3
movsd  -0x38(%rbp),%xmm2
mulsd  -0x40(%rbp),%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x24(%rbp)
jmp    1241 <mp_mul_f_cmuladd+0x111>
movl   $0x0,-0x28(%rbp)
cmpl   $0x8,-0x28(%rbp)
jge    13c6 <mp_mul_f_cmuladd+0x296>
mov    -0x98(%rbp),%rdi
call   1050 <getc@plt>
mov    %al,%dl
mov    -0x88(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0xa0(%rbp),%rdi
call   1050 <getc@plt>
mov    %al,%dl
mov    -0x90(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1376 <mp_mul_f_cmuladd+0x246>
movsd  -0x80(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x60(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x98(%rbp),%rdi
call   1030 <fclose@plt>
mov    -0xa0(%rbp),%rdi
call   1030 <fclose@plt>
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<mp_mul_f_cmuladd>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %rcx,%rbp
mov    %rsi,%r14
mov    %edi,%r13d
lea    0xec2(%rip),%r12
mov    %rdx,%rdi
mov    %r12,%rsi
call   1040 <fopen@plt>
mov    %rax,%r15
mov    %rbp,%rdi
mov    %r12,%rsi
call   1040 <fopen@plt>
mov    %rax,%rbp
xor    %ebx,%ebx
nopl   0x0(%rax)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x20(%rsp,%rbx,1)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,(%rsp,%rbx,1)
add    $0x1,%rbx
cmp    $0x18,%rbx
jne    1170 <mp_mul_f_cmuladd+0x40>
movsd  (%rsp),%xmm0
movsd  %xmm0,(%r14)
movupd 0x8(%rsp),%xmm0
movupd 0x8(%r14),%xmm1
movupd 0x28(%rsp),%xmm2
mulpd  %xmm1,%xmm2
addpd  %xmm0,%xmm2
movupd %xmm2,0x8(%r14)
mov    %r13d,%ebx
cmp    $0x4,%r13d
jl     13af <mp_mul_f_cmuladd+0x27f>
mov    %ebx,%r12d
mov    $0x3,%r13d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x20(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x21(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x1(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x22(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x2(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x23(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x3(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x24(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x4(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x25(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x5(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x26(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x6(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x27(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x7(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x28(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x8(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x29(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x9(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x2a(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0xa(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x2b(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0xb(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x2c(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0xc(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x2d(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0xd(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x2e(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0xe(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x2f(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0xf(%rsp)
movapd 0x20(%rsp),%xmm0
movsd  (%r14,%r13,8),%xmm1
movsd  0x8(%r14,%r13,8),%xmm2
movapd %xmm2,%xmm3
xorpd  0xc82(%rip),%xmm3
unpcklpd %xmm2,%xmm3
unpcklpd %xmm1,%xmm1
mulpd  %xmm0,%xmm1
shufpd $0x1,%xmm0,%xmm0
mulpd  %xmm3,%xmm0
addpd  %xmm0,%xmm1
addpd  (%rsp),%xmm1
movupd %xmm1,(%r14,%r13,8)
add    $0x2,%r13
cmp    %r12,%r13
jb     11e0 <mp_mul_f_cmuladd+0xb0>
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x20(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x21(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x1(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x22(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x2(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x23(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x3(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x24(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x4(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x25(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x5(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x26(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x6(%rsp)
mov    %r15,%rdi
call   1050 <getc@plt>
mov    %al,0x27(%rsp)
mov    %rbp,%rdi
call   1050 <getc@plt>
mov    %al,0x7(%rsp)
movslq %ebx,%rax
movsd  0x8(%r14,%rax,8),%xmm0
mulsd  0x20(%rsp),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,0x8(%r14,%rax,8)
mov    %r15,%rdi
call   1030 <fclose@plt>
mov    %rbp,%rdi
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <fclose@plt>

LLVM IR (-O0):
@.str = [3 x i8] c"rb\00"
define void @mp_mul_f_cmuladd(i32 %0, double* %1, i8* %2, i8* %3) {
%5 = alloca i32
%6 = alloca double*
%7 = alloca i8*
%8 = alloca i8*
%9 = alloca i32
%10 = alloca i32
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca [3 x double]
%16 = alloca [3 x double]
%17 = alloca i8*
%18 = alloca i8*
%19 = alloca i32*
%20 = alloca i32*
store i32 %0, i32* %5
store double* %1, double** %6
store i8* %2, i8** %7
store i8* %3, i8** %8
%21 = getelementptr inbounds [3 x double], [3 x double]* %15, i64 0, i64 0
%22 = bitcast double* %21 to i8*
store i8* %22, i8** %17
%23 = getelementptr inbounds [3 x double], [3 x double]* %16, i64 0, i64 0
%24 = bitcast double* %23 to i8*
store i8* %24, i8** %18
%25 = load i8*, i8** %7
%26 = call i32* @fopen(i8* %25, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0))
store i32* %26, i32** %19
%27 = load i8*, i8** %8
%28 = call i32* @fopen(i8* %27, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0))
store i32* %28, i32** %20
store i32 0, i32* %10
br label %29
29: ; preds = %47, %4
%30 = load i32, i32* %10
%31 = icmp slt i32 %30, 24
br i1 %31, label %32, label %50
32: ; preds = %29
%33 = load i32*, i32** %19
%34 = call i64 @getc(i32* %33)
%35 = trunc i64 %34 to i8
%36 = load i8*, i8** %17
%37 = load i32, i32* %10
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds i8, i8* %36, i64 %38
store i8 %35, i8* %39
%40 = load i32*, i32** %20
%41 = call i64 @getc(i32* %40)
%42 = trunc i64 %41 to i8
%43 = load i8*, i8** %18
%44 = load i32, i32* %10
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds i8, i8* %43, i64 %45
store i8 %42, i8* %46
br label %47
47: ; preds = %32
%48 = load i32, i32* %10
%49 = add nsw i32 %48, 1
store i32 %49, i32* %10
br label %29
50: ; preds = %29
%51 = getelementptr inbounds [3 x double], [3 x double]* %16, i64 0, i64 0
%52 = load double, double* %51
%53 = load double*, double** %6
%54 = getelementptr inbounds double, double* %53, i64 0
store double %52, double* %54
%55 = getelementptr inbounds [3 x double], [3 x double]* %16, i64 0, i64 1
%56 = load double, double* %55
%57 = load double*, double** %6
%58 = getelementptr inbounds double, double* %57, i64 1
%59 = load double, double* %58
%60 = getelementptr inbounds [3 x double], [3 x double]* %15, i64 0, i64 1
%61 = load double, double* %60
%62 = call double @llvm.fmuladd.f64(double %59, double %61, double %56)
%63 = load double*, double** %6
%64 = getelementptr inbounds double, double* %63, i64 1
store double %62, double* %64
%65 = getelementptr inbounds [3 x double], [3 x double]* %16, i64 0, i64 2
%66 = load double, double* %65
%67 = load double*, double** %6
%68 = getelementptr inbounds double, double* %67, i64 2
%69 = load double, double* %68
%70 = getelementptr inbounds [3 x double], [3 x double]* %15, i64 0, i64 2
%71 = load double, double* %70
%72 = call double @llvm.fmuladd.f64(double %69, double %71, double %66)
%73 = load double*, double** %6
%74 = getelementptr inbounds double, double* %73, i64 2
store double %72, double* %74
store i32 3, i32* %9
br label %75
75: ; preds = %145, %50
%76 = load i32, i32* %9
%77 = load i32, i32* %5
%78 = icmp slt i32 %76, %77
br i1 %78, label %79, label %148
79: ; preds = %75
store i32 0, i32* %10
br label %80
80: ; preds = %98, %79
%81 = load i32, i32* %10
%82 = icmp slt i32 %81, 16
br i1 %82, label %83, label %101
83: ; preds = %80
%84 = load i32*, i32** %19
%85 = call i64 @getc(i32* %84)
%86 = trunc i64 %85 to i8
%87 = load i8*, i8** %17
%88 = load i32, i32* %10
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds i8, i8* %87, i64 %89
store i8 %86, i8* %90
%91 = load i32*, i32** %20
%92 = call i64 @getc(i32* %91)
%93 = trunc i64 %92 to i8
%94 = load i8*, i8** %18
%95 = load i32, i32* %10
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds i8, i8* %94, i64 %96
store i8 %93, i8* %97
br label %98
98: ; preds = %83
%99 = load i32, i32* %10
%100 = add nsw i32 %99, 1
store i32 %100, i32* %10
br label %80
101: ; preds = %80
%102 = load double*, double** %6
%103 = load i32, i32* %9
%104 = sext i32 %103 to i64
%105 = getelementptr inbounds double, double* %102, i64 %104
%106 = load double, double* %105
store double %106, double* %11
%107 = load double*, double** %6
%108 = load i32, i32* %9
%109 = add nsw i32 %108, 1
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds double, double* %107, i64 %110
%112 = load double, double* %111
store double %112, double* %12
%113 = getelementptr inbounds [3 x double], [3 x double]* %15, i64 0, i64 0
%114 = load double, double* %113
store double %114, double* %13
%115 = getelementptr inbounds [3 x double], [3 x double]* %15, i64 0, i64 1
%116 = load double, double* %115
store double %116, double* %14
%117 = getelementptr inbounds [3 x double], [3 x double]* %16, i64 0, i64 0
%118 = load double, double* %117
%119 = load double, double* %11
%120 = load double, double* %13
%121 = load double, double* %12
%122 = load double, double* %14
%123 = fmul double %121, %122
%124 = fneg double %123
%125 = call double @llvm.fmuladd.f64(double %119, double %120, double %124)
%126 = fadd double %118, %125
%127 = load double*, double** %6
%128 = load i32, i32* %9
%129 = sext i32 %128 to i64
%130 = getelementptr inbounds double, double* %127, i64 %129
store double %126, double* %130
%131 = getelementptr inbounds [3 x double], [3 x double]* %16, i64 0, i64 1
%132 = load double, double* %131
%133 = load double, double* %11
%134 = load double, double* %14
%135 = load double, double* %12
%136 = load double, double* %13
%137 = fmul double %135, %136
%138 = call double @llvm.fmuladd.f64(double %133, double %134, double %137)
%139 = fadd double %132, %138
%140 = load double*, double** %6
%141 = load i32, i32* %9
%142 = add nsw i32 %141, 1
%143 = sext i32 %142 to i64
%144 = getelementptr inbounds double, double* %140, i64 %143
store double %139, double* %144
br label %145
145: ; preds = %101
%146 = load i32, i32* %9
%147 = add nsw i32 %146, 2
store i32 %147, i32* %9
br label %75
148: ; preds = %75
store i32 0, i32* %10
br label %149
149: ; preds = %167, %148
%150 = load i32, i32* %10
%151 = icmp slt i32 %150, 8
br i1 %151, label %152, label %170
152: ; preds = %149
%153 = load i32*, i32** %19
%154 = call i64 @getc(i32* %153)
%155 = trunc i64 %154 to i8
%156 = load i8*, i8** %17
%157 = load i32, i32* %10
%158 = sext i32 %157 to i64
%159 = getelementptr inbounds i8, i8* %156, i64 %158
store i8 %155, i8* %159
%160 = load i32*, i32** %20
%161 = call i64 @getc(i32* %160)
%162 = trunc i64 %161 to i8
%163 = load i8*, i8** %18
%164 = load i32, i32* %10
%165 = sext i32 %164 to i64
%166 = getelementptr inbounds i8, i8* %163, i64 %165
store i8 %162, i8* %166
br label %167
167: ; preds = %152
%168 = load i32, i32* %10
%169 = add nsw i32 %168, 1
store i32 %169, i32* %10
br label %149
170: ; preds = %149
%171 = getelementptr inbounds [3 x double], [3 x double]* %16, i64 0, i64 0
%172 = load double, double* %171
%173 = load double*, double** %6
%174 = load i32, i32* %5
%175 = add nsw i32 %174, 1
%176 = sext i32 %175 to i64
%177 = getelementptr inbounds double, double* %173, i64 %176
%178 = load double, double* %177
%179 = getelementptr inbounds [3 x double], [3 x double]* %15, i64 0, i64 0
%180 = load double, double* %179
%181 = call double @llvm.fmuladd.f64(double %178, double %180, double %172)
%182 = load double*, double** %6
%183 = load i32, i32* %5
%184 = add nsw i32 %183, 1
%185 = sext i32 %184 to i64
%186 = getelementptr inbounds double, double* %182, i64 %185
store double %181, double* %186
%187 = load i32*, i32** %19
%188 = call i32 @fclose(i32* %187)
%189 = load i32*, i32** %20
%190 = call i32 @fclose(i32* %189)
ret void
}

LLVM IR (-O2):
@.str = [3 x i8] c"rb\00"
define void @mp_mul_f_cmuladd(i32 %0, double* nocapture %1, i8* nocapture %2, i8* nocapture %3) {
%5 = alloca [3 x double]
%6 = alloca [3 x double]
%7 = bitcast [3 x double]* %5 to i8*
%8 = bitcast [3 x double]* %6 to i8*
%9 = tail call i32* @fopen(i8* %2, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0))
%10 = tail call i32* @fopen(i8* %3, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0))
br label %11
11: ; preds = %4, %11
%12 = phi i64 [ 0, %4 ], [ %19, %11 ]
%13 = tail call i64 @getc(i32* %9)
%14 = trunc i64 %13 to i8
%15 = getelementptr inbounds i8, i8* %7, i64 %12
store i8 %14, i8* %15
%16 = tail call i64 @getc(i32* %10)
%17 = trunc i64 %16 to i8
%18 = getelementptr inbounds i8, i8* %8, i64 %12
store i8 %17, i8* %18
%19 = add nuw nsw i64 %12, 1
%20 = icmp eq i64 %19, 24
br i1 %20, label %21, label %11
21: ; preds = %11
%22 = getelementptr inbounds [3 x double], [3 x double]* %6, i64 0, i64 0
%23 = load double, double* %22
store double %23, double* %1
%24 = getelementptr inbounds [3 x double], [3 x double]* %6, i64 0, i64 1
%25 = getelementptr inbounds double, double* %1, i64 1
%26 = getelementptr inbounds [3 x double], [3 x double]* %5, i64 0, i64 1
%27 = bitcast double* %24 to <2 x double>*
%28 = load <2 x double>, <2 x double>* %27
%29 = bitcast double* %25 to <2 x double>*
%30 = load <2 x double>, <2 x double>* %29
%31 = bitcast double* %26 to <2 x double>*
%32 = load <2 x double>, <2 x double>* %31
%33 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %30, <2 x double> %32, <2 x double> %28)
%34 = bitcast double* %25 to <2 x double>*
store <2 x double> %33, <2 x double>* %34
%35 = icmp sgt i32 %0, 3
br i1 %35, label %36, label %157
36: ; preds = %21
%37 = zext i32 %0 to i64
%38 = getelementptr inbounds i8, i8* %7, i64 1
%39 = getelementptr inbounds i8, i8* %8, i64 1
%40 = getelementptr inbounds i8, i8* %7, i64 2
%41 = getelementptr inbounds i8, i8* %8, i64 2
%42 = getelementptr inbounds i8, i8* %7, i64 3
%43 = getelementptr inbounds i8, i8* %8, i64 3
%44 = getelementptr inbounds i8, i8* %7, i64 4
%45 = getelementptr inbounds i8, i8* %8, i64 4
%46 = getelementptr inbounds i8, i8* %7, i64 5
%47 = getelementptr inbounds i8, i8* %8, i64 5
%48 = getelementptr inbounds i8, i8* %7, i64 6
%49 = getelementptr inbounds i8, i8* %8, i64 6
%50 = getelementptr inbounds i8, i8* %7, i64 7
%51 = getelementptr inbounds i8, i8* %8, i64 7
%52 = getelementptr inbounds [3 x double], [3 x double]* %5, i64 0, i64 1
%53 = bitcast double* %52 to i8*
%54 = getelementptr inbounds [3 x double], [3 x double]* %6, i64 0, i64 1
%55 = bitcast double* %54 to i8*
%56 = getelementptr inbounds i8, i8* %7, i64 9
%57 = getelementptr inbounds i8, i8* %8, i64 9
%58 = getelementptr inbounds i8, i8* %7, i64 10
%59 = getelementptr inbounds i8, i8* %8, i64 10
%60 = getelementptr inbounds i8, i8* %7, i64 11
%61 = getelementptr inbounds i8, i8* %8, i64 11
%62 = getelementptr inbounds i8, i8* %7, i64 12
%63 = getelementptr inbounds i8, i8* %8, i64 12
%64 = getelementptr inbounds i8, i8* %7, i64 13
%65 = getelementptr inbounds i8, i8* %8, i64 13
%66 = getelementptr inbounds i8, i8* %7, i64 14
%67 = getelementptr inbounds i8, i8* %8, i64 14
%68 = getelementptr inbounds i8, i8* %7, i64 15
%69 = getelementptr inbounds i8, i8* %8, i64 15
%70 = bitcast [3 x double]* %5 to <2 x double>*
%71 = bitcast [3 x double]* %6 to <2 x double>*
br label %72
72: ; preds = %36, %72
%73 = phi i64 [ 3, %36 ], [ %155, %72 ]
%74 = tail call i64 @getc(i32* %9)
%75 = trunc i64 %74 to i8
store i8 %75, i8* %7
%76 = tail call i64 @getc(i32* %10)
%77 = trunc i64 %76 to i8
store i8 %77, i8* %8
%78 = tail call i64 @getc(i32* %9)
%79 = trunc i64 %78 to i8
store i8 %79, i8* %38
%80 = tail call i64 @getc(i32* %10)
%81 = trunc i64 %80 to i8
store i8 %81, i8* %39
%82 = tail call i64 @getc(i32* %9)
%83 = trunc i64 %82 to i8
store i8 %83, i8* %40
%84 = tail call i64 @getc(i32* %10)
%85 = trunc i64 %84 to i8
store i8 %85, i8* %41
%86 = tail call i64 @getc(i32* %9)
%87 = trunc i64 %86 to i8
store i8 %87, i8* %42
%88 = tail call i64 @getc(i32* %10)
%89 = trunc i64 %88 to i8
store i8 %89, i8* %43
%90 = tail call i64 @getc(i32* %9)
%91 = trunc i64 %90 to i8
store i8 %91, i8* %44
%92 = tail call i64 @getc(i32* %10)
%93 = trunc i64 %92 to i8
store i8 %93, i8* %45
%94 = tail call i64 @getc(i32* %9)
%95 = trunc i64 %94 to i8
store i8 %95, i8* %46
%96 = tail call i64 @getc(i32* %10)
%97 = trunc i64 %96 to i8
store i8 %97, i8* %47
%98 = tail call i64 @getc(i32* %9)
%99 = trunc i64 %98 to i8
store i8 %99, i8* %48
%100 = tail call i64 @getc(i32* %10)
%101 = trunc i64 %100 to i8
store i8 %101, i8* %49
%102 = tail call i64 @getc(i32* %9)
%103 = trunc i64 %102 to i8
store i8 %103, i8* %50
%104 = tail call i64 @getc(i32* %10)
%105 = trunc i64 %104 to i8
store i8 %105, i8* %51
%106 = tail call i64 @getc(i32* %9)
%107 = trunc i64 %106 to i8
store i8 %107, i8* %53
%108 = tail call i64 @getc(i32* %10)
%109 = trunc i64 %108 to i8
store i8 %109, i8* %55
%110 = tail call i64 @getc(i32* %9)
%111 = trunc i64 %110 to i8
store i8 %111, i8* %56
%112 = tail call i64 @getc(i32* %10)
%113 = trunc i64 %112 to i8
store i8 %113, i8* %57
%114 = tail call i64 @getc(i32* %9)
%115 = trunc i64 %114 to i8
store i8 %115, i8* %58
%116 = tail call i64 @getc(i32* %10)
%117 = trunc i64 %116 to i8
store i8 %117, i8* %59
%118 = tail call i64 @getc(i32* %9)
%119 = trunc i64 %118 to i8
store i8 %119, i8* %60
%120 = tail call i64 @getc(i32* %10)
%121 = trunc i64 %120 to i8
store i8 %121, i8* %61
%122 = tail call i64 @getc(i32* %9)
%123 = trunc i64 %122 to i8
store i8 %123, i8* %62
%124 = tail call i64 @getc(i32* %10)
%125 = trunc i64 %124 to i8
store i8 %125, i8* %63
%126 = tail call i64 @getc(i32* %9)
%127 = trunc i64 %126 to i8
store i8 %127, i8* %64
%128 = tail call i64 @getc(i32* %10)
%129 = trunc i64 %128 to i8
store i8 %129, i8* %65
%130 = tail call i64 @getc(i32* %9)
%131 = trunc i64 %130 to i8
store i8 %131, i8* %66
%132 = tail call i64 @getc(i32* %10)
%133 = trunc i64 %132 to i8
store i8 %133, i8* %67
%134 = tail call i64 @getc(i32* %9)
%135 = trunc i64 %134 to i8
store i8 %135, i8* %68
%136 = tail call i64 @getc(i32* %10)
%137 = trunc i64 %136 to i8
store i8 %137, i8* %69
%138 = getelementptr inbounds double, double* %1, i64 %73
%139 = load double, double* %138
%140 = add nuw nsw i64 %73, 1
%141 = getelementptr inbounds double, double* %1, i64 %140
%142 = load double, double* %141
%143 = load <2 x double>, <2 x double>* %70
%144 = fneg double %142
%145 = load <2 x double>, <2 x double>* %71
%146 = insertelement <2 x double> poison, double %142, i64 0
%147 = insertelement <2 x double> %146, double %144, i64 1
%148 = fmul <2 x double> %143, %147
%149 = shufflevector <2 x double> %148, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%150 = insertelement <2 x double> poison, double %139, i64 0
%151 = shufflevector <2 x double> %150, <2 x double> poison, <2 x i32> zeroinitializer
%152 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %151, <2 x double> %143, <2 x double> %149)
%153 = fadd <2 x double> %145, %152
%154 = bitcast double* %138 to <2 x double>*
store <2 x double> %153, <2 x double>* %154
%155 = add nuw nsw i64 %73, 2
%156 = icmp ult i64 %155, %37
br i1 %156, label %72, label %157
157: ; preds = %72, %21
%158 = getelementptr inbounds [3 x double], [3 x double]* %5, i64 0, i64 0
%159 = tail call i64 @getc(i32* %9)
%160 = trunc i64 %159 to i8
store i8 %160, i8* %7
%161 = tail call i64 @getc(i32* %10)
%162 = trunc i64 %161 to i8
store i8 %162, i8* %8
%163 = tail call i64 @getc(i32* %9)
%164 = trunc i64 %163 to i8
%165 = getelementptr inbounds i8, i8* %7, i64 1
store i8 %164, i8* %165
%166 = tail call i64 @getc(i32* %10)
%167 = trunc i64 %166 to i8
%168 = getelementptr inbounds i8, i8* %8, i64 1
store i8 %167, i8* %168
%169 = tail call i64 @getc(i32* %9)
%170 = trunc i64 %169 to i8
%171 = getelementptr inbounds i8, i8* %7, i64 2
store i8 %170, i8* %171
%172 = tail call i64 @getc(i32* %10)
%173 = trunc i64 %172 to i8
%174 = getelementptr inbounds i8, i8* %8, i64 2
store i8 %173, i8* %174
%175 = tail call i64 @getc(i32* %9)
%176 = trunc i64 %175 to i8
%177 = getelementptr inbounds i8, i8* %7, i64 3
store i8 %176, i8* %177
%178 = tail call i64 @getc(i32* %10)
%179 = trunc i64 %178 to i8
%180 = getelementptr inbounds i8, i8* %8, i64 3
store i8 %179, i8* %180
%181 = tail call i64 @getc(i32* %9)
%182 = trunc i64 %181 to i8
%183 = getelementptr inbounds i8, i8* %7, i64 4
store i8 %182, i8* %183
%184 = tail call i64 @getc(i32* %10)
%185 = trunc i64 %184 to i8
%186 = getelementptr inbounds i8, i8* %8, i64 4
store i8 %185, i8* %186
%187 = tail call i64 @getc(i32* %9)
%188 = trunc i64 %187 to i8
%189 = getelementptr inbounds i8, i8* %7, i64 5
store i8 %188, i8* %189
%190 = tail call i64 @getc(i32* %10)
%191 = trunc i64 %190 to i8
%192 = getelementptr inbounds i8, i8* %8, i64 5
store i8 %191, i8* %192
%193 = tail call i64 @getc(i32* %9)
%194 = trunc i64 %193 to i8
%195 = getelementptr inbounds i8, i8* %7, i64 6
store i8 %194, i8* %195
%196 = tail call i64 @getc(i32* %10)
%197 = trunc i64 %196 to i8
%198 = getelementptr inbounds i8, i8* %8, i64 6
store i8 %197, i8* %198
%199 = tail call i64 @getc(i32* %9)
%200 = trunc i64 %199 to i8
%201 = getelementptr inbounds i8, i8* %7, i64 7
store i8 %200, i8* %201
%202 = tail call i64 @getc(i32* %10)
%203 = trunc i64 %202 to i8
%204 = getelementptr inbounds i8, i8* %8, i64 7
store i8 %203, i8* %204
%205 = load double, double* %22
%206 = add nsw i32 %0, 1
%207 = sext i32 %206 to i64
%208 = getelementptr inbounds double, double* %1, i64 %207
%209 = load double, double* %208
%210 = load double, double* %158
%211 = tail call double @llvm.fmuladd.f64(double %209, double %210, double %205)
store double %211, double* %208
%212 = tail call i32 @fclose(i32* %9)
%213 = tail call i32 @fclose(i32* %10)
ret void
}