#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 */
float cosf(float);
float sinf(float);
int ax2omf(float *ax, float *om, int p) {
  float c, s, omc, q, om_temp[9];
  int i, j;
  float epsijk;
  epsijk = (float)p;
  c = cosf(ax[3]);
  s = sinf(ax[3]);
  omc = 1.0f - c;
  for (i = 0; i < 3; i++) {
    *(om + i + 3 * i) = ax[i] * ax[i] * omc + c;
  }
  q = omc * ax[0] * ax[1];
  *(om + 3) = q + s * ax[2];
  *(om + 1) = q - s * ax[2];
  q = omc * ax[1] * ax[2];
  *(om + 7) = q + s * ax[0];
  *(om + 5) = q - s * ax[0];
  q = omc * ax[2] * ax[0];
  *(om + 2) = q + s * ax[1];
  *(om + 6) = q - s * ax[1];
  if (p < 0) {
    for (i = 0; i < 9; i++) {
      om_temp[i] = om[i];
    }
    for (i = 0; i < 3; i++) {
      for (j = 0; j < 3; j++) {
        om[i + 3 * j] = om_temp[3 * i + j];
      }
    }
  }
  return 1;
}


x86 assembly (clang, -O0):
<ax2omf>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
cvtsi2ssl -0x14(%rbp),%xmm0
movss  %xmm0,-0x5c(%rbp)
mov    -0x8(%rbp),%rax
movss  0xc(%rax),%xmm0
call   1040 <cosf@plt>
movss  %xmm0,-0x18(%rbp)
mov    -0x8(%rbp),%rax
movss  0xc(%rax),%xmm0
call   1030 <sinf@plt>
movss  %xmm0,-0x1c(%rbp)
movss  0xe95(%rip),%xmm0
subss  -0x18(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
movl   $0x0,-0x54(%rbp)
cmpl   $0x3,-0x54(%rbp)
jge    11db <ax2omf+0xbb>
mov    -0x8(%rbp),%rax
movslq -0x54(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x54(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  -0x20(%rbp),%xmm2
movss  -0x18(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
imul   $0x3,-0x54(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    117c <ax2omf+0x5c>
movss  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mulss  (%rax),%xmm0
mov    -0x8(%rbp),%rax
mulss  0x4(%rax),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm1
movss  -0x1c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x8(%rax),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0xc(%rax)
movss  -0x24(%rbp),%xmm1
movss  -0x1c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x8(%rax),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0x4(%rax)
movss  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mulss  0x4(%rax),%xmm0
mov    -0x8(%rbp),%rax
mulss  0x8(%rax),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm1
movss  -0x1c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0x1c(%rax)
movss  -0x24(%rbp),%xmm1
movss  -0x1c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0x14(%rax)
movss  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mulss  0x8(%rax),%xmm0
mov    -0x8(%rbp),%rax
mulss  (%rax),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm1
movss  -0x1c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x4(%rax),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0x8(%rax)
movss  -0x24(%rbp),%xmm1
movss  -0x1c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movss  0x4(%rax),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movss  %xmm0,0x18(%rax)
cmpl   $0x0,-0x14(%rbp)
jge    13d6 <ax2omf+0x2b6>
movl   $0x0,-0x54(%rbp)
cmpl   $0x9,-0x54(%rbp)
jge    136a <ax2omf+0x24a>
mov    -0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movslq -0x54(%rbp),%rax
movss  %xmm0,-0x50(%rbp,%rax,4)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    133b <ax2omf+0x21b>
movl   $0x0,-0x54(%rbp)
cmpl   $0x3,-0x54(%rbp)
jge    13d1 <ax2omf+0x2b1>
movl   $0x0,-0x58(%rbp)
cmpl   $0x3,-0x58(%rbp)
jge    13be <ax2omf+0x29e>
imul   $0x3,-0x54(%rbp),%eax
add    -0x58(%rbp),%eax
cltq
movss  -0x50(%rbp,%rax,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x54(%rbp),%ecx
imul   $0x3,-0x58(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1382 <ax2omf+0x262>
jmp    13c3 <ax2omf+0x2a3>
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    1371 <ax2omf+0x251>
jmp    13d6 <ax2omf+0x2b6>
mov    $0x1,%eax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ax2omf>:
push   %rbp
push   %r14
push   %rbx
sub    $0x10,%rsp
mov    %edx,%r14d
mov    %rsi,%rbx
mov    %rdi,%rbp
movss  0xc(%rdi),%xmm0
call   1040 <cosf@plt>
movss  %xmm0,0xc(%rsp)
movss  0xc(%rbp),%xmm0
call   1030 <sinf@plt>
movss  0xead(%rip),%xmm1
movss  0xc(%rsp),%xmm2
subss  %xmm2,%xmm1
movss  0x0(%rbp),%xmm8
mulss  %xmm8,%xmm8
mulss  %xmm1,%xmm8
addss  %xmm2,%xmm8
movss  %xmm8,(%rbx)
movss  0x4(%rbp),%xmm3
mulss  %xmm3,%xmm3
mulss  %xmm1,%xmm3
addss  %xmm2,%xmm3
movss  %xmm3,0x10(%rbx)
movss  0x8(%rbp),%xmm4
mulss  %xmm4,%xmm4
mulss  %xmm1,%xmm4
addss  %xmm2,%xmm4
movss  %xmm4,0x20(%rbx)
movss  0x0(%rbp),%xmm5
mulss  %xmm1,%xmm5
mulss  0x4(%rbp),%xmm5
movss  0x8(%rbp),%xmm9
mulss  %xmm0,%xmm9
addss  %xmm5,%xmm9
movss  %xmm9,0xc(%rbx)
movss  0x8(%rbp),%xmm2
mulss  %xmm0,%xmm2
subss  %xmm2,%xmm5
movss  %xmm5,0x4(%rbx)
movss  0x4(%rbp),%xmm7
mulss  %xmm1,%xmm7
mulss  0x8(%rbp),%xmm7
movss  0x0(%rbp),%xmm2
mulss  %xmm0,%xmm2
addss  %xmm7,%xmm2
movss  %xmm2,0x1c(%rbx)
movss  0x0(%rbp),%xmm6
mulss  %xmm0,%xmm6
subss  %xmm6,%xmm7
movss  %xmm7,0x14(%rbx)
mulss  0x8(%rbp),%xmm1
mulss  0x0(%rbp),%xmm1
movss  0x4(%rbp),%xmm6
mulss  %xmm0,%xmm6
addss  %xmm1,%xmm6
movss  %xmm6,0x8(%rbx)
mulss  0x4(%rbp),%xmm0
subss  %xmm0,%xmm1
movss  %xmm1,0x18(%rbx)
test   %r14d,%r14d
jns    1268 <ax2omf+0x148>
movss  %xmm8,(%rbx)
movss  %xmm5,0xc(%rbx)
movss  %xmm6,0x18(%rbx)
movss  %xmm9,0x4(%rbx)
movss  %xmm3,0x10(%rbx)
movss  %xmm7,0x1c(%rbx)
movss  %xmm1,0x8(%rbx)
movss  %xmm2,0x14(%rbx)
movss  %xmm4,0x20(%rbx)
mov    $0x1,%eax
add    $0x10,%rsp
pop    %rbx
pop    %r14
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @ax2omf(float* nocapture %0, float* nocapture %1, i32 %2) {
%4 = getelementptr inbounds float, float* %0, i64 3
%5 = load float, float* %4
%6 = tail call float @cosf(float %5)
%7 = load float, float* %4
%8 = tail call float @sinf(float %7)
%9 = fsub float 1.000000e+00, %6
%10 = load float, float* %0
%11 = fmul float %10, %10
%12 = tail call float @llvm.fmuladd.f32(float %11, float %9, float %6)
store float %12, float* %1
%13 = getelementptr inbounds float, float* %0, i64 1
%14 = load float, float* %13
%15 = fmul float %14, %14
%16 = tail call float @llvm.fmuladd.f32(float %15, float %9, float %6)
%17 = getelementptr inbounds float, float* %1, i64 1
%18 = getelementptr inbounds float, float* %1, i64 4
store float %16, float* %18
%19 = getelementptr inbounds float, float* %0, i64 2
%20 = load float, float* %19
%21 = fmul float %20, %20
%22 = tail call float @llvm.fmuladd.f32(float %21, float %9, float %6)
%23 = getelementptr inbounds float, float* %1, i64 2
%24 = getelementptr inbounds float, float* %1, i64 8
store float %22, float* %24
%25 = load float, float* %0
%26 = fmul float %9, %25
%27 = load float, float* %13
%28 = fmul float %26, %27
%29 = load float, float* %19
%30 = tail call float @llvm.fmuladd.f32(float %8, float %29, float %28)
%31 = getelementptr inbounds float, float* %1, i64 3
store float %30, float* %31
%32 = load float, float* %19
%33 = fneg float %8
%34 = tail call float @llvm.fmuladd.f32(float %33, float %32, float %28)
store float %34, float* %17
%35 = load float, float* %13
%36 = fmul float %9, %35
%37 = load float, float* %19
%38 = fmul float %36, %37
%39 = load float, float* %0
%40 = tail call float @llvm.fmuladd.f32(float %8, float %39, float %38)
%41 = getelementptr inbounds float, float* %1, i64 7
store float %40, float* %41
%42 = load float, float* %0
%43 = tail call float @llvm.fmuladd.f32(float %33, float %42, float %38)
%44 = getelementptr inbounds float, float* %1, i64 5
store float %43, float* %44
%45 = load float, float* %19
%46 = fmul float %9, %45
%47 = load float, float* %0
%48 = fmul float %46, %47
%49 = load float, float* %13
%50 = tail call float @llvm.fmuladd.f32(float %8, float %49, float %48)
store float %50, float* %23
%51 = load float, float* %13
%52 = tail call float @llvm.fmuladd.f32(float %33, float %51, float %48)
%53 = getelementptr inbounds float, float* %1, i64 6
store float %52, float* %53
%54 = icmp slt i32 %2, 0
br i1 %54, label %55, label %58
55: ; preds = %3
%56 = getelementptr inbounds float, float* %1, i64 4
%57 = getelementptr inbounds float, float* %1, i64 8
store float %12, float* %1
store float %34, float* %31
store float %50, float* %53
store float %30, float* %17
store float %16, float* %56
store float %43, float* %41
store float %52, float* %23
store float %40, float* %44
store float %22, float* %57
br label %58
58: ; preds = %55, %3
ret i32 1
}