#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 atan2(double, double);
double cos(double);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *);
double sin(double);
double sqrt(double);
int /*<<< orphan*/ stderr;
double tan(double);
void fault_params_to_okada_form(double sx1, double sy1, double sx2, double sy2,
                                double dip, double D, double bd, double *strike,
                                double *L, double *W, double *ofx, double *ofy,
                                double *ofxe, double *ofye, double *tfx,
                                double *tfy, double *tfxe, double *tfye) {
  if (D < bd) {
    fprintf(stderr, "Burial depth is greater than locking depth!\n");
    fprintf(stderr, "Fault width will be negative!\n");
  }
  if (D == bd) {
    fprintf(stderr, "Burial depth is equal to locking depth!\n");
    fprintf(stderr, "Fault width will be zero!\n");
  }
  *strike = atan2(sy1 - sy2, sx1 - sx2) + M_PI;
  *L = sqrt((sx2 - sx1) * (sx2 - sx1) + (sy2 - sy1) * (sy2 - sy1));
  *W = (D - bd) / sin(dip);
  *ofx = sx1 + D / tan(dip) + sin(*strike);
  *ofy = sy1 - D / tan(dip) + cos(*strike);
  *ofxe = sx2 + D / tan(dip) + sin(*strike);
  *ofye = sy2 - D / tan(dip) + cos(*strike);
  *tfx = sx1 + bd / tan(dip) + sin(*strike);
  *tfy = sy1 - bd / tan(dip) + cos(*strike);
  *tfxe = sx2 + bd / tan(dip) + sin(*strike);
  *tfye = sy2 - bd / tan(dip) + cos(*strike);
}


x86 assembly (clang, -O0):
<fault_params_to_okada_form>:
push   %rbp
mov    %rsp,%rbp
sub    $0x130,%rsp
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  %xmm2,-0x18(%rbp)
movsd  %xmm3,-0x20(%rbp)
movsd  %xmm4,-0x28(%rbp)
movsd  %xmm5,-0x30(%rbp)
movsd  %xmm6,-0x38(%rbp)
mov    %rdi,-0x40(%rbp)
mov    %rsi,-0x48(%rbp)
mov    %rdx,-0x50(%rbp)
mov    %rcx,-0x58(%rbp)
mov    %r8,-0x60(%rbp)
mov    %r9,-0x68(%rbp)
movsd  -0x30(%rbp),%xmm1
movsd  -0x38(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11f8 <fault_params_to_okada_form+0x98>
mov    0x2e03(%rip),%rax
mov    (%rax),%edi
lea    0xe22(%rip),%rsi
call   1040 <fprintf@plt>
mov    0x2dee(%rip),%rax
mov    (%rax),%edi
lea    0xe3a(%rip),%rsi
call   1040 <fprintf@plt>
movsd  -0x30(%rbp),%xmm0
ucomisd -0x38(%rbp),%xmm0
jne    1238 <fault_params_to_okada_form+0xd8>
jp     1238 <fault_params_to_okada_form+0xd8>
mov    0x2dc3(%rip),%rax
mov    (%rax),%edi
lea    0xe2e(%rip),%rsi
call   1040 <fprintf@plt>
mov    0x2dae(%rip),%rax
mov    (%rax),%edi
lea    0xe42(%rip),%rsi
call   1040 <fprintf@plt>
movsd  -0x10(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  -0x8(%rbp),%xmm1
subsd  -0x18(%rbp),%xmm1
call   1030 <atan2@plt>
mov    0x2d90(%rip),%rax
addsd  (%rax),%xmm0
mov    -0x40(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x18(%rbp),%xmm0
subsd  -0x8(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm2
subsd  -0x8(%rbp),%xmm2
movsd  -0x20(%rbp),%xmm1
subsd  -0x10(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm3
subsd  -0x10(%rbp),%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1080 <sqrt@plt>
mov    -0x48(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x30(%rbp),%xmm0
subsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x130(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
call   1050 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x130(%rbp),%xmm0
ff
divsd  %xmm1,%xmm0
mov    -0x50(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x128(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
call   1060 <tan@plt>
movsd  -0x128(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x120(%rbp),%xmm0
ff
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x118(%rbp)
ff
mov    -0x40(%rbp),%rax
movsd  (%rax),%xmm0
call   1050 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x118(%rbp),%xmm0
ff
addsd  %xmm1,%xmm0
mov    -0x58(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x108(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x110(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
call   1060 <tan@plt>
movsd  -0x110(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x108(%rbp),%xmm0
ff
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
mov    -0x40(%rbp),%rax
movsd  (%rax),%xmm0
call   1070 <cos@plt>
movaps %xmm0,%xmm1
movsd  -0x100(%rbp),%xmm0
ff
addsd  %xmm1,%xmm0
mov    -0x60(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
call   1060 <tan@plt>
movsd  -0xf8(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0xf0(%rbp),%xmm0
ff
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
mov    -0x40(%rbp),%rax
movsd  (%rax),%xmm0
call   1050 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0xe8(%rbp),%xmm0
ff
addsd  %xmm1,%xmm0
mov    -0x68(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
call   1060 <tan@plt>
movsd  -0xe0(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0xd8(%rbp),%xmm0
ff
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
mov    -0x40(%rbp),%rax
movsd  (%rax),%xmm0
call   1070 <cos@plt>
movaps %xmm0,%xmm1
movsd  -0xd0(%rbp),%xmm0
ff
addsd  %xmm1,%xmm0
mov    0x10(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
call   1060 <tan@plt>
movsd  -0xc8(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0xc0(%rbp),%xmm0
ff
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
mov    -0x40(%rbp),%rax
movsd  (%rax),%xmm0
call   1050 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0xb8(%rbp),%xmm0
ff
addsd  %xmm1,%xmm0
mov    0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
call   1060 <tan@plt>
movsd  -0xb0(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0xa8(%rbp),%xmm0
ff
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
mov    -0x40(%rbp),%rax
movsd  (%rax),%xmm0
call   1070 <cos@plt>
movaps %xmm0,%xmm1
movsd  -0xa0(%rbp),%xmm0
ff
addsd  %xmm1,%xmm0
mov    0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
call   1060 <tan@plt>
movsd  -0x98(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x90(%rbp),%xmm0
ff
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x40(%rbp),%rax
movsd  (%rax),%xmm0
call   1050 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x88(%rbp),%xmm0
ff
addsd  %xmm1,%xmm0
mov    0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x28(%rbp),%xmm0
call   1060 <tan@plt>
movsd  -0x80(%rbp),%xmm1
movaps %xmm0,%xmm2
movsd  -0x78(%rbp),%xmm0
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x40(%rbp),%rax
movsd  (%rax),%xmm0
call   1070 <cos@plt>
movaps %xmm0,%xmm1
movsd  -0x70(%rbp),%xmm0
addsd  %xmm1,%xmm0
mov    0x30(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x130,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<fault_params_to_okada_form>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %r9,0x50(%rsp)
mov    %r8,0x48(%rsp)
mov    %rcx,%r12
mov    %rdx,%r13
mov    %rsi,%r15
mov    %rdi,%rbx
movsd  %xmm4,0x10(%rsp)
movsd  %xmm3,0x28(%rsp)
movsd  %xmm2,0x40(%rsp)
movsd  %xmm1,0x18(%rsp)
movsd  %xmm0,0x38(%rsp)
movsd  %xmm6,0x20(%rsp)
movsd  %xmm5,0x8(%rsp)
ucomisd %xmm5,%xmm6
jbe    11d9 <fault_params_to_okada_form+0x79>
mov    0x2e1d(%rip),%r14
mov    (%r14),%edi
lea    0xe3b(%rip),%rsi
call   1040 <fprintf@plt>
mov    (%r14),%edi
lea    0xe59(%rip),%rsi
call   1040 <fprintf@plt>
movsd  0x8(%rsp),%xmm0
ucomisd 0x20(%rsp),%xmm0
jne    120e <fault_params_to_okada_form+0xae>
jp     120e <fault_params_to_okada_form+0xae>
mov    0x2de8(%rip),%rbp
mov    0x0(%rbp),%edi
lea    0xe52(%rip),%rsi
call   1040 <fprintf@plt>
mov    0x0(%rbp),%edi
lea    0xe6c(%rip),%rsi
call   1040 <fprintf@plt>
mov    0xb0(%rsp),%r14
mov    0xa8(%rsp),%rbp
movsd  0x18(%rsp),%xmm0
subsd  0x28(%rsp),%xmm0
movsd  0x38(%rsp),%xmm1
subsd  0x40(%rsp),%xmm1
call   1030 <atan2@plt>
mov    0x2da6(%rip),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rbx)
movsd  0x40(%rsp),%xmm0
subsd  0x38(%rsp),%xmm0
movsd  0x28(%rsp),%xmm1
subsd  0x18(%rsp),%xmm1
mulsd  %xmm1,%xmm1
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     127e <fault_params_to_okada_form+0x11e>
sqrtsd %xmm0,%xmm0
jmp    1283 <fault_params_to_okada_form+0x123>
call   1080 <sqrt@plt>
movsd  %xmm0,(%r15)
movsd  0x8(%rsp),%xmm0
subsd  0x20(%rsp),%xmm0
movsd  %xmm0,0x30(%rsp)
movsd  0x10(%rsp),%xmm0
call   1050 <sin@plt>
movsd  0x30(%rsp),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%r13)
movsd  0x10(%rsp),%xmm0
call   1060 <tan@plt>
movsd  0x8(%rsp),%xmm1
divsd  %xmm0,%xmm1
addsd  0x38(%rsp),%xmm1
movsd  %xmm1,0x30(%rsp)
movsd  (%rbx),%xmm0
call   1050 <sin@plt>
addsd  0x30(%rsp),%xmm0
movsd  %xmm0,(%r12)
movsd  0x10(%rsp),%xmm0
call   1060 <tan@plt>
movsd  0x8(%rsp),%xmm1
divsd  %xmm0,%xmm1
movsd  0x18(%rsp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x30(%rsp)
movsd  (%rbx),%xmm0
call   1070 <cos@plt>
addsd  0x30(%rsp),%xmm0
mov    0x48(%rsp),%rax
movsd  %xmm0,(%rax)
movsd  0x10(%rsp),%xmm0
call   1060 <tan@plt>
movsd  0x8(%rsp),%xmm1
divsd  %xmm0,%xmm1
addsd  0x40(%rsp),%xmm1
movsd  %xmm1,0x48(%rsp)
movsd  (%rbx),%xmm0
call   1050 <sin@plt>
addsd  0x48(%rsp),%xmm0
mov    0x50(%rsp),%rax
movsd  %xmm0,(%rax)
movsd  0x10(%rsp),%xmm0
call   1060 <tan@plt>
movsd  0x8(%rsp),%xmm1
divsd  %xmm0,%xmm1
movsd  0x28(%rsp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  (%rbx),%xmm0
call   1070 <cos@plt>
addsd  0x8(%rsp),%xmm0
mov    0x90(%rsp),%rax
movsd  %xmm0,(%rax)
movsd  0x10(%rsp),%xmm0
call   1060 <tan@plt>
movsd  0x20(%rsp),%xmm1
divsd  %xmm0,%xmm1
addsd  0x38(%rsp),%xmm1
movsd  %xmm1,0x8(%rsp)
movsd  (%rbx),%xmm0
call   1050 <sin@plt>
addsd  0x8(%rsp),%xmm0
mov    0x98(%rsp),%rax
movsd  %xmm0,(%rax)
movsd  0x10(%rsp),%xmm0
call   1060 <tan@plt>
movsd  0x20(%rsp),%xmm1
divsd  %xmm0,%xmm1
movsd  0x18(%rsp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x18(%rsp)
movsd  (%rbx),%xmm0
call   1070 <cos@plt>
addsd  0x18(%rsp),%xmm0
mov    0xa0(%rsp),%rax
movsd  %xmm0,(%rax)
movsd  0x10(%rsp),%xmm0
call   1060 <tan@plt>
movsd  0x20(%rsp),%xmm1
divsd  %xmm0,%xmm1
addsd  0x40(%rsp),%xmm1
movsd  %xmm1,0x18(%rsp)
movsd  (%rbx),%xmm0
call   1050 <sin@plt>
addsd  0x18(%rsp),%xmm0
movsd  %xmm0,0x0(%rbp)
movsd  0x10(%rsp),%xmm0
call   1060 <tan@plt>
movsd  0x20(%rsp),%xmm1
divsd  %xmm0,%xmm1
movsd  0x28(%rsp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,0x28(%rsp)
movsd  (%rbx),%xmm0
call   1070 <cos@plt>
addsd  0x28(%rsp),%xmm0
movsd  %xmm0,(%r14)
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@stderr = global i32 0
@.str = [45 x i8] c"Burial depth is greater than locking depth!\0A\00"
@.str.1 = [31 x i8] c"Fault width will be negative!\0A\00"
@.str.2 = [41 x i8] c"Burial depth is equal to locking depth!\0A\00"
@.str.3 = [27 x i8] c"Fault width will be zero!\0A\00"
@M_PI = global double 0.000000e+00
define void @fault_params_to_okada_form(double %0, double %1, double %2, double %3, double %4, double %5, double %6, double* nocapture %7, double* nocapture writeonly %8, double* nocapture writeonly %9, double* nocapture writeonly %10, double* nocapture writeonly %11, double* nocapture writeonly %12, double* nocapture writeonly %13, double* nocapture writeonly %14, double* nocapture writeonly %15, double* nocapture writeonly %16, double* nocapture writeonly %17) {
%19 = fcmp olt double %5, %6
br i1 %19, label %20, label %25
20: ; preds = %18
%21 = load i32, i32* @stderr
%22 = tail call i32 @fprintf(i32 %21, i8* getelementptr inbounds ([45 x i8], [45 x i8]* @.str, i64 0, i64 0))
%23 = load i32, i32* @stderr
%24 = tail call i32 @fprintf(i32 %23, i8* getelementptr inbounds ([31 x i8], [31 x i8]* @.str.1, i64 0, i64 0))
br label %25
25: ; preds = %20, %18
%26 = fcmp oeq double %5, %6
br i1 %26, label %27, label %32
27: ; preds = %25
%28 = load i32, i32* @stderr
%29 = tail call i32 @fprintf(i32 %28, i8* getelementptr inbounds ([41 x i8], [41 x i8]* @.str.2, i64 0, i64 0))
%30 = load i32, i32* @stderr
%31 = tail call i32 @fprintf(i32 %30, i8* getelementptr inbounds ([27 x i8], [27 x i8]* @.str.3, i64 0, i64 0))
br label %32
32: ; preds = %27, %25
%33 = fsub double %1, %3
%34 = fsub double %0, %2
%35 = tail call double @atan2(double %33, double %34)
%36 = load double, double* @M_PI
%37 = fadd double %35, %36
store double %37, double* %7
%38 = fsub double %2, %0
%39 = fsub double %3, %1
%40 = fmul double %39, %39
%41 = tail call double @llvm.fmuladd.f64(double %38, double %38, double %40)
%42 = tail call double @sqrt(double %41)
store double %42, double* %8
%43 = fsub double %5, %6
%44 = tail call double @sin(double %4)
%45 = fdiv double %43, %44
store double %45, double* %9
%46 = tail call double @tan(double %4)
%47 = fdiv double %5, %46
%48 = fadd double %47, %0
%49 = load double, double* %7
%50 = tail call double @sin(double %49)
%51 = fadd double %48, %50
store double %51, double* %10
%52 = tail call double @tan(double %4)
%53 = fdiv double %5, %52
%54 = fsub double %1, %53
%55 = load double, double* %7
%56 = tail call double @cos(double %55)
%57 = fadd double %54, %56
store double %57, double* %11
%58 = tail call double @tan(double %4)
%59 = fdiv double %5, %58
%60 = fadd double %59, %2
%61 = load double, double* %7
%62 = tail call double @sin(double %61)
%63 = fadd double %60, %62
store double %63, double* %12
%64 = tail call double @tan(double %4)
%65 = fdiv double %5, %64
%66 = fsub double %3, %65
%67 = load double, double* %7
%68 = tail call double @cos(double %67)
%69 = fadd double %66, %68
store double %69, double* %13
%70 = tail call double @tan(double %4)
%71 = fdiv double %6, %70
%72 = fadd double %71, %0
%73 = load double, double* %7
%74 = tail call double @sin(double %73)
%75 = fadd double %72, %74
store double %75, double* %14
%76 = tail call double @tan(double %4)
%77 = fdiv double %6, %76
%78 = fsub double %1, %77
%79 = load double, double* %7
%80 = tail call double @cos(double %79)
%81 = fadd double %78, %80
store double %81, double* %15
%82 = tail call double @tan(double %4)
%83 = fdiv double %6, %82
%84 = fadd double %83, %2
%85 = load double, double* %7
%86 = tail call double @sin(double %85)
%87 = fadd double %84, %86
store double %87, double* %16
%88 = tail call double @tan(double %4)
%89 = fdiv double %6, %88
%90 = fsub double %3, %89
%91 = load double, double* %7
%92 = tail call double @cos(double %91)
%93 = fadd double %90, %92
store double %93, double* %17
ret void
}