#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 */
float exp(int);
int /*<<< orphan*/ zapit(float *, int);
int gen_brune_stf(float *slip, float *t0, float *stf, int nt, float *dt,
                  float *z0) {
  int it, nstf;
  float t95, tend, sfac, tfac;
  float sum;
  zapit(stf, nt);
  t95 = 1.745 * exp(1.0) * (*t0);
  tend = 3.0 * t95;
  nstf = (int)((tend) / (*dt) + 0.5);
  if (nstf > nt)
    nstf = nt;
  if (nstf == 0)
    return (0);
  sfac = (*slip) / (*t0);
  tfac = (*dt) / (*t0);
  for (it = 0; it < nstf; it++)
    stf[it] = sfac * (it * tfac) * exp(-it * tfac);
  sum = 0.0;
  for (it = 0; it < nstf; it++)
    sum = sum + (*dt) * stf[it];
  if (sum <= 0.0)
    return (0);
  return (nstf);
}


x86 assembly (clang, -O0):
<gen_brune_stf>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %ecx,-0x24(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
mov    -0x20(%rbp),%rdi
mov    -0x24(%rbp),%esi
call   1040 <zapit@plt>
mov    $0x1,%edi
call   1030 <exp@plt>
cvtss2sd %xmm0,%xmm1
movsd  0xeaf(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x44(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe79(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x48(%rbp),%xmm0
mov    -0x30(%rbp),%rax
divss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe4b(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x24(%rbp),%eax
jle    11d2 <gen_brune_stf+0xb2>
mov    -0x24(%rbp),%eax
mov    %eax,-0x40(%rbp)
cmpl   $0x0,-0x40(%rbp)
jne    11e8 <gen_brune_stf+0xc8>
movl   $0x0,-0x4(%rbp)
jmp    12f3 <gen_brune_stf+0x1d3>
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x18(%rbp),%rax
divss  (%rax),%xmm0
movss  %xmm0,-0x4c(%rbp)
mov    -0x30(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x18(%rbp),%rax
divss  (%rax),%xmm0
movss  %xmm0,-0x50(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    127b <gen_brune_stf+0x15b>
movss  -0x4c(%rbp),%xmm0
cvtsi2ssl -0x3c(%rbp),%xmm1
mulss  -0x50(%rbp),%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x58(%rbp)
xor    %eax,%eax
sub    -0x3c(%rbp),%eax
cvtsi2ss %eax,%xmm0
mulss  -0x50(%rbp),%xmm0
cvttss2si %xmm0,%edi
call   1030 <exp@plt>
movaps %xmm0,%xmm1
movss  -0x58(%rbp),%xmm0
mulss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1219 <gen_brune_stf+0xf9>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x54(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    12cb <gen_brune_stf+0x1ab>
movss  -0x54(%rbp),%xmm1
mov    -0x30(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x54(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    128a <gen_brune_stf+0x16a>
movss  -0x54(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jb     12ed <gen_brune_stf+0x1cd>
movl   $0x0,-0x4(%rbp)
jmp    12f3 <gen_brune_stf+0x1d3>
mov    -0x40(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<gen_brune_stf>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %r8,%r14
mov    %ecx,%ebp
mov    %rdx,%r13
mov    %rsi,%r12
mov    %rdi,%r15
mov    %rdx,%rdi
mov    %ecx,%esi
call   1040 <zapit@plt>
mov    $0x1,%edi
call   1030 <exp@plt>
cvtss2sd %xmm0,%xmm1
mulsd  0xea4(%rip),%xmm1
movss  (%r12),%xmm0
cvtss2sd %xmm0,%xmm2
mulsd  %xmm1,%xmm2
xorps  %xmm1,%xmm1
cvtsd2ss %xmm2,%xmm1
mulss  0xe97(%rip),%xmm1
movss  (%r14),%xmm2
divss  %xmm2,%xmm1
cvtss2sd %xmm1,%xmm1
addsd  0xe7a(%rip),%xmm1
cvttsd2si %xmm1,%r12d
cmp    %ebp,%r12d
cmovg  %ebp,%r12d
xor    %eax,%eax
test   %r12d,%r12d
je     12d2 <gen_brune_stf+0x1b2>
jle    12d2 <gen_brune_stf+0x1b2>
movss  (%r15),%xmm1
divss  %xmm0,%xmm1
movss  %xmm1,0xc(%rsp)
divss  %xmm0,%xmm2
mov    %r12d,%r15d
xor    %ebp,%ebp
xor    %ebx,%ebx
movss  %xmm2,0x10(%rsp)
nopl   0x0(%rax,%rax,1)
xorps  %xmm1,%xmm1
cvtsi2ss %ebx,%xmm1
mulss  %xmm2,%xmm1
xorps  %xmm0,%xmm0
cvtsi2ss %ebp,%xmm0
mulss  0xc(%rsp),%xmm1
movss  %xmm1,0x14(%rsp)
mulss  %xmm2,%xmm0
cvttss2si %xmm0,%edi
call   1030 <exp@plt>
movss  0x10(%rsp),%xmm2
mulss  0x14(%rsp),%xmm0
movss  %xmm0,0x0(%r13,%rbx,4)
add    $0x1,%rbx
add    $0xffffffff,%ebp
cmp    %rbx,%r15
jne    11d0 <gen_brune_stf+0xb0>
test   %r12d,%r12d
jle    123b <gen_brune_stf+0x11b>
movss  (%r14),%xmm0
lea    -0x1(%r15),%rcx
mov    %r15d,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
jae    1242 <gen_brune_stf+0x122>
xorps  %xmm1,%xmm1
xor    %ecx,%ecx
jmp    1295 <gen_brune_stf+0x175>
xor    %eax,%eax
jmp    12d2 <gen_brune_stf+0x1b2>
and    $0xfffffffc,%r15d
xorps  %xmm1,%xmm1
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movss  0x0(%r13,%rcx,4),%xmm2
mulss  %xmm0,%xmm2
addss  %xmm1,%xmm2
movss  0x4(%r13,%rcx,4),%xmm1
mulss  %xmm0,%xmm1
addss  %xmm2,%xmm1
movss  0x8(%r13,%rcx,4),%xmm2
mulss  %xmm0,%xmm2
addss  %xmm1,%xmm2
movss  0xc(%r13,%rcx,4),%xmm1
mulss  %xmm0,%xmm1
addss  %xmm2,%xmm1
add    $0x4,%rcx
cmp    %rcx,%r15
jne    1250 <gen_brune_stf+0x130>
test   %rax,%rax
je     12c6 <gen_brune_stf+0x1a6>
lea    0x0(,%rcx,4),%rcx
add    %r13,%rcx
xor    %edx,%edx
nopw   0x0(%rax,%rax,1)
movss  (%rcx,%rdx,4),%xmm2
mulss  %xmm0,%xmm2
addss  %xmm2,%xmm1
add    $0x1,%rdx
cmp    %rdx,%rax
jne    12b0 <gen_brune_stf+0x190>
xor    %eax,%eax
xorps  %xmm0,%xmm0
ucomiss %xmm1,%xmm0
cmovb  %r12d,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @gen_brune_stf(float* %0, float* %1, float* %2, i32 %3, float* %4, float* %5) {
%7 = alloca i32
%8 = alloca float*
%9 = alloca float*
%10 = alloca float*
%11 = alloca i32
%12 = alloca float*
%13 = alloca float*
%14 = alloca i32
%15 = alloca i32
%16 = alloca float
%17 = alloca float
%18 = alloca float
%19 = alloca float
%20 = alloca float
store float* %0, float** %8
store float* %1, float** %9
store float* %2, float** %10
store i32 %3, i32* %11
store float* %4, float** %12
store float* %5, float** %13
%21 = load float*, float** %10
%22 = load i32, i32* %11
%23 = call i32 @zapit(float* %21, i32 %22)
%24 = call float @exp(i32 1)
%25 = fpext float %24 to double
%26 = fmul double 1.745000e+00, %25
%27 = load float*, float** %9
%28 = load float, float* %27
%29 = fpext float %28 to double
%30 = fmul double %26, %29
%31 = fptrunc double %30 to float
store float %31, float* %16
%32 = load float, float* %16
%33 = fpext float %32 to double
%34 = fmul double 3.000000e+00, %33
%35 = fptrunc double %34 to float
store float %35, float* %17
%36 = load float, float* %17
%37 = load float*, float** %12
%38 = load float, float* %37
%39 = fdiv float %36, %38
%40 = fpext float %39 to double
%41 = fadd double %40, 5.000000e-01
%42 = fptosi double %41 to i32
store i32 %42, i32* %15
%43 = load i32, i32* %15
%44 = load i32, i32* %11
%45 = icmp sgt i32 %43, %44
br i1 %45, label %46, label %48
46: ; preds = %6
%47 = load i32, i32* %11
store i32 %47, i32* %15
br label %48
48: ; preds = %46, %6
%49 = load i32, i32* %15
%50 = icmp eq i32 %49, 0
br i1 %50, label %51, label %52
51: ; preds = %48
store i32 0, i32* %7
br label %114
52: ; preds = %48
%53 = load float*, float** %8
%54 = load float, float* %53
%55 = load float*, float** %9
%56 = load float, float* %55
%57 = fdiv float %54, %56
store float %57, float* %18
%58 = load float*, float** %12
%59 = load float, float* %58
%60 = load float*, float** %9
%61 = load float, float* %60
%62 = fdiv float %59, %61
store float %62, float* %19
store i32 0, i32* %14
br label %63
63: ; preds = %86, %52
%64 = load i32, i32* %14
%65 = load i32, i32* %15
%66 = icmp slt i32 %64, %65
br i1 %66, label %67, label %89
67: ; preds = %63
%68 = load float, float* %18
%69 = load i32, i32* %14
%70 = sitofp i32 %69 to float
%71 = load float, float* %19
%72 = fmul float %70, %71
%73 = fmul float %68, %72
%74 = load i32, i32* %14
%75 = sub nsw i32 0, %74
%76 = sitofp i32 %75 to float
%77 = load float, float* %19
%78 = fmul float %76, %77
%79 = fptosi float %78 to i32
%80 = call float @exp(i32 %79)
%81 = fmul float %73, %80
%82 = load float*, float** %10
%83 = load i32, i32* %14
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds float, float* %82, i64 %84
store float %81, float* %85
br label %86
86: ; preds = %67
%87 = load i32, i32* %14
%88 = add nsw i32 %87, 1
store i32 %88, i32* %14
br label %63
89: ; preds = %63
store float 0.000000e+00, float* %20
store i32 0, i32* %14
br label %90
90: ; preds = %104, %89
%91 = load i32, i32* %14
%92 = load i32, i32* %15
%93 = icmp slt i32 %91, %92
br i1 %93, label %94, label %107
94: ; preds = %90
%95 = load float, float* %20
%96 = load float*, float** %12
%97 = load float, float* %96
%98 = load float*, float** %10
%99 = load i32, i32* %14
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds float, float* %98, i64 %100
%102 = load float, float* %101
%103 = call float @llvm.fmuladd.f32(float %97, float %102, float %95)
store float %103, float* %20
br label %104
104: ; preds = %94
%105 = load i32, i32* %14
%106 = add nsw i32 %105, 1
store i32 %106, i32* %14
br label %90
107: ; preds = %90
%108 = load float, float* %20
%109 = fpext float %108 to double
%110 = fcmp ole double %109, 0.000000e+00
br i1 %110, label %111, label %112
111: ; preds = %107
store i32 0, i32* %7
br label %114
112: ; preds = %107
%113 = load i32, i32* %15
store i32 %113, i32* %7
br label %114
114: ; preds = %112, %111, %51
%115 = load i32, i32* %7
ret i32 %115
}

LLVM IR (-O2):
define i32 @gen_brune_stf(float* nocapture %0, float* nocapture %1, float* %2, i32 %3, float* nocapture %4, float* nocapture readnone %5) {
%7 = tail call i32 @zapit(float* %2, i32 %3)
%8 = tail call float @exp(i32 1)
%9 = fpext float %8 to double
%10 = fmul double %9, 1.745000e+00
%11 = load float, float* %1
%12 = fpext float %11 to double
%13 = fmul double %10, %12
%14 = fptrunc double %13 to float
%15 = fmul float %14, 3.000000e+00
%16 = load float, float* %4
%17 = fdiv float %15, %16
%18 = fpext float %17 to double
%19 = fadd double %18, 5.000000e-01
%20 = fptosi double %19 to i32
%21 = icmp sgt i32 %20, %3
%22 = select i1 %21, i32 %3, i32 %20
%23 = icmp eq i32 %22, 0
br i1 %23, label %96, label %24
24: ; preds = %6
%25 = load float, float* %0
%26 = fdiv float %25, %11
%27 = fdiv float %16, %11
%28 = icmp sgt i32 %22, 0
br i1 %28, label %29, label %96
29: ; preds = %24
%30 = zext i32 %22 to i64
br label %39
31: ; preds = %39
br i1 %28, label %32, label %96
32: ; preds = %31
%33 = load float, float* %4
%34 = add nsw i64 %30, -1
%35 = and i64 %30, 3
%36 = icmp ult i64 %34, 3
br i1 %36, label %77, label %37
37: ; preds = %32
%38 = and i64 %30, 4294967292
br label %55
39: ; preds = %29, %39
%40 = phi i64 [ 0, %29 ], [ %53, %39 ]
%41 = trunc i64 %40 to i32
%42 = sitofp i32 %41 to float
%43 = fmul float %27, %42
%44 = fmul float %26, %43
%45 = trunc i64 %40 to i32
%46 = sub i32 0, %45
%47 = sitofp i32 %46 to float
%48 = fmul float %27, %47
%49 = fptosi float %48 to i32
%50 = tail call float @exp(i32 %49)
%51 = fmul float %44, %50
%52 = getelementptr inbounds float, float* %2, i64 %40
store float %51, float* %52
%53 = add nuw nsw i64 %40, 1
%54 = icmp eq i64 %53, %30
br i1 %54, label %31, label %39
55: ; preds = %55, %37
%56 = phi i64 [ 0, %37 ], [ %74, %55 ]
%57 = phi float [ 0.000000e+00, %37 ], [ %73, %55 ]
%58 = phi i64 [ 0, %37 ], [ %75, %55 ]
%59 = getelementptr inbounds float, float* %2, i64 %56
%60 = load float, float* %59
%61 = tail call float @llvm.fmuladd.f32(float %33, float %60, float %57)
%62 = or i64 %56, 1
%63 = getelementptr inbounds float, float* %2, i64 %62
%64 = load float, float* %63
%65 = tail call float @llvm.fmuladd.f32(float %33, float %64, float %61)
%66 = or i64 %56, 2
%67 = getelementptr inbounds float, float* %2, i64 %66
%68 = load float, float* %67
%69 = tail call float @llvm.fmuladd.f32(float %33, float %68, float %65)
%70 = or i64 %56, 3
%71 = getelementptr inbounds float, float* %2, i64 %70
%72 = load float, float* %71
%73 = tail call float @llvm.fmuladd.f32(float %33, float %72, float %69)
%74 = add nuw nsw i64 %56, 4
%75 = add i64 %58, 4
%76 = icmp eq i64 %75, %38
br i1 %76, label %77, label %55
77: ; preds = %55, %32
%78 = phi float [ undef, %32 ], [ %73, %55 ]
%79 = phi i64 [ 0, %32 ], [ %74, %55 ]
%80 = phi float [ 0.000000e+00, %32 ], [ %73, %55 ]
%81 = icmp eq i64 %35, 0
br i1 %81, label %92, label %82
82: ; preds = %77, %82
%83 = phi i64 [ %89, %82 ], [ %79, %77 ]
%84 = phi float [ %88, %82 ], [ %80, %77 ]
%85 = phi i64 [ %90, %82 ], [ 0, %77 ]
%86 = getelementptr inbounds float, float* %2, i64 %83
%87 = load float, float* %86
%88 = tail call float @llvm.fmuladd.f32(float %33, float %87, float %84)
%89 = add nuw nsw i64 %83, 1
%90 = add i64 %85, 1
%91 = icmp eq i64 %90, %35
br i1 %91, label %92, label %82
92: ; preds = %82, %77
%93 = phi float [ %78, %77 ], [ %88, %82 ]
%94 = fcmp ugt float %93, 0.000000e+00
%95 = select i1 %94, i32 %22, i32 0
br label %96
96: ; preds = %24, %31, %92, %6
%97 = phi i32 [ 0, %6 ], [ 0, %31 ], [ %95, %92 ], [ 0, %24 ]
ret i32 %97
}