#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 */
int /*<<< orphan*/ free_vector(float *, int, unsigned long);
float *vector(int, unsigned long);
void tridag(float a[], float b[], float c[], float r[], float u[],
            unsigned long n) {
  unsigned long j;
  float bet, *gam;
  gam = vector(1, n);
  u[1] = r[1] / (bet = b[1]);
  for (j = 2; j <= n; j++) {
    gam[j] = c[j - 1] / bet;
    bet = b[j] - a[j] * gam[j];
    u[j] = (r[j] - a[j] * u[j - 1]) / bet;
  }
  for (j = (n - 1); j >= 1; j--)
    u[j] -= gam[j + 1] * u[j + 1];
  free_vector(gam, 1, n);
}


x86 assembly (clang, -O0):
<tridag>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x30(%rbp),%rsi
mov    $0x1,%edi
call   1040 <vector@plt>
mov    %rax,-0x48(%rbp)
mov    -0x20(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm1
movss  %xmm1,-0x3c(%rbp)
divss  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movss  %xmm0,0x4(%rax)
movq   $0x2,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x30(%rbp),%rax
ja     1253 <tridag+0x133>
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    $0x1,%rcx
movss  (%rax,%rcx,4),%xmm0
divss  -0x3c(%rbp),%xmm0
mov    -0x48(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x48(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    $0x1,%rcx
movss  (%rax,%rcx,4),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
divss  -0x3c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    117e <tridag+0x5e>
mov    -0x30(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x38(%rbp)
cmpq   $0x1,-0x38(%rbp)
jb     12bf <tridag+0x19f>
mov    -0x48(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  0x4(%rax,%rcx,4),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  0x4(%rax,%rcx,4),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movd   %xmm0,%edx
xor    $0x80000000,%edx
movd   %edx,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x38(%rbp)
jmp    125f <tridag+0x13f>
mov    -0x48(%rbp),%rdi
mov    -0x30(%rbp),%rdx
mov    $0x1,%esi
call   1030 <free_vector@plt>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<tridag>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %r9,%r14
mov    %r8,%rbx
mov    %rcx,%r15
mov    %rdx,%r12
mov    %rsi,%r13
mov    %rdi,%rbp
mov    $0x1,%edi
mov    %r9,%rsi
call   1040 <vector@plt>
movss  0x4(%r15),%xmm1
movss  0x4(%r13),%xmm0
divss  %xmm0,%xmm1
movss  %xmm1,0x4(%rbx)
cmp    $0x2,%r14
jb     123a <tridag+0x11a>
lea    0x8(%rax),%rcx
lea    (%rbx,%r14,4),%rdx
add    $0x4,%rdx
cmp    %rdx,%rcx
jae    11da <tridag+0xba>
lea    (%rax,%r14,4),%rsi
add    $0x4,%rsi
lea    0x4(%rbx),%rdx
cmp    %rsi,%rdx
jae    11da <tridag+0xba>
mov    $0x1,%esi
movss  (%r12,%rsi,4),%xmm1
divss  %xmm0,%xmm1
movss  %xmm1,-0x4(%rcx,%rsi,4)
movss  0x4(%r13,%rsi,4),%xmm0
movss  0x4(%rbp,%rsi,4),%xmm2
mulss  %xmm2,%xmm1
subss  %xmm1,%xmm0
movss  0x4(%r15,%rsi,4),%xmm1
mulss  -0x4(%rdx,%rsi,4),%xmm2
subss  %xmm2,%xmm1
divss  %xmm0,%xmm1
movss  %xmm1,(%rdx,%rsi,4)
add    $0x1,%rsi
cmp    %rsi,%r14
jne    1190 <tridag+0x70>
jmp    123a <tridag+0x11a>
movss  0x4(%rbx),%xmm1
mov    $0x1,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%r12,%rdx,4),%xmm2
divss  %xmm0,%xmm2
movss  %xmm2,-0x4(%rcx,%rdx,4)
movss  0x4(%r13,%rdx,4),%xmm0
movss  0x4(%rbp,%rdx,4),%xmm3
mulss  %xmm3,%xmm2
subss  %xmm2,%xmm0
movss  0x4(%r15,%rdx,4),%xmm2
mulss  %xmm3,%xmm1
subss  %xmm1,%xmm2
divss  %xmm0,%xmm2
movss  %xmm2,0x4(%rbx,%rdx,4)
add    $0x1,%rdx
movaps %xmm2,%xmm1
cmp    %rdx,%r14
jne    11f0 <tridag+0xd0>
mov    %r14,%rcx
add    $0xffffffffffffffff,%rcx
je     12c9 <tridag+0x1a9>
movss  (%rbx,%r14,4),%xmm0
test   $0x1,%cl
jne    1260 <tridag+0x140>
mov    %rcx,%rdx
mov    %r14,%rcx
cmp    $0x2,%r14
jne    1285 <tridag+0x165>
jmp    12c9 <tridag+0x1a9>
mulss  (%rax,%r14,4),%xmm0
movss  -0x4(%rbx,%r14,4),%xmm1
subss  %xmm0,%xmm1
movss  %xmm1,-0x4(%rbx,%r14,4)
lea    -0x2(%r14),%rdx
movaps %xmm1,%xmm0
cmp    $0x2,%r14
je     12c9 <tridag+0x1a9>
add    $0xffffffffffffffff,%rdx
nopl   0x0(%rax)
movss  (%rbx,%rdx,4),%xmm1
movss  0x4(%rbx,%rdx,4),%xmm2
mulss  (%rax,%rcx,4),%xmm0
subss  %xmm0,%xmm2
movss  %xmm2,0x4(%rbx,%rdx,4)
mulss  0x4(%rax,%rdx,4),%xmm2
subss  %xmm2,%xmm1
movss  %xmm1,(%rbx,%rdx,4)
mov    %rdx,%rcx
add    $0xfffffffffffffffe,%rdx
movaps %xmm1,%xmm0
cmp    $0xffffffffffffffff,%rdx
jne    1290 <tridag+0x170>
mov    %rax,%rdi
mov    $0x1,%esi
mov    %r14,%rdx
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free_vector@plt>

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

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