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

/* Variables and functions */
int /*<<< orphan*/ printf(char *, int, float);
int /*<<< orphan*/ stdlib_strided_smskrsqrt(int, float *, int, int *, int,
                                            float *, int);
int main() {
  float x[] = {0.0, 4.0, 9.0, 12.0, 24.0, 64.0, 81.0, 101.0};
  uint8_t m[] = {0, 0, 1, 0, 1, 1, 0, 0};
  float 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 strideM = 2;
  int64_t strideY = 2;
  stdlib_strided_smskrsqrt(N, x, strideX, m, strideM, y, strideY);
  for (int i = 0; i < 8; i++) {
    printf("y[ %i ] = %f\n", i, y[i]);
  }
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
movl   $0x0,-0x4(%rbp)
mov    0xeb7(%rip),%rax
mov    %rax,-0x30(%rbp)
mov    0xeb4(%rip),%rax
mov    %rax,-0x28(%rbp)
mov    0xeb1(%rip),%rax
mov    %rax,-0x20(%rbp)
mov    0xeae(%rip),%rax
mov    %rax,-0x18(%rbp)
mov    0xeab(%rip),%rax
mov    %rax,-0x50(%rbp)
mov    0xea8(%rip),%rax
mov    %rax,-0x48(%rbp)
mov    0xea5(%rip),%rax
mov    %rax,-0x40(%rbp)
mov    0xea2(%rip),%rax
mov    %rax,-0x38(%rbp)
lea    -0x70(%rbp),%rdi
xor    %esi,%esi
mov    $0x20,%edx
call   1050 <memset@plt>
movl   $0x4,-0x74(%rbp)
movl   $0x2,-0x78(%rbp)
movl   $0x2,-0x7c(%rbp)
movl   $0x2,-0x80(%rbp)
mov    -0x74(%rbp),%edi
lea    -0x30(%rbp),%rsi
mov    -0x78(%rbp),%edx
lea    -0x50(%rbp),%rcx
mov    -0x7c(%rbp),%r8d
lea    -0x70(%rbp),%r9
mov    -0x80(%rbp),%eax
mov    %eax,(%rsp)
call   1030 <stdlib_strided_smskrsqrt@plt>
movl   $0x0,-0x84(%rbp)
cmpl   $0x8,-0x84(%rbp)
jge    1231 <main+0x101>
mov    -0x84(%rbp),%esi
movslq -0x84(%rbp),%rax
movss  -0x70(%rbp,%rax,4),%xmm0
lea    0xe28(%rip),%rdi
call   1040 <printf@plt>
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    11f1 <main+0xc1>
mov    -0x4(%rbp),%eax
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbx
sub    $0x70,%rsp
movaps 0xee4(%rip),%xmm0
movaps %xmm0,0x60(%rsp)
movaps 0xec8(%rip),%xmm0
movaps %xmm0,0x50(%rsp)
movaps 0xeec(%rip),%xmm0
movaps %xmm0,0x40(%rsp)
movaps 0xed0(%rip),%xmm0
movaps %xmm0,0x30(%rsp)
xorps  %xmm0,%xmm0
movaps %xmm0,0x20(%rsp)
movaps %xmm0,0x10(%rsp)
movl   $0x2,(%rsp)
lea    0x50(%rsp),%rsi
lea    0x30(%rsp),%rcx
lea    0x10(%rsp),%r9
mov    $0x4,%edi
mov    $0x2,%edx
mov    $0x2,%r8d
call   1030 <stdlib_strided_smskrsqrt@plt>
movss  0x10(%rsp),%xmm0
lea    0xea6(%rip),%rbx
mov    %rbx,%rdi
xor    %esi,%esi
call   1040 <printf@plt>
movss  0x14(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x1,%esi
call   1040 <printf@plt>
movss  0x18(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x2,%esi
call   1040 <printf@plt>
movss  0x1c(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x3,%esi
call   1040 <printf@plt>
movss  0x20(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x4,%esi
call   1040 <printf@plt>
movss  0x24(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x5,%esi
call   1040 <printf@plt>
movss  0x28(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x6,%esi
call   1040 <printf@plt>
movss  0x2c(%rsp),%xmm0
mov    %rbx,%rdi
mov    $0x7,%esi
call   1040 <printf@plt>
xor    %eax,%eax
add    $0x70,%rsp
pop    %rbx
ret

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

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