#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 scalar_t__ int64_t;

/* Variables and functions */
double fabs(double const);
scalar_t__ stdlib_base_is_nan(double);
double stdlib_strided_dnanmaxabs(const int64_t N, const double *X,
                                 const int64_t stride) {
  double max;
  int64_t ix;
  int64_t i;
  double v;
  if (N <= 0) {
    return 0.0 / 0.0;
  }
  if (N == 1 || stride == 0) {
    return fabs(X[0]);
  }
  if (stride < 0) {
    ix = (1 - N) * stride;
  } else {
    ix = 0;
  }
  for (i = 0; i < N; i++) {
    v = X[ix];
    if (v == v) {
      break;
    }
    ix += stride;
  }
  if (i == N) {
    return 0.0 / 0.0;
  }
  max = fabs(v);
  i += 1;
  for (; i < N; i++) {
    ix += stride;
    v = fabs(X[ix]);
    if (stdlib_base_is_nan(v)) {
      continue;
    }
    if (v > max) {
      max = v;
    }
  }
  return max;
}


x86 assembly (clang, -O0):
<stdlib_strided_dnanmaxabs>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
cmpq   $0x0,-0x10(%rbp)
jg     1141 <stdlib_strided_dnanmaxabs+0x31>
movsd  0xed9(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    12c9 <stdlib_strided_dnanmaxabs+0x1b9>
cmpq   $0x1,-0x10(%rbp)
je     1157 <stdlib_strided_dnanmaxabs+0x47>
cmpq   $0x0,-0x20(%rbp)
jne    1174 <stdlib_strided_dnanmaxabs+0x64>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movaps 0xe9a(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    12c9 <stdlib_strided_dnanmaxabs+0x1b9>
cmpq   $0x0,-0x20(%rbp)
jge    1196 <stdlib_strided_dnanmaxabs+0x86>
mov    $0x1,%eax
sub    -0x10(%rbp),%rax
imul   -0x20(%rbp),%rax
mov    %rax,-0x30(%rbp)
jmp    119e <stdlib_strided_dnanmaxabs+0x8e>
movq   $0x0,-0x30(%rbp)
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    11fe <stdlib_strided_dnanmaxabs+0xee>
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
ucomisd -0x40(%rbp),%xmm0
jne    11e1 <stdlib_strided_dnanmaxabs+0xd1>
jp     11e1 <stdlib_strided_dnanmaxabs+0xd1>
jmp    11fe <stdlib_strided_dnanmaxabs+0xee>
mov    -0x20(%rbp),%rax
add    -0x30(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11a6 <stdlib_strided_dnanmaxabs+0x96>
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jne    121e <stdlib_strided_dnanmaxabs+0x10e>
movsd  0xdfc(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    12c9 <stdlib_strided_dnanmaxabs+0x1b9>
movsd  -0x40(%rbp),%xmm0
movaps 0xdd6(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    12bf <stdlib_strided_dnanmaxabs+0x1af>
mov    -0x20(%rbp),%rax
add    -0x30(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xd93(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
call   1030 <stdlib_base_is_nan@plt>
cmp    $0x0,%rax
je     128f <stdlib_strided_dnanmaxabs+0x17f>
jmp    12ae <stdlib_strided_dnanmaxabs+0x19e>
movsd  -0x40(%rbp),%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    12a9 <stdlib_strided_dnanmaxabs+0x199>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    12ae <stdlib_strided_dnanmaxabs+0x19e>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    123f <stdlib_strided_dnanmaxabs+0x12f>
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<stdlib_strided_dnanmaxabs>:
push   %r15
push   %r14
push   %rbx
sub    $0x20,%rsp
movsd  0xedf(%rip),%xmm0
test   %rdi,%rdi
jle    118a <stdlib_strided_dnanmaxabs+0x7a>
cmp    $0x1,%rdi
je     1198 <stdlib_strided_dnanmaxabs+0x88>
mov    %rdx,%r14
test   %rdx,%rdx
je     1198 <stdlib_strided_dnanmaxabs+0x88>
movsd  %xmm0,0x8(%rsp)
mov    $0x1,%eax
mov    $0x1,%ecx
sub    %rdi,%rcx
imul   %r14,%rcx
mov    %r14,%rdx
sar    $0x3f,%rdx
and    %rcx,%rdx
lea    -0x1(%rdi),%rbx
lea    (%rdx,%r14,1),%rcx
shl    $0x3,%r14
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rsi,%rdx,8),%xmm0
ucomisd %xmm0,%xmm0
jnp    11ad <stdlib_strided_dnanmaxabs+0x9d>
add    %r14,%rsi
add    $0x1,%rax
add    $0xffffffffffffffff,%rbx
jb     1170 <stdlib_strided_dnanmaxabs+0x60>
jmp    1190 <stdlib_strided_dnanmaxabs+0x80>
movsd  %xmm0,0x8(%rsp)
movsd  0x8(%rsp),%xmm0
jmp    11a3 <stdlib_strided_dnanmaxabs+0x93>
movsd  (%rsi),%xmm0
andps  0xe6d(%rip),%xmm0
add    $0x20,%rsp
pop    %rbx
pop    %r14
pop    %r15
ret
cmp    $0xffffffffffffffff,%rbx
je     1190 <stdlib_strided_dnanmaxabs+0x80>
andpd  0xe55(%rip),%xmm0
movsd  %xmm0,0x8(%rsp)
cmp    %rdi,%rax
jge    1190 <stdlib_strided_dnanmaxabs+0x80>
lea    (%rsi,%rcx,8),%r15
jmp    11d9 <stdlib_strided_dnanmaxabs+0xc9>
nopl   0x0(%rax)
add    %r14,%r15
add    $0xffffffffffffffff,%rbx
je     1190 <stdlib_strided_dnanmaxabs+0x80>
movsd  (%r15),%xmm0
andpd  0xe2a(%rip),%xmm0
movapd %xmm0,0x10(%rsp)
call   1030 <stdlib_base_is_nan@plt>
test   %rax,%rax
jne    11d0 <stdlib_strided_dnanmaxabs+0xc0>
movapd 0x10(%rsp),%xmm0
maxsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
jmp    11d0 <stdlib_strided_dnanmaxabs+0xc0>

LLVM IR (-O0):
define double @stdlib_strided_dnanmaxabs(i64 %0, double* %1, i64 %2) {
%4 = alloca double
%5 = alloca i64
%6 = alloca double*
%7 = alloca i64
%8 = alloca double
%9 = alloca i64
%10 = alloca i64
%11 = alloca double
store i64 %0, i64* %5
store double* %1, double** %6
store i64 %2, i64* %7
%12 = load i64, i64* %5
%13 = icmp sle i64 %12, 0
br i1 %13, label %14, label %15
14: ; preds = %3
store double 0x7FF8000000000000, double* %4
br label %95
15: ; preds = %3
%16 = load i64, i64* %5
%17 = icmp eq i64 %16, 1
br i1 %17, label %21, label %18
18: ; preds = %15
%19 = load i64, i64* %7
%20 = icmp eq i64 %19, 0
br i1 %20, label %21, label %26
21: ; preds = %18, %15
%22 = load double*, double** %6
%23 = getelementptr inbounds double, double* %22, i64 0
%24 = load double, double* %23
%25 = call double @llvm.fabs.f64(double %24)
store double %25, double* %4
br label %95
26: ; preds = %18
%27 = load i64, i64* %7
%28 = icmp slt i64 %27, 0
br i1 %28, label %29, label %34
29: ; preds = %26
%30 = load i64, i64* %5
%31 = sub nsw i64 1, %30
%32 = load i64, i64* %7
%33 = mul nsw i64 %31, %32
store i64 %33, i64* %9
br label %35
34: ; preds = %26
store i64 0, i64* %9
br label %35
35: ; preds = %34, %29
store i64 0, i64* %10
br label %36
36: ; preds = %53, %35
%37 = load i64, i64* %10
%38 = load i64, i64* %5
%39 = icmp slt i64 %37, %38
br i1 %39, label %40, label %56
40: ; preds = %36
%41 = load double*, double** %6
%42 = load i64, i64* %9
%43 = getelementptr inbounds double, double* %41, i64 %42
%44 = load double, double* %43
store double %44, double* %11
%45 = load double, double* %11
%46 = load double, double* %11
%47 = fcmp oeq double %45, %46
br i1 %47, label %48, label %49
48: ; preds = %40
br label %56
49: ; preds = %40
%50 = load i64, i64* %7
%51 = load i64, i64* %9
%52 = add nsw i64 %51, %50
store i64 %52, i64* %9
br label %53
53: ; preds = %49
%54 = load i64, i64* %10
%55 = add nsw i64 %54, 1
store i64 %55, i64* %10
br label %36
56: ; preds = %48, %36
%57 = load i64, i64* %10
%58 = load i64, i64* %5
%59 = icmp eq i64 %57, %58
br i1 %59, label %60, label %61
60: ; preds = %56
store double 0x7FF8000000000000, double* %4
br label %95
61: ; preds = %56
%62 = load double, double* %11
%63 = call double @llvm.fabs.f64(double %62)
store double %63, double* %8
%64 = load i64, i64* %10
%65 = add nsw i64 %64, 1
store i64 %65, i64* %10
br label %66
66: ; preds = %90, %61
%67 = load i64, i64* %10
%68 = load i64, i64* %5
%69 = icmp slt i64 %67, %68
br i1 %69, label %70, label %93
70: ; preds = %66
%71 = load i64, i64* %7
%72 = load i64, i64* %9
%73 = add nsw i64 %72, %71
store i64 %73, i64* %9
%74 = load double*, double** %6
%75 = load i64, i64* %9
%76 = getelementptr inbounds double, double* %74, i64 %75
%77 = load double, double* %76
%78 = call double @llvm.fabs.f64(double %77)
store double %78, double* %11
%79 = load double, double* %11
%80 = call i64 @stdlib_base_is_nan(double %79)
%81 = icmp ne i64 %80, 0
br i1 %81, label %82, label %83
82: ; preds = %70
br label %90
83: ; preds = %70
%84 = load double, double* %11
%85 = load double, double* %8
%86 = fcmp ogt double %84, %85
br i1 %86, label %87, label %89
87: ; preds = %83
%88 = load double, double* %11
store double %88, double* %8
br label %89
89: ; preds = %87, %83
br label %90
90: ; preds = %89, %82
%91 = load i64, i64* %10
%92 = add nsw i64 %91, 1
store i64 %92, i64* %10
br label %66
93: ; preds = %66
%94 = load double, double* %8
store double %94, double* %4
br label %95
95: ; preds = %93, %60, %21, %14
%96 = load double, double* %4
ret double %96
}

LLVM IR (-O2):
define double @stdlib_strided_dnanmaxabs(i64 %0, double* nocapture %1, i64 %2) {
%4 = icmp slt i64 %0, 1
br i1 %4, label %48, label %5
5: ; preds = %3
%6 = icmp eq i64 %0, 1
%7 = icmp eq i64 %2, 0
%8 = or i1 %6, %7
br i1 %8, label %9, label %12
9: ; preds = %5
%10 = load double, double* %1
%11 = tail call double @llvm.fabs.f64(double %10)
br label %48
12: ; preds = %5
%13 = icmp slt i64 %2, 0
%14 = sub nsw i64 1, %0
%15 = mul nsw i64 %14, %2
%16 = select i1 %13, i64 %15, i64 0
br label %17
17: ; preds = %12, %23
%18 = phi i64 [ %25, %23 ], [ 0, %12 ]
%19 = phi i64 [ %24, %23 ], [ %16, %12 ]
%20 = getelementptr inbounds double, double* %1, i64 %19
%21 = load double, double* %20
%22 = fcmp ord double %21, 0.000000e+00
br i1 %22, label %27, label %23
23: ; preds = %17
%24 = add nsw i64 %19, %2
%25 = add nuw nsw i64 %18, 1
%26 = icmp eq i64 %25, %0
br i1 %26, label %48, label %17
27: ; preds = %17
%28 = icmp eq i64 %18, %0
br i1 %28, label %48, label %29
29: ; preds = %27
%30 = tail call double @llvm.fabs.f64(double %21)
%31 = add nuw nsw i64 %18, 1
%32 = icmp slt i64 %31, %0
br i1 %32, label %33, label %48
33: ; preds = %29, %33
%34 = phi i64 [ %46, %33 ], [ %31, %29 ]
%35 = phi i64 [ %37, %33 ], [ %19, %29 ]
%36 = phi double [ %45, %33 ], [ %30, %29 ]
%37 = add nsw i64 %35, %2
%38 = getelementptr inbounds double, double* %1, i64 %37
%39 = load double, double* %38
%40 = tail call double @llvm.fabs.f64(double %39)
%41 = tail call i64 @stdlib_base_is_nan(double %40)
%42 = icmp eq i64 %41, 0
%43 = fcmp ogt double %40, %36
%44 = select i1 %42, i1 %43, i1 false
%45 = select i1 %44, double %40, double %36
%46 = add i64 %34, 1
%47 = icmp eq i64 %46, %0
br i1 %47, label %48, label %33
48: ; preds = %23, %33, %29, %27, %3, %9
%49 = phi double [ %11, %9 ], [ 0x7FF8000000000000, %3 ], [ 0x7FF8000000000000, %27 ], [ %30, %29 ], [ %45, %33 ], [ 0x7FF8000000000000, %23 ]
ret double %49
}