#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 M_PI;
double abs(double);
double cos(double);
int /*<<< orphan*/ fault_params_to_okada_form(double, double, double, double,
                                              double, double, double, double *,
                                              double *, double *, double *,
                                              double *, double *, double *,
                                              double *, double *, double *,
                                              double *);
double sin(double);
void getM(double *fi, int nmod, double mu, double *M) {
  int j;
  double *fault;
  double *model;
  double fx1, fy1;
  double fx2, fy2;
  double dip;
  double D;
  double Pr;
  double bd;
  double ss;
  double ds;
  double ts;
  double strike;
  double L, W;
  double ofx, ofy;
  double ofxe, ofye;
  double tfx, tfy;
  double tfxe, tfye;
  double length;
  double width;
  double depth;
  double str;
  double east;
  double north;
  for (j = 0; j < nmod; j++) {
    fault = &fi[11 * j];
    model = &M[10 * j];
    fx1 = fault[0];
    fy1 = fault[1];
    fx2 = fault[2];
    fy2 = fault[3];
    dip = fault[4];
    D = fault[5];
    Pr = fault[6];
    bd = fault[7];
    ss = fault[8];
    ds = fault[9];
    ts = fault[10];
    fault_params_to_okada_form(fx1, fy1, fx2, fy2, dip * M_PI / 180.0, D, bd,
                               &strike, &L, &W, &ofx, &ofy, &ofxe, &ofye, &tfx,
                               &tfy, &tfxe, &tfye);
    length = L;
    width = W;
    if (dip < 0) {
      east = (L / 2) * cos(strike) + tfx;
      north = (L / 2) * sin(strike) + tfy;
      depth = bd;
      width = abs(width);
      ss = -1 * ss;
      ds = -1 * ds;
      ts = -1 * ts;
    } else {
      east = (L / 2) * cos(strike) + ofx;
      north = (L / 2) * sin(strike) + ofy;
      depth = D;
    }
    str = (180.0 / M_PI) * ((M_PI / 2) - strike);
    model[0] = length;
    model[1] = width;
    model[2] = depth;
    model[3] = dip;
    model[4] = str;
    model[5] = east;
    model[6] = north;
    model[7] = ss;
    model[8] = ds;
    model[9] = ts;
  }
}


x86 assembly (clang, -O0):
<getM>:
push   %rbp
mov    %rsp,%rbp
push   %r14
push   %rbx
sub    $0x160,%rsp
mov    %rdi,-0x18(%rbp)
mov    %esi,-0x1c(%rbp)
movsd  %xmm0,-0x28(%rbp)
mov    %rdx,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    15b8 <getM+0x478>
mov    -0x18(%rbp),%rax
imul   $0xb,-0x34(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x30(%rbp),%rax
imul   $0xa,-0x34(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
mov    -0x40(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x40(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x40(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x40(%rbp),%rax
movsd  0x18(%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x40(%rbp),%rax
movsd  0x20(%rax),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x40(%rbp),%rax
movsd  0x28(%rax),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x40(%rbp),%rax
movsd  0x30(%rax),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x40(%rbp),%rax
movsd  0x38(%rax),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x40(%rbp),%rax
movsd  0x40(%rax),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x40(%rbp),%rax
movsd  0x48(%rax),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x40(%rbp),%rax
movsd  0x50(%rax),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm1
movsd  -0x60(%rbp),%xmm2
movsd  -0x68(%rbp),%xmm3
movsd  -0x70(%rbp),%xmm4
mov    0x2d86(%rip),%rax
mulsd  (%rax),%xmm4
movsd  0xd92(%rip),%xmm5
divsd  %xmm5,%xmm4
movsd  -0x78(%rbp),%xmm5
movsd  -0x88(%rbp),%xmm6
ff
lea    -0xa8(%rbp),%rdi
lea    -0xb0(%rbp),%rsi
lea    -0xb8(%rbp),%rdx
lea    -0xc0(%rbp),%rcx
lea    -0xc8(%rbp),%r8
lea    -0xd0(%rbp),%r9
lea    -0xd8(%rbp),%r14
lea    -0xe0(%rbp),%rbx
lea    -0xe8(%rbp),%r11
lea    -0xf0(%rbp),%r10
lea    -0xf8(%rbp),%rax
mov    %r14,(%rsp)
mov    %rbx,0x8(%rsp)
mov    %r11,0x10(%rsp)
mov    %r10,0x18(%rsp)
mov    %rax,0x20(%rsp)
call   1050 <fault_params_to_okada_form@plt>
movsd  -0xb0(%rbp),%xmm0
ff
movsd  %xmm0,-0x100(%rbp)
ff
movsd  -0xb8(%rbp),%xmm0
movsd  %xmm0,-0x108(%rbp)
ff
xorps  %xmm0,%xmm0
ucomisd -0x70(%rbp),%xmm0
jbe    1421 <getM+0x2e1>
movsd  -0xb0(%rbp),%xmm0
ff
movsd  0xce1(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x138(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
call   1060 <cos@plt>
movaps %xmm0,%xmm2
movsd  -0x138(%rbp),%xmm0
ff
movsd  -0xe0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
movsd  0xc95(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x130(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
call   1040 <sin@plt>
movaps %xmm0,%xmm2
movsd  -0x130(%rbp),%xmm0
ff
movsd  -0xe8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x128(%rbp)
ff
movsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x110(%rbp)
ff
movsd  -0x108(%rbp),%xmm0
ff
call   1030 <abs@plt>
movsd  %xmm0,-0x108(%rbp)
ff
movsd  0xc34(%rip),%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0x90(%rbp)
ff
movsd  0xc1c(%rip),%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  0xc04(%rip),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0xa0(%rbp)
ff
jmp    14c6 <getM+0x386>
movsd  -0xb0(%rbp),%xmm0
ff
movsd  0xbd7(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x148(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
call   1060 <cos@plt>
movaps %xmm0,%xmm2
movsd  -0x148(%rbp),%xmm0
ff
movsd  -0xc0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
movsd  0xb8b(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x140(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
call   1040 <sin@plt>
movaps %xmm0,%xmm2
movsd  -0x140(%rbp),%xmm0
ff
movsd  -0xc8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x128(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x110(%rbp)
ff
movsd  0xb32(%rip),%xmm0
mov    0x2b13(%rip),%rax
divsd  (%rax),%xmm0
mov    0x2b08(%rip),%rax
movsd  (%rax),%xmm1
movsd  0xb1c(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  -0xa8(%rbp),%xmm1
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x118(%rbp)
ff
movsd  -0x100(%rbp),%xmm0
ff
mov    -0x48(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x108(%rbp),%xmm0
ff
mov    -0x48(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movsd  -0x110(%rbp),%xmm0
ff
mov    -0x48(%rbp),%rax
movsd  %xmm0,0x10(%rax)
movsd  -0x70(%rbp),%xmm0
mov    -0x48(%rbp),%rax
movsd  %xmm0,0x18(%rax)
movsd  -0x118(%rbp),%xmm0
ff
mov    -0x48(%rbp),%rax
movsd  %xmm0,0x20(%rax)
movsd  -0x120(%rbp),%xmm0
ff
mov    -0x48(%rbp),%rax
movsd  %xmm0,0x28(%rax)
movsd  -0x128(%rbp),%xmm0
ff
mov    -0x48(%rbp),%rax
movsd  %xmm0,0x30(%rax)
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x48(%rbp),%rax
movsd  %xmm0,0x38(%rax)
movsd  -0x98(%rbp),%xmm0
ff
mov    -0x48(%rbp),%rax
movsd  %xmm0,0x40(%rax)
movsd  -0xa0(%rbp),%xmm0
ff
mov    -0x48(%rbp),%rax
movsd  %xmm0,0x48(%rax)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1165 <getM+0x25>
add    $0x160,%rsp
pop    %rbx
pop    %r14
pop    %rbp
ret

x86 assembly (clang, -O2):
<getM>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xf8,%rsp
test   %esi,%esi
jle    143d <getM+0x2fd>
mov    %rdx,%r14
mov    %rdi,%rbx
mov    %esi,%eax
add    %rax,%rax
lea    (%rax,%rax,4),%r13
xor    %r12d,%r12d
mov    0x2e76(%rip),%r15
movsd  0xe85(%rip),%xmm8
jmp    1276 <getM+0x136>
movsd  0x8(%rsp),%xmm2
mulsd  0xe7a(%rip),%xmm2
movapd %xmm2,0x50(%rsp)
movsd  (%rsp),%xmm0
movsd  0x80(%rsp),%xmm2
movaps %xmm2,0x10(%rsp)
call   1040 <sin@plt>
movapd 0x10(%rsp),%xmm7
movapd 0xa0(%rsp),%xmm2
unpcklpd 0x50(%rsp),%xmm2
movapd 0x90(%rsp),%xmm1
unpcklpd %xmm0,%xmm1
mulpd  %xmm2,%xmm1
movhpd 0x78(%rsp),%xmm7
addpd  %xmm1,%xmm7
movsd  0x40(%rsp),%xmm4
movsd  0xe15(%rip),%xmm8
movapd 0xb0(%rsp),%xmm2
movapd 0xc0(%rsp),%xmm3
movsd  0x28(%rsp),%xmm6
movsd  (%r15),%xmm0
movapd %xmm8,%xmm1
divsd  %xmm0,%xmm1
mulsd  0xdef(%rip),%xmm0
subsd  (%rsp),%xmm0
mulsd  %xmm1,%xmm0
movsd  0x30(%rsp),%xmm1
movsd  %xmm1,(%r14,%rbp,8)
movsd  %xmm6,0x8(%r14,%rbp,8)
movsd  %xmm4,0x10(%r14,%rbp,8)
movsd  0x38(%rsp),%xmm1
movsd  %xmm1,0x18(%r14,%rbp,8)
movsd  %xmm0,0x20(%r14,%rbp,8)
movupd %xmm7,0x28(%r14,%rbp,8)
movupd %xmm3,0x38(%r14,%rbp,8)
movsd  %xmm2,0x48(%r14,%rbp,8)
add    $0xa,%r12
add    $0x58,%rbx
cmp    %r12,%r13
je     143d <getM+0x2fd>
mov    %r12d,%ebp
and    $0xfffffffe,%ebp
movsd  (%rbx),%xmm0
movsd  0x8(%rbx),%xmm1
movsd  0x10(%rbx),%xmm2
movsd  0x18(%rbx),%xmm3
movsd  0x20(%rbx),%xmm7
movsd  %xmm7,0x38(%rsp)
movsd  0x28(%rbx),%xmm5
movsd  0x38(%rbx),%xmm6
movups 0x40(%rbx),%xmm4
movaps %xmm4,0xc0(%rsp)
movsd  0x50(%rbx),%xmm4
movaps %xmm4,0xb0(%rsp)
movsd  (%r15),%xmm4
mulsd  %xmm7,%xmm4
divsd  %xmm8,%xmm4
sub    $0x8,%rsp
movsd  %xmm5,0x48(%rsp)
movsd  %xmm6,0x58(%rsp)
lea    0x8(%rsp),%rdi
lea    0x10(%rsp),%rsi
lea    0x90(%rsp),%rdx
lea    0x88(%rsp),%rcx
lea    0x80(%rsp),%r8
lea    0xf8(%rsp),%r9
lea    0xe0(%rsp),%rax
push   %rax
lea    0xf0(%rsp),%rax
push   %rax
lea    0x80(%rsp),%rax
push   %rax
lea    0x90(%rsp),%rax
push   %rax
lea    0x110(%rsp),%rax
push   %rax
call   1050 <fault_params_to_okada_form@plt>
add    $0x30,%rsp
movsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x30(%rsp)
movsd  0x88(%rsp),%xmm0
movsd  %xmm0,0x28(%rsp)
movsd  (%rsp),%xmm0
call   1060 <cos@plt>
movapd %xmm0,%xmm1
xorpd  %xmm0,%xmm0
ucomisd 0x38(%rsp),%xmm0
movsd  0x30(%rsp),%xmm0
mulsd  0xc8c(%rip),%xmm0
movapd %xmm0,0xa0(%rsp)
movapd %xmm1,0x90(%rsp)
jbe    1180 <getM+0x40>
movsd  0x8(%rsp),%xmm2
mulsd  0xc66(%rip),%xmm2
movapd %xmm2,0x40(%rsp)
movsd  (%rsp),%xmm0
movsd  0x70(%rsp),%xmm2
movaps %xmm2,0x10(%rsp)
call   1040 <sin@plt>
movapd 0xa0(%rsp),%xmm2
unpcklpd 0x40(%rsp),%xmm2
movapd 0x90(%rsp),%xmm1
unpcklpd %xmm0,%xmm1
mulpd  %xmm2,%xmm1
movapd 0x10(%rsp),%xmm0
movhpd 0x68(%rsp),%xmm0
addpd  %xmm1,%xmm0
movapd %xmm0,0x10(%rsp)
movsd  0x28(%rsp),%xmm0
call   1030 <abs@plt>
movapd 0x10(%rsp),%xmm7
movaps %xmm0,%xmm6
movapd 0xc11(%rip),%xmm0
movapd 0xc0(%rsp),%xmm3
xorpd  %xmm0,%xmm3
movapd 0xb0(%rsp),%xmm2
xorpd  %xmm0,%xmm2
movsd  0xbce(%rip),%xmm8
movsd  0x50(%rsp),%xmm4
jmp    1203 <getM+0xc3>
add    $0xf8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@M_PI = global double 0.000000e+00
define void @getM(double* nocapture %0, i32 %1, double %2, double* nocapture writeonly %3) {
%5 = alloca double
%6 = alloca double
%7 = alloca double
%8 = alloca double
%9 = alloca double
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = bitcast double* %5 to i8*
%17 = bitcast double* %6 to i8*
%18 = bitcast double* %7 to i8*
%19 = bitcast double* %8 to i8*
%20 = bitcast double* %9 to i8*
%21 = bitcast double* %10 to i8*
%22 = bitcast double* %11 to i8*
%23 = bitcast double* %12 to i8*
%24 = bitcast double* %13 to i8*
%25 = bitcast double* %14 to i8*
%26 = bitcast double* %15 to i8*
%27 = icmp sgt i32 %1, 0
br i1 %27, label %28, label %119
28: ; preds = %4
%29 = zext i32 %1 to i64
br label %30
30: ; preds = %28, %96
%31 = phi i64 [ 0, %28 ], [ %117, %96 ]
%32 = mul nuw nsw i64 %31, 11
%33 = getelementptr inbounds double, double* %0, i64 %32
%34 = mul i64 %31, 10
%35 = and i64 %34, 4294967294
%36 = getelementptr inbounds double, double* %3, i64 %35
%37 = load double, double* %33
%38 = getelementptr inbounds double, double* %33, i64 1
%39 = load double, double* %38
%40 = getelementptr inbounds double, double* %33, i64 2
%41 = load double, double* %40
%42 = getelementptr inbounds double, double* %33, i64 3
%43 = load double, double* %42
%44 = getelementptr inbounds double, double* %33, i64 4
%45 = load double, double* %44
%46 = getelementptr inbounds double, double* %33, i64 5
%47 = load double, double* %46
%48 = getelementptr inbounds double, double* %33, i64 7
%49 = load double, double* %48
%50 = getelementptr inbounds double, double* %33, i64 8
%51 = bitcast double* %50 to <2 x double>*
%52 = load <2 x double>, <2 x double>* %51
%53 = getelementptr inbounds double, double* %33, i64 10
%54 = load double, double* %53
%55 = load double, double* @M_PI
%56 = fmul double %45, %55
%57 = fdiv double %56, 1.800000e+02
%58 = call i32 @fault_params_to_okada_form(double %37, double %39, double %41, double %43, double %57, double %47, double %49, double* %5, double* %6, double* %7, double* %8, double* %9, double* %10, double* %11, double* %12, double* %13, double* %14, double* %15)
%59 = load double, double* %6
%60 = load double, double* %7
%61 = fcmp olt double %45, 0.000000e+00
%62 = fmul double %59, 5.000000e-01
%63 = load double, double* %5
%64 = call double @cos(double %63)
br i1 %61, label %65, label %82
65: ; preds = %30
%66 = load double, double* %12
%67 = load double, double* %6
%68 = fmul double %67, 5.000000e-01
%69 = load double, double* %5
%70 = call double @sin(double %69)
%71 = load double, double* %13
%72 = insertelement <2 x double> poison, double %62, i64 0
%73 = insertelement <2 x double> %72, double %68, i64 1
%74 = insertelement <2 x double> poison, double %64, i64 0
%75 = insertelement <2 x double> %74, double %70, i64 1
%76 = insertelement <2 x double> poison, double %66, i64 0
%77 = insertelement <2 x double> %76, double %71, i64 1
%78 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %73, <2 x double> %75, <2 x double> %77)
%79 = call double @abs(double %60)
%80 = fneg <2 x double> %52
%81 = fneg double %54
br label %96
82: ; preds = %30
%83 = load double, double* %8
%84 = load double, double* %6
%85 = fmul double %84, 5.000000e-01
%86 = load double, double* %5
%87 = call double @sin(double %86)
%88 = load double, double* %9
%89 = insertelement <2 x double> poison, double %62, i64 0
%90 = insertelement <2 x double> %89, double %85, i64 1
%91 = insertelement <2 x double> poison, double %64, i64 0
%92 = insertelement <2 x double> %91, double %87, i64 1
%93 = insertelement <2 x double> poison, double %83, i64 0
%94 = insertelement <2 x double> %93, double %88, i64 1
%95 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %90, <2 x double> %92, <2 x double> %94)
br label %96
96: ; preds = %82, %65
%97 = phi double [ %81, %65 ], [ %54, %82 ]
%98 = phi double [ %79, %65 ], [ %60, %82 ]
%99 = phi double [ %49, %65 ], [ %47, %82 ]
%100 = phi <2 x double> [ %78, %65 ], [ %95, %82 ]
%101 = phi <2 x double> [ %80, %65 ], [ %52, %82 ]
%102 = load double, double* @M_PI
%103 = fdiv double 1.800000e+02, %102
%104 = fmul double %102, 5.000000e-01
%105 = load double, double* %5
%106 = fsub double %104, %105
%107 = fmul double %103, %106
store double %59, double* %36
%108 = getelementptr inbounds double, double* %36, i64 1
store double %98, double* %108
%109 = getelementptr inbounds double, double* %36, i64 2
store double %99, double* %109
%110 = getelementptr inbounds double, double* %36, i64 3
store double %45, double* %110
%111 = getelementptr inbounds double, double* %36, i64 4
store double %107, double* %111
%112 = getelementptr inbounds double, double* %36, i64 5
%113 = bitcast double* %112 to <2 x double>*
store <2 x double> %100, <2 x double>* %113
%114 = getelementptr inbounds double, double* %36, i64 7
%115 = bitcast double* %114 to <2 x double>*
store <2 x double> %101, <2 x double>* %115
%116 = getelementptr inbounds double, double* %36, i64 9
store double %97, double* %116
%117 = add nuw nsw i64 %31, 1
%118 = icmp eq i64 %117, %29
br i1 %118, label %119, label %30
119: ; preds = %96, %4
ret void
}