int wavezopt(restwav, nwave, swave, szopt, wave, zopt, stat)
double restwav;
int nwave, swave, szopt;
const double wave[];
double zopt[];
int stat[];
{
  double r;
  register int iwave, *statp;
  register const double *wavep;
  register double *zoptp;
  if (restwav == 0.0) {
    return 2;
  }
  r = 1.0 / restwav;
  wavep = wave;
  zoptp = zopt;
  statp = stat;
  for (iwave = 0; iwave < nwave; iwave++) {
    *zoptp = r * (*wavep) - 1.0;
    *(statp++) = 0;
    wavep += swave;
    zoptp += szopt;
  }
  return 0;
}


x86 assembly (clang, -O0):
<wavezopt>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
mov    %edi,-0x14(%rbp)
mov    %esi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1142 <wavezopt+0x42>
jp     1142 <wavezopt+0x42>
movl   $0x2,-0x4(%rbp)
jmp    11f8 <wavezopt+0xf8>
movsd  0xeb6(%rip),%xmm0
divsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x28(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x30(%rbp),%rax
mov    %rax,-0x60(%rbp)
mov    -0x38(%rbp),%rax
mov    %rax,-0x50(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11f1 <wavezopt+0xf1>
movsd  -0x40(%rbp),%xmm0
mov    -0x58(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xe70(%rip),%xmm1
addsd  %xmm1,%xmm0
mov    -0x60(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x50(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x50(%rbp)
movl   $0x0,(%rax)
mov    -0x18(%rbp),%ecx
mov    -0x58(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x1c(%rbp),%ecx
mov    -0x60(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x60(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1173 <wavezopt+0x73>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<wavezopt>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
mov    $0x2,%eax
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1122 <wavezopt+0x22>
jp     1122 <wavezopt+0x22>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %eax,%eax
test   %edi,%edi
jle    1119 <wavezopt+0x19>
movsd  0xed0(%rip),%xmm1
divsd  %xmm0,%xmm1
movslq %esi,%r14
movslq %edx,%r15
lea    -0x1(%rdi),%edx
cmp    $0x3,%edx
jb     1221 <wavezopt+0x121>
mov    %edi,%r12d
and    $0xfffffffc,%r12d
mov    %r15,%r10
shl    $0x5,%r10
mov    %r14,%r11
shl    $0x5,%r11
lea    0x0(,%r14,8),%rsi
lea    0x0(,%r15,8),%rdx
movsd  0xe95(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r8,%rbx
mov    %rcx,%rbp
movsd  (%rcx),%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
movsd  %xmm2,(%r8)
movl   $0x0,(%r9)
add    %rsi,%rcx
movsd  0x0(%rbp,%rsi,1),%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
add    %rdx,%r8
movsd  %xmm2,(%rbx,%rdx,1)
movl   $0x0,0x4(%r9)
movsd  (%rsi,%rcx,1),%xmm2
mulsd  %xmm1,%xmm2
add    %rsi,%rcx
addsd  %xmm0,%xmm2
movsd  %xmm2,(%rdx,%r8,1)
add    %rdx,%r8
movl   $0x0,0x8(%r9)
movsd  (%rsi,%rcx,1),%xmm2
mulsd  %xmm1,%xmm2
add    %rsi,%rcx
addsd  %xmm0,%xmm2
movsd  %xmm2,(%rdx,%r8,1)
add    %rdx,%r8
movl   $0x0,0xc(%r9)
add    $0x10,%r9
add    %rsi,%rcx
add    %rdx,%r8
add    $0xfffffffc,%r12d
jne    1180 <wavezopt+0x80>
add    %r10,%rbx
add    %r11,%rbp
mov    %rbx,%r8
mov    %rbp,%rcx
test   $0x3,%dil
je     1119 <wavezopt+0x19>
shl    $0x3,%r15
shl    $0x3,%r14
and    $0x3,%edi
xor    %eax,%eax
movsd  0xdc8(%rip),%xmm0
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rcx),%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
movsd  %xmm2,(%r8)
movl   $0x0,(%r9,%rdx,4)
add    %r15,%r8
add    %r14,%rcx
add    $0x1,%rdx
cmp    %edx,%edi
jne    1250 <wavezopt+0x150>
jmp    1119 <wavezopt+0x19>

LLVM IR (-O0):
define i32 @wavezopt(double %0, i32 %1, i32 %2, i32 %3, double* %4, double* %5, i32* %6) {
%8 = alloca i32
%9 = alloca double
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca double*
%14 = alloca double*
%15 = alloca i32*
%16 = alloca double
%17 = alloca i32
%18 = alloca i32*
%19 = alloca double*
%20 = alloca double*
store double %0, double* %9
store i32 %1, i32* %10
store i32 %2, i32* %11
store i32 %3, i32* %12
store double* %4, double** %13
store double* %5, double** %14
store i32* %6, i32** %15
%21 = load double, double* %9
%22 = fcmp oeq double %21, 0.000000e+00
br i1 %22, label %23, label %24
23: ; preds = %7
store i32 2, i32* %8
br label %54
24: ; preds = %7
%25 = load double, double* %9
%26 = fdiv double 1.000000e+00, %25
store double %26, double* %16
%27 = load double*, double** %13
store double* %27, double** %19
%28 = load double*, double** %14
store double* %28, double** %20
%29 = load i32*, i32** %15
store i32* %29, i32** %18
store i32 0, i32* %17
br label %30
30: ; preds = %50, %24
%31 = load i32, i32* %17
%32 = load i32, i32* %10
%33 = icmp slt i32 %31, %32
br i1 %33, label %34, label %53
34: ; preds = %30
%35 = load double, double* %16
%36 = load double*, double** %19
%37 = load double, double* %36
%38 = call double @llvm.fmuladd.f64(double %35, double %37, double -1.000000e+00)
%39 = load double*, double** %20
store double %38, double* %39
%40 = load i32*, i32** %18
%41 = getelementptr inbounds i32, i32* %40, i32 1
store i32* %41, i32** %18
store i32 0, i32* %40
%42 = load i32, i32* %11
%43 = load double*, double** %19
%44 = sext i32 %42 to i64
%45 = getelementptr inbounds double, double* %43, i64 %44
store double* %45, double** %19
%46 = load i32, i32* %12
%47 = load double*, double** %20
%48 = sext i32 %46 to i64
%49 = getelementptr inbounds double, double* %47, i64 %48
store double* %49, double** %20
br label %50
50: ; preds = %34
%51 = load i32, i32* %17
%52 = add nsw i32 %51, 1
store i32 %52, i32* %17
br label %30
53: ; preds = %30
store i32 0, i32* %8
br label %54
54: ; preds = %53, %23
%55 = load i32, i32* %8
ret i32 %55
}

LLVM IR (-O2):
define i32 @wavezopt(double %0, i32 %1, i32 %2, i32 %3, double* nocapture %4, double* nocapture writeonly %5, i32* nocapture writeonly %6) {
%8 = fcmp oeq double %0, 0.000000e+00
br i1 %8, label %64, label %9
9: ; preds = %7
%10 = fdiv double 1.000000e+00, %0
%11 = icmp sgt i32 %1, 0
br i1 %11, label %12, label %64
12: ; preds = %9
%13 = sext i32 %2 to i64
%14 = sext i32 %3 to i64
%15 = add i32 %1, -1
%16 = and i32 %1, 3
%17 = icmp ult i32 %15, 3
br i1 %17, label %47, label %18
18: ; preds = %12
%19 = and i32 %1, -4
br label %20
20: ; preds = %20, %18
%21 = phi double* [ %5, %18 ], [ %44, %20 ]
%22 = phi double* [ %4, %18 ], [ %43, %20 ]
%23 = phi i32* [ %6, %18 ], [ %42, %20 ]
%24 = phi i32 [ 0, %18 ], [ %45, %20 ]
%25 = load double, double* %22
%26 = tail call double @llvm.fmuladd.f64(double %10, double %25, double -1.000000e+00)
store double %26, double* %21
%27 = getelementptr inbounds i32, i32* %23, i64 1
store i32 0, i32* %23
%28 = getelementptr inbounds double, double* %22, i64 %13
%29 = getelementptr inbounds double, double* %21, i64 %14
%30 = load double, double* %28
%31 = tail call double @llvm.fmuladd.f64(double %10, double %30, double -1.000000e+00)
store double %31, double* %29
%32 = getelementptr inbounds i32, i32* %23, i64 2
store i32 0, i32* %27
%33 = getelementptr inbounds double, double* %28, i64 %13
%34 = getelementptr inbounds double, double* %29, i64 %14
%35 = load double, double* %33
%36 = tail call double @llvm.fmuladd.f64(double %10, double %35, double -1.000000e+00)
store double %36, double* %34
%37 = getelementptr inbounds i32, i32* %23, i64 3
store i32 0, i32* %32
%38 = getelementptr inbounds double, double* %33, i64 %13
%39 = getelementptr inbounds double, double* %34, i64 %14
%40 = load double, double* %38
%41 = tail call double @llvm.fmuladd.f64(double %10, double %40, double -1.000000e+00)
store double %41, double* %39
%42 = getelementptr inbounds i32, i32* %23, i64 4
store i32 0, i32* %37
%43 = getelementptr inbounds double, double* %38, i64 %13
%44 = getelementptr inbounds double, double* %39, i64 %14
%45 = add i32 %24, 4
%46 = icmp eq i32 %45, %19
br i1 %46, label %47, label %20
47: ; preds = %20, %12
%48 = phi double* [ %5, %12 ], [ %44, %20 ]
%49 = phi double* [ %4, %12 ], [ %43, %20 ]
%50 = phi i32* [ %6, %12 ], [ %42, %20 ]
%51 = icmp eq i32 %16, 0
br i1 %51, label %64, label %52
52: ; preds = %47, %52
%53 = phi double* [ %61, %52 ], [ %48, %47 ]
%54 = phi double* [ %60, %52 ], [ %49, %47 ]
%55 = phi i32* [ %59, %52 ], [ %50, %47 ]
%56 = phi i32 [ %62, %52 ], [ 0, %47 ]
%57 = load double, double* %54
%58 = tail call double @llvm.fmuladd.f64(double %10, double %57, double -1.000000e+00)
store double %58, double* %53
%59 = getelementptr inbounds i32, i32* %55, i64 1
store i32 0, i32* %55
%60 = getelementptr inbounds double, double* %54, i64 %13
%61 = getelementptr inbounds double, double* %53, i64 %14
%62 = add i32 %56, 1
%63 = icmp eq i32 %62, %16
br i1 %63, label %64, label %52
64: ; preds = %47, %52, %9, %7
%65 = phi i32 [ 2, %7 ], [ 0, %9 ], [ 0, %52 ], [ 0, %47 ]
ret i32 %65
}