#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__ calloc(int, int);
int /*<<< orphan*/ free(double *);
void otrma(double *c, double *a, double *b, int n) {
  double z, *q0, *p, *s, *t;
  int i, j, k;
  q0 = (double *)calloc(n, sizeof(double));
  for (i = 0; i < n; ++i, ++c) {
    for (j = 0, t = b; j < n; ++j) {
      for (k = 0, s = a + i * n, z = 0.; k < n; ++k)
        z += *t++ * *s++;
      q0[j] = z;
    }
    for (j = 0, p = c, t = a; j < n; ++j, p += n) {
      for (k = 0, s = q0, z = 0.; k < n; ++k)
        z += *t++ * *s++;
      *p = z;
    }
  }
  free(q0);
}


x86 assembly (clang, -O0):
<otrma>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    -0x1c(%rbp),%edi
mov    $0x8,%esi
call   1040 <calloc@plt>
mov    %rax,-0x30(%rbp)
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12eb <otrma+0x1cb>
movl   $0x0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1210 <otrma+0xf0>
movl   $0x0,-0x54(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x54(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11f0 <otrma+0xd0>
mov    -0x48(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x48(%rbp)
movsd  (%rax),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    119e <otrma+0x7e>
movsd  -0x28(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    116a <otrma+0x4a>
movl   $0x0,-0x50(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12cc <otrma+0x1ac>
movl   $0x0,-0x54(%rbp)
mov    -0x30(%rbp),%rax
mov    %rax,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x54(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    129c <otrma+0x17c>
mov    -0x48(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x48(%rbp)
movsd  (%rax),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    124a <otrma+0x12a>
movsd  -0x28(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
mov    -0x1c(%rbp),%ecx
mov    -0x38(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
jmp    1227 <otrma+0x107>
jmp    12d1 <otrma+0x1b1>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x8(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x8(%rbp)
jmp    114f <otrma+0x2f>
mov    -0x30(%rbp),%rdi
call   1030 <free@plt>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<otrma>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %ecx,%r13d
mov    %rdx,0x10(%rsp)
mov    %rsi,(%rsp)
mov    %rdi,%r12
mov    %ecx,%edi
mov    $0x8,%esi
call   1040 <calloc@plt>
mov    %rax,%rdi
test   %r13d,%r13d
jle    135e <otrma+0x23e>
mov    %r13d,%r14d
lea    -0x1(%r13),%r10d
lea    0x0(,%r14,8),%rax
mov    %rax,0x8(%rsp)
mov    %r14d,%r11d
and    $0xfffffffc,%r11d
mov    %r13d,%ebx
and    $0x3,%ebx
shl    $0x3,%rbx
xor    %r9d,%r9d
mov    (%rsp),%rsi
jmp    11a6 <otrma+0x86>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r9
add    $0x8,%r12
add    0x8(%rsp),%rsi
cmp    %r14,%r9
je     135e <otrma+0x23e>
mov    %r9,%rax
imul   %r14,%rax
mov    (%rsp),%rcx
lea    (%rcx,%rax,8),%r8
xor    %edx,%edx
mov    0x10(%rsp),%rbp
jmp    11d2 <otrma+0xb2>
xchg   %ax,%ax
movsd  %xmm0,(%rdi,%rdx,8)
add    $0x1,%rdx
cmp    %r14,%rdx
je     1290 <otrma+0x170>
xorpd  %xmm0,%xmm0
cmp    $0x3,%r10d
jae    11f0 <otrma+0xd0>
mov    %r8,%rax
test   $0x3,%r13b
jne    125f <otrma+0x13f>
jmp    11c0 <otrma+0xa0>
nopw   0x0(%rax,%rax,1)
xorpd  %xmm0,%xmm0
xor    %eax,%eax
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rax,8),%xmm1
mulsd  (%rsi,%rax,8),%xmm1
movsd  0x8(%rbp,%rax,8),%xmm2
addsd  %xmm0,%xmm1
mulsd  0x8(%rsi,%rax,8),%xmm2
addsd  %xmm1,%xmm2
movsd  0x10(%rbp,%rax,8),%xmm1
mulsd  0x10(%rsi,%rax,8),%xmm1
addsd  %xmm2,%xmm1
movsd  0x18(%rbp,%rax,8),%xmm0
mulsd  0x18(%rsi,%rax,8),%xmm0
addsd  %xmm1,%xmm0
add    $0xffffffffffffffe0,%rcx
add    $0x4,%rax
cmp    %eax,%r11d
jne    1200 <otrma+0xe0>
sub    %rcx,%rbp
mov    %rsi,%rax
sub    %rcx,%rax
test   $0x3,%r13b
je     11c0 <otrma+0xa0>
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rcx,1),%xmm1
mulsd  (%rax,%rcx,1),%xmm1
addsd  %xmm1,%xmm0
add    $0x8,%rcx
cmp    %ecx,%ebx
jne    1270 <otrma+0x150>
add    %rcx,%rbp
jmp    11c0 <otrma+0xa0>
nop
xor    %r15d,%r15d
mov    (%rsp),%rcx
mov    %r12,%rbp
jmp    12b7 <otrma+0x197>
nopl   0x0(%rax)
movsd  %xmm0,0x0(%rbp)
add    $0x1,%r15d
lea    0x0(%rbp,%r14,8),%rbp
cmp    %r13d,%r15d
je     1190 <otrma+0x70>
xorpd  %xmm0,%xmm0
cmp    $0x3,%r10d
jae    12d0 <otrma+0x1b0>
mov    %rdi,%rax
test   $0x3,%r13b
jne    133e <otrma+0x21e>
jmp    12a0 <otrma+0x180>
nopl   0x0(%rax)
xor    %eax,%eax
xor    %r8d,%r8d
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rcx,%rax,8),%xmm1
mulsd  (%rdi,%rax,8),%xmm1
movsd  0x8(%rcx,%rax,8),%xmm2
addsd  %xmm0,%xmm1
mulsd  0x8(%rdi,%rax,8),%xmm2
addsd  %xmm1,%xmm2
movsd  0x10(%rcx,%rax,8),%xmm1
mulsd  0x10(%rdi,%rax,8),%xmm1
addsd  %xmm2,%xmm1
movsd  0x18(%rcx,%rax,8),%xmm0
mulsd  0x18(%rdi,%rax,8),%xmm0
addsd  %xmm1,%xmm0
add    $0xffffffffffffffe0,%r8
add    $0x4,%rax
cmp    %eax,%r11d
jne    12e0 <otrma+0x1c0>
sub    %r8,%rcx
mov    %rdi,%rax
sub    %r8,%rax
test   $0x3,%r13b
je     12a0 <otrma+0x180>
xor    %edx,%edx
movsd  (%rcx,%rdx,1),%xmm1
mulsd  (%rax,%rdx,1),%xmm1
addsd  %xmm1,%xmm0
add    $0x8,%rdx
cmp    %edx,%ebx
jne    1340 <otrma+0x220>
add    %rdx,%rcx
jmp    12a0 <otrma+0x180>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

LLVM IR (-O2):
define void @otrma(double* nocapture writeonly %0, double* nocapture %1, double* nocapture %2, i32 %3) {
%5 = tail call i64 @calloc(i32 %3, i32 8)
%6 = inttoptr i64 %5 to double*
%7 = icmp sgt i32 %3, 0
br i1 %7, label %8, label %138
8: ; preds = %4
%9 = zext i32 %3 to i64
%10 = zext i32 %3 to i64
%11 = add i32 %3, -1
%12 = and i32 %3, 3
%13 = icmp ult i32 %11, 3
%14 = and i32 %3, -4
%15 = icmp eq i32 %12, 0
%16 = and i32 %3, 3
%17 = icmp ult i32 %11, 3
%18 = and i32 %3, -4
%19 = icmp eq i32 %16, 0
br label %20
20: ; preds = %134, %8
%21 = phi i64 [ 0, %8 ], [ %135, %134 ]
%22 = phi double* [ %0, %8 ], [ %136, %134 ]
%23 = mul nsw i64 %21, %9
%24 = getelementptr inbounds double, double* %1, i64 %23
br label %25
25: ; preds = %73, %20
%26 = phi i64 [ 0, %20 ], [ %77, %73 ]
%27 = phi double* [ %2, %20 ], [ %74, %73 ]
br i1 %13, label %55, label %28
28: ; preds = %25, %28
%29 = phi double* [ %48, %28 ], [ %27, %25 ]
%30 = phi double* [ %50, %28 ], [ %24, %25 ]
%31 = phi double [ %52, %28 ], [ 0.000000e+00, %25 ]
%32 = phi i32 [ %53, %28 ], [ 0, %25 ]
%33 = getelementptr inbounds double, double* %29, i64 1
%34 = load double, double* %29
%35 = getelementptr inbounds double, double* %30, i64 1
%36 = load double, double* %30
%37 = tail call double @llvm.fmuladd.f64(double %34, double %36, double %31)
%38 = getelementptr inbounds double, double* %29, i64 2
%39 = load double, double* %33
%40 = getelementptr inbounds double, double* %30, i64 2
%41 = load double, double* %35
%42 = tail call double @llvm.fmuladd.f64(double %39, double %41, double %37)
%43 = getelementptr inbounds double, double* %29, i64 3
%44 = load double, double* %38
%45 = getelementptr inbounds double, double* %30, i64 3
%46 = load double, double* %40
%47 = tail call double @llvm.fmuladd.f64(double %44, double %46, double %42)
%48 = getelementptr inbounds double, double* %29, i64 4
%49 = load double, double* %43
%50 = getelementptr inbounds double, double* %30, i64 4
%51 = load double, double* %45
%52 = tail call double @llvm.fmuladd.f64(double %49, double %51, double %47)
%53 = add i32 %32, 4
%54 = icmp eq i32 %53, %14
br i1 %54, label %55, label %28
55: ; preds = %28, %25
%56 = phi double* [ undef, %25 ], [ %48, %28 ]
%57 = phi double [ undef, %25 ], [ %52, %28 ]
%58 = phi double* [ %27, %25 ], [ %48, %28 ]
%59 = phi double* [ %24, %25 ], [ %50, %28 ]
%60 = phi double [ 0.000000e+00, %25 ], [ %52, %28 ]
br i1 %15, label %73, label %61
61: ; preds = %55, %61
%62 = phi double* [ %66, %61 ], [ %58, %55 ]
%63 = phi double* [ %68, %61 ], [ %59, %55 ]
%64 = phi double [ %70, %61 ], [ %60, %55 ]
%65 = phi i32 [ %71, %61 ], [ 0, %55 ]
%66 = getelementptr inbounds double, double* %62, i64 1
%67 = load double, double* %62
%68 = getelementptr inbounds double, double* %63, i64 1
%69 = load double, double* %63
%70 = tail call double @llvm.fmuladd.f64(double %67, double %69, double %64)
%71 = add i32 %65, 1
%72 = icmp eq i32 %71, %12
br i1 %72, label %73, label %61
73: ; preds = %61, %55
%74 = phi double* [ %56, %55 ], [ %66, %61 ]
%75 = phi double [ %57, %55 ], [ %70, %61 ]
%76 = getelementptr inbounds double, double* %6, i64 %26
store double %75, double* %76
%77 = add nuw nsw i64 %26, 1
%78 = icmp eq i64 %77, %10
br i1 %78, label %79, label %25
79: ; preds = %73, %128
%80 = phi i32 [ %131, %128 ], [ 0, %73 ]
%81 = phi double* [ %129, %128 ], [ %1, %73 ]
%82 = phi double* [ %132, %128 ], [ %22, %73 ]
br i1 %17, label %110, label %83
83: ; preds = %79, %83
%84 = phi double* [ %103, %83 ], [ %81, %79 ]
%85 = phi double* [ %105, %83 ], [ %6, %79 ]
%86 = phi double [ %107, %83 ], [ 0.000000e+00, %79 ]
%87 = phi i32 [ %108, %83 ], [ 0, %79 ]
%88 = getelementptr inbounds double, double* %84, i64 1
%89 = load double, double* %84
%90 = getelementptr inbounds double, double* %85, i64 1
%91 = load double, double* %85
%92 = tail call double @llvm.fmuladd.f64(double %89, double %91, double %86)
%93 = getelementptr inbounds double, double* %84, i64 2
%94 = load double, double* %88
%95 = getelementptr inbounds double, double* %85, i64 2
%96 = load double, double* %90
%97 = tail call double @llvm.fmuladd.f64(double %94, double %96, double %92)
%98 = getelementptr inbounds double, double* %84, i64 3
%99 = load double, double* %93
%100 = getelementptr inbounds double, double* %85, i64 3
%101 = load double, double* %95
%102 = tail call double @llvm.fmuladd.f64(double %99, double %101, double %97)
%103 = getelementptr inbounds double, double* %84, i64 4
%104 = load double, double* %98
%105 = getelementptr inbounds double, double* %85, i64 4
%106 = load double, double* %100
%107 = tail call double @llvm.fmuladd.f64(double %104, double %106, double %102)
%108 = add i32 %87, 4
%109 = icmp eq i32 %108, %18
br i1 %109, label %110, label %83
110: ; preds = %83, %79
%111 = phi double* [ undef, %79 ], [ %103, %83 ]
%112 = phi double [ undef, %79 ], [ %107, %83 ]
%113 = phi double* [ %81, %79 ], [ %103, %83 ]
%114 = phi double* [ %6, %79 ], [ %105, %83 ]
%115 = phi double [ 0.000000e+00, %79 ], [ %107, %83 ]
br i1 %19, label %128, label %116
116: ; preds = %110, %116
%117 = phi double* [ %121, %116 ], [ %113, %110 ]
%118 = phi double* [ %123, %116 ], [ %114, %110 ]
%119 = phi double [ %125, %116 ], [ %115, %110 ]
%120 = phi i32 [ %126, %116 ], [ 0, %110 ]
%121 = getelementptr inbounds double, double* %117, i64 1
%122 = load double, double* %117
%123 = getelementptr inbounds double, double* %118, i64 1
%124 = load double, double* %118
%125 = tail call double @llvm.fmuladd.f64(double %122, double %124, double %119)
%126 = add i32 %120, 1
%127 = icmp eq i32 %126, %16
br i1 %127, label %128, label %116
128: ; preds = %116, %110
%129 = phi double* [ %111, %110 ], [ %121, %116 ]
%130 = phi double [ %112, %110 ], [ %125, %116 ]
store double %130, double* %82
%131 = add nuw nsw i32 %80, 1
%132 = getelementptr inbounds double, double* %82, i64 %9
%133 = icmp eq i32 %131, %3
br i1 %133, label %134, label %79
134: ; preds = %128
%135 = add nuw nsw i64 %21, 1
%136 = getelementptr inbounds double, double* %22, i64 1
%137 = icmp eq i64 %135, %10
br i1 %137, label %138, label %20
138: ; preds = %134, %4
%139 = tail call i32 @free(double* %6)
ret void
}