#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 fabs(double);
void ddrand(int n, double *x, int incx, int *seeds) {
  int ix, xix;
  if (n < 1)
    return;
  for (ix = 1; ix <= 1 + (n - 1) * incx; ix = ix + incx) {
    seeds[1] = 171 * (seeds[1] % 177) - 2 * (seeds[1] / 177);
    seeds[2] = 172 * (seeds[2] % 176) - 35 * (seeds[2] / 176);
    seeds[3] = 170 * (seeds[3] % 178) - 63 * (seeds[3] / 178);
    if (seeds[1] < 0)
      seeds[1] = seeds[1] + 30269;
    if (seeds[2] < 0)
      seeds[2] = seeds[2] + 30307;
    if (seeds[3] < 0)
      seeds[3] = seeds[3] + 30323;
    x[ix] = ((double)seeds[1]) / 30269.0 + ((double)seeds[2]) / 30307.0 +
            ((double)seeds[3]) / 30323.0;
    xix = (int)x[ix];
    x[ix] = fabs(x[ix] - xix);
  }
}


x86 assembly (clang, -O0):
<ddrand>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
cmpl   $0x1,-0x4(%rbp)
jge    1121 <ddrand+0x21>
jmp    12d6 <ddrand+0x1d6>
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x14(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     12d6 <ddrand+0x1d6>
mov    -0x20(%rbp),%rax
mov    0x4(%rax),%eax
mov    $0xb1,%ecx
cltd
idiv   %ecx
imul   $0xab,%edx,%ecx
mov    -0x20(%rbp),%rax
mov    0x4(%rax),%eax
mov    $0xb1,%esi
cltd
idiv   %esi
shl    $0x1,%eax
sub    %eax,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,0x4(%rax)
mov    -0x20(%rbp),%rax
mov    0x8(%rax),%eax
mov    $0xb0,%ecx
cltd
idiv   %ecx
imul   $0xac,%edx,%ecx
mov    -0x20(%rbp),%rax
mov    0x8(%rax),%eax
mov    $0xb0,%esi
cltd
idiv   %esi
imul   $0x23,%eax,%eax
sub    %eax,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,0x8(%rax)
mov    -0x20(%rbp),%rax
mov    0xc(%rax),%eax
mov    $0xb2,%ecx
cltd
idiv   %ecx
imul   $0xaa,%edx,%ecx
mov    -0x20(%rbp),%rax
mov    0xc(%rax),%eax
mov    $0xb2,%esi
cltd
idiv   %esi
imul   $0x3f,%eax,%eax
sub    %eax,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,0xc(%rax)
mov    -0x20(%rbp),%rax
cmpl   $0x0,0x4(%rax)
jge    11f2 <ddrand+0xf2>
mov    -0x20(%rbp),%rax
mov    0x4(%rax),%ecx
add    $0x763d,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,0x4(%rax)
mov    -0x20(%rbp),%rax
cmpl   $0x0,0x8(%rax)
jge    1214 <ddrand+0x114>
mov    -0x20(%rbp),%rax
mov    0x8(%rax),%ecx
add    $0x7663,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,0x8(%rax)
mov    -0x20(%rbp),%rax
cmpl   $0x0,0xc(%rax)
jge    1236 <ddrand+0x136>
mov    -0x20(%rbp),%rax
mov    0xc(%rax),%ecx
add    $0x7673,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,0xc(%rax)
mov    -0x20(%rbp),%rax
cvtsi2sdl 0x4(%rax),%xmm0
movsd  0xdd9(%rip),%xmm1
divsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
cvtsi2sdl 0x8(%rax),%xmm1
movsd  0xdbc(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
cvtsi2sdl 0xc(%rax),%xmm1
movsd  0xd9b(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
cvttsd2si (%rax,%rcx,8),%eax
mov    %eax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x28(%rbp),%xmm1
subsd  %xmm1,%xmm0
movaps 0xd49(%rip),%xmm1
pand   %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    -0x14(%rbp),%eax
mov    %eax,-0x24(%rbp)
jmp    1128 <ddrand+0x28>
pop    %rbp
ret

x86 assembly (clang, -O2):
<ddrand>:
test   %edi,%edi
jle    129b <ddrand+0x19b>
add    $0xffffffff,%edi
imul   %edx,%edi
test   %edi,%edi
js     129b <ddrand+0x19b>
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
add    $0x1,%edi
mov    0x4(%rcx),%eax
mov    0x8(%rcx),%r10d
mov    0xc(%rcx),%r9d
movslq %edx,%r8
mov    $0x1,%r11d
movsd  0xec3(%rip),%xmm0
movsd  0xec3(%rip),%xmm1
movsd  0xec3(%rip),%xmm2
movapd 0xecb(%rip),%xmm3
cs nopw 0x0(%rax,%rax,1)
nop
movslq %eax,%rbx
imul   $0xffffffffb92143fb,%rbx,%rdx
shr    $0x20,%rdx
add    %ebx,%edx
mov    %edx,%ebp
shr    $0x1f,%ebp
imul   $0x46debc05,%rbx,%rax
sar    $0x7,%edx
shr    $0x20,%rax
sub    %ebx,%eax
mov    %ebx,%r14d
mov    %eax,%r12d
shr    $0x1f,%r12d
add    %ebp,%edx
sar    $0x7,%eax
movslq %r10d,%r15
imul   $0x2e8ba2e9,%r15,%rbp
mov    %rbp,%rbx
add    %r12d,%eax
shr    $0x3f,%rbx
shr    $0x20,%rbp
sar    $0x5,%ebp
imul   $0xb1,%edx,%edx
add    %ebx,%ebp
imul   $0xb0,%ebp,%ebx
sub    %ebx,%r15d
sub    %edx,%r14d
imul   $0xffffffdd,%ebp,%r10d
movslq %r9d,%rbp
imul   $0xffffffffb81702e1,%rbp,%rdx
imul   $0xac,%r15d,%r9d
shr    $0x20,%rdx
add    %ebp,%edx
mov    %edx,%ebx
shr    $0x1f,%ebx
imul   $0xab,%r14d,%r14d
sar    $0x7,%edx
add    %ebx,%edx
imul   $0xb2,%edx,%r15d
lea    (%r14,%rax,2),%r12d
sub    %r15d,%ebp
mov    %edx,%ebx
shl    $0x6,%ebx
sub    %ebx,%edx
imul   $0xaa,%ebp,%r15d
test   %r12d,%r12d
lea    0x763d(%r14,%rax,2),%eax
cmovns %r12d,%eax
lea    (%r10,%r9,1),%ebp
test   %ebp,%ebp
lea    0x7663(%r10,%r9,1),%r10d
lea    (%rdx,%r15,1),%ebx
cmovns %ebp,%r10d
lea    (%rdx,%r15,1),%r9d
add    $0x7673,%r9d
test   %ebx,%ebx
cmovns %ebx,%r9d
xorps  %xmm4,%xmm4
cvtsi2sd %eax,%xmm4
cvtsi2sd %r10d,%xmm5
divsd  %xmm0,%xmm4
divsd  %xmm1,%xmm5
cvtsi2sd %r9d,%xmm6
addsd  %xmm4,%xmm5
divsd  %xmm2,%xmm6
addsd  %xmm5,%xmm6
cvttpd2dq %xmm6,%xmm4
cvtdq2pd %xmm4,%xmm4
subsd  %xmm4,%xmm6
andpd  %xmm3,%xmm6
movlpd %xmm6,(%rsi,%r11,8)
add    %r8,%r11
cmp    %rdi,%r11
jle    1160 <ddrand+0x60>
mov    %eax,0x4(%rcx)
mov    %r10d,0x8(%rcx)
mov    %r9d,0xc(%rcx)
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @ddrand(i32 %0, double* %1, i32 %2, i32* %3) {
%5 = alloca i32
%6 = alloca double*
%7 = alloca i32
%8 = alloca i32*
%9 = alloca i32
%10 = alloca i32
store i32 %0, i32* %5
store double* %1, double** %6
store i32 %2, i32* %7
store i32* %3, i32** %8
%11 = load i32, i32* %5
%12 = icmp slt i32 %11, 1
br i1 %12, label %13, label %14
13: ; preds = %4
br label %143
14: ; preds = %4
store i32 1, i32* %9
br label %15
15: ; preds = %139, %14
%16 = load i32, i32* %9
%17 = load i32, i32* %5
%18 = sub nsw i32 %17, 1
%19 = load i32, i32* %7
%20 = mul nsw i32 %18, %19
%21 = add nsw i32 1, %20
%22 = icmp sle i32 %16, %21
br i1 %22, label %23, label %143
23: ; preds = %15
%24 = load i32*, i32** %8
%25 = getelementptr inbounds i32, i32* %24, i64 1
%26 = load i32, i32* %25
%27 = srem i32 %26, 177
%28 = mul nsw i32 171, %27
%29 = load i32*, i32** %8
%30 = getelementptr inbounds i32, i32* %29, i64 1
%31 = load i32, i32* %30
%32 = sdiv i32 %31, 177
%33 = mul nsw i32 2, %32
%34 = sub nsw i32 %28, %33
%35 = load i32*, i32** %8
%36 = getelementptr inbounds i32, i32* %35, i64 1
store i32 %34, i32* %36
%37 = load i32*, i32** %8
%38 = getelementptr inbounds i32, i32* %37, i64 2
%39 = load i32, i32* %38
%40 = srem i32 %39, 176
%41 = mul nsw i32 172, %40
%42 = load i32*, i32** %8
%43 = getelementptr inbounds i32, i32* %42, i64 2
%44 = load i32, i32* %43
%45 = sdiv i32 %44, 176
%46 = mul nsw i32 35, %45
%47 = sub nsw i32 %41, %46
%48 = load i32*, i32** %8
%49 = getelementptr inbounds i32, i32* %48, i64 2
store i32 %47, i32* %49
%50 = load i32*, i32** %8
%51 = getelementptr inbounds i32, i32* %50, i64 3
%52 = load i32, i32* %51
%53 = srem i32 %52, 178
%54 = mul nsw i32 170, %53
%55 = load i32*, i32** %8
%56 = getelementptr inbounds i32, i32* %55, i64 3
%57 = load i32, i32* %56
%58 = sdiv i32 %57, 178
%59 = mul nsw i32 63, %58
%60 = sub nsw i32 %54, %59
%61 = load i32*, i32** %8
%62 = getelementptr inbounds i32, i32* %61, i64 3
store i32 %60, i32* %62
%63 = load i32*, i32** %8
%64 = getelementptr inbounds i32, i32* %63, i64 1
%65 = load i32, i32* %64
%66 = icmp slt i32 %65, 0
br i1 %66, label %67, label %74
67: ; preds = %23
%68 = load i32*, i32** %8
%69 = getelementptr inbounds i32, i32* %68, i64 1
%70 = load i32, i32* %69
%71 = add nsw i32 %70, 30269
%72 = load i32*, i32** %8
%73 = getelementptr inbounds i32, i32* %72, i64 1
store i32 %71, i32* %73
br label %74
74: ; preds = %67, %23
%75 = load i32*, i32** %8
%76 = getelementptr inbounds i32, i32* %75, i64 2
%77 = load i32, i32* %76
%78 = icmp slt i32 %77, 0
br i1 %78, label %79, label %86
79: ; preds = %74
%80 = load i32*, i32** %8
%81 = getelementptr inbounds i32, i32* %80, i64 2
%82 = load i32, i32* %81
%83 = add nsw i32 %82, 30307
%84 = load i32*, i32** %8
%85 = getelementptr inbounds i32, i32* %84, i64 2
store i32 %83, i32* %85
br label %86
86: ; preds = %79, %74
%87 = load i32*, i32** %8
%88 = getelementptr inbounds i32, i32* %87, i64 3
%89 = load i32, i32* %88
%90 = icmp slt i32 %89, 0
br i1 %90, label %91, label %98
91: ; preds = %86
%92 = load i32*, i32** %8
%93 = getelementptr inbounds i32, i32* %92, i64 3
%94 = load i32, i32* %93
%95 = add nsw i32 %94, 30323
%96 = load i32*, i32** %8
%97 = getelementptr inbounds i32, i32* %96, i64 3
store i32 %95, i32* %97
br label %98
98: ; preds = %91, %86
%99 = load i32*, i32** %8
%100 = getelementptr inbounds i32, i32* %99, i64 1
%101 = load i32, i32* %100
%102 = sitofp i32 %101 to double
%103 = fdiv double %102, 3.026900e+04
%104 = load i32*, i32** %8
%105 = getelementptr inbounds i32, i32* %104, i64 2
%106 = load i32, i32* %105
%107 = sitofp i32 %106 to double
%108 = fdiv double %107, 3.030700e+04
%109 = fadd double %103, %108
%110 = load i32*, i32** %8
%111 = getelementptr inbounds i32, i32* %110, i64 3
%112 = load i32, i32* %111
%113 = sitofp i32 %112 to double
%114 = fdiv double %113, 3.032300e+04
%115 = fadd double %109, %114
%116 = load double*, double** %6
%117 = load i32, i32* %9
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds double, double* %116, i64 %118
store double %115, double* %119
%120 = load double*, double** %6
%121 = load i32, i32* %9
%122 = sext i32 %121 to i64
%123 = getelementptr inbounds double, double* %120, i64 %122
%124 = load double, double* %123
%125 = fptosi double %124 to i32
store i32 %125, i32* %10
%126 = load double*, double** %6
%127 = load i32, i32* %9
%128 = sext i32 %127 to i64
%129 = getelementptr inbounds double, double* %126, i64 %128
%130 = load double, double* %129
%131 = load i32, i32* %10
%132 = sitofp i32 %131 to double
%133 = fsub double %130, %132
%134 = call double @llvm.fabs.f64(double %133)
%135 = load double*, double** %6
%136 = load i32, i32* %9
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds double, double* %135, i64 %137
store double %134, double* %138
br label %139
139: ; preds = %98
%140 = load i32, i32* %9
%141 = load i32, i32* %7
%142 = add nsw i32 %140, %141
store i32 %142, i32* %9
br label %15
143: ; preds = %13, %15
ret void
}

LLVM IR (-O2):
define void @ddrand(i32 %0, double* nocapture writeonly %1, i32 %2, i32* nocapture %3) {
%5 = icmp slt i32 %0, 1
br i1 %5, label %65, label %6
6: ; preds = %4
%7 = add nsw i32 %0, -1
%8 = mul nsw i32 %7, %2
%9 = icmp slt i32 %8, 0
br i1 %9, label %65, label %10
10: ; preds = %6
%11 = add nuw nsw i32 %8, 1
%12 = getelementptr inbounds i32, i32* %3, i64 1
%13 = getelementptr inbounds i32, i32* %3, i64 2
%14 = getelementptr inbounds i32, i32* %3, i64 3
%15 = load i32, i32* %12
%16 = load i32, i32* %13
%17 = load i32, i32* %14
%18 = sext i32 %2 to i64
%19 = zext i32 %11 to i64
br label %20
20: ; preds = %10, %20
%21 = phi i64 [ 1, %10 ], [ %62, %20 ]
%22 = phi i32 [ %17, %10 ], [ %48, %20 ]
%23 = phi i32 [ %16, %10 ], [ %45, %20 ]
%24 = phi i32 [ %15, %10 ], [ %42, %20 ]
%25 = srem i32 %24, 177
%26 = mul nsw i32 %25, 171
%27 = sdiv i32 %24, -177
%28 = shl nsw i32 %27, 1
%29 = add i32 %28, %26
%30 = srem i32 %23, 176
%31 = mul nsw i32 %30, 172
%32 = sdiv i32 %23, 176
%33 = mul nsw i32 %32, -35
%34 = add nsw i32 %33, %31
%35 = srem i32 %22, 178
%36 = mul nsw i32 %35, 170
%37 = sdiv i32 %22, 178
%38 = mul nsw i32 %37, -63
%39 = add nsw i32 %38, %36
%40 = icmp slt i32 %29, 0
%41 = add nsw i32 %29, 30269
%42 = select i1 %40, i32 %41, i32 %29
%43 = icmp slt i32 %34, 0
%44 = add nsw i32 %34, 30307
%45 = select i1 %43, i32 %44, i32 %34
%46 = icmp slt i32 %39, 0
%47 = add nsw i32 %39, 30323
%48 = select i1 %46, i32 %47, i32 %39
%49 = sitofp i32 %42 to double
%50 = fdiv double %49, 3.026900e+04
%51 = sitofp i32 %45 to double
%52 = fdiv double %51, 3.030700e+04
%53 = fadd double %50, %52
%54 = sitofp i32 %48 to double
%55 = fdiv double %54, 3.032300e+04
%56 = fadd double %53, %55
%57 = getelementptr inbounds double, double* %1, i64 %21
%58 = fptosi double %56 to i32
%59 = sitofp i32 %58 to double
%60 = fsub double %56, %59
%61 = tail call double @llvm.fabs.f64(double %60)
store double %61, double* %57
%62 = add i64 %21, %18
%63 = icmp sgt i64 %62, %19
br i1 %63, label %64, label %20
64: ; preds = %20
store i32 %42, i32* %12
store i32 %45, i32* %13
store i32 %48, i32* %14
br label %65
65: ; preds = %6, %64, %4
ret void
}