#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 */
typedef int int64_t;

/* Variables and functions */
int /*<<< orphan*/ printf(char *, int, double);
int /*<<< orphan*/ stdlib_strided_drsqrt(int, double *, int, double *, int);
int main() {
  double x[] = {0.0, 4.0, 9.0, 12.0, 24.0, 64.0, 81.0, 101.0};
  double y[] = {0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0};
  int64_t N = 4;
  int64_t strideX = 2;
  int64_t strideY = 2;
  stdlib_strided_drsqrt(N, x, strideX, y, strideY);
  for (int i = 0; i < 8; i++) {
    printf("y[ %i ] = %lf\n", i, y[i]);
  }
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
movl   $0x0,-0x4(%rbp)
lea    -0x50(%rbp),%rdi
lea    0xea3(%rip),%rsi
mov    $0x40,%edx
call   1060 <memcpy@plt>
lea    -0x90(%rbp),%rdi
xor    %esi,%esi
mov    $0x40,%edx
call   1050 <memset@plt>
movl   $0x4,-0x94(%rbp)
movl   $0x2,-0x98(%rbp)
movl   $0x2,-0x9c(%rbp)
mov    -0x94(%rbp),%edi
lea    -0x50(%rbp),%rsi
mov    -0x98(%rbp),%edx
lea    -0x90(%rbp),%rcx
mov    -0x9c(%rbp),%r8d
call   1030 <stdlib_strided_drsqrt@plt>
movl   $0x0,-0xa0(%rbp)
cmpl   $0x8,-0xa0(%rbp)
jge    1208 <main+0xc8>
mov    -0xa0(%rbp),%esi
movslq -0xa0(%rbp),%rax
movsd  -0x90(%rbp,%rax,8),%xmm0
ff ff
lea    0xe51(%rip),%rdi
call   1040 <printf@plt>
mov    -0xa0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa0(%rbp)
jmp    11c5 <main+0x85>
mov    -0x4(%rbp),%eax
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbx
sub    $0x80,%rsp
movaps 0xf01(%rip),%xmm0
movaps %xmm0,0x70(%rsp)
movaps 0xee5(%rip),%xmm0
movaps %xmm0,0x60(%rsp)
movaps 0xec9(%rip),%xmm0
movaps %xmm0,0x50(%rsp)
movaps 0xead(%rip),%xmm0
movaps %xmm0,0x40(%rsp)
xorps  %xmm0,%xmm0
movaps %xmm0,0x30(%rsp)
movaps %xmm0,0x20(%rsp)
movaps %xmm0,0x10(%rsp)
movaps %xmm0,(%rsp)
lea    0x40(%rsp),%rsi
mov    %rsp,%rcx
mov    $0x4,%edi
mov    $0x2,%edx
mov    $0x2,%r8d
call   1030 <stdlib_strided_drsqrt@plt>
movsd  (%rsp),%xmm0
lea    0xea9(%rip),%rbx
mov    %rbx,%rdi
xor    %esi,%esi
call   1040 <printf@plt>
movsd  0x8(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%esi
call   1040 <printf@plt>
movsd  0x10(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x2,%esi
call   1040 <printf@plt>
movsd  0x18(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x3,%esi
call   1040 <printf@plt>
movsd  0x20(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x4,%esi
call   1040 <printf@plt>
movsd  0x28(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x5,%esi
call   1040 <printf@plt>
movsd  0x30(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x6,%esi
call   1040 <printf@plt>
movsd  0x38(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x7,%esi
call   1040 <printf@plt>
xor    %eax,%eax
add    $0x80,%rsp
pop    %rbx
ret

LLVM IR (-O0):
@__const.main.x = [8 x double] [double 0.000000e+00, double 4.000000e+00, double 9.000000e+00, double 1.200000e+01, double 2.400000e+01, double 6.400000e+01, double 8.100000e+01, double 1.010000e+02]
@.str = [15 x i8] c"y[ %i ] = %lf\0A\00"
define i32 @main() {
%1 = alloca i32
%2 = alloca [8 x double]
%3 = alloca [8 x double]
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
store i32 0, i32* %1
%8 = bitcast [8 x double]* %2 to i8*
%9 = bitcast [8 x double]* %3 to i8*
store i32 4, i32* %4
store i32 2, i32* %5
store i32 2, i32* %6
%10 = load i32, i32* %4
%11 = getelementptr inbounds [8 x double], [8 x double]* %2, i64 0, i64 0
%12 = load i32, i32* %5
%13 = getelementptr inbounds [8 x double], [8 x double]* %3, i64 0, i64 0
%14 = load i32, i32* %6
%15 = call i32 @stdlib_strided_drsqrt(i32 %10, double* %11, i32 %12, double* %13, i32 %14)
store i32 0, i32* %7
br label %16
16: ; preds = %26, %0
%17 = load i32, i32* %7
%18 = icmp slt i32 %17, 8
br i1 %18, label %19, label %29
19: ; preds = %16
%20 = load i32, i32* %7
%21 = load i32, i32* %7
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds [8 x double], [8 x double]* %3, i64 0, i64 %22
%24 = load double, double* %23
%25 = call i32 @printf(i8* getelementptr inbounds ([15 x i8], [15 x i8]* @.str, i64 0, i64 0), i32 %20, double %24)
br label %26
26: ; preds = %19
%27 = load i32, i32* %7
%28 = add nsw i32 %27, 1
store i32 %28, i32* %7
br label %16
29: ; preds = %16
%30 = load i32, i32* %1
ret i32 %30
}

LLVM IR (-O2):
@__const.main.x = [8 x double] [double 0.000000e+00, double 4.000000e+00, double 9.000000e+00, double 1.200000e+01, double 2.400000e+01, double 6.400000e+01, double 8.100000e+01, double 1.010000e+02]
@.str = [15 x i8] c"y[ %i ] = %lf\0A\00"
define i32 @main() {
%1 = alloca [8 x double]
%2 = alloca [8 x double]
%3 = bitcast [8 x double]* %1 to i8*
%4 = bitcast [8 x double]* %2 to i8*
%5 = getelementptr inbounds [8 x double], [8 x double]* %1, i64 0, i64 0
%6 = getelementptr inbounds [8 x double], [8 x double]* %2, i64 0, i64 0
%7 = call i32 @stdlib_strided_drsqrt(i32 4, double* %5, i32 2, double* %6, i32 2)
%8 = load double, double* %6
%9 = call i32 @printf(i8* getelementptr inbounds ([15 x i8], [15 x i8]* @.str, i64 0, i64 0), i32 0, double %8)
%10 = getelementptr inbounds [8 x double], [8 x double]* %2, i64 0, i64 1
%11 = load double, double* %10
%12 = call i32 @printf(i8* getelementptr inbounds ([15 x i8], [15 x i8]* @.str, i64 0, i64 0), i32 1, double %11)
%13 = getelementptr inbounds [8 x double], [8 x double]* %2, i64 0, i64 2
%14 = load double, double* %13
%15 = call i32 @printf(i8* getelementptr inbounds ([15 x i8], [15 x i8]* @.str, i64 0, i64 0), i32 2, double %14)
%16 = getelementptr inbounds [8 x double], [8 x double]* %2, i64 0, i64 3
%17 = load double, double* %16
%18 = call i32 @printf(i8* getelementptr inbounds ([15 x i8], [15 x i8]* @.str, i64 0, i64 0), i32 3, double %17)
%19 = getelementptr inbounds [8 x double], [8 x double]* %2, i64 0, i64 4
%20 = load double, double* %19
%21 = call i32 @printf(i8* getelementptr inbounds ([15 x i8], [15 x i8]* @.str, i64 0, i64 0), i32 4, double %20)
%22 = getelementptr inbounds [8 x double], [8 x double]* %2, i64 0, i64 5
%23 = load double, double* %22
%24 = call i32 @printf(i8* getelementptr inbounds ([15 x i8], [15 x i8]* @.str, i64 0, i64 0), i32 5, double %23)
%25 = getelementptr inbounds [8 x double], [8 x double]* %2, i64 0, i64 6
%26 = load double, double* %25
%27 = call i32 @printf(i8* getelementptr inbounds ([15 x i8], [15 x i8]* @.str, i64 0, i64 0), i32 6, double %26)
%28 = getelementptr inbounds [8 x double], [8 x double]* %2, i64 0, i64 7
%29 = load double, double* %28
%30 = call i32 @printf(i8* getelementptr inbounds ([15 x i8], [15 x i8]* @.str, i64 0, i64 0), i32 7, double %29)
ret i32 0
}