#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 *_glp_lib_xcalloc(int, int);
int /*<<< orphan*/ _glp_lib_xfree(double *);
void _glp_mat_adat_numeric(int m, int n, int P_per[], int A_ptr[], int A_ind[],
                           double A_val[], double D_diag[], int S_ptr[],
                           int S_ind[], double S_val[], double S_diag[]) {
  int i, j, t, ii, jj, tt, beg, end, beg1, end1, k;
  double sum, *work;
  work = _glp_lib_xcalloc(1 + n, sizeof(double));
  for (j = 1; j <= n; j++)
    work[j] = 0.0;
  for (ii = 1; ii <= m; ii++) {
    i = P_per[ii];
    beg = A_ptr[i], end = A_ptr[i + 1];
    for (t = beg; t < end; t++)
      work[A_ind[t]] = A_val[t];
    beg = S_ptr[ii], end = S_ptr[ii + 1];
    for (t = beg; t < end; t++) {
      jj = S_ind[t];
      j = P_per[jj];
      sum = 0.0;
      beg1 = A_ptr[j], end1 = A_ptr[j + 1];
      for (tt = beg1; tt < end1; tt++) {
        k = A_ind[tt];
        sum += work[k] * D_diag[k] * A_val[tt];
      }
      S_val[t] = sum;
    }
    sum = 0.0;
    beg = A_ptr[i], end = A_ptr[i + 1];
    for (t = beg; t < end; t++) {
      k = A_ind[t];
      sum += A_val[t] * D_diag[k] * A_val[t];
      work[k] = 0.0;
    }
    S_diag[ii] = sum;
  }
  _glp_lib_xfree(work);
  return;
}


x86 assembly (clang, -O0):
<_glp_mat_adat_numeric>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    -0x8(%rbp),%edi
add    $0x1,%edi
mov    $0x8,%esi
call   1030 <_glp_lib_xcalloc@plt>
mov    %rax,-0x68(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jg     1197 <_glp_mat_adat_numeric+0x77>
mov    -0x68(%rbp),%rax
movslq -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    116d <_glp_mat_adat_numeric+0x4d>
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     13da <_glp_mat_adat_numeric+0x2ba>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x44(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x48(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x48(%rbp),%eax
jge    121b <_glp_mat_adat_numeric+0xfb>
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x68(%rbp),%rax
mov    -0x20(%rbp),%rcx
movslq -0x34(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11df <_glp_mat_adat_numeric+0xbf>
mov    0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x44(%rbp)
mov    0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x48(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x48(%rbp),%eax
jge    131a <_glp_mat_adat_numeric+0x1fa>
mov    0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x50(%rbp)
mov    -0x4c(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x50(%rbp),%eax
jge    12fa <_glp_mat_adat_numeric+0x1da>
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x54(%rbp)
mov    -0x68(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1299 <_glp_mat_adat_numeric+0x179>
movsd  -0x60(%rbp),%xmm0
mov    0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1242 <_glp_mat_adat_numeric+0x122>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x44(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x48(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x48(%rbp),%eax
jge    13ba <_glp_mat_adat_numeric+0x29a>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x54(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x68(%rbp),%rax
movslq -0x54(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1349 <_glp_mat_adat_numeric+0x229>
movsd  -0x60(%rbp),%xmm0
mov    0x30(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    119e <_glp_mat_adat_numeric+0x7e>
mov    -0x68(%rbp),%rdi
call   1040 <_glp_lib_xfree@plt>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<_glp_mat_adat_numeric>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r9,%r14
mov    %r8,%rbx
mov    %rcx,%r15
mov    %rdx,%r12
mov    %esi,%ebp
mov    %rdi,(%rsp)
mov    0x80(%rsp),%rax
mov    %rax,0x10(%rsp)
mov    0x78(%rsp),%rax
mov    %rax,0x20(%rsp)
mov    0x70(%rsp),%rax
mov    %rax,0x18(%rsp)
mov    0x68(%rsp),%rax
mov    %rax,0x8(%rsp)
mov    0x60(%rsp),%r13
lea    0x1(%rbp),%edi
mov    $0x8,%esi
call   1030 <_glp_lib_xcalloc@plt>
mov    %rbp,%rcx
mov    %rax,%rbp
test   %ecx,%ecx
jle    11a8 <_glp_mat_adat_numeric+0x78>
lea    0x8(%rbp),%rdi
mov    %ecx,%edx
shl    $0x3,%rdx
xor    %esi,%esi
call   1050 <memset@plt>
cmpl   $0x0,(%rsp)
jle    1416 <_glp_mat_adat_numeric+0x2e6>
mov    (%rsp),%rax
add    $0x1,%eax
mov    %rax,(%rsp)
mov    $0x1,%esi
xorpd  %xmm0,%xmm0
jmp    11e5 <_glp_mat_adat_numeric+0xb5>
nopl   0x0(%rax,%rax,1)
mov    0x10(%rsp),%rax
movsd  %xmm1,(%rax,%r8,8)
cmp    (%rsp),%rsi
je     1416 <_glp_mat_adat_numeric+0x2e6>
mov    %rsi,%r8
movslq (%r12,%rsi,4),%rax
movslq (%r15,%rax,4),%r9
movslq 0x4(%r15,%rax,4),%r10
cmp    %r10d,%r9d
jge    128f <_glp_mat_adat_numeric+0x15f>
mov    %r10d,%ecx
sub    %r9d,%ecx
mov    %r9,%rdx
not    %rdx
add    %r10,%rdx
mov    %r9,%rax
and    $0x3,%rcx
je     123a <_glp_mat_adat_numeric+0x10a>
mov    %r9,%rax
nopl   0x0(%rax)
movsd  (%r14,%rax,8),%xmm1
movslq (%rbx,%rax,4),%rsi
movsd  %xmm1,0x0(%rbp,%rsi,8)
add    $0x1,%rax
add    $0xffffffffffffffff,%rcx
jne    1220 <_glp_mat_adat_numeric+0xf0>
cmp    $0x3,%rdx
jb     128f <_glp_mat_adat_numeric+0x15f>
movsd  (%r14,%rax,8),%xmm1
movslq (%rbx,%rax,4),%rcx
movsd  %xmm1,0x0(%rbp,%rcx,8)
movsd  0x8(%r14,%rax,8),%xmm1
movslq 0x4(%rbx,%rax,4),%rcx
movsd  %xmm1,0x0(%rbp,%rcx,8)
movsd  0x10(%r14,%rax,8),%xmm1
movslq 0x8(%rbx,%rax,4),%rcx
movsd  %xmm1,0x0(%rbp,%rcx,8)
movsd  0x18(%r14,%rax,8),%xmm1
movslq 0xc(%rbx,%rax,4),%rcx
movsd  %xmm1,0x0(%rbp,%rcx,8)
add    $0x4,%rax
cmp    %rax,%r10
jne    1240 <_glp_mat_adat_numeric+0x110>
mov    0x8(%rsp),%rax
movslq (%rax,%r8,4),%rsi
movslq 0x4(%rax,%r8,4),%r11
cmp    %r11d,%esi
jl     1377 <_glp_mat_adat_numeric+0x247>
lea    0x1(%r8),%rsi
xorpd  %xmm1,%xmm1
cmp    %r10d,%r9d
jge    11d0 <_glp_mat_adat_numeric+0xa0>
mov    %r10d,%ecx
sub    %r9d,%ecx
xorpd  %xmm1,%xmm1
mov    %r9,%rax
test   $0x1,%cl
je     12f3 <_glp_mat_adat_numeric+0x1c3>
movslq (%rbx,%r9,4),%rax
movsd  (%r14,%r9,8),%xmm2
movsd  0x0(%r13,%rax,8),%xmm1
mulsd  %xmm2,%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movq   $0x0,0x0(%rbp,%rax,8)
lea    0x1(%r9),%rax
not    %r9
add    %r10,%r9
je     11d0 <_glp_mat_adat_numeric+0xa0>
nop
movslq (%rbx,%rax,4),%rcx
movsd  (%r14,%rax,8),%xmm2
movsd  0x0(%r13,%rcx,8),%xmm3
mulsd  %xmm2,%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm1,%xmm3
movq   $0x0,0x0(%rbp,%rcx,8)
movslq 0x4(%rbx,%rax,4),%rcx
movsd  0x8(%r14,%rax,8),%xmm2
movsd  0x0(%r13,%rcx,8),%xmm1
mulsd  %xmm2,%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm3,%xmm1
movq   $0x0,0x0(%rbp,%rcx,8)
add    $0x2,%rax
cmp    %rax,%r10
jne    1300 <_glp_mat_adat_numeric+0x1d0>
jmp    11d0 <_glp_mat_adat_numeric+0xa0>
nopl   0x0(%rax)
mov    0x20(%rsp),%rax
movsd  %xmm1,(%rax,%rsi,8)
add    $0x1,%rsi
cmp    %r11,%rsi
je     12a6 <_glp_mat_adat_numeric+0x176>
mov    0x18(%rsp),%rax
movslq (%rax,%rsi,4),%rax
movslq (%r12,%rax,4),%rax
movslq (%r15,%rax,4),%rcx
movslq 0x4(%r15,%rax,4),%rax
xorpd  %xmm1,%xmm1
cmp    %eax,%ecx
jge    1360 <_glp_mat_adat_numeric+0x230>
mov    %eax,%edi
sub    %ecx,%edi
xorpd  %xmm1,%xmm1
mov    %rcx,%rdx
test   $0x1,%dil
je     13c5 <_glp_mat_adat_numeric+0x295>
movslq (%rbx,%rcx,4),%rdx
movsd  0x0(%rbp,%rdx,8),%xmm1
mulsd  0x0(%r13,%rdx,8),%xmm1
mulsd  (%r14,%rcx,8),%xmm1
addsd  %xmm0,%xmm1
lea    0x1(%rcx),%rdx
not    %rcx
add    %rax,%rcx
je     1360 <_glp_mat_adat_numeric+0x230>
nopl   (%rax)
movslq (%rbx,%rdx,4),%rcx
movsd  0x0(%rbp,%rcx,8),%xmm2
mulsd  0x0(%r13,%rcx,8),%xmm2
mulsd  (%r14,%rdx,8),%xmm2
addsd  %xmm1,%xmm2
movslq 0x4(%rbx,%rdx,4),%rcx
movsd  0x0(%rbp,%rcx,8),%xmm1
mulsd  0x0(%r13,%rcx,8),%xmm1
mulsd  0x8(%r14,%rdx,8),%xmm1
addsd  %xmm2,%xmm1
add    $0x2,%rdx
cmp    %rdx,%rax
jne    13d0 <_glp_mat_adat_numeric+0x2a0>
jmp    1360 <_glp_mat_adat_numeric+0x230>
mov    %rbp,%rdi
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <_glp_lib_xfree@plt>

LLVM IR (-O0):
define void @_glp_mat_adat_numeric(i32 %0, i32 %1, i32* %2, i32* %3, i32* %4, double* %5, double* %6, i32* %7, i32* %8, double* %9, double* %10) {
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32*
%15 = alloca i32*
%16 = alloca i32*
%17 = alloca double*
%18 = alloca double*
%19 = alloca i32*
%20 = alloca i32*
%21 = alloca double*
%22 = alloca double*
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca i32
%27 = alloca i32
%28 = alloca i32
%29 = alloca i32
%30 = alloca i32
%31 = alloca i32
%32 = alloca i32
%33 = alloca i32
%34 = alloca double
%35 = alloca double*
store i32 %0, i32* %12
store i32 %1, i32* %13
store i32* %2, i32** %14
store i32* %3, i32** %15
store i32* %4, i32** %16
store double* %5, double** %17
store double* %6, double** %18
store i32* %7, i32** %19
store i32* %8, i32** %20
store double* %9, double** %21
store double* %10, double** %22
%36 = load i32, i32* %13
%37 = add nsw i32 1, %36
%38 = call double* @_glp_lib_xcalloc(i32 %37, i32 8)
store double* %38, double** %35
store i32 1, i32* %24
br label %39
39: ; preds = %48, %11
%40 = load i32, i32* %24
%41 = load i32, i32* %13
%42 = icmp sle i32 %40, %41
br i1 %42, label %43, label %51
43: ; preds = %39
%44 = load double*, double** %35
%45 = load i32, i32* %24
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds double, double* %44, i64 %46
store double 0.000000e+00, double* %47
br label %48
48: ; preds = %43
%49 = load i32, i32* %24
%50 = add nsw i32 %49, 1
store i32 %50, i32* %24
br label %39
51: ; preds = %39
store i32 1, i32* %26
br label %52
52: ; preds = %229, %51
%53 = load i32, i32* %26
%54 = load i32, i32* %12
%55 = icmp sle i32 %53, %54
br i1 %55, label %56, label %232
56: ; preds = %52
%57 = load i32*, i32** %14
%58 = load i32, i32* %26
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds i32, i32* %57, i64 %59
%61 = load i32, i32* %60
store i32 %61, i32* %23
%62 = load i32*, i32** %15
%63 = load i32, i32* %23
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds i32, i32* %62, i64 %64
%66 = load i32, i32* %65
store i32 %66, i32* %29
%67 = load i32*, i32** %15
%68 = load i32, i32* %23
%69 = add nsw i32 %68, 1
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds i32, i32* %67, i64 %70
%72 = load i32, i32* %71
store i32 %72, i32* %30
%73 = load i32, i32* %29
store i32 %73, i32* %25
br label %74
74: ; preds = %92, %56
%75 = load i32, i32* %25
%76 = load i32, i32* %30
%77 = icmp slt i32 %75, %76
br i1 %77, label %78, label %95
78: ; preds = %74
%79 = load double*, double** %17
%80 = load i32, i32* %25
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %79, i64 %81
%83 = load double, double* %82
%84 = load double*, double** %35
%85 = load i32*, i32** %16
%86 = load i32, i32* %25
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds i32, i32* %85, i64 %87
%89 = load i32, i32* %88
%90 = sext i32 %89 to i64
%91 = getelementptr inbounds double, double* %84, i64 %90
store double %83, double* %91
br label %92
92: ; preds = %78
%93 = load i32, i32* %25
%94 = add nsw i32 %93, 1
store i32 %94, i32* %25
br label %74
95: ; preds = %74
%96 = load i32*, i32** %19
%97 = load i32, i32* %26
%98 = sext i32 %97 to i64
%99 = getelementptr inbounds i32, i32* %96, i64 %98
%100 = load i32, i32* %99
store i32 %100, i32* %29
%101 = load i32*, i32** %19
%102 = load i32, i32* %26
%103 = add nsw i32 %102, 1
%104 = sext i32 %103 to i64
%105 = getelementptr inbounds i32, i32* %101, i64 %104
%106 = load i32, i32* %105
store i32 %106, i32* %30
%107 = load i32, i32* %29
store i32 %107, i32* %25
br label %108
108: ; preds = %172, %95
%109 = load i32, i32* %25
%110 = load i32, i32* %30
%111 = icmp slt i32 %109, %110
br i1 %111, label %112, label %175
112: ; preds = %108
%113 = load i32*, i32** %20
%114 = load i32, i32* %25
%115 = sext i32 %114 to i64
%116 = getelementptr inbounds i32, i32* %113, i64 %115
%117 = load i32, i32* %116
store i32 %117, i32* %27
%118 = load i32*, i32** %14
%119 = load i32, i32* %27
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds i32, i32* %118, i64 %120
%122 = load i32, i32* %121
store i32 %122, i32* %24
store double 0.000000e+00, double* %34
%123 = load i32*, i32** %15
%124 = load i32, i32* %24
%125 = sext i32 %124 to i64
%126 = getelementptr inbounds i32, i32* %123, i64 %125
%127 = load i32, i32* %126
store i32 %127, i32* %31
%128 = load i32*, i32** %15
%129 = load i32, i32* %24
%130 = add nsw i32 %129, 1
%131 = sext i32 %130 to i64
%132 = getelementptr inbounds i32, i32* %128, i64 %131
%133 = load i32, i32* %132
store i32 %133, i32* %32
%134 = load i32, i32* %31
store i32 %134, i32* %28
br label %135
135: ; preds = %163, %112
%136 = load i32, i32* %28
%137 = load i32, i32* %32
%138 = icmp slt i32 %136, %137
br i1 %138, label %139, label %166
139: ; preds = %135
%140 = load i32*, i32** %16
%141 = load i32, i32* %28
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds i32, i32* %140, i64 %142
%144 = load i32, i32* %143
store i32 %144, i32* %33
%145 = load double*, double** %35
%146 = load i32, i32* %33
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds double, double* %145, i64 %147
%149 = load double, double* %148
%150 = load double*, double** %18
%151 = load i32, i32* %33
%152 = sext i32 %151 to i64
%153 = getelementptr inbounds double, double* %150, i64 %152
%154 = load double, double* %153
%155 = fmul double %149, %154
%156 = load double*, double** %17
%157 = load i32, i32* %28
%158 = sext i32 %157 to i64
%159 = getelementptr inbounds double, double* %156, i64 %158
%160 = load double, double* %159
%161 = load double, double* %34
%162 = call double @llvm.fmuladd.f64(double %155, double %160, double %161)
store double %162, double* %34
br label %163
163: ; preds = %139
%164 = load i32, i32* %28
%165 = add nsw i32 %164, 1
store i32 %165, i32* %28
br label %135
166: ; preds = %135
%167 = load double, double* %34
%168 = load double*, double** %21
%169 = load i32, i32* %25
%170 = sext i32 %169 to i64
%171 = getelementptr inbounds double, double* %168, i64 %170
store double %167, double* %171
br label %172
172: ; preds = %166
%173 = load i32, i32* %25
%174 = add nsw i32 %173, 1
store i32 %174, i32* %25
br label %108
175: ; preds = %108
store double 0.000000e+00, double* %34
%176 = load i32*, i32** %15
%177 = load i32, i32* %23
%178 = sext i32 %177 to i64
%179 = getelementptr inbounds i32, i32* %176, i64 %178
%180 = load i32, i32* %179
store i32 %180, i32* %29
%181 = load i32*, i32** %15
%182 = load i32, i32* %23
%183 = add nsw i32 %182, 1
%184 = sext i32 %183 to i64
%185 = getelementptr inbounds i32, i32* %181, i64 %184
%186 = load i32, i32* %185
store i32 %186, i32* %30
%187 = load i32, i32* %29
store i32 %187, i32* %25
br label %188
188: ; preds = %220, %175
%189 = load i32, i32* %25
%190 = load i32, i32* %30
%191 = icmp slt i32 %189, %190
br i1 %191, label %192, label %223
192: ; preds = %188
%193 = load i32*, i32** %16
%194 = load i32, i32* %25
%195 = sext i32 %194 to i64
%196 = getelementptr inbounds i32, i32* %193, i64 %195
%197 = load i32, i32* %196
store i32 %197, i32* %33
%198 = load double*, double** %17
%199 = load i32, i32* %25
%200 = sext i32 %199 to i64
%201 = getelementptr inbounds double, double* %198, i64 %200
%202 = load double, double* %201
%203 = load double*, double** %18
%204 = load i32, i32* %33
%205 = sext i32 %204 to i64
%206 = getelementptr inbounds double, double* %203, i64 %205
%207 = load double, double* %206
%208 = fmul double %202, %207
%209 = load double*, double** %17
%210 = load i32, i32* %25
%211 = sext i32 %210 to i64
%212 = getelementptr inbounds double, double* %209, i64 %211
%213 = load double, double* %212
%214 = load double, double* %34
%215 = call double @llvm.fmuladd.f64(double %208, double %213, double %214)
store double %215, double* %34
%216 = load double*, double** %35
%217 = load i32, i32* %33
%218 = sext i32 %217 to i64
%219 = getelementptr inbounds double, double* %216, i64 %218
store double 0.000000e+00, double* %219
br label %220
220: ; preds = %192
%221 = load i32, i32* %25
%222 = add nsw i32 %221, 1
store i32 %222, i32* %25
br label %188
223: ; preds = %188
%224 = load double, double* %34
%225 = load double*, double** %22
%226 = load i32, i32* %26
%227 = sext i32 %226 to i64
%228 = getelementptr inbounds double, double* %225, i64 %227
store double %224, double* %228
br label %229
229: ; preds = %223
%230 = load i32, i32* %26
%231 = add nsw i32 %230, 1
store i32 %231, i32* %26
br label %52
232: ; preds = %52
%233 = load double*, double** %35
%234 = call i32 @_glp_lib_xfree(double* %233)
ret void
}

LLVM IR (-O2):
define void @_glp_mat_adat_numeric(i32 %0, i32 %1, i32* nocapture %2, i32* nocapture %3, i32* nocapture %4, double* nocapture %5, double* nocapture %6, i32* nocapture %7, i32* nocapture %8, double* nocapture writeonly %9, double* nocapture writeonly %10) {
%12 = add nsw i32 %1, 1
%13 = tail call double* @_glp_lib_xcalloc(i32 %12, i32 8)
%14 = icmp slt i32 %1, 1
br i1 %14, label %20, label %15
15: ; preds = %11
%16 = getelementptr double, double* %13, i64 1
%17 = bitcast double* %16 to i8*
%18 = zext i32 %1 to i64
%19 = shl nuw nsw i64 %18, 3
br label %20
20: ; preds = %15, %11
%21 = icmp slt i32 %0, 1
br i1 %21, label %231, label %22
22: ; preds = %20
%23 = add nuw i32 %0, 1
%24 = zext i32 %23 to i64
br label %25
25: ; preds = %22, %227
%26 = phi i64 [ 1, %22 ], [ %94, %227 ]
%27 = getelementptr inbounds i32, i32* %2, i64 %26
%28 = load i32, i32* %27
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds i32, i32* %3, i64 %29
%31 = load i32, i32* %30
%32 = add nsw i32 %28, 1
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds i32, i32* %3, i64 %33
%35 = load i32, i32* %34
%36 = icmp slt i32 %31, %35
br i1 %36, label %37, label %91
37: ; preds = %25
%38 = sext i32 %31 to i64
%39 = sext i32 %35 to i64
%40 = sub nsw i64 %39, %38
%41 = xor i64 %38, -1
%42 = add nsw i64 %41, %39
%43 = and i64 %40, 3
%44 = icmp eq i64 %43, 0
br i1 %44, label %57, label %45
45: ; preds = %37, %45
%46 = phi i64 [ %54, %45 ], [ %38, %37 ]
%47 = phi i64 [ %55, %45 ], [ 0, %37 ]
%48 = getelementptr inbounds double, double* %5, i64 %46
%49 = load double, double* %48
%50 = getelementptr inbounds i32, i32* %4, i64 %46
%51 = load i32, i32* %50
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds double, double* %13, i64 %52
store double %49, double* %53
%54 = add nsw i64 %46, 1
%55 = add i64 %47, 1
%56 = icmp eq i64 %55, %43
br i1 %56, label %57, label %45
57: ; preds = %45, %37
%58 = phi i64 [ %38, %37 ], [ %54, %45 ]
%59 = icmp ult i64 %42, 3
br i1 %59, label %91, label %60
60: ; preds = %57, %60
%61 = phi i64 [ %89, %60 ], [ %58, %57 ]
%62 = getelementptr inbounds double, double* %5, i64 %61
%63 = load double, double* %62
%64 = getelementptr inbounds i32, i32* %4, i64 %61
%65 = load i32, i32* %64
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds double, double* %13, i64 %66
store double %63, double* %67
%68 = add nsw i64 %61, 1
%69 = getelementptr inbounds double, double* %5, i64 %68
%70 = load double, double* %69
%71 = getelementptr inbounds i32, i32* %4, i64 %68
%72 = load i32, i32* %71
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %13, i64 %73
store double %70, double* %74
%75 = add nsw i64 %61, 2
%76 = getelementptr inbounds double, double* %5, i64 %75
%77 = load double, double* %76
%78 = getelementptr inbounds i32, i32* %4, i64 %75
%79 = load i32, i32* %78
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds double, double* %13, i64 %80
store double %77, double* %81
%82 = add nsw i64 %61, 3
%83 = getelementptr inbounds double, double* %5, i64 %82
%84 = load double, double* %83
%85 = getelementptr inbounds i32, i32* %4, i64 %82
%86 = load i32, i32* %85
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds double, double* %13, i64 %87
store double %84, double* %88
%89 = add nsw i64 %61, 4
%90 = icmp eq i64 %89, %39
br i1 %90, label %91, label %60
91: ; preds = %57, %60, %25
%92 = getelementptr inbounds i32, i32* %7, i64 %26
%93 = load i32, i32* %92
%94 = add nuw nsw i64 %26, 1
%95 = getelementptr inbounds i32, i32* %7, i64 %94
%96 = load i32, i32* %95
%97 = icmp slt i32 %93, %96
br i1 %97, label %98, label %101
98: ; preds = %91
%99 = sext i32 %93 to i64
%100 = sext i32 %96 to i64
br label %127
101: ; preds = %196, %91
br i1 %36, label %102, label %227
102: ; preds = %101
%103 = sext i32 %31 to i64
%104 = sext i32 %35 to i64
%105 = sub nsw i64 %104, %103
%106 = xor i64 %103, -1
%107 = and i64 %105, 1
%108 = icmp eq i64 %107, 0
br i1 %108, label %121, label %109
109: ; preds = %102
%110 = getelementptr inbounds i32, i32* %4, i64 %103
%111 = load i32, i32* %110
%112 = getelementptr inbounds double, double* %5, i64 %103
%113 = load double, double* %112
%114 = sext i32 %111 to i64
%115 = getelementptr inbounds double, double* %6, i64 %114
%116 = load double, double* %115
%117 = fmul double %113, %116
%118 = tail call double @llvm.fmuladd.f64(double %117, double %113, double 0.000000e+00)
%119 = getelementptr inbounds double, double* %13, i64 %114
store double 0.000000e+00, double* %119
%120 = add nsw i64 %103, 1
br label %121
121: ; preds = %109, %102
%122 = phi double [ undef, %102 ], [ %118, %109 ]
%123 = phi i64 [ %103, %102 ], [ %120, %109 ]
%124 = phi double [ 0.000000e+00, %102 ], [ %118, %109 ]
%125 = sub nsw i64 0, %104
%126 = icmp eq i64 %106, %125
br i1 %126, label %227, label %201
127: ; preds = %98, %196
%128 = phi i64 [ %99, %98 ], [ %199, %196 ]
%129 = getelementptr inbounds i32, i32* %8, i64 %128
%130 = load i32, i32* %129
%131 = sext i32 %130 to i64
%132 = getelementptr inbounds i32, i32* %2, i64 %131
%133 = load i32, i32* %132
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds i32, i32* %3, i64 %134
%136 = load i32, i32* %135
%137 = add nsw i32 %133, 1
%138 = sext i32 %137 to i64
%139 = getelementptr inbounds i32, i32* %3, i64 %138
%140 = load i32, i32* %139
%141 = icmp slt i32 %136, %140
br i1 %141, label %142, label %196
142: ; preds = %127
%143 = sext i32 %136 to i64
%144 = sext i32 %140 to i64
%145 = sub nsw i64 %144, %143
%146 = xor i64 %143, -1
%147 = and i64 %145, 1
%148 = icmp eq i64 %147, 0
br i1 %148, label %162, label %149
149: ; preds = %142
%150 = getelementptr inbounds i32, i32* %4, i64 %143
%151 = load i32, i32* %150
%152 = sext i32 %151 to i64
%153 = getelementptr inbounds double, double* %13, i64 %152
%154 = load double, double* %153
%155 = getelementptr inbounds double, double* %6, i64 %152
%156 = load double, double* %155
%157 = fmul double %154, %156
%158 = getelementptr inbounds double, double* %5, i64 %143
%159 = load double, double* %158
%160 = tail call double @llvm.fmuladd.f64(double %157, double %159, double 0.000000e+00)
%161 = add nsw i64 %143, 1
br label %162
162: ; preds = %149, %142
%163 = phi double [ undef, %142 ], [ %160, %149 ]
%164 = phi i64 [ %143, %142 ], [ %161, %149 ]
%165 = phi double [ 0.000000e+00, %142 ], [ %160, %149 ]
%166 = sub nsw i64 0, %144
%167 = icmp eq i64 %146, %166
br i1 %167, label %196, label %168
168: ; preds = %162, %168
%169 = phi i64 [ %194, %168 ], [ %164, %162 ]
%170 = phi double [ %193, %168 ], [ %165, %162 ]
%171 = getelementptr inbounds i32, i32* %4, i64 %169
%172 = load i32, i32* %171
%173 = sext i32 %172 to i64
%174 = getelementptr inbounds double, double* %13, i64 %173
%175 = load double, double* %174
%176 = getelementptr inbounds double, double* %6, i64 %173
%177 = load double, double* %176
%178 = fmul double %175, %177
%179 = getelementptr inbounds double, double* %5, i64 %169
%180 = load double, double* %179
%181 = tail call double @llvm.fmuladd.f64(double %178, double %180, double %170)
%182 = add nsw i64 %169, 1
%183 = getelementptr inbounds i32, i32* %4, i64 %182
%184 = load i32, i32* %183
%185 = sext i32 %184 to i64
%186 = getelementptr inbounds double, double* %13, i64 %185
%187 = load double, double* %186
%188 = getelementptr inbounds double, double* %6, i64 %185
%189 = load double, double* %188
%190 = fmul double %187, %189
%191 = getelementptr inbounds double, double* %5, i64 %182
%192 = load double, double* %191
%193 = tail call double @llvm.fmuladd.f64(double %190, double %192, double %181)
%194 = add nsw i64 %169, 2
%195 = icmp eq i64 %194, %144
br i1 %195, label %196, label %168
196: ; preds = %162, %168, %127
%197 = phi double [ 0.000000e+00, %127 ], [ %163, %162 ], [ %193, %168 ]
%198 = getelementptr inbounds double, double* %9, i64 %128
store double %197, double* %198
%199 = add nsw i64 %128, 1
%200 = icmp eq i64 %199, %100
br i1 %200, label %101, label %127
201: ; preds = %121, %201
%202 = phi i64 [ %225, %201 ], [ %123, %121 ]
%203 = phi double [ %223, %201 ], [ %124, %121 ]
%204 = getelementptr inbounds i32, i32* %4, i64 %202
%205 = load i32, i32* %204
%206 = getelementptr inbounds double, double* %5, i64 %202
%207 = load double, double* %206
%208 = sext i32 %205 to i64
%209 = getelementptr inbounds double, double* %6, i64 %208
%210 = load double, double* %209
%211 = fmul double %207, %210
%212 = tail call double @llvm.fmuladd.f64(double %211, double %207, double %203)
%213 = getelementptr inbounds double, double* %13, i64 %208
store double 0.000000e+00, double* %213
%214 = add nsw i64 %202, 1
%215 = getelementptr inbounds i32, i32* %4, i64 %214
%216 = load i32, i32* %215
%217 = getelementptr inbounds double, double* %5, i64 %214
%218 = load double, double* %217
%219 = sext i32 %216 to i64
%220 = getelementptr inbounds double, double* %6, i64 %219
%221 = load double, double* %220
%222 = fmul double %218, %221
%223 = tail call double @llvm.fmuladd.f64(double %222, double %218, double %212)
%224 = getelementptr inbounds double, double* %13, i64 %219
store double 0.000000e+00, double* %224
%225 = add nsw i64 %202, 2
%226 = icmp eq i64 %225, %104
br i1 %226, label %227, label %201
227: ; preds = %121, %201, %101
%228 = phi double [ 0.000000e+00, %101 ], [ %122, %121 ], [ %223, %201 ]
%229 = getelementptr inbounds double, double* %10, i64 %26
store double %228, double* %229
%230 = icmp eq i64 %94, %24
br i1 %230, label %231, label %25
231: ; preds = %227, %20
%232 = tail call i32 @_glp_lib_xfree(double* %13)
ret void
}