void ZuN(double *Z, double *u, double *T1, double *T2, int len) {
  Z[0 + 0 * len * 7] = -1;
  Z[0 + 6 * len * 7] = (u[0]) * T1[0] + (u[1]) * T1[3] + T1[6];
  Z[1 + 1 * len * 7] = -1;
  Z[1 + 7 * len * 7] = (u[0]) * T1[0] + (u[1]) * T1[3] + T1[6];
  Z[2 + 2 * len * 7] = -1;
  Z[2 + 6 * len * 7] = -((u[0]) * T1[0] + (u[1]) * T1[3] + T1[6]) *
                       ((u[5]) * T2[0] + (u[6]) * T2[3] + T2[6]);
  Z[2 + 7 * len * 7] = -((u[0]) * T1[0] + (u[1]) * T1[3] + T1[6]) *
                       ((u[5]) * T2[1] + (u[6]) * T2[4] + T2[7]);
  Z[3 + 3 * len * 7] = -1;
  Z[3 + 6 * len * 7] = (u[0]) * T1[1] + (u[1]) * T1[4] + T1[7];
  Z[4 + 4 * len * 7] = -1;
  Z[4 + 7 * len * 7] = (u[0]) * T1[1] + (u[1]) * T1[4] + T1[7];
  Z[5 + 5 * len * 7] = -1;
  Z[5 + 6 * len * 7] = -((u[0]) * T1[1] + (u[1]) * T1[4] + T1[7]) *
                       ((u[5]) * T2[0] + (u[6]) * T2[3] + T2[6]);
  Z[5 + 7 * len * 7] = -((u[0]) * T1[1] + (u[1]) * T1[4] + T1[7]) *
                       ((u[5]) * T2[1] + (u[6]) * T2[4] + T2[7]);
  Z[6 + 8 * len * 7] = -1;
  Z[6 + 6 * len * 7] = -((u[5]) * T2[0] + (u[6]) * T2[3] + T2[6]);
  Z[6 + 7 * len * 7] = -((u[5]) * T2[1] + (u[6]) * T2[4] + T2[7]);
}


x86 assembly (clang, -O0):
<ZuN>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    -0x8(%rbp),%rax
imul   $0x0,-0x24(%rbp),%ecx
imul   $0x7,%ecx,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  0xecf(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x18(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
addsd  0x30(%rax),%xmm0
mov    -0x8(%rbp),%rax
imul   $0x6,-0x24(%rbp),%ecx
imul   $0x7,%ecx,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x0,%ecx
imul   $0x7,%ecx,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  0xe66(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x18(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
addsd  0x30(%rax),%xmm0
mov    -0x8(%rbp),%rax
imul   $0x7,-0x24(%rbp),%ecx
imul   $0x7,%ecx,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
imul   $0x7,%ecx,%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  0xdfd(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x18(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
addsd  0x30(%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm1
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm3
mov    -0x10(%rbp),%rax
movsd  0x30(%rax),%xmm2
mov    -0x20(%rbp),%rax
mulsd  0x18(%rax),%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mov    -0x20(%rbp),%rax
addsd  0x30(%rax),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
imul   $0x6,-0x24(%rbp),%ecx
imul   $0x7,%ecx,%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x18(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
addsd  0x30(%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm1
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm3
mov    -0x10(%rbp),%rax
movsd  0x30(%rax),%xmm2
mov    -0x20(%rbp),%rax
mulsd  0x20(%rax),%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mov    -0x20(%rbp),%rax
addsd  0x38(%rax),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
imul   $0x7,-0x24(%rbp),%ecx
imul   $0x7,%ecx,%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
imul   $0x3,-0x24(%rbp),%ecx
imul   $0x7,%ecx,%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
movsd  0xcae(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm2
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x20(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
addsd  0x38(%rax),%xmm0
mov    -0x8(%rbp),%rax
imul   $0x6,-0x24(%rbp),%ecx
imul   $0x7,%ecx,%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x2,%ecx
imul   $0x7,%ecx,%ecx
add    $0x4,%ecx
movslq %ecx,%rcx
movsd  0xc44(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm2
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x20(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
addsd  0x38(%rax),%xmm0
mov    -0x8(%rbp),%rax
imul   $0x7,-0x24(%rbp),%ecx
imul   $0x7,%ecx,%ecx
add    $0x4,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
imul   $0x5,-0x24(%rbp),%ecx
imul   $0x7,%ecx,%ecx
add    $0x5,%ecx
movslq %ecx,%rcx
movsd  0xbdc(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm2
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x20(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
addsd  0x38(%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm1
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm3
mov    -0x10(%rbp),%rax
movsd  0x30(%rax),%xmm2
mov    -0x20(%rbp),%rax
mulsd  0x18(%rax),%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mov    -0x20(%rbp),%rax
addsd  0x30(%rax),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
imul   $0x6,-0x24(%rbp),%ecx
imul   $0x7,%ecx,%ecx
add    $0x5,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm2
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x20(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
addsd  0x38(%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm1
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm3
mov    -0x10(%rbp),%rax
movsd  0x30(%rax),%xmm2
mov    -0x20(%rbp),%rax
mulsd  0x20(%rax),%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mov    -0x20(%rbp),%rax
addsd  0x38(%rax),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
imul   $0x7,-0x24(%rbp),%ecx
imul   $0x7,%ecx,%ecx
add    $0x5,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x3,%ecx
imul   $0x7,%ecx,%ecx
add    $0x6,%ecx
movslq %ecx,%rcx
movsd  0xa89(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm0
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x10(%rbp),%rax
movsd  0x30(%rax),%xmm1
mov    -0x20(%rbp),%rax
mulsd  0x18(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
addsd  0x30(%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x8(%rbp),%rax
imul   $0x6,-0x24(%rbp),%ecx
imul   $0x7,%ecx,%ecx
add    $0x6,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm0
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm2
mov    -0x10(%rbp),%rax
movsd  0x30(%rax),%xmm1
mov    -0x20(%rbp),%rax
mulsd  0x20(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
addsd  0x38(%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x8(%rbp),%rax
imul   $0x7,-0x24(%rbp),%ecx
imul   $0x7,%ecx,%ecx
add    $0x6,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
pop    %rbp
ret

x86 assembly (clang, -O2):
<ZuN>:
movabs $0xbff0000000000000,%r9
mov    %r9,(%rdi)
movsd  0x8(%rsi),%xmm0
mulsd  0x18(%rdx),%xmm0
movsd  (%rsi),%xmm1
mulsd  (%rdx),%xmm1
addsd  %xmm0,%xmm1
addsd  0x30(%rdx),%xmm1
imul   $0x2a,%r8d,%eax
movslq %eax,%r10
movsd  %xmm1,(%rdi,%r10,8)
lea    0x0(,%r8,8),%eax
sub    %r8d,%eax
cltq
mov    %r9,0x8(%rdi,%rax,8)
movsd  0x8(%rsi),%xmm0
mulsd  0x18(%rdx),%xmm0
movsd  (%rsi),%xmm1
mulsd  (%rdx),%xmm1
addsd  %xmm0,%xmm1
addsd  0x30(%rdx),%xmm1
imul   $0x31,%r8d,%eax
movslq %eax,%r11
movsd  %xmm1,0x8(%rdi,%r11,8)
mov    %r8d,%eax
shl    $0x4,%eax
sub    %r8d,%eax
sub    %r8d,%eax
cltq
mov    %r9,0x10(%rdi,%rax,8)
movsd  (%rsi),%xmm1
movsd  0x8(%rsi),%xmm0
mulsd  0x18(%rdx),%xmm0
mulsd  (%rdx),%xmm1
addsd  %xmm0,%xmm1
addsd  0x30(%rdx),%xmm1
movapd 0xe5a(%rip),%xmm0
xorpd  %xmm0,%xmm1
movsd  0x28(%rsi),%xmm2
movsd  0x30(%rsi),%xmm3
mulsd  0x18(%rcx),%xmm3
mulsd  (%rcx),%xmm2
addsd  %xmm3,%xmm2
addsd  0x30(%rcx),%xmm2
mulsd  %xmm1,%xmm2
movsd  %xmm2,0x10(%rdi,%r10,8)
movsd  0x8(%rsi),%xmm1
mulsd  0x18(%rdx),%xmm1
movsd  (%rsi),%xmm2
mulsd  (%rdx),%xmm2
addsd  %xmm1,%xmm2
addsd  0x30(%rdx),%xmm2
movsd  0x28(%rsi),%xmm1
movsd  0x30(%rsi),%xmm3
mulsd  0x20(%rcx),%xmm3
mulsd  0x8(%rcx),%xmm1
xorpd  %xmm0,%xmm2
addsd  %xmm3,%xmm1
addsd  0x38(%rcx),%xmm1
mulsd  %xmm2,%xmm1
movsd  %xmm1,0x10(%rdi,%r11,8)
lea    (%r8,%r8,4),%eax
lea    (%r8,%rax,4),%eax
cltq
mov    %r9,0x18(%rdi,%rax,8)
movsd  (%rsi),%xmm1
movsd  0x8(%rsi),%xmm2
mulsd  0x20(%rdx),%xmm2
mulsd  0x8(%rdx),%xmm1
addsd  %xmm2,%xmm1
addsd  0x38(%rdx),%xmm1
movsd  %xmm1,0x18(%rdi,%r10,8)
lea    (%r8,%r8,8),%eax
lea    (%rax,%rax,2),%eax
add    %r8d,%eax
cltq
mov    %r9,0x20(%rdi,%rax,8)
movsd  0x8(%rsi),%xmm1
mulsd  0x20(%rdx),%xmm1
movsd  (%rsi),%xmm2
mulsd  0x8(%rdx),%xmm2
addsd  %xmm1,%xmm2
addsd  0x38(%rdx),%xmm2
movsd  %xmm2,0x20(%rdi,%r11,8)
imul   $0x23,%r8d,%eax
cltq
mov    %r9,0x28(%rdi,%rax,8)
movsd  (%rsi),%xmm1
movsd  0x8(%rsi),%xmm2
mulsd  0x20(%rdx),%xmm2
mulsd  0x8(%rdx),%xmm1
addsd  %xmm2,%xmm1
addsd  0x38(%rdx),%xmm1
xorpd  %xmm0,%xmm1
movsd  0x28(%rsi),%xmm2
movsd  0x30(%rsi),%xmm3
mulsd  0x18(%rcx),%xmm3
mulsd  (%rcx),%xmm2
addsd  %xmm3,%xmm2
addsd  0x30(%rcx),%xmm2
mulsd  %xmm1,%xmm2
movsd  %xmm2,0x28(%rdi,%r10,8)
movsd  0x8(%rsi),%xmm1
mulsd  0x20(%rdx),%xmm1
movsd  (%rsi),%xmm2
mulsd  0x8(%rdx),%xmm2
addsd  %xmm1,%xmm2
addsd  0x38(%rdx),%xmm2
movsd  0x28(%rsi),%xmm1
movsd  0x30(%rsi),%xmm3
mulsd  0x20(%rcx),%xmm3
mulsd  0x8(%rcx),%xmm1
xorpd  %xmm0,%xmm2
addsd  %xmm3,%xmm1
addsd  0x38(%rcx),%xmm1
mulsd  %xmm2,%xmm1
movsd  %xmm1,0x28(%rdi,%r11,8)
imul   $0x38,%r8d,%eax
or     $0x6,%eax
cltq
mov    %r9,(%rdi,%rax,8)
movsd  0x28(%rsi),%xmm1
movsd  0x30(%rsi),%xmm2
mulsd  0x18(%rcx),%xmm2
mulsd  (%rcx),%xmm1
addsd  %xmm2,%xmm1
addsd  0x30(%rcx),%xmm1
xorpd  %xmm0,%xmm1
movlpd %xmm1,0x30(%rdi,%r10,8)
movsd  0x30(%rsi),%xmm1
mulsd  0x20(%rcx),%xmm1
movsd  0x28(%rsi),%xmm2
mulsd  0x8(%rcx),%xmm2
addsd  %xmm1,%xmm2
addsd  0x38(%rcx),%xmm2
xorpd  %xmm0,%xmm2
movlpd %xmm2,0x30(%rdi,%r11,8)
ret

LLVM IR (-O0):
define void @ZuN(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
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
%11 = load double*, double** %6
%12 = load i32, i32* %10
%13 = mul nsw i32 0, %12
%14 = mul nsw i32 %13, 7
%15 = add nsw i32 0, %14
%16 = sext i32 %15 to i64
%17 = getelementptr inbounds double, double* %11, i64 %16
store double -1.000000e+00, double* %17
%18 = load double*, double** %7
%19 = getelementptr inbounds double, double* %18, i64 0
%20 = load double, double* %19
%21 = load double*, double** %8
%22 = getelementptr inbounds double, double* %21, i64 0
%23 = load double, double* %22
%24 = load double*, double** %7
%25 = getelementptr inbounds double, double* %24, i64 1
%26 = load double, double* %25
%27 = load double*, double** %8
%28 = getelementptr inbounds double, double* %27, i64 3
%29 = load double, double* %28
%30 = fmul double %26, %29
%31 = call double @llvm.fmuladd.f64(double %20, double %23, double %30)
%32 = load double*, double** %8
%33 = getelementptr inbounds double, double* %32, i64 6
%34 = load double, double* %33
%35 = fadd double %31, %34
%36 = load double*, double** %6
%37 = load i32, i32* %10
%38 = mul nsw i32 6, %37
%39 = mul nsw i32 %38, 7
%40 = add nsw i32 0, %39
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds double, double* %36, i64 %41
store double %35, double* %42
%43 = load double*, double** %6
%44 = load i32, i32* %10
%45 = mul nsw i32 1, %44
%46 = mul nsw i32 %45, 7
%47 = add nsw i32 1, %46
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds double, double* %43, i64 %48
store double -1.000000e+00, double* %49
%50 = load double*, double** %7
%51 = getelementptr inbounds double, double* %50, i64 0
%52 = load double, double* %51
%53 = load double*, double** %8
%54 = getelementptr inbounds double, double* %53, i64 0
%55 = load double, double* %54
%56 = load double*, double** %7
%57 = getelementptr inbounds double, double* %56, i64 1
%58 = load double, double* %57
%59 = load double*, double** %8
%60 = getelementptr inbounds double, double* %59, i64 3
%61 = load double, double* %60
%62 = fmul double %58, %61
%63 = call double @llvm.fmuladd.f64(double %52, double %55, double %62)
%64 = load double*, double** %8
%65 = getelementptr inbounds double, double* %64, i64 6
%66 = load double, double* %65
%67 = fadd double %63, %66
%68 = load double*, double** %6
%69 = load i32, i32* %10
%70 = mul nsw i32 7, %69
%71 = mul nsw i32 %70, 7
%72 = add nsw i32 1, %71
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %68, i64 %73
store double %67, double* %74
%75 = load double*, double** %6
%76 = load i32, i32* %10
%77 = mul nsw i32 2, %76
%78 = mul nsw i32 %77, 7
%79 = add nsw i32 2, %78
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds double, double* %75, i64 %80
store double -1.000000e+00, double* %81
%82 = load double*, double** %7
%83 = getelementptr inbounds double, double* %82, i64 0
%84 = load double, double* %83
%85 = load double*, double** %8
%86 = getelementptr inbounds double, double* %85, i64 0
%87 = load double, double* %86
%88 = load double*, double** %7
%89 = getelementptr inbounds double, double* %88, i64 1
%90 = load double, double* %89
%91 = load double*, double** %8
%92 = getelementptr inbounds double, double* %91, i64 3
%93 = load double, double* %92
%94 = fmul double %90, %93
%95 = call double @llvm.fmuladd.f64(double %84, double %87, double %94)
%96 = load double*, double** %8
%97 = getelementptr inbounds double, double* %96, i64 6
%98 = load double, double* %97
%99 = fadd double %95, %98
%100 = fneg double %99
%101 = load double*, double** %7
%102 = getelementptr inbounds double, double* %101, i64 5
%103 = load double, double* %102
%104 = load double*, double** %9
%105 = getelementptr inbounds double, double* %104, i64 0
%106 = load double, double* %105
%107 = load double*, double** %7
%108 = getelementptr inbounds double, double* %107, i64 6
%109 = load double, double* %108
%110 = load double*, double** %9
%111 = getelementptr inbounds double, double* %110, i64 3
%112 = load double, double* %111
%113 = fmul double %109, %112
%114 = call double @llvm.fmuladd.f64(double %103, double %106, double %113)
%115 = load double*, double** %9
%116 = getelementptr inbounds double, double* %115, i64 6
%117 = load double, double* %116
%118 = fadd double %114, %117
%119 = fmul double %100, %118
%120 = load double*, double** %6
%121 = load i32, i32* %10
%122 = mul nsw i32 6, %121
%123 = mul nsw i32 %122, 7
%124 = add nsw i32 2, %123
%125 = sext i32 %124 to i64
%126 = getelementptr inbounds double, double* %120, i64 %125
store double %119, double* %126
%127 = load double*, double** %7
%128 = getelementptr inbounds double, double* %127, i64 0
%129 = load double, double* %128
%130 = load double*, double** %8
%131 = getelementptr inbounds double, double* %130, i64 0
%132 = load double, double* %131
%133 = load double*, double** %7
%134 = getelementptr inbounds double, double* %133, i64 1
%135 = load double, double* %134
%136 = load double*, double** %8
%137 = getelementptr inbounds double, double* %136, i64 3
%138 = load double, double* %137
%139 = fmul double %135, %138
%140 = call double @llvm.fmuladd.f64(double %129, double %132, double %139)
%141 = load double*, double** %8
%142 = getelementptr inbounds double, double* %141, i64 6
%143 = load double, double* %142
%144 = fadd double %140, %143
%145 = fneg double %144
%146 = load double*, double** %7
%147 = getelementptr inbounds double, double* %146, i64 5
%148 = load double, double* %147
%149 = load double*, double** %9
%150 = getelementptr inbounds double, double* %149, i64 1
%151 = load double, double* %150
%152 = load double*, double** %7
%153 = getelementptr inbounds double, double* %152, i64 6
%154 = load double, double* %153
%155 = load double*, double** %9
%156 = getelementptr inbounds double, double* %155, i64 4
%157 = load double, double* %156
%158 = fmul double %154, %157
%159 = call double @llvm.fmuladd.f64(double %148, double %151, double %158)
%160 = load double*, double** %9
%161 = getelementptr inbounds double, double* %160, i64 7
%162 = load double, double* %161
%163 = fadd double %159, %162
%164 = fmul double %145, %163
%165 = load double*, double** %6
%166 = load i32, i32* %10
%167 = mul nsw i32 7, %166
%168 = mul nsw i32 %167, 7
%169 = add nsw i32 2, %168
%170 = sext i32 %169 to i64
%171 = getelementptr inbounds double, double* %165, i64 %170
store double %164, double* %171
%172 = load double*, double** %6
%173 = load i32, i32* %10
%174 = mul nsw i32 3, %173
%175 = mul nsw i32 %174, 7
%176 = add nsw i32 3, %175
%177 = sext i32 %176 to i64
%178 = getelementptr inbounds double, double* %172, i64 %177
store double -1.000000e+00, double* %178
%179 = load double*, double** %7
%180 = getelementptr inbounds double, double* %179, i64 0
%181 = load double, double* %180
%182 = load double*, double** %8
%183 = getelementptr inbounds double, double* %182, i64 1
%184 = load double, double* %183
%185 = load double*, double** %7
%186 = getelementptr inbounds double, double* %185, i64 1
%187 = load double, double* %186
%188 = load double*, double** %8
%189 = getelementptr inbounds double, double* %188, i64 4
%190 = load double, double* %189
%191 = fmul double %187, %190
%192 = call double @llvm.fmuladd.f64(double %181, double %184, double %191)
%193 = load double*, double** %8
%194 = getelementptr inbounds double, double* %193, i64 7
%195 = load double, double* %194
%196 = fadd double %192, %195
%197 = load double*, double** %6
%198 = load i32, i32* %10
%199 = mul nsw i32 6, %198
%200 = mul nsw i32 %199, 7
%201 = add nsw i32 3, %200
%202 = sext i32 %201 to i64
%203 = getelementptr inbounds double, double* %197, i64 %202
store double %196, double* %203
%204 = load double*, double** %6
%205 = load i32, i32* %10
%206 = mul nsw i32 4, %205
%207 = mul nsw i32 %206, 7
%208 = add nsw i32 4, %207
%209 = sext i32 %208 to i64
%210 = getelementptr inbounds double, double* %204, i64 %209
store double -1.000000e+00, double* %210
%211 = load double*, double** %7
%212 = getelementptr inbounds double, double* %211, i64 0
%213 = load double, double* %212
%214 = load double*, double** %8
%215 = getelementptr inbounds double, double* %214, i64 1
%216 = load double, double* %215
%217 = load double*, double** %7
%218 = getelementptr inbounds double, double* %217, i64 1
%219 = load double, double* %218
%220 = load double*, double** %8
%221 = getelementptr inbounds double, double* %220, i64 4
%222 = load double, double* %221
%223 = fmul double %219, %222
%224 = call double @llvm.fmuladd.f64(double %213, double %216, double %223)
%225 = load double*, double** %8
%226 = getelementptr inbounds double, double* %225, i64 7
%227 = load double, double* %226
%228 = fadd double %224, %227
%229 = load double*, double** %6
%230 = load i32, i32* %10
%231 = mul nsw i32 7, %230
%232 = mul nsw i32 %231, 7
%233 = add nsw i32 4, %232
%234 = sext i32 %233 to i64
%235 = getelementptr inbounds double, double* %229, i64 %234
store double %228, double* %235
%236 = load double*, double** %6
%237 = load i32, i32* %10
%238 = mul nsw i32 5, %237
%239 = mul nsw i32 %238, 7
%240 = add nsw i32 5, %239
%241 = sext i32 %240 to i64
%242 = getelementptr inbounds double, double* %236, i64 %241
store double -1.000000e+00, double* %242
%243 = load double*, double** %7
%244 = getelementptr inbounds double, double* %243, i64 0
%245 = load double, double* %244
%246 = load double*, double** %8
%247 = getelementptr inbounds double, double* %246, i64 1
%248 = load double, double* %247
%249 = load double*, double** %7
%250 = getelementptr inbounds double, double* %249, i64 1
%251 = load double, double* %250
%252 = load double*, double** %8
%253 = getelementptr inbounds double, double* %252, i64 4
%254 = load double, double* %253
%255 = fmul double %251, %254
%256 = call double @llvm.fmuladd.f64(double %245, double %248, double %255)
%257 = load double*, double** %8
%258 = getelementptr inbounds double, double* %257, i64 7
%259 = load double, double* %258
%260 = fadd double %256, %259
%261 = fneg double %260
%262 = load double*, double** %7
%263 = getelementptr inbounds double, double* %262, i64 5
%264 = load double, double* %263
%265 = load double*, double** %9
%266 = getelementptr inbounds double, double* %265, i64 0
%267 = load double, double* %266
%268 = load double*, double** %7
%269 = getelementptr inbounds double, double* %268, i64 6
%270 = load double, double* %269
%271 = load double*, double** %9
%272 = getelementptr inbounds double, double* %271, i64 3
%273 = load double, double* %272
%274 = fmul double %270, %273
%275 = call double @llvm.fmuladd.f64(double %264, double %267, double %274)
%276 = load double*, double** %9
%277 = getelementptr inbounds double, double* %276, i64 6
%278 = load double, double* %277
%279 = fadd double %275, %278
%280 = fmul double %261, %279
%281 = load double*, double** %6
%282 = load i32, i32* %10
%283 = mul nsw i32 6, %282
%284 = mul nsw i32 %283, 7
%285 = add nsw i32 5, %284
%286 = sext i32 %285 to i64
%287 = getelementptr inbounds double, double* %281, i64 %286
store double %280, double* %287
%288 = load double*, double** %7
%289 = getelementptr inbounds double, double* %288, i64 0
%290 = load double, double* %289
%291 = load double*, double** %8
%292 = getelementptr inbounds double, double* %291, i64 1
%293 = load double, double* %292
%294 = load double*, double** %7
%295 = getelementptr inbounds double, double* %294, i64 1
%296 = load double, double* %295
%297 = load double*, double** %8
%298 = getelementptr inbounds double, double* %297, i64 4
%299 = load double, double* %298
%300 = fmul double %296, %299
%301 = call double @llvm.fmuladd.f64(double %290, double %293, double %300)
%302 = load double*, double** %8
%303 = getelementptr inbounds double, double* %302, i64 7
%304 = load double, double* %303
%305 = fadd double %301, %304
%306 = fneg double %305
%307 = load double*, double** %7
%308 = getelementptr inbounds double, double* %307, i64 5
%309 = load double, double* %308
%310 = load double*, double** %9
%311 = getelementptr inbounds double, double* %310, i64 1
%312 = load double, double* %311
%313 = load double*, double** %7
%314 = getelementptr inbounds double, double* %313, i64 6
%315 = load double, double* %314
%316 = load double*, double** %9
%317 = getelementptr inbounds double, double* %316, i64 4
%318 = load double, double* %317
%319 = fmul double %315, %318
%320 = call double @llvm.fmuladd.f64(double %309, double %312, double %319)
%321 = load double*, double** %9
%322 = getelementptr inbounds double, double* %321, i64 7
%323 = load double, double* %322
%324 = fadd double %320, %323
%325 = fmul double %306, %324
%326 = load double*, double** %6
%327 = load i32, i32* %10
%328 = mul nsw i32 7, %327
%329 = mul nsw i32 %328, 7
%330 = add nsw i32 5, %329
%331 = sext i32 %330 to i64
%332 = getelementptr inbounds double, double* %326, i64 %331
store double %325, double* %332
%333 = load double*, double** %6
%334 = load i32, i32* %10
%335 = mul nsw i32 8, %334
%336 = mul nsw i32 %335, 7
%337 = add nsw i32 6, %336
%338 = sext i32 %337 to i64
%339 = getelementptr inbounds double, double* %333, i64 %338
store double -1.000000e+00, double* %339
%340 = load double*, double** %7
%341 = getelementptr inbounds double, double* %340, i64 5
%342 = load double, double* %341
%343 = load double*, double** %9
%344 = getelementptr inbounds double, double* %343, i64 0
%345 = load double, double* %344
%346 = load double*, double** %7
%347 = getelementptr inbounds double, double* %346, i64 6
%348 = load double, double* %347
%349 = load double*, double** %9
%350 = getelementptr inbounds double, double* %349, i64 3
%351 = load double, double* %350
%352 = fmul double %348, %351
%353 = call double @llvm.fmuladd.f64(double %342, double %345, double %352)
%354 = load double*, double** %9
%355 = getelementptr inbounds double, double* %354, i64 6
%356 = load double, double* %355
%357 = fadd double %353, %356
%358 = fneg double %357
%359 = load double*, double** %6
%360 = load i32, i32* %10
%361 = mul nsw i32 6, %360
%362 = mul nsw i32 %361, 7
%363 = add nsw i32 6, %362
%364 = sext i32 %363 to i64
%365 = getelementptr inbounds double, double* %359, i64 %364
store double %358, double* %365
%366 = load double*, double** %7
%367 = getelementptr inbounds double, double* %366, i64 5
%368 = load double, double* %367
%369 = load double*, double** %9
%370 = getelementptr inbounds double, double* %369, i64 1
%371 = load double, double* %370
%372 = load double*, double** %7
%373 = getelementptr inbounds double, double* %372, i64 6
%374 = load double, double* %373
%375 = load double*, double** %9
%376 = getelementptr inbounds double, double* %375, i64 4
%377 = load double, double* %376
%378 = fmul double %374, %377
%379 = call double @llvm.fmuladd.f64(double %368, double %371, double %378)
%380 = load double*, double** %9
%381 = getelementptr inbounds double, double* %380, i64 7
%382 = load double, double* %381
%383 = fadd double %379, %382
%384 = fneg double %383
%385 = load double*, double** %6
%386 = load i32, i32* %10
%387 = mul nsw i32 7, %386
%388 = mul nsw i32 %387, 7
%389 = add nsw i32 6, %388
%390 = sext i32 %389 to i64
%391 = getelementptr inbounds double, double* %385, i64 %390
store double %384, double* %391
ret void
}

LLVM IR (-O2):
define void @ZuN(double* nocapture writeonly %0, double* nocapture %1, double* nocapture %2, double* nocapture %3, i32 %4) {
store double -1.000000e+00, double* %0
%6 = load double, double* %1
%7 = load double, double* %2
%8 = getelementptr inbounds double, double* %1, i64 1
%9 = load double, double* %8
%10 = getelementptr inbounds double, double* %2, i64 3
%11 = load double, double* %10
%12 = fmul double %9, %11
%13 = tail call double @llvm.fmuladd.f64(double %6, double %7, double %12)
%14 = getelementptr inbounds double, double* %2, i64 6
%15 = load double, double* %14
%16 = fadd double %15, %13
%17 = mul i32 %4, 42
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds double, double* %0, i64 %18
store double %16, double* %19
%20 = mul nsw i32 %4, 7
%21 = add nsw i32 %20, 1
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds double, double* %0, i64 %22
store double -1.000000e+00, double* %23
%24 = load double, double* %1
%25 = load double, double* %2
%26 = load double, double* %8
%27 = load double, double* %10
%28 = fmul double %26, %27
%29 = tail call double @llvm.fmuladd.f64(double %24, double %25, double %28)
%30 = load double, double* %14
%31 = fadd double %30, %29
%32 = mul i32 %4, 49
%33 = add nsw i32 %32, 1
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds double, double* %0, i64 %34
store double %31, double* %35
%36 = mul nsw i32 %4, 14
%37 = add nsw i32 %36, 2
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds double, double* %0, i64 %38
store double -1.000000e+00, double* %39
%40 = load double, double* %1
%41 = load double, double* %2
%42 = load double, double* %8
%43 = load double, double* %10
%44 = fmul double %42, %43
%45 = tail call double @llvm.fmuladd.f64(double %40, double %41, double %44)
%46 = load double, double* %14
%47 = fadd double %46, %45
%48 = fneg double %47
%49 = getelementptr inbounds double, double* %1, i64 5
%50 = load double, double* %49
%51 = load double, double* %3
%52 = getelementptr inbounds double, double* %1, i64 6
%53 = load double, double* %52
%54 = getelementptr inbounds double, double* %3, i64 3
%55 = load double, double* %54
%56 = fmul double %53, %55
%57 = tail call double @llvm.fmuladd.f64(double %50, double %51, double %56)
%58 = getelementptr inbounds double, double* %3, i64 6
%59 = load double, double* %58
%60 = fadd double %59, %57
%61 = fmul double %60, %48
%62 = add nsw i32 %17, 2
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds double, double* %0, i64 %63
store double %61, double* %64
%65 = load double, double* %1
%66 = load double, double* %2
%67 = load double, double* %8
%68 = load double, double* %10
%69 = fmul double %67, %68
%70 = tail call double @llvm.fmuladd.f64(double %65, double %66, double %69)
%71 = load double, double* %14
%72 = fadd double %71, %70
%73 = fneg double %72
%74 = load double, double* %49
%75 = getelementptr inbounds double, double* %3, i64 1
%76 = load double, double* %75
%77 = load double, double* %52
%78 = getelementptr inbounds double, double* %3, i64 4
%79 = load double, double* %78
%80 = fmul double %77, %79
%81 = tail call double @llvm.fmuladd.f64(double %74, double %76, double %80)
%82 = getelementptr inbounds double, double* %3, i64 7
%83 = load double, double* %82
%84 = fadd double %83, %81
%85 = fmul double %84, %73
%86 = add nsw i32 %32, 2
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds double, double* %0, i64 %87
store double %85, double* %88
%89 = mul i32 %4, 21
%90 = add nsw i32 %89, 3
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds double, double* %0, i64 %91
store double -1.000000e+00, double* %92
%93 = load double, double* %1
%94 = getelementptr inbounds double, double* %2, i64 1
%95 = load double, double* %94
%96 = load double, double* %8
%97 = getelementptr inbounds double, double* %2, i64 4
%98 = load double, double* %97
%99 = fmul double %96, %98
%100 = tail call double @llvm.fmuladd.f64(double %93, double %95, double %99)
%101 = getelementptr inbounds double, double* %2, i64 7
%102 = load double, double* %101
%103 = fadd double %102, %100
%104 = add nsw i32 %17, 3
%105 = sext i32 %104 to i64
%106 = getelementptr inbounds double, double* %0, i64 %105
store double %103, double* %106
%107 = mul nsw i32 %4, 28
%108 = add nsw i32 %107, 4
%109 = sext i32 %108 to i64
%110 = getelementptr inbounds double, double* %0, i64 %109
store double -1.000000e+00, double* %110
%111 = load double, double* %1
%112 = load double, double* %94
%113 = load double, double* %8
%114 = load double, double* %97
%115 = fmul double %113, %114
%116 = tail call double @llvm.fmuladd.f64(double %111, double %112, double %115)
%117 = load double, double* %101
%118 = fadd double %117, %116
%119 = add nsw i32 %32, 4
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds double, double* %0, i64 %120
store double %118, double* %121
%122 = mul i32 %4, 35
%123 = add nsw i32 %122, 5
%124 = sext i32 %123 to i64
%125 = getelementptr inbounds double, double* %0, i64 %124
store double -1.000000e+00, double* %125
%126 = load double, double* %1
%127 = load double, double* %94
%128 = load double, double* %8
%129 = load double, double* %97
%130 = fmul double %128, %129
%131 = tail call double @llvm.fmuladd.f64(double %126, double %127, double %130)
%132 = load double, double* %101
%133 = fadd double %132, %131
%134 = fneg double %133
%135 = load double, double* %49
%136 = load double, double* %3
%137 = load double, double* %52
%138 = load double, double* %54
%139 = fmul double %137, %138
%140 = tail call double @llvm.fmuladd.f64(double %135, double %136, double %139)
%141 = load double, double* %58
%142 = fadd double %141, %140
%143 = fmul double %142, %134
%144 = add nsw i32 %17, 5
%145 = sext i32 %144 to i64
%146 = getelementptr inbounds double, double* %0, i64 %145
store double %143, double* %146
%147 = load double, double* %1
%148 = load double, double* %94
%149 = load double, double* %8
%150 = load double, double* %97
%151 = fmul double %149, %150
%152 = tail call double @llvm.fmuladd.f64(double %147, double %148, double %151)
%153 = load double, double* %101
%154 = fadd double %153, %152
%155 = fneg double %154
%156 = load double, double* %49
%157 = load double, double* %75
%158 = load double, double* %52
%159 = load double, double* %78
%160 = fmul double %158, %159
%161 = tail call double @llvm.fmuladd.f64(double %156, double %157, double %160)
%162 = load double, double* %82
%163 = fadd double %162, %161
%164 = fmul double %163, %155
%165 = add nsw i32 %32, 5
%166 = sext i32 %165 to i64
%167 = getelementptr inbounds double, double* %0, i64 %166
store double %164, double* %167
%168 = mul nsw i32 %4, 56
%169 = or i32 %168, 6
%170 = sext i32 %169 to i64
%171 = getelementptr inbounds double, double* %0, i64 %170
store double -1.000000e+00, double* %171
%172 = load double, double* %49
%173 = load double, double* %3
%174 = load double, double* %52
%175 = load double, double* %54
%176 = fmul double %174, %175
%177 = tail call double @llvm.fmuladd.f64(double %172, double %173, double %176)
%178 = load double, double* %58
%179 = fadd double %178, %177
%180 = fneg double %179
%181 = add nsw i32 %17, 6
%182 = sext i32 %181 to i64
%183 = getelementptr inbounds double, double* %0, i64 %182
store double %180, double* %183
%184 = load double, double* %49
%185 = load double, double* %75
%186 = load double, double* %52
%187 = load double, double* %78
%188 = fmul double %186, %187
%189 = tail call double @llvm.fmuladd.f64(double %184, double %185, double %188)
%190 = load double, double* %82
%191 = fadd double %190, %189
%192 = fneg double %191
%193 = add nsw i32 %32, 6
%194 = sext i32 %193 to i64
%195 = getelementptr inbounds double, double* %0, i64 %194
store double %192, double* %195
ret void
}