#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 */
double *GaussProjCal(double, double, double);
double *LMLE(double *, double *, int);
int /*<<< orphan*/ free(double *);
double *gaosifansuan(double, double, double);
double m(int);
scalar_t__ malloc(int);
double *zuobiao(double *L, double *B, int *RSSI, int n, double L0) {
  double xy[10][2];
  double *length = (double *)malloc(n * sizeof(double));
  double *XY = (double *)malloc(2 * sizeof(double));
  double X, Y;
  double *BL = (double *)malloc(2 * sizeof(double));
  double *xyl = (double *)malloc((2 * n) * sizeof(double));
  for (int i = 0; i < n; i++) {
    double *a = GaussProjCal(L[i], B[i], L0);
    xy[i][0] = a[0];
    xy[i][1] = a[1];
    free(a);
    a = NULL;
    length[i] = m(RSSI[i]);
  }
  for (int i = 0, j = 0; i < 2 * n, j < n; i++, j++) {
    xyl[i] = xy[j][0];
    xyl[i + 1] = xy[j][1];
    i = i + 1;
  }
  XY = LMLE(xyl, length, n);
  X = XY[0];
  Y = XY[1];
  BL = gaosifansuan(X, Y, L0);
  free(length);
  length = NULL;
  free(XY);
  XY = NULL;
  free(xyl);
  xyl = NULL;
  return BL;
}


x86 assembly (clang, -O0):
<zuobiao>:
push   %rbp
mov    %rsp,%rbp
sub    $0x120,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movsd  %xmm0,-0x28(%rbp)
movslq -0x1c(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0xd8(%rbp)
mov    $0x10,%edi
call   1060 <malloc@plt>
mov    %rax,-0xe0(%rbp)
mov    $0x10,%edi
call   1060 <malloc@plt>
mov    %rax,-0xf8(%rbp)
mov    -0x1c(%rbp),%eax
shl    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0x100(%rbp)
movl   $0x0,-0x104(%rbp)
mov    -0x104(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12b6 <zuobiao+0x156>
mov    -0x8(%rbp),%rax
movslq -0x104(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x104(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x28(%rbp),%xmm2
call   1050 <GaussProjCal@plt>
mov    %rax,-0x110(%rbp)
mov    -0x110(%rbp),%rax
movsd  (%rax),%xmm0
movslq -0x104(%rbp),%rcx
lea    -0xd0(%rbp),%rax
shl    $0x4,%rcx
add    %rcx,%rax
movsd  %xmm0,(%rax)
mov    -0x110(%rbp),%rax
movsd  0x8(%rax),%xmm0
movslq -0x104(%rbp),%rcx
lea    -0xd0(%rbp),%rax
shl    $0x4,%rcx
add    %rcx,%rax
movsd  %xmm0,0x8(%rax)
mov    -0x110(%rbp),%rdi
call   1030 <free@plt>
movq   $0x0,-0x110(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x104(%rbp),%rcx
mov    (%rax,%rcx,4),%edi
call   1070 <m@plt>
mov    -0xd8(%rbp),%rax
movslq -0x104(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x104(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x104(%rbp)
jmp    11db <zuobiao+0x7b>
movl   $0x0,-0x114(%rbp)
movl   $0x0,-0x118(%rbp)
mov    -0x118(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1369 <zuobiao+0x209>
movslq -0x118(%rbp),%rcx
lea    -0xd0(%rbp),%rax
shl    $0x4,%rcx
add    %rcx,%rax
movsd  (%rax),%xmm0
mov    -0x100(%rbp),%rax
movslq -0x114(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movslq -0x118(%rbp),%rcx
lea    -0xd0(%rbp),%rax
shl    $0x4,%rcx
add    %rcx,%rax
movsd  0x8(%rax),%xmm0
mov    -0x100(%rbp),%rax
mov    -0x114(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x114(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x114(%rbp)
mov    -0x114(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x114(%rbp)
mov    -0x118(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x118(%rbp)
jmp    12ca <zuobiao+0x16a>
mov    -0x100(%rbp),%rdi
mov    -0xd8(%rbp),%rsi
mov    -0x1c(%rbp),%edx
call   1080 <LMLE@plt>
mov    %rax,-0xe0(%rbp)
mov    -0xe0(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
mov    -0xe0(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  -0xe8(%rbp),%xmm0
ff
movsd  -0xf0(%rbp),%xmm1
ff
movsd  -0x28(%rbp),%xmm2
call   1040 <gaosifansuan@plt>
mov    %rax,-0xf8(%rbp)
mov    -0xd8(%rbp),%rdi
call   1030 <free@plt>
movq   $0x0,-0xd8(%rbp)
mov    -0xe0(%rbp),%rdi
call   1030 <free@plt>
movq   $0x0,-0xe0(%rbp)
mov    -0x100(%rbp),%rdi
call   1030 <free@plt>
movq   $0x0,-0x100(%rbp)
mov    -0xf8(%rbp),%rax
add    $0x120,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<zuobiao>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xc8,%rsp
movsd  %xmm0,(%rsp)
mov    %ecx,%r12d
mov    %rdx,%r14
mov    %rsi,%rbx
mov    %rdi,%rbp
lea    0x0(,%r12,8),%edi
call   1060 <malloc@plt>
mov    %rax,%r15
mov    $0x10,%edi
call   1060 <malloc@plt>
mov    $0x10,%edi
call   1060 <malloc@plt>
mov    %r12d,%edi
shl    $0x4,%edi
call   1060 <malloc@plt>
mov    %rax,%r13
test   %r12d,%r12d
jle    1370 <zuobiao+0x210>
mov    %r13,0x10(%rsp)
mov    %r12,0x18(%rsp)
mov    %r12d,%eax
mov    %rax,0x8(%rsp)
lea    0x0(,%rax,4),%r13
xor    %r12d,%r12d
nopw   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%r12,2),%xmm0
movsd  (%rbx,%r12,2),%xmm1
movsd  (%rsp),%xmm2
call   1050 <GaussProjCal@plt>
movups (%rax),%xmm0
movaps %xmm0,0x20(%rsp,%r12,4)
mov    %rax,%rdi
call   1030 <free@plt>
mov    (%r14,%r12,1),%edi
call   1070 <m@plt>
movsd  %xmm0,(%r15,%r12,2)
add    $0x4,%r12
cmp    %r12,%r13
jne    11e0 <zuobiao+0x80>
mov    0x18(%rsp),%r12
test   %r12d,%r12d
mov    0x10(%rsp),%r13
mov    0x8(%rsp),%rbp
jle    1370 <zuobiao+0x210>
cmp    $0x1,%r12d
jne    1247 <zuobiao+0xe7>
xor    %ecx,%ecx
xor    %eax,%eax
jmp    1324 <zuobiao+0x1c4>
mov    %ebp,%eax
and    $0xfffffffe,%eax
lea    -0x2(%rax),%rdx
mov    %rdx,%rsi
shr    %rsi
add    $0x1,%rsi
mov    %esi,%ecx
and    $0x3,%ecx
cmp    $0x6,%rdx
jae    1269 <zuobiao+0x109>
xor    %edx,%edx
jmp    12e5 <zuobiao+0x185>
and    $0xfffffffffffffffc,%rsi
mov    $0x30,%edi
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movaps -0x40(%rsp,%rdi,2),%xmm0
movaps -0x30(%rsp,%rdi,2),%xmm1
movups %xmm1,-0x50(%r13,%rdi,2)
movups %xmm0,-0x60(%r13,%rdi,2)
movaps -0x20(%rsp,%rdi,2),%xmm0
movaps -0x10(%rsp,%rdi,2),%xmm1
movups %xmm1,-0x30(%r13,%rdi,2)
movups %xmm0,-0x40(%r13,%rdi,2)
movaps (%rsp,%rdi,2),%xmm0
movaps 0x10(%rsp,%rdi,2),%xmm1
movups %xmm1,-0x10(%r13,%rdi,2)
movups %xmm0,-0x20(%r13,%rdi,2)
movaps 0x20(%rsp,%rdi,2),%xmm0
movaps 0x30(%rsp,%rdi,2),%xmm1
movups %xmm1,0x10(%r13,%rdi,2)
movups %xmm0,0x0(%r13,%rdi,2)
add    $0x8,%rdx
add    $0x40,%rdi
add    $0xfffffffffffffffc,%rsi
jne    1280 <zuobiao+0x120>
test   %rcx,%rcx
je     131b <zuobiao+0x1bb>
shl    $0x4,%rdx
lea    (%rsp,%rdx,1),%rsi
add    $0x20,%rsi
add    %r13,%rdx
shl    $0x5,%rcx
xor    %edi,%edi
nop
movaps (%rsi,%rdi,1),%xmm0
movaps 0x10(%rsi,%rdi,1),%xmm1
movups %xmm1,0x10(%rdx,%rdi,1)
movups %xmm0,(%rdx,%rdi,1)
add    $0x20,%rdi
cmp    %rdi,%rcx
jne    1300 <zuobiao+0x1a0>
cmp    %rbp,%rax
je     1370 <zuobiao+0x210>
lea    (%rax,%rax,1),%rcx
lea    0x0(,%rcx,8),%rcx
add    %r13,%rcx
sub    %rax,%rbp
shl    $0x4,%rax
add    %rsp,%rax
add    $0x20,%rax
mov    $0x8,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  -0x8(%rax,%rdx,1),%xmm0
movsd  %xmm0,-0x8(%rcx,%rdx,1)
movsd  (%rax,%rdx,1),%xmm0
movsd  %xmm0,(%rcx,%rdx,1)
add    $0x10,%rdx
add    $0xffffffffffffffff,%rbp
jne    1350 <zuobiao+0x1f0>
mov    %r13,%rdi
mov    %r15,%rsi
mov    %r12d,%edx
call   1080 <LMLE@plt>
mov    %rax,%rbx
movsd  (%rax),%xmm0
movsd  0x8(%rax),%xmm1
movsd  (%rsp),%xmm2
call   1040 <gaosifansuan@plt>
mov    %rax,%rbp
mov    %r15,%rdi
call   1030 <free@plt>
mov    %rbx,%rdi
call   1030 <free@plt>
mov    %r13,%rdi
call   1030 <free@plt>
mov    %rbp,%rax
add    $0xc8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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