#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 maximo(double *, int);
double normaEuclideana(double *, int);
int /*<<< orphan*/ normaliza(double *, int, double);
double calculaMu(double x[], int n) {
  int i;
  double max;
  double gama, sigma;
  max = maximo(x, n);
  if (max == 0)
    gama = 0;
  else {
    normaliza(x, n, max);
    sigma = normaEuclideana(x, n);
    if (x[0] < 0)
      sigma *= (-1);
    x[0] += sigma;
    gama = 1 / (sigma * x[0]);
    for (i = 1; i < n; i++)
      x[i] /= gama;
    x[0] = sigma;
    gama *= max;
  }
  return gama;
}


x86 assembly (clang, -O0):
<calculaMu>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
call   1040 <maximo@plt>
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1175 <calculaMu+0x45>
jp     1175 <calculaMu+0x45>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    1241 <calculaMu+0x111>
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
movsd  -0x18(%rbp),%xmm0
call   1030 <normaliza@plt>
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
call   1050 <normaEuclideana@plt>
movsd  %xmm0,-0x28(%rbp)
mov    -0x8(%rbp),%rax
xorps  %xmm0,%xmm0
ucomisd (%rax),%xmm0
jbe    11ba <calculaMu+0x8a>
movsd  0xe50(%rip),%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
movsd  -0x28(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mulsd  (%rax),%xmm1
movsd  0xe28(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1225 <calculaMu+0xf5>
movsd  -0x20(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    11f0 <calculaMu+0xc0>
movsd  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x18(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<calculaMu>:
push   %rbp
push   %rbx
push   %rax
mov    %esi,%ebp
mov    %rdi,%rbx
call   1040 <maximo@plt>
movapd %xmm0,%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
jne    1151 <calculaMu+0x21>
jnp    1278 <calculaMu+0x148>
mov    %rbx,%rdi
mov    %ebp,%esi
movsd  %xmm1,(%rsp)
movapd %xmm1,%xmm0
call   1030 <normaliza@plt>
mov    %rbx,%rdi
mov    %ebp,%esi
call   1050 <normaEuclideana@plt>
movsd  (%rbx),%xmm2
movapd 0xe86(%rip),%xmm3
xorpd  %xmm0,%xmm3
xorpd  %xmm4,%xmm4
movapd %xmm2,%xmm1
cmpltsd %xmm4,%xmm1
andpd  %xmm1,%xmm3
andnpd %xmm0,%xmm1
orpd   %xmm3,%xmm1
addsd  %xmm1,%xmm2
mulsd  %xmm1,%xmm2
movsd  0xe69(%rip),%xmm0
divsd  %xmm2,%xmm0
cmp    $0x2,%ebp
jl     1267 <calculaMu+0x137>
mov    %ebp,%eax
lea    -0x1(%rax),%rcx
mov    $0x1,%edx
cmp    $0x2,%rcx
jb     1250 <calculaMu+0x120>
mov    %rcx,%rdx
and    $0xfffffffffffffffe,%rdx
movapd %xmm0,%xmm2
unpcklpd %xmm0,%xmm2
lea    -0x2(%rdx),%rsi
mov    %rsi,%rbp
shr    %rbp
add    $0x1,%rbp
test   %rsi,%rsi
je     127f <calculaMu+0x14f>
mov    %rbp,%rsi
and    $0xfffffffffffffffe,%rsi
mov    $0x1,%edi
nopl   0x0(%rax,%rax,1)
movupd (%rbx,%rdi,8),%xmm3
movupd 0x10(%rbx,%rdi,8),%xmm4
divpd  %xmm2,%xmm3
movupd %xmm3,(%rbx,%rdi,8)
divpd  %xmm2,%xmm4
movupd %xmm4,0x10(%rbx,%rdi,8)
add    $0x4,%rdi
add    $0xfffffffffffffffe,%rsi
jne    1200 <calculaMu+0xd0>
test   $0x1,%bpl
je     123c <calculaMu+0x10c>
movupd (%rbx,%rdi,8),%xmm3
divpd  %xmm2,%xmm3
movupd %xmm3,(%rbx,%rdi,8)
cmp    %rdx,%rcx
je     1267 <calculaMu+0x137>
or     $0x1,%rcx
mov    %rcx,%rdx
nopl   0x0(%rax,%rax,1)
movsd  (%rbx,%rdx,8),%xmm2
divsd  %xmm0,%xmm2
movsd  %xmm2,(%rbx,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1250 <calculaMu+0x120>
movsd  %xmm1,(%rbx)
movsd  (%rsp),%xmm1
mulsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
add    $0x8,%rsp
pop    %rbx
pop    %rbp
ret
mov    $0x1,%edi
test   $0x1,%bpl
jne    122e <calculaMu+0xfe>
jmp    123c <calculaMu+0x10c>

LLVM IR (-O0):
define double @calculaMu(double* %0, i32 %1) {
%3 = alloca double*
%4 = alloca i32
%5 = alloca i32
%6 = alloca double
%7 = alloca double
%8 = alloca double
store double* %0, double** %3
store i32 %1, i32* %4
%9 = load double*, double** %3
%10 = load i32, i32* %4
%11 = call double @maximo(double* %9, i32 %10)
store double %11, double* %6
%12 = load double, double* %6
%13 = fcmp oeq double %12, 0.000000e+00
br i1 %13, label %14, label %15
14: ; preds = %2
store double 0.000000e+00, double* %7
br label %64
15: ; preds = %2
%16 = load double*, double** %3
%17 = load i32, i32* %4
%18 = load double, double* %6
%19 = call i32 @normaliza(double* %16, i32 %17, double %18)
%20 = load double*, double** %3
%21 = load i32, i32* %4
%22 = call double @normaEuclideana(double* %20, i32 %21)
store double %22, double* %8
%23 = load double*, double** %3
%24 = getelementptr inbounds double, double* %23, i64 0
%25 = load double, double* %24
%26 = fcmp olt double %25, 0.000000e+00
br i1 %26, label %27, label %30
27: ; preds = %15
%28 = load double, double* %8
%29 = fmul double %28, -1.000000e+00
store double %29, double* %8
br label %30
30: ; preds = %27, %15
%31 = load double, double* %8
%32 = load double*, double** %3
%33 = getelementptr inbounds double, double* %32, i64 0
%34 = load double, double* %33
%35 = fadd double %34, %31
store double %35, double* %33
%36 = load double, double* %8
%37 = load double*, double** %3
%38 = getelementptr inbounds double, double* %37, i64 0
%39 = load double, double* %38
%40 = fmul double %36, %39
%41 = fdiv double 1.000000e+00, %40
store double %41, double* %7
store i32 1, i32* %5
br label %42
42: ; preds = %54, %30
%43 = load i32, i32* %5
%44 = load i32, i32* %4
%45 = icmp slt i32 %43, %44
br i1 %45, label %46, label %57
46: ; preds = %42
%47 = load double, double* %7
%48 = load double*, double** %3
%49 = load i32, i32* %5
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds double, double* %48, i64 %50
%52 = load double, double* %51
%53 = fdiv double %52, %47
store double %53, double* %51
br label %54
54: ; preds = %46
%55 = load i32, i32* %5
%56 = add nsw i32 %55, 1
store i32 %56, i32* %5
br label %42
57: ; preds = %42
%58 = load double, double* %8
%59 = load double*, double** %3
%60 = getelementptr inbounds double, double* %59, i64 0
store double %58, double* %60
%61 = load double, double* %6
%62 = load double, double* %7
%63 = fmul double %62, %61
store double %63, double* %7
br label %64
64: ; preds = %57, %14
%65 = load double, double* %7
ret double %65
}

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