void get_face_2d_utopia(int iX, int iY, int nX, int nY, int X_nbr, int Y_nbr,
                        int nrm, int tid, double delta_t, double u_face,
                        double u_face_t, double *sc_face, double *sc_old) {
  int nbrid;
  int offset = 1;
  int iwind_t = 1;
  double tmp;
  double tmp_2;
  double sc_face_tmp;
  double wind = 1.;
  double wind_t = 1.;
  double crnt;
  double crnt_t;
  crnt = u_face * delta_t * 1.;
  if (crnt * nrm > 0.)
    offset = 0;
  nbrid = (iY + (Y_nbr)) * nX + iX + (X_nbr);
  sc_face_tmp = 0.5 * (*(sc_old + nbrid) + *(sc_old + tid) -
                       crnt * nrm * (*(sc_old + nbrid) - *(sc_old + tid)));
  tmp = (1. - crnt * crnt) * 0.1666666666666666;
  nbrid = (iY + (offset)*Y_nbr) * nX + iX + (offset)*X_nbr;
  sc_face_tmp += *(sc_old + nbrid) * 2. * tmp;
  nbrid = (iY + (offset + nrm) * Y_nbr) * nX + iX + (offset + nrm) * X_nbr;
  sc_face_tmp -= *(sc_old + nbrid) * tmp;
  nbrid = (iY + (offset - nrm) * Y_nbr) * nX + iX + (offset - nrm) * X_nbr;
  sc_face_tmp -= *(sc_old + nbrid) * tmp;
  nbrid = (iY + (offset)*Y_nbr) * nX + iX + (offset)*X_nbr;
  sc_face_tmp -= *(sc_old + nbrid) * 0.5 * u_face_t * delta_t * 1. * u_face_t *
                 delta_t * 1.;
  nbrid = (iY + (offset)*Y_nbr + (1 - Y_nbr * nrm)) * nX + iX + (offset)*X_nbr +
          (1 - X_nbr * nrm);
  sc_face_tmp += *(sc_old + nbrid) * u_face_t * delta_t * 1. *
                 (-0.25 + 0.25 * u_face_t * delta_t * 1.);
  nbrid = (iY + (offset)*Y_nbr - (1 - Y_nbr * nrm)) * nX + iX + (offset)*X_nbr -
          (1 - X_nbr * nrm);
  sc_face_tmp += *(sc_old + nbrid) * u_face_t * delta_t * 1. *
                 (0.25 + 0.25 * u_face_t * delta_t * 1.);
  if (crnt < 0.)
    wind = -1.;
  crnt_t = u_face_t * delta_t * 1.;
  if (crnt_t < 0.) {
    wind_t = -1.;
    iwind_t = -1;
  }
  tmp_2 = 0.25 * crnt_t * (wind - crnt);
  sc_face_tmp += *(sc_old + tid) * wind_t * nrm * tmp_2;
  nbrid = (iY + (Y_nbr)) * nX + iX + (X_nbr);
  sc_face_tmp -= *(sc_old + nbrid) * wind_t * nrm * tmp_2;
  nbrid = (iY - iwind_t * (1 - Y_nbr * nrm)) * nX + iX -
          iwind_t * (1 - X_nbr * nrm);
  sc_face_tmp -= *(sc_old + nbrid) * wind_t * nrm * tmp_2;
  nbrid = (iY + (Y_nbr)-iwind_t * (1 - Y_nbr * nrm)) * nX + iX +
          (X_nbr)-iwind_t * (1 - X_nbr * nrm);
  sc_face_tmp += *(sc_old + nbrid) * wind_t * nrm * tmp_2;
  *sc_face = sc_face_tmp;
}


x86 assembly (clang, -O0):
<get_face_2d_utopia>:
push   %rbp
mov    %rsp,%rbp
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    %r8d,-0x14(%rbp)
mov    %r9d,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
movsd  %xmm1,-0x28(%rbp)
movsd  %xmm2,-0x30(%rbp)
movl   $0x1,-0x38(%rbp)
movl   $0x1,-0x3c(%rbp)
movsd  0xeb5(%rip),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  0xea8(%rip),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x28(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  0xe91(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
cvtsi2sdl 0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    119a <get_face_2d_utopia+0x9a>
movl   $0x0,-0x38(%rbp)
mov    -0x8(%rbp),%eax
add    -0x18(%rbp),%eax
imul   -0xc(%rbp),%eax
add    -0x4(%rbp),%eax
add    -0x14(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x28(%rbp),%rax
movslq 0x18(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  -0x70(%rbp),%xmm1
cvtsi2sdl 0x10(%rbp),%xmm2
mulsd  %xmm2,%xmm1
mov    0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x28(%rbp),%rax
movslq 0x18(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  0xe02(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x70(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xdb4(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  0xdd0(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x8(%rbp),%eax
mov    -0x38(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    %ecx,%eax
imul   -0xc(%rbp),%eax
add    -0x4(%rbp),%eax
mov    -0x38(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x34(%rbp)
mov    0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  0xd90(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movsd  -0x48(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%eax
mov    -0x38(%rbp),%ecx
add    0x10(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    %ecx,%eax
imul   -0xc(%rbp),%eax
add    -0x4(%rbp),%eax
mov    -0x38(%rbp),%ecx
add    0x10(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x34(%rbp)
mov    0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x48(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%eax
mov    -0x38(%rbp),%ecx
sub    0x10(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    %ecx,%eax
imul   -0xc(%rbp),%eax
add    -0x4(%rbp),%eax
mov    -0x38(%rbp),%ecx
sub    0x10(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x34(%rbp)
mov    0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x48(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%eax
mov    -0x38(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    %ecx,%eax
imul   -0xc(%rbp),%eax
add    -0x4(%rbp),%eax
mov    -0x38(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x34(%rbp)
mov    0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  0xc7d(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
mulsd  -0x30(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  0xc4e(%rip),%xmm1
mulsd  %xmm1,%xmm0
mulsd  -0x30(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%eax
mov    -0x38(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    %ecx,%eax
mov    -0x18(%rbp),%edx
imul   0x10(%rbp),%edx
mov    $0x1,%ecx
sub    %edx,%ecx
add    %ecx,%eax
imul   -0xc(%rbp),%eax
add    -0x4(%rbp),%eax
mov    -0x38(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    %ecx,%eax
mov    -0x14(%rbp),%edx
imul   0x10(%rbp),%edx
mov    $0x1,%ecx
sub    %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x34(%rbp)
mov    0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x30(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  0xbbd(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xbb9(%rip),%xmm2
mulsd  -0x30(%rbp),%xmm2
mulsd  -0x20(%rbp),%xmm2
movsd  0xbaf(%rip),%xmm1
addsd  %xmm1,%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%eax
mov    -0x38(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    %ecx,%eax
mov    -0x18(%rbp),%edx
imul   0x10(%rbp),%edx
mov    $0x1,%ecx
sub    %edx,%ecx
sub    %ecx,%eax
imul   -0xc(%rbp),%eax
add    -0x4(%rbp),%eax
mov    -0x38(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    %ecx,%eax
mov    -0x14(%rbp),%edx
imul   0x10(%rbp),%edx
mov    $0x1,%ecx
sub    %edx,%ecx
sub    %ecx,%eax
mov    %eax,-0x34(%rbp)
mov    0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x30(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  0xb2b(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xb27(%rip),%xmm2
mulsd  -0x30(%rbp),%xmm2
mulsd  -0x20(%rbp),%xmm2
movsd  0xb15(%rip),%xmm1
addsd  %xmm1,%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x70(%rbp),%xmm0
jbe    1524 <get_face_2d_utopia+0x424>
movsd  0xb11(%rip),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x30(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  0xaca(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x78(%rbp),%xmm0
jbe    1561 <get_face_2d_utopia+0x461>
movsd  0xadb(%rip),%xmm0
movsd  %xmm0,-0x68(%rbp)
movl   $0xffffffff,-0x3c(%rbp)
movsd  0xa9f(%rip),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm1
subsd  -0x70(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    0x28(%rbp),%rax
movslq 0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x68(%rbp),%xmm0
cvtsi2sdl 0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x50(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%eax
add    -0x18(%rbp),%eax
imul   -0xc(%rbp),%eax
add    -0x4(%rbp),%eax
add    -0x14(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x68(%rbp),%xmm0
cvtsi2sdl 0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x50(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%eax
mov    -0x3c(%rbp),%ecx
mov    -0x18(%rbp),%esi
imul   0x10(%rbp),%esi
mov    $0x1,%edx
sub    %esi,%edx
imul   %edx,%ecx
sub    %ecx,%eax
imul   -0xc(%rbp),%eax
add    -0x4(%rbp),%eax
mov    -0x3c(%rbp),%ecx
mov    -0x14(%rbp),%esi
imul   0x10(%rbp),%esi
mov    $0x1,%edx
sub    %esi,%edx
imul   %edx,%ecx
sub    %ecx,%eax
mov    %eax,-0x34(%rbp)
mov    0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x68(%rbp),%xmm0
cvtsi2sdl 0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x50(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x8(%rbp),%eax
add    -0x18(%rbp),%eax
mov    -0x3c(%rbp),%ecx
mov    -0x18(%rbp),%esi
imul   0x10(%rbp),%esi
mov    $0x1,%edx
sub    %esi,%edx
imul   %edx,%ecx
sub    %ecx,%eax
imul   -0xc(%rbp),%eax
add    -0x4(%rbp),%eax
add    -0x14(%rbp),%eax
mov    -0x3c(%rbp),%ecx
mov    -0x14(%rbp),%esi
imul   0x10(%rbp),%esi
mov    $0x1,%edx
sub    %esi,%edx
imul   %edx,%ecx
sub    %ecx,%eax
mov    %eax,-0x34(%rbp)
mov    0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x68(%rbp),%xmm0
cvtsi2sdl 0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x50(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
mov    0x20(%rbp),%rax
movsd  %xmm0,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<get_face_2d_utopia>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %edx,%r10d
mov    0x38(%rsp),%r12d
mulsd  %xmm0,%xmm1
cvtsi2sd %r12d,%xmm9
movapd %xmm1,%xmm10
mulsd  %xmm9,%xmm10
lea    (%r9,%rsi,1),%r14d
mov    %r14d,%r15d
imul   %edx,%r15d
lea    (%r8,%rdi,1),%eax
mov    %eax,-0x4(%rsp)
add    %eax,%r15d
xor    %edx,%edx
xorpd  %xmm8,%xmm8
xor    %ebx,%ebx
ucomisd %xmm8,%xmm10
mov    $0x0,%eax
cmovbe %r9d,%eax
mov    0x50(%rsp),%r11
cmovbe %r8d,%edx
setbe  %bl
add    %esi,%eax
mov    %eax,%ebp
imul   %r10d,%ebp
add    %edi,%edx
add    %edx,%ebp
movslq %ebp,%rbp
movsd  (%r11,%rbp,8),%xmm15
lea    (%rbx,%r12,1),%r13d
mov    %r13d,%ebp
imul   %r9d,%ebp
add    %esi,%ebp
imul   %r10d,%ebp
imul   %r8d,%r13d
add    %edi,%r13d
add    %ebp,%r13d
sub    %r12d,%ebx
mov    %ebx,%ebp
imul   %r9d,%ebp
add    %esi,%ebp
imul   %r10d,%ebp
imul   %r8d,%ebx
add    %edi,%ebx
add    %ebp,%ebx
imul   %r12d,%r9d
mov    $0x1,%ebp
sub    %r9d,%ebp
lea    (%rax,%rbp,1),%r9d
imul   %r10d,%r9d
imul   %r8d,%r12d
mov    %edx,%ecx
sub    %r12d,%ecx
lea    (%rcx,%r9,1),%r8d
add    $0x1,%r8d
lea    (%r12,%rdx,1),%ecx
add    $0xffffffff,%ecx
sub    %ebp,%eax
imul   %r10d,%eax
add    %eax,%ecx
movslq %ecx,%rax
movsd  (%r11,%rax,8),%xmm14
movapd %xmm1,%xmm12
movapd %xmm1,%xmm7
movsd  0xe0c(%rip),%xmm11
cmpltsd %xmm8,%xmm7
movsd  0xe15(%rip),%xmm13
movapd %xmm7,%xmm4
andpd  %xmm13,%xmm4
andnpd %xmm11,%xmm7
orpd   %xmm4,%xmm7
subsd  %xmm1,%xmm7
movslq %r15d,%rcx
mulsd  %xmm1,%xmm12
movapd %xmm0,%xmm4
mulsd  %xmm2,%xmm4
movapd %xmm4,%xmm1
cmpltsd %xmm8,%xmm1
andpd  %xmm1,%xmm13
andnpd %xmm11,%xmm1
subsd  %xmm12,%xmm11
movslq %r8d,%rax
movsd  0xdc2(%rip),%xmm12
movsd  (%r11,%rax,8),%xmm3
xor    %eax,%eax
ucomisd %xmm4,%xmm8
mulsd  %xmm2,%xmm14
mulsd  %xmm2,%xmm3
movsd  0xdbb(%rip),%xmm8
movapd %xmm15,%xmm5
mulsd  %xmm8,%xmm5
mulsd  %xmm2,%xmm5
mulsd  %xmm0,%xmm5
mulsd  %xmm2,%xmm5
mulsd  %xmm12,%xmm2
mulsd  %xmm0,%xmm2
orpd   %xmm13,%xmm1
mulsd  %xmm0,%xmm14
movapd %xmm2,%xmm6
addsd  %xmm12,%xmm6
mulsd  %xmm14,%xmm6
movsd  (%r11,%rcx,8),%xmm13
movslq 0x40(%rsp),%rcx
mulsd  %xmm0,%xmm3
addsd  0xd65(%rip),%xmm2
mulsd  %xmm3,%xmm2
movsd  (%r11,%rcx,8),%xmm3
mulsd  %xmm0,%xmm5
movapd %xmm13,%xmm0
subsd  %xmm3,%xmm0
mulsd  %xmm10,%xmm0
movapd %xmm1,%xmm10
mulsd  %xmm13,%xmm10
addsd  %xmm3,%xmm13
subsd  %xmm0,%xmm13
mulsd  0xd14(%rip),%xmm11
mulsd  %xmm8,%xmm13
movslq %r13d,%rcx
addsd  %xmm15,%xmm15
mulsd  %xmm11,%xmm15
addsd  %xmm13,%xmm15
movsd  (%r11,%rcx,8),%xmm0
mulsd  %xmm11,%xmm0
subsd  %xmm0,%xmm15
movslq %ebx,%rcx
mulsd  (%r11,%rcx,8),%xmm11
subsd  %xmm11,%xmm15
subsd  %xmm5,%xmm15
addsd  %xmm2,%xmm15
addsd  %xmm6,%xmm15
mulsd  %xmm12,%xmm4
mulsd  %xmm4,%xmm7
mulsd  %xmm1,%xmm3
mulsd  %xmm9,%xmm3
mulsd  %xmm7,%xmm3
addsd  %xmm15,%xmm3
lea    -0x1(%r12),%ecx
setbe  %al
add    %eax,%eax
add    $0xffffffff,%eax
imul   %eax,%ebp
imul   %eax,%ecx
sub    %ebp,%esi
imul   %r10d,%esi
add    %ecx,%edi
add    %esi,%edi
mulsd  %xmm9,%xmm10
mulsd  %xmm7,%xmm10
subsd  %xmm10,%xmm3
movslq %edi,%rax
movsd  (%r11,%rax,8),%xmm0
sub    %ebp,%r14d
imul   %r10d,%r14d
add    -0x4(%rsp),%ecx
add    %r14d,%ecx
mulsd  %xmm1,%xmm0
mulsd  %xmm9,%xmm0
mulsd  %xmm7,%xmm0
subsd  %xmm0,%xmm3
movslq %ecx,%rax
movsd  (%r11,%rax,8),%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm9,%xmm0
mulsd  %xmm7,%xmm0
addsd  %xmm3,%xmm0
mov    0x48(%rsp),%rax
movsd  %xmm0,(%rax)
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @get_face_2d_utopia(i32 %0, i32 %1, i32 %2, i32 %3, i32 %4, i32 %5, i32 %6, i32 %7, double %8, double %9, double %10, double* nocapture writeonly %11, double* nocapture %12) {
%14 = fmul double %8, %9
%15 = sitofp i32 %6 to double
%16 = fmul double %14, %15
%17 = fcmp ule double %16, 0.000000e+00
%18 = zext i1 %17 to i32
%19 = add nsw i32 %5, %1
%20 = mul nsw i32 %19, %2
%21 = add i32 %4, %0
%22 = add i32 %21, %20
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double, double* %12, i64 %23
%25 = load double, double* %24
%26 = sext i32 %7 to i64
%27 = getelementptr inbounds double, double* %12, i64 %26
%28 = load double, double* %27
%29 = fneg double %14
%30 = tail call double @llvm.fmuladd.f64(double %29, double %14, double 1.000000e+00)
%31 = fmul double %30, 0x3FC5555555555553
%32 = select i1 %17, i32 %5, i32 0
%33 = add nsw i32 %32, %1
%34 = mul nsw i32 %33, %2
%35 = select i1 %17, i32 %4, i32 0
%36 = add i32 %35, %0
%37 = add i32 %36, %34
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds double, double* %12, i64 %38
%40 = load double, double* %39
%41 = add nsw i32 %18, %6
%42 = mul nsw i32 %41, %5
%43 = add nsw i32 %42, %1
%44 = mul nsw i32 %43, %2
%45 = mul nsw i32 %41, %4
%46 = add i32 %45, %0
%47 = add i32 %46, %44
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds double, double* %12, i64 %48
%50 = load double, double* %49
%51 = sub nsw i32 %18, %6
%52 = mul nsw i32 %51, %5
%53 = add nsw i32 %52, %1
%54 = mul nsw i32 %53, %2
%55 = mul nsw i32 %51, %4
%56 = add i32 %55, %0
%57 = add i32 %56, %54
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double, double* %12, i64 %58
%60 = load double, double* %59
%61 = mul nsw i32 %6, %5
%62 = sub nsw i32 1, %61
%63 = add nsw i32 %33, %62
%64 = mul nsw i32 %63, %2
%65 = mul nsw i32 %6, %4
%66 = add i32 %65, -1
%67 = sub nsw i32 1, %65
%68 = add i32 %36, %67
%69 = add i32 %68, %64
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds double, double* %12, i64 %70
%72 = load double, double* %71
%73 = fmul double %10, 2.500000e-01
%74 = fmul double %73, %8
%75 = sub nsw i32 %33, %62
%76 = mul nsw i32 %75, %2
%77 = add i32 %36, %66
%78 = add i32 %77, %76
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds double, double* %12, i64 %79
%81 = load double, double* %80
%82 = fcmp olt double %14, 0.000000e+00
%83 = select i1 %82, double -1.000000e+00, double 1.000000e+00
%84 = fmul double %8, %10
%85 = fcmp olt double %84, 0.000000e+00
%86 = select i1 %85, i32 -1, i32 1
%87 = select i1 %85, double -1.000000e+00, double 1.000000e+00
%88 = fmul double %81, %10
%89 = fmul double %88, %8
%90 = fadd double %74, 2.500000e-01
%91 = fmul double %72, %10
%92 = fmul double %91, %8
%93 = fadd double %74, -2.500000e-01
%94 = fmul double %40, 5.000000e-01
%95 = fmul double %94, %10
%96 = fmul double %95, %8
%97 = fmul double %96, %10
%98 = fmul double %97, %8
%99 = fneg double %60
%100 = fneg double %50
%101 = fmul double %40, 2.000000e+00
%102 = fneg double %16
%103 = fsub double %25, %28
%104 = fadd double %25, %28
%105 = tail call double @llvm.fmuladd.f64(double %102, double %103, double %104)
%106 = fmul double %105, 5.000000e-01
%107 = tail call double @llvm.fmuladd.f64(double %101, double %31, double %106)
%108 = tail call double @llvm.fmuladd.f64(double %100, double %31, double %107)
%109 = tail call double @llvm.fmuladd.f64(double %99, double %31, double %108)
%110 = fsub double %109, %98
%111 = tail call double @llvm.fmuladd.f64(double %92, double %93, double %110)
%112 = tail call double @llvm.fmuladd.f64(double %89, double %90, double %111)
%113 = fmul double %84, 2.500000e-01
%114 = fsub double %83, %14
%115 = fmul double %113, %114
%116 = fmul double %87, %28
%117 = fmul double %116, %15
%118 = tail call double @llvm.fmuladd.f64(double %117, double %115, double %112)
%119 = fneg double %25
%120 = fmul double %87, %119
%121 = fmul double %120, %15
%122 = tail call double @llvm.fmuladd.f64(double %121, double %115, double %118)
%123 = mul nsw i32 %86, %62
%124 = sub nsw i32 %1, %123
%125 = mul nsw i32 %124, %2
%126 = mul i32 %66, %86
%127 = add i32 %126, %0
%128 = add i32 %127, %125
%129 = sext i32 %128 to i64
%130 = getelementptr inbounds double, double* %12, i64 %129
%131 = load double, double* %130
%132 = fneg double %131
%133 = fmul double %87, %132
%134 = fmul double %133, %15
%135 = tail call double @llvm.fmuladd.f64(double %134, double %115, double %122)
%136 = sub nsw i32 %19, %123
%137 = mul nsw i32 %136, %2
%138 = add i32 %21, %126
%139 = add i32 %138, %137
%140 = sext i32 %139 to i64
%141 = getelementptr inbounds double, double* %12, i64 %140
%142 = load double, double* %141
%143 = fmul double %87, %142
%144 = fmul double %143, %15
%145 = tail call double @llvm.fmuladd.f64(double %144, double %115, double %135)
store double %145, double* %11
ret void
}