void LoadXY2Matrix(double *mat, int n, double x, double y) {
  int i, N, nm1;
  double ax;
  N = n + 1;
  nm1 = (n - 1) * N;
  for (i = n - 1, ax = 1.; i >= 0; i--, ax *= x) {
    mat[nm1 + i] += ax;
    mat[i * N + n] += ax * y;
  }
  for (i = n - 2; i >= 0; i--, ax *= x)
    mat[i * N] += ax;
}


x86 assembly (clang, -O0):
<LoadXY2Matrix>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
imul   -0x28(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
movsd  0xec4(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
cmpl   $0x0,-0x24(%rbp)
jl     11b1 <LoadXY2Matrix+0xb1>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x28(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
movsd  -0x18(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    1141 <LoadXY2Matrix+0x41>
mov    -0xc(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jl     11fe <LoadXY2Matrix+0xfe>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x28(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
movsd  -0x18(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    11ba <LoadXY2Matrix+0xba>
pop    %rbp
ret

x86 assembly (clang, -O2):
<LoadXY2Matrix>:
test   %esi,%esi
jle    120d <LoadXY2Matrix+0x10d>
lea    0x1(%rsi),%r8d
lea    -0x1(%rsi),%eax
imul   %r8d,%eax
mov    %esi,%r10d
cltq
mov    %r8d,%r9d
lea    0x1(%r10),%rdx
lea    (%rdi,%rax,8),%rax
lea    -0x1(%r10),%rcx
imul   %r9,%rcx
add    %r10,%rcx
lea    (%rdi,%rcx,8),%rcx
shl    $0x3,%r9
neg    %r9
movsd  0xebf(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  -0x10(%rax,%rdx,8),%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,-0x10(%rax,%rdx,8)
movapd %xmm2,%xmm3
mulsd  %xmm1,%xmm3
addsd  (%rcx),%xmm3
movsd  %xmm3,(%rcx)
mulsd  %xmm0,%xmm2
add    $0xffffffffffffffff,%rdx
add    %r9,%rcx
cmp    $0x1,%rdx
ja     1150 <LoadXY2Matrix+0x50>
cmp    $0x2,%esi
jl     120d <LoadXY2Matrix+0x10d>
add    $0xfffffffe,%esi
movslq %r8d,%r8
mov    %rsi,%rax
test   $0x1,%sil
jne    11b3 <LoadXY2Matrix+0xb3>
imul   %r8,%rax
movsd  (%rdi,%rax,8),%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,(%rdi,%rax,8)
lea    -0x1(%rsi),%rax
mulsd  %xmm0,%xmm2
test   %rsi,%rsi
je     120d <LoadXY2Matrix+0x10d>
mov    %rax,%rcx
imul   %r8,%rcx
mov    %r8,%rsi
shl    $0x4,%rsi
neg    %rsi
lea    0x1(%rax),%rdx
add    $0xffffffffffffffff,%rax
imul   %r8,%rax
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rdi,%rcx,8),%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,(%rdi,%rcx,8)
mulsd  %xmm0,%xmm2
movsd  (%rdi,%rax,8),%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,(%rdi,%rax,8)
mulsd  %xmm0,%xmm2
add    %rsi,%rdi
add    $0xfffffffffffffffe,%rdx
jne    11e0 <LoadXY2Matrix+0xe0>
ret

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

LLVM IR (-O2):
define void @LoadXY2Matrix(double* nocapture %0, i32 %1, double %2, double %3) {
%5 = add nsw i32 %1, 1
%6 = icmp sgt i32 %1, 0
br i1 %6, label %7, label %63
7: ; preds = %4
%8 = add nsw i32 %1, -1
%9 = mul nsw i32 %8, %5
%10 = zext i32 %1 to i64
%11 = sext i32 %9 to i64
%12 = zext i32 %5 to i64
br label %13
13: ; preds = %7, %13
%14 = phi i64 [ %10, %7 ], [ %16, %13 ]
%15 = phi double [ 1.000000e+00, %7 ], [ %26, %13 ]
%16 = add nsw i64 %14, -1
%17 = add nsw i64 %16, %11
%18 = getelementptr inbounds double, double* %0, i64 %17
%19 = load double, double* %18
%20 = fadd double %15, %19
store double %20, double* %18
%21 = mul nsw i64 %16, %12
%22 = add nsw i64 %21, %10
%23 = getelementptr inbounds double, double* %0, i64 %22
%24 = load double, double* %23
%25 = tail call double @llvm.fmuladd.f64(double %15, double %3, double %24)
store double %25, double* %23
%26 = fmul double %15, %2
%27 = icmp ugt i64 %14, 1
br i1 %27, label %13, label %28
28: ; preds = %13
%29 = icmp sgt i32 %1, 1
br i1 %29, label %30, label %63
30: ; preds = %28
%31 = add nsw i32 %1, -2
%32 = zext i32 %31 to i64
%33 = sext i32 %5 to i64
%34 = and i64 %32, 1
%35 = icmp eq i64 %34, 0
br i1 %35, label %36, label %43
36: ; preds = %30
%37 = mul nsw i64 %32, %33
%38 = getelementptr inbounds double, double* %0, i64 %37
%39 = load double, double* %38
%40 = fadd double %26, %39
store double %40, double* %38
%41 = add nsw i64 %32, -1
%42 = fmul double %26, %2
br label %43
43: ; preds = %36, %30
%44 = phi i64 [ %32, %30 ], [ %41, %36 ]
%45 = phi double [ %26, %30 ], [ %42, %36 ]
%46 = icmp eq i32 %31, 0
br i1 %46, label %63, label %47
47: ; preds = %43, %47
%48 = phi i64 [ %60, %47 ], [ %44, %43 ]
%49 = phi double [ %61, %47 ], [ %45, %43 ]
%50 = mul nsw i64 %48, %33
%51 = getelementptr inbounds double, double* %0, i64 %50
%52 = load double, double* %51
%53 = fadd double %49, %52
store double %53, double* %51
%54 = add nsw i64 %48, -1
%55 = fmul double %49, %2
%56 = mul nsw i64 %54, %33
%57 = getelementptr inbounds double, double* %0, i64 %56
%58 = load double, double* %57
%59 = fadd double %55, %58
store double %59, double* %57
%60 = add nsw i64 %48, -2
%61 = fmul double %55, %2
%62 = icmp eq i64 %54, 0
br i1 %62, label %63, label %47
63: ; preds = %43, %47, %4, %28
ret void
}