#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 Inner_Product(double *, double *, int);
int /*<<< orphan*/ Multiply_Sym_Matrix_by_Vector_ut(double *, double *,
                                                    double *, int);
void Update_H(double *H, double *dx, double *dg, double *Hg, int n) {
  double *Hdg;
  double dxtdg;
  double dgtHdg;
  double *pH;
  int i, j;
  dxtdg = Inner_Product(dx, dg, n);
  Multiply_Sym_Matrix_by_Vector_ut(Hg, H, dg, n);
  dgtHdg = Inner_Product(dg, Hg, n);
  dxtdg = 1.0 / dxtdg;
  dgtHdg = 1.0 / dgtHdg;
  pH = H;
  for (i = 0; i < n; i++)
    for (j = i; j < n; j++)
      *pH++ += dxtdg * dx[i] * dx[j] - dgtHdg * Hg[i] * Hg[j];
}


x86 assembly (clang, -O0):
<Update_H>:
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    %r8d,-0x24(%rbp)
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    -0x24(%rbp),%edx
call   1030 <Inner_Product@plt>
movsd  %xmm0,-0x38(%rbp)
mov    -0x20(%rbp),%rdi
mov    -0x8(%rbp),%rsi
mov    -0x18(%rbp),%rdx
mov    -0x24(%rbp),%ecx
call   1040 <Multiply_Sym_Matrix_by_Vector_ut@plt>
mov    -0x18(%rbp),%rdi
mov    -0x20(%rbp),%rsi
mov    -0x24(%rbp),%edx
call   1030 <Inner_Product@plt>
movsd  %xmm0,-0x40(%rbp)
movsd  0xe7e(%rip),%xmm0
divsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xe6c(%rip),%xmm0
divsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x48(%rbp)
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1260 <Update_H+0x140>
mov    -0x4c(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    124d <Update_H+0x12d>
movsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x4c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x50(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x48(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x48(%rbp)
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    11bf <Update_H+0x9f>
jmp    1252 <Update_H+0x132>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11ad <Update_H+0x8d>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<Update_H>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r8d,%r12d
mov    %rcx,%r14
mov    %rdx,%rbp
mov    %rsi,%rbx
mov    %rdi,%r15
mov    %rsi,%rdi
mov    %rdx,%rsi
mov    %r8d,%edx
call   1030 <Inner_Product@plt>
movsd  %xmm0,(%rsp)
mov    %r14,%rdi
mov    %r15,%rsi
mov    %rbp,%rdx
mov    %r12d,%ecx
call   1040 <Multiply_Sym_Matrix_by_Vector_ut@plt>
mov    %rbp,%rdi
mov    %r14,%rsi
mov    %r12d,%edx
call   1030 <Inner_Product@plt>
test   %r12d,%r12d
jle    13d1 <Update_H+0x2b1>
movsd  0xe80(%rip),%xmm1
divsd  (%rsp),%xmm1
movsd  0xe7b(%rip),%xmm2
divsd  %xmm0,%xmm2
mov    %r12d,%eax
lea    (%rbx,%rax,8),%rcx
mov    %rcx,0x18(%rsp)
lea    (%r14,%rax,8),%rcx
mov    %rcx,0x10(%rsp)
mov    %rax,%rcx
neg    %rcx
mov    %rcx,0x20(%rsp)
xor    %r11d,%r11d
mov    %r14,%rdi
mov    %rbx,%rbp
mov    %rax,%r8
jmp    11d8 <Update_H+0xb8>
nop
add    $0xffffffffffffffff,%r8
add    $0x8,%rbp
add    $0x8,%rdi
mov    %r13,%r15
cmp    %rax,%r11
je     13d1 <Update_H+0x2b1>
mov    %r11,%rdx
mov    %rax,%rsi
sub    %r11,%rsi
add    $0x1,%r11
cmp    $0x2,%rsi
jae    1210 <Update_H+0xf0>
mov    %rdx,%r10
mov    %r15,%r13
mov    %eax,%ecx
sub    %r10d,%ecx
test   $0x1,%cl
jne    1300 <Update_H+0x1e0>
mov    %r10,%rsi
mov    %r13,%rcx
jmp    133a <Update_H+0x21a>
nopw   0x0(%rax,%rax,1)
mov    %r8,(%rsp)
lea    (%rbx,%rdx,8),%r12
lea    (%r14,%rdx,8),%r9
lea    (%r14,%r11,8),%r10
mov    %r11,0x8(%rsp)
lea    (%rbx,%r11,8),%rcx
lea    (%r15,%rsi,8),%r11
cmp    0x18(%rsp),%r15
setb   %r8b
cmp    %r11,%r12
setb   %r13b
cmp    %rcx,%r15
setb   %cl
cmp    0x10(%rsp),%r15
setb   %r12b
cmp    %r11,%r9
setb   %r9b
cmp    %r10,%r15
setb   %r10b
test   %r13b,%r8b
jne    13bd <Update_H+0x29d>
and    %r13b,%cl
mov    (%rsp),%r8
mov    0x8(%rsp),%r11
jne    11eb <Update_H+0xcb>
and    %r9b,%r12b
jne    11eb <Update_H+0xcb>
and    %r9b,%r10b
jne    11eb <Update_H+0xcb>
mov    %r8,%r12
and    $0xfffffffffffffffe,%r12
mov    %rsi,%r9
and    $0xfffffffffffffffe,%r9
lea    (%rdx,%r9,1),%r10
lea    (%r15,%r9,8),%r13
movsd  (%rbx,%rdx,8),%xmm0
mulsd  %xmm1,%xmm0
unpcklpd %xmm0,%xmm0
movsd  (%r14,%rdx,8),%xmm3
mulsd  %xmm2,%xmm3
unpcklpd %xmm3,%xmm3
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movupd 0x0(%rbp,%rcx,8),%xmm4
mulpd  %xmm0,%xmm4
movupd (%rdi,%rcx,8),%xmm5
mulpd  %xmm3,%xmm5
addpd  %xmm4,%xmm5
movupd (%r15,%rcx,8),%xmm4
addpd  %xmm5,%xmm4
movupd %xmm4,(%r15,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%r12
jne    12c0 <Update_H+0x1a0>
cmp    %r9,%rsi
je     11c0 <Update_H+0xa0>
jmp    11f1 <Update_H+0xd1>
xchg   %ax,%ax
movsd  (%rbx,%rdx,8),%xmm0
mulsd  %xmm1,%xmm0
movsd  (%r14,%rdx,8),%xmm3
mulsd  %xmm2,%xmm3
mulsd  (%r14,%r10,8),%xmm3
mulsd  (%rbx,%r10,8),%xmm0
addsd  %xmm3,%xmm0
lea    0x8(%r13),%rcx
addsd  0x0(%r13),%xmm0
movsd  %xmm0,0x0(%r13)
lea    0x1(%r10),%rsi
mov    %rcx,%r13
not    %r10
cmp    0x20(%rsp),%r10
je     11c0 <Update_H+0xa0>
nopl   0x0(%rax,%rax,1)
movsd  (%rbx,%rdx,8),%xmm0
mulsd  %xmm1,%xmm0
movsd  (%r14,%rdx,8),%xmm3
mulsd  %xmm2,%xmm3
mulsd  (%r14,%rsi,8),%xmm3
mulsd  (%rbx,%rsi,8),%xmm0
addsd  %xmm3,%xmm0
addsd  (%rcx),%xmm0
movsd  %xmm0,(%rcx)
movsd  (%rbx,%rdx,8),%xmm0
mulsd  %xmm1,%xmm0
movsd  (%r14,%rdx,8),%xmm3
mulsd  %xmm2,%xmm3
mulsd  0x8(%r14,%rsi,8),%xmm3
mulsd  0x8(%rbx,%rsi,8),%xmm0
addsd  %xmm3,%xmm0
addsd  0x8(%rcx),%xmm0
movsd  %xmm0,0x8(%rcx)
add    $0x10,%rcx
add    $0x2,%rsi
cmp    %rsi,%rax
jne    1350 <Update_H+0x230>
mov    %rcx,%r13
jmp    11c0 <Update_H+0xa0>
mov    %rdx,%r10
mov    %r15,%r13
mov    (%rsp),%r8
mov    0x8(%rsp),%r11
jmp    11f1 <Update_H+0xd1>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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