void simp3(float **a, int i1, int k1, int ip, int kp) {
  int kk, ii;
  float piv;
  piv = 1.0 / a[ip + 1][kp + 1];
  for (ii = 1; ii <= i1 + 1; ii++)
    if (ii - 1 != ip) {
      a[ii][kp + 1] *= piv;
      for (kk = 1; kk <= k1 + 1; kk++)
        if (kk - 1 != kp)
          a[ii][kk] -= a[ip + 1][kk] * a[ii][kp + 1];
    }
  for (kk = 1; kk <= k1 + 1; kk++)
    if (kk - 1 != kp)
      a[ip + 1][kk] *= -piv;
  a[ip + 1][kp + 1] = piv;
}


x86 assembly (clang, -O0):
<simp3>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x18(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xec0(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     124e <simp3+0x14e>
mov    -0x20(%rbp),%eax
sub    $0x1,%eax
cmp    -0x14(%rbp),%eax
je     123b <simp3+0x13b>
movss  -0x24(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x18(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     1236 <simp3+0x136>
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
cmp    -0x18(%rbp),%eax
je     1223 <simp3+0x123>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x18(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x1c(%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)
jmp    1228 <simp3+0x128>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    119f <simp3+0x9f>
jmp    123b <simp3+0x13b>
jmp    1240 <simp3+0x140>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1154 <simp3+0x54>
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     12b9 <simp3+0x1b9>
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
cmp    -0x18(%rbp),%eax
je     12a6 <simp3+0x1a6>
movss  -0x24(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x1c(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    12ab <simp3+0x1ab>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1255 <simp3+0x155>
movss  -0x24(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x18(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
pop    %rbp
ret

x86 assembly (clang, -O2):
<simp3>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
movslq %ecx,%rax
mov    0x8(%rdi,%rax,8),%r14
movslq %r8d,%rax
lea    0x1(%rax),%r8
movss  0xee0(%rip),%xmm0
divss  0x4(%r14,%rax,4),%xmm0
mov    %eax,%r9d
test   %esi,%esi
js     121a <simp3+0x11a>
xor    %r10d,%r10d
test   %edx,%edx
cmovg  %edx,%r10d
mov    %ecx,%r11d
add    $0x2,%esi
add    $0x1,%r10d
mov    %r9,%r15
neg    %r15
mov    %r10d,%r12d
and    $0xfffffffe,%r12d
neg    %r12
mov    $0x1,%r13d
jmp    116d <simp3+0x6d>
nopl   (%rax)
add    $0x1,%r13
cmp    %rsi,%r13
je     121a <simp3+0x11a>
lea    -0x1(%r13),%rax
cmp    %r11,%rax
je     1160 <simp3+0x60>
mov    (%rdi,%r13,8),%rcx
movss  (%rcx,%r8,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,(%rcx,%r8,4)
test   %edx,%edx
js     1160 <simp3+0x60>
mov    $0x1,%ebx
jne    11d2 <simp3+0xd2>
test   $0x1,%r10b
je     1160 <simp3+0x60>
lea    -0x1(%rbx),%rax
cmp    %r9,%rax
je     1160 <simp3+0x60>
movss  (%r14,%rbx,4),%xmm1
movss  (%rcx,%rbx,4),%xmm2
mulss  (%rcx,%r8,4),%xmm1
subss  %xmm1,%xmm2
movss  %xmm2,(%rcx,%rbx,4)
jmp    1160 <simp3+0x60>
lea    (%r12,%rbx,1),%rax
add    $0x2,%rax
add    $0x2,%rbx
cmp    $0x1,%rax
je     1195 <simp3+0x95>
lea    (%r15,%rbx,1),%rax
cmp    $0x1,%rax
je     11f6 <simp3+0xf6>
movss  (%r14,%rbx,4),%xmm1
movss  (%rcx,%rbx,4),%xmm2
mulss  (%rcx,%r8,4),%xmm1
subss  %xmm1,%xmm2
movss  %xmm2,(%rcx,%rbx,4)
test   %rax,%rax
je     11c0 <simp3+0xc0>
movss  0x4(%r14,%rbx,4),%xmm1
movss  0x4(%rcx,%rbx,4),%xmm2
mulss  (%rcx,%r8,4),%xmm1
subss  %xmm1,%xmm2
movss  %xmm2,0x4(%rcx,%rbx,4)
jmp    11c0 <simp3+0xc0>
test   %edx,%edx
js     12b3 <simp3+0x1b3>
movaps 0xde7(%rip),%xmm1
xorps  %xmm0,%xmm1
lea    0x1(%rdx),%eax
test   %edx,%edx
je     1295 <simp3+0x195>
mov    %r9,%rdx
neg    %rdx
mov    %eax,%esi
and    $0xfffffffe,%esi
neg    %rsi
mov    $0x1,%ecx
jmp    1262 <simp3+0x162>
nopl   0x0(%rax,%rax,1)
lea    (%rsi,%rcx,1),%rdi
add    $0x2,%rdi
add    $0x2,%rcx
cmp    $0x1,%rdi
je     129a <simp3+0x19a>
lea    (%rdx,%rcx,1),%rdi
cmp    $0x1,%rdi
je     127c <simp3+0x17c>
movss  (%r14,%rcx,4),%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,(%r14,%rcx,4)
test   %rdi,%rdi
je     1250 <simp3+0x150>
movss  0x4(%r14,%rcx,4),%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,0x4(%r14,%rcx,4)
jmp    1250 <simp3+0x150>
mov    $0x1,%ecx
test   $0x1,%al
je     12b3 <simp3+0x1b3>
lea    -0x1(%rcx),%rax
cmp    %r9,%rax
je     12b3 <simp3+0x1b3>
mulss  (%r14,%rcx,4),%xmm1
movss  %xmm1,(%r14,%rcx,4)
movss  %xmm0,(%r14,%r8,4)
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

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

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