int freqvrad(restfrq, nfreq, sfreq, svrad, freq, vrad, stat)
double restfrq;
int nfreq, sfreq, svrad;
const double freq[];
double vrad[];
int stat[];
{
  double r;
  register int ifreq, *statp;
  register const double *freqp;
  register double *vradp;
  if (restfrq == 0.0) {
    return 2;
  }
  r = 2.99792458e8 / restfrq;
  freqp = freq;
  vradp = vrad;
  statp = stat;
  for (ifreq = 0; ifreq < nfreq; ifreq++) {
    *vradp = r * (restfrq - *freqp);
    *(statp++) = 0;
    freqp += sfreq;
    vradp += svrad;
  }
  return 0;
}


x86 assembly (clang, -O0):
<freqvrad>:
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 <freqvrad+0x42>
jp     1142 <freqvrad+0x42>
movl   $0x2,-0x4(%rbp)
jmp    11f1 <freqvrad+0xf1>
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    11ea <freqvrad+0xea>
movsd  -0x40(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm1
mov    -0x58(%rbp),%rax
subsd  (%rax),%xmm1
mulsd  %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 <freqvrad+0x73>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<freqvrad>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
mov    $0x2,%eax
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1122 <freqvrad+0x22>
jp     1122 <freqvrad+0x22>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %eax,%eax
test   %edi,%edi
jle    1119 <freqvrad+0x19>
movsd  0xed0(%rip),%xmm1
divsd  %xmm0,%xmm1
movslq %esi,%r14
movslq %edx,%r15
lea    -0x1(%rdi),%edx
cmp    $0x3,%edx
jb     1211 <freqvrad+0x111>
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
nopl   0x0(%rax,%rax,1)
mov    %r8,%rbx
mov    %rcx,%rbp
movapd %xmm0,%xmm2
subsd  (%rcx),%xmm2
mulsd  %xmm1,%xmm2
movsd  %xmm2,(%r8)
movl   $0x0,(%r9)
add    %rsi,%rcx
movapd %xmm0,%xmm2
subsd  0x0(%rbp,%rsi,1),%xmm2
mulsd  %xmm1,%xmm2
add    %rdx,%r8
movsd  %xmm2,(%rbx,%rdx,1)
movl   $0x0,0x4(%r9)
movapd %xmm0,%xmm2
subsd  (%rsi,%rcx,1),%xmm2
add    %rsi,%rcx
mulsd  %xmm1,%xmm2
movsd  %xmm2,(%rdx,%r8,1)
add    %rdx,%r8
movl   $0x0,0x8(%r9)
movapd %xmm0,%xmm2
subsd  (%rsi,%rcx,1),%xmm2
add    %rsi,%rcx
mulsd  %xmm1,%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    1170 <freqvrad+0x70>
add    %r10,%rbx
add    %r11,%rbp
mov    %rbx,%r8
mov    %rbp,%rcx
test   $0x3,%dil
je     1119 <freqvrad+0x19>
shl    $0x3,%r15
shl    $0x3,%r14
and    $0x3,%edi
xor    %eax,%eax
xor    %edx,%edx
nopw   0x0(%rax,%rax,1)
movapd %xmm0,%xmm2
subsd  (%rcx),%xmm2
mulsd  %xmm1,%xmm2
movsd  %xmm2,(%r8)
movl   $0x0,(%r9,%rdx,4)
add    %r15,%r8
add    %r14,%rcx
add    $0x1,%rdx
cmp    %edx,%edi
jne    1230 <freqvrad+0x130>
jmp    1119 <freqvrad+0x19>

LLVM IR (-O0):
define i32 @freqvrad(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 %56
24: ; preds = %7
%25 = load double, double* %9
%26 = fdiv double 0x41B1DE784A000000, %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 = %52, %24
%31 = load i32, i32* %17
%32 = load i32, i32* %10
%33 = icmp slt i32 %31, %32
br i1 %33, label %34, label %55
34: ; preds = %30
%35 = load double, double* %16
%36 = load double, double* %9
%37 = load double*, double** %19
%38 = load double, double* %37
%39 = fsub double %36, %38
%40 = fmul double %35, %39
%41 = load double*, double** %20
store double %40, double* %41
%42 = load i32*, i32** %18
%43 = getelementptr inbounds i32, i32* %42, i32 1
store i32* %43, i32** %18
store i32 0, i32* %42
%44 = load i32, i32* %11
%45 = load double*, double** %19
%46 = sext i32 %44 to i64
%47 = getelementptr inbounds double, double* %45, i64 %46
store double* %47, double** %19
%48 = load i32, i32* %12
%49 = load double*, double** %20
%50 = sext i32 %48 to i64
%51 = getelementptr inbounds double, double* %49, i64 %50
store double* %51, double** %20
br label %52
52: ; preds = %34
%53 = load i32, i32* %17
%54 = add nsw i32 %53, 1
store i32 %54, i32* %17
br label %30
55: ; preds = %30
store i32 0, i32* %8
br label %56
56: ; preds = %55, %23
%57 = load i32, i32* %8
ret i32 %57
}

LLVM IR (-O2):
define i32 @freqvrad(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 %69, label %9
9: ; preds = %7
%10 = fdiv double 0x41B1DE784A000000, %0
%11 = icmp sgt i32 %1, 0
br i1 %11, label %12, label %69
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 %51, label %18
18: ; preds = %12
%19 = and i32 %1, -4
br label %20
20: ; preds = %20, %18
%21 = phi double* [ %5, %18 ], [ %48, %20 ]
%22 = phi double* [ %4, %18 ], [ %47, %20 ]
%23 = phi i32* [ %6, %18 ], [ %46, %20 ]
%24 = phi i32 [ 0, %18 ], [ %49, %20 ]
%25 = load double, double* %22
%26 = fsub double %0, %25
%27 = fmul double %10, %26
store double %27, double* %21
%28 = getelementptr inbounds i32, i32* %23, i64 1
store i32 0, i32* %23
%29 = getelementptr inbounds double, double* %22, i64 %13
%30 = getelementptr inbounds double, double* %21, i64 %14
%31 = load double, double* %29
%32 = fsub double %0, %31
%33 = fmul double %10, %32
store double %33, double* %30
%34 = getelementptr inbounds i32, i32* %23, i64 2
store i32 0, i32* %28
%35 = getelementptr inbounds double, double* %29, i64 %13
%36 = getelementptr inbounds double, double* %30, i64 %14
%37 = load double, double* %35
%38 = fsub double %0, %37
%39 = fmul double %10, %38
store double %39, double* %36
%40 = getelementptr inbounds i32, i32* %23, i64 3
store i32 0, i32* %34
%41 = getelementptr inbounds double, double* %35, i64 %13
%42 = getelementptr inbounds double, double* %36, i64 %14
%43 = load double, double* %41
%44 = fsub double %0, %43
%45 = fmul double %10, %44
store double %45, double* %42
%46 = getelementptr inbounds i32, i32* %23, i64 4
store i32 0, i32* %40
%47 = getelementptr inbounds double, double* %41, i64 %13
%48 = getelementptr inbounds double, double* %42, i64 %14
%49 = add i32 %24, 4
%50 = icmp eq i32 %49, %19
br i1 %50, label %51, label %20
51: ; preds = %20, %12
%52 = phi double* [ %5, %12 ], [ %48, %20 ]
%53 = phi double* [ %4, %12 ], [ %47, %20 ]
%54 = phi i32* [ %6, %12 ], [ %46, %20 ]
%55 = icmp eq i32 %16, 0
br i1 %55, label %69, label %56
56: ; preds = %51, %56
%57 = phi double* [ %66, %56 ], [ %52, %51 ]
%58 = phi double* [ %65, %56 ], [ %53, %51 ]
%59 = phi i32* [ %64, %56 ], [ %54, %51 ]
%60 = phi i32 [ %67, %56 ], [ 0, %51 ]
%61 = load double, double* %58
%62 = fsub double %0, %61
%63 = fmul double %10, %62
store double %63, double* %57
%64 = getelementptr inbounds i32, i32* %59, i64 1
store i32 0, i32* %59
%65 = getelementptr inbounds double, double* %58, i64 %13
%66 = getelementptr inbounds double, double* %57, i64 %14
%67 = add i32 %60, 1
%68 = icmp eq i32 %67, %16
br i1 %68, label %69, label %56
69: ; preds = %51, %56, %9, %7
%70 = phi i32 [ 2, %7 ], [ 0, %9 ], [ 0, %56 ], [ 0, %51 ]
ret i32 %70
}