#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*/ copyCol(double **, double *, int, int);
double sqrt(double);
double calculaUSigmaGama(double A[1000][1000], double norma[1000],
                         double u[1000], double max, int index, int n, int m) {
  int i;
  double gama, sigma, aux;
  sigma = sqrt(norma[index]);
  if (sigma < 1E-5)
    return 1E-5;
  copyCol(A, u, index, n);
  if (u[index] < 0)
    sigma *= (-1);
  u[index] += sigma;
  gama = (u[index] / sigma);
  aux = u[index];
  for (i = index; i < n; i++)
    u[i] /= aux;
  if (sigma > 0)
    sigma *= -1;
  return gama;
}


x86 assembly (clang, -O0):
<calculaUSigmaGama>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
mov    %ecx,-0x2c(%rbp)
mov    %r8d,-0x30(%rbp)
mov    %r9d,-0x34(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x48(%rbp)
movsd  0xe9d(%rip),%xmm0
ucomisd -0x48(%rbp),%xmm0
jbe    1180 <calculaUSigmaGama+0x60>
movsd  0xe8a(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1264 <calculaUSigmaGama+0x144>
mov    -0x10(%rbp),%rdi
mov    -0x20(%rbp),%rsi
mov    -0x2c(%rbp),%edx
mov    -0x30(%rbp),%ecx
call   1040 <copyCol@plt>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
xorps  %xmm0,%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jbe    11bb <calculaUSigmaGama+0x9b>
movsd  0xe57(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x30(%rbp),%eax
jge    1236 <calculaUSigmaGama+0x116>
movsd  -0x50(%rbp),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1201 <calculaUSigmaGama+0xe1>
movsd  -0x48(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    125a <calculaUSigmaGama+0x13a>
movsd  0xdb8(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<calculaUSigmaGama>:
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x18,%rsp
mov    %r8d,%r15d
mov    %rdx,%rbx
mov    %rdi,%rbp
movslq %ecx,%r14
movsd  (%rsi,%r14,8),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     114f <calculaUSigmaGama+0x2f>
xorps  %xmm1,%xmm1
sqrtsd %xmm0,%xmm1
jmp    1158 <calculaUSigmaGama+0x38>
call   1030 <sqrt@plt>
movapd %xmm0,%xmm1
movsd  0xea0(%rip),%xmm0
ucomisd %xmm1,%xmm0
ja     12a7 <calculaUSigmaGama+0x187>
mov    %rbp,%rdi
mov    %rbx,%rsi
mov    %r14d,%edx
mov    %r15d,%ecx
movapd %xmm1,(%rsp)
call   1040 <copyCol@plt>
movsd  (%rbx,%r14,8),%xmm1
movapd 0xe82(%rip),%xmm0
movapd (%rsp),%xmm4
xorpd  %xmm4,%xmm0
xorpd  %xmm2,%xmm2
movapd %xmm1,%xmm3
cmpltsd %xmm2,%xmm3
andpd  %xmm3,%xmm0
andnpd %xmm4,%xmm3
orpd   %xmm0,%xmm3
addsd  %xmm3,%xmm1
movsd  %xmm1,(%rbx,%r14,8)
movapd %xmm1,%xmm0
divsd  %xmm3,%xmm0
cmp    %r15d,%r14d
jge    12a7 <calculaUSigmaGama+0x187>
movslq %r15d,%rax
movapd %xmm1,%xmm2
divsd  %xmm1,%xmm2
movsd  %xmm2,(%rbx,%r14,8)
lea    0x1(%r14),%rcx
cmp    %rax,%rcx
je     12a7 <calculaUSigmaGama+0x187>
mov    %r14,%rdx
not    %rdx
add    %rax,%rdx
cmp    $0x2,%rdx
jb     1290 <calculaUSigmaGama+0x170>
mov    %rdx,%r8
and    $0xfffffffffffffffe,%r8
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
lea    -0x2(%r8),%rsi
mov    %rsi,%r9
shr    %r9
add    $0x1,%r9
test   %rsi,%rsi
je     12b2 <calculaUSigmaGama+0x192>
mov    %r9,%rdi
and    $0xfffffffffffffffe,%rdi
lea    (%rbx,%r14,8),%rsi
add    $0x18,%rsi
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd -0x10(%rsi,%rbp,8),%xmm3
movupd (%rsi,%rbp,8),%xmm4
divpd  %xmm2,%xmm3
movupd %xmm3,-0x10(%rsi,%rbp,8)
divpd  %xmm2,%xmm4
movupd %xmm4,(%rsi,%rbp,8)
add    $0x4,%rbp
add    $0xfffffffffffffffe,%rdi
jne    1240 <calculaUSigmaGama+0x120>
test   $0x1,%r9b
je     127f <calculaUSigmaGama+0x15f>
add    %rcx,%rbp
movupd (%rbx,%rbp,8),%xmm3
divpd  %xmm2,%xmm3
movupd %xmm3,(%rbx,%rbp,8)
cmp    %r8,%rdx
je     12a7 <calculaUSigmaGama+0x187>
add    %r8,%rcx
nopw   0x0(%rax,%rax,1)
movsd  (%rbx,%rcx,8),%xmm2
divsd  %xmm1,%xmm2
movsd  %xmm2,(%rbx,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1290 <calculaUSigmaGama+0x170>
add    $0x18,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %ebp,%ebp
test   $0x1,%r9b
jne    126e <calculaUSigmaGama+0x14e>
jmp    127f <calculaUSigmaGama+0x15f>

LLVM IR (-O0):
define double @calculaUSigmaGama([1000 x double]* %0, double* %1, double* %2, double %3, i32 %4, i32 %5, i32 %6) {
%8 = alloca double
%9 = alloca [1000 x double]*
%10 = alloca double*
%11 = alloca double*
%12 = alloca double
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca double
%18 = alloca double
%19 = alloca double
store [1000 x double]* %0, [1000 x double]** %9
store double* %1, double** %10
store double* %2, double** %11
store double %3, double* %12
store i32 %4, i32* %13
store i32 %5, i32* %14
store i32 %6, i32* %15
%20 = load double*, double** %10
%21 = load i32, i32* %13
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds double, double* %20, i64 %22
%24 = load double, double* %23
%25 = call double @sqrt(double %24)
store double %25, double* %18
%26 = load double, double* %18
%27 = fcmp olt double %26, 1.000000e-05
br i1 %27, label %28, label %29
28: ; preds = %7
store double 1.000000e-05, double* %8
br label %89
29: ; preds = %7
%30 = load [1000 x double]*, [1000 x double]** %9
%31 = bitcast [1000 x double]* %30 to double**
%32 = load double*, double** %11
%33 = load i32, i32* %13
%34 = load i32, i32* %14
%35 = call i32 @copyCol(double** %31, double* %32, i32 %33, i32 %34)
%36 = load double*, double** %11
%37 = load i32, i32* %13
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds double, double* %36, i64 %38
%40 = load double, double* %39
%41 = fcmp olt double %40, 0.000000e+00
br i1 %41, label %42, label %45
42: ; preds = %29
%43 = load double, double* %18
%44 = fmul double %43, -1.000000e+00
store double %44, double* %18
br label %45
45: ; preds = %42, %29
%46 = load double, double* %18
%47 = load double*, double** %11
%48 = load i32, i32* %13
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %47, i64 %49
%51 = load double, double* %50
%52 = fadd double %51, %46
store double %52, double* %50
%53 = load double*, double** %11
%54 = load i32, i32* %13
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %53, i64 %55
%57 = load double, double* %56
%58 = load double, double* %18
%59 = fdiv double %57, %58
store double %59, double* %17
%60 = load double*, double** %11
%61 = load i32, i32* %13
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double, double* %60, i64 %62
%64 = load double, double* %63
store double %64, double* %19
%65 = load i32, i32* %13
store i32 %65, i32* %16
br label %66
66: ; preds = %78, %45
%67 = load i32, i32* %16
%68 = load i32, i32* %14
%69 = icmp slt i32 %67, %68
br i1 %69, label %70, label %81
70: ; preds = %66
%71 = load double, double* %19
%72 = load double*, double** %11
%73 = load i32, i32* %16
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds double, double* %72, i64 %74
%76 = load double, double* %75
%77 = fdiv double %76, %71
store double %77, double* %75
br label %78
78: ; preds = %70
%79 = load i32, i32* %16
%80 = add nsw i32 %79, 1
store i32 %80, i32* %16
br label %66
81: ; preds = %66
%82 = load double, double* %18
%83 = fcmp ogt double %82, 0.000000e+00
br i1 %83, label %84, label %87
84: ; preds = %81
%85 = load double, double* %18
%86 = fmul double %85, -1.000000e+00
store double %86, double* %18
br label %87
87: ; preds = %84, %81
%88 = load double, double* %17
store double %88, double* %8
br label %89
89: ; preds = %87, %28
%90 = load double, double* %8
ret double %90
}

LLVM IR (-O2):
define double @calculaUSigmaGama([1000 x double]* %0, double* nocapture %1, double* %2, double %3, i32 %4, i32 %5, i32 %6) {
%8 = sext i32 %4 to i64
%9 = getelementptr inbounds double, double* %1, i64 %8
%10 = load double, double* %9
%11 = tail call double @sqrt(double %10)
%12 = fcmp olt double %11, 1.000000e-05
br i1 %12, label %87, label %13
13: ; preds = %7
%14 = bitcast [1000 x double]* %0 to double**
%15 = tail call i32 @copyCol(double** %14, double* %2, i32 %4, i32 %5)
%16 = getelementptr inbounds double, double* %2, i64 %8
%17 = load double, double* %16
%18 = fcmp olt double %17, 0.000000e+00
%19 = fneg double %11
%20 = select i1 %18, double %19, double %11
%21 = fadd double %17, %20
store double %21, double* %16
%22 = fdiv double %21, %20
%23 = icmp slt i32 %4, %5
br i1 %23, label %24, label %87
24: ; preds = %13
%25 = sext i32 %5 to i64
%26 = getelementptr inbounds double, double* %2, i64 %8
%27 = fdiv double %21, %21
store double %27, double* %26
%28 = add nsw i64 %8, 1
%29 = icmp eq i64 %28, %25
br i1 %29, label %87, label %30
30: ; preds = %24
%31 = xor i64 %8, -1
%32 = add nsw i64 %31, %25
%33 = icmp ult i64 %32, 2
br i1 %33, label %77, label %34
34: ; preds = %30
%35 = and i64 %32, -2
%36 = add nsw i64 %28, %35
%37 = insertelement <2 x double> poison, double %21, i64 0
%38 = shufflevector <2 x double> %37, <2 x double> poison, <2 x i32> zeroinitializer
%39 = add nsw i64 %35, -2
%40 = lshr exact i64 %39, 1
%41 = add nuw i64 %40, 1
%42 = and i64 %41, 1
%43 = icmp eq i64 %39, 0
br i1 %43, label %65, label %44
44: ; preds = %34
%45 = and i64 %41, -2
br label %46
46: ; preds = %46, %44
%47 = phi i64 [ 0, %44 ], [ %62, %46 ]
%48 = phi i64 [ 0, %44 ], [ %63, %46 ]
%49 = add i64 %28, %47
%50 = getelementptr inbounds double, double* %2, i64 %49
%51 = bitcast double* %50 to <2 x double>*
%52 = load <2 x double>, <2 x double>* %51
%53 = fdiv <2 x double> %52, %38
%54 = bitcast double* %50 to <2 x double>*
store <2 x double> %53, <2 x double>* %54
%55 = or i64 %47, 2
%56 = add i64 %28, %55
%57 = getelementptr inbounds double, double* %2, i64 %56
%58 = bitcast double* %57 to <2 x double>*
%59 = load <2 x double>, <2 x double>* %58
%60 = fdiv <2 x double> %59, %38
%61 = bitcast double* %57 to <2 x double>*
store <2 x double> %60, <2 x double>* %61
%62 = add nuw i64 %47, 4
%63 = add nuw i64 %48, 2
%64 = icmp eq i64 %63, %45
br i1 %64, label %65, label %46
65: ; preds = %46, %34
%66 = phi i64 [ 0, %34 ], [ %62, %46 ]
%67 = icmp eq i64 %42, 0
br i1 %67, label %75, label %68
68: ; preds = %65
%69 = add i64 %28, %66
%70 = getelementptr inbounds double, double* %2, i64 %69
%71 = bitcast double* %70 to <2 x double>*
%72 = load <2 x double>, <2 x double>* %71
%73 = fdiv <2 x double> %72, %38
%74 = bitcast double* %70 to <2 x double>*
store <2 x double> %73, <2 x double>* %74
br label %75
75: ; preds = %65, %68
%76 = icmp eq i64 %32, %35
br i1 %76, label %87, label %77
77: ; preds = %30, %75
%78 = phi i64 [ %28, %30 ], [ %36, %75 ]
br label %79
79: ; preds = %77, %79
%80 = phi i64 [ %85, %79 ], [ %78, %77 ]
%81 = getelementptr inbounds double, double* %2, i64 %80
%82 = load double, double* %81
%83 = getelementptr inbounds double, double* %2, i64 %80
%84 = fdiv double %82, %21
store double %84, double* %83
%85 = add nsw i64 %80, 1
%86 = icmp eq i64 %85, %25
br i1 %86, label %87, label %79
87: ; preds = %79, %24, %75, %13, %7
%88 = phi double [ 1.000000e-05, %7 ], [ %22, %13 ], [ %22, %75 ], [ %22, %24 ], [ %22, %79 ]
ret double %88
}