#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 */

/* Variables and functions */
scalar_t__ frexp(float, int *);
void stbiw__linear_to_rgbe(unsigned char *rgbe, float *linear) {
  int exponent;
  float maxcomp =
      ((linear[0]) > (((linear[1]) > (linear[2]) ? (linear[1]) : (linear[2])))
           ? (linear[0])
           : (((linear[1]) > (linear[2]) ? (linear[1]) : (linear[2]))));
  if (maxcomp < 1e-32) {
    rgbe[0] = rgbe[1] = rgbe[2] = rgbe[3] = 0;
  } else {
    float normalize = (float)frexp(maxcomp, &exponent) * 256.0f / maxcomp;
    rgbe[0] = (unsigned char)(linear[0] * normalize);
    rgbe[1] = (unsigned char)(linear[1] * normalize);
    rgbe[2] = (unsigned char)(linear[2] * normalize);
    rgbe[3] = (unsigned char)(exponent + 128);
  }
}


x86 assembly (clang, -O0):
<stbiw__linear_to_rgbe>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    -0x10(%rbp),%rax
ucomiss 0x8(%rax),%xmm0
jbe    1157 <stbiw__linear_to_rgbe+0x47>
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm0
movss  %xmm0,-0x24(%rbp)
jmp    1165 <stbiw__linear_to_rgbe+0x55>
mov    -0x10(%rbp),%rax
movss  0x8(%rax),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x20(%rbp),%xmm0
movss  -0x24(%rbp),%xmm1
ucomiss %xmm1,%xmm0
jbe    118a <stbiw__linear_to_rgbe+0x7a>
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x28(%rbp)
jmp    11cc <stbiw__linear_to_rgbe+0xbc>
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    -0x10(%rbp),%rax
ucomiss 0x8(%rax),%xmm0
jbe    11b4 <stbiw__linear_to_rgbe+0xa4>
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm0
movss  %xmm0,-0x2c(%rbp)
jmp    11c2 <stbiw__linear_to_rgbe+0xb2>
mov    -0x10(%rbp),%rax
movss  0x8(%rax),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x18(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe19(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1215 <stbiw__linear_to_rgbe+0x105>
mov    -0x8(%rbp),%rax
movb   $0x0,0x3(%rax)
mov    -0x8(%rbp),%rax
movb   $0x0,0x2(%rax)
mov    -0x8(%rbp),%rax
movb   $0x0,0x1(%rax)
mov    -0x8(%rbp),%rax
movb   $0x0,(%rax)
jmp    12ae <stbiw__linear_to_rgbe+0x19e>
movss  -0x18(%rbp),%xmm0
lea    -0x14(%rbp),%rdi
call   1030 <frexp@plt>
cvtsi2ss %rax,%xmm0
movss  0xdd8(%rip),%xmm1
mulss  %xmm1,%xmm0
movss  -0x18(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x1c(%rbp)
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
movss  -0x1c(%rbp),%xmm1
mulss  %xmm1,%xmm0
cvttss2si %xmm0,%eax
mov    %al,%cl
mov    -0x8(%rbp),%rax
mov    %cl,(%rax)
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm0
movss  -0x1c(%rbp),%xmm1
mulss  %xmm1,%xmm0
cvttss2si %xmm0,%eax
mov    %al,%cl
mov    -0x8(%rbp),%rax
mov    %cl,0x1(%rax)
mov    -0x10(%rbp),%rax
movss  0x8(%rax),%xmm0
movss  -0x1c(%rbp),%xmm1
mulss  %xmm1,%xmm0
cvttss2si %xmm0,%eax
mov    %al,%cl
mov    -0x8(%rbp),%rax
mov    %cl,0x2(%rax)
mov    -0x14(%rbp),%eax
add    $0x80,%eax
mov    %al,%cl
mov    -0x8(%rbp),%rax
mov    %cl,0x3(%rax)
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<stbiw__linear_to_rgbe>:
push   %r14
push   %rbx
push   %rax
mov    %rdi,%r14
movss  (%rsi),%xmm0
movss  0x4(%rsi),%xmm1
maxss  0x8(%rsi),%xmm1
maxss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm2
movsd  0xecb(%rip),%xmm1
ucomisd %xmm2,%xmm1
jbe    1147 <stbiw__linear_to_rgbe+0x37>
movw   $0x0,(%r14)
xor    %eax,%eax
xor    %ecx,%ecx
jmp    119d <stbiw__linear_to_rgbe+0x8d>
mov    %rsi,%rbx
lea    0x4(%rsp),%rdi
movss  %xmm0,(%rsp)
call   1030 <frexp@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %rax,%xmm0
mulss  0xe9f(%rip),%xmm0
divss  (%rsp),%xmm0
movss  (%rbx),%xmm1
mulss  %xmm0,%xmm1
cvttss2si %xmm1,%eax
mov    %al,(%r14)
movss  0x4(%rbx),%xmm1
mulss  %xmm0,%xmm1
cvttss2si %xmm1,%eax
mov    %al,0x1(%r14)
mulss  0x8(%rbx),%xmm0
cvttss2si %xmm0,%ecx
mov    0x4(%rsp),%al
xor    $0x80,%al
mov    %al,0x3(%r14)
mov    %cl,0x2(%r14)
add    $0x8,%rsp
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @stbiw__linear_to_rgbe(i8* %0, float* %1) {
%3 = alloca i8*
%4 = alloca float*
%5 = alloca i32
%6 = alloca float
%7 = alloca float
store i8* %0, i8** %3
store float* %1, float** %4
%8 = load float*, float** %4
%9 = getelementptr inbounds float, float* %8, i64 0
%10 = load float, float* %9
%11 = load float*, float** %4
%12 = getelementptr inbounds float, float* %11, i64 1
%13 = load float, float* %12
%14 = load float*, float** %4
%15 = getelementptr inbounds float, float* %14, i64 2
%16 = load float, float* %15
%17 = fcmp ogt float %13, %16
br i1 %17, label %18, label %22
18: ; preds = %2
%19 = load float*, float** %4
%20 = getelementptr inbounds float, float* %19, i64 1
%21 = load float, float* %20
br label %26
22: ; preds = %2
%23 = load float*, float** %4
%24 = getelementptr inbounds float, float* %23, i64 2
%25 = load float, float* %24
br label %26
26: ; preds = %22, %18
%27 = phi float [ %21, %18 ], [ %25, %22 ]
%28 = fcmp ogt float %10, %27
br i1 %28, label %29, label %33
29: ; preds = %26
%30 = load float*, float** %4
%31 = getelementptr inbounds float, float* %30, i64 0
%32 = load float, float* %31
br label %51
33: ; preds = %26
%34 = load float*, float** %4
%35 = getelementptr inbounds float, float* %34, i64 1
%36 = load float, float* %35
%37 = load float*, float** %4
%38 = getelementptr inbounds float, float* %37, i64 2
%39 = load float, float* %38
%40 = fcmp ogt float %36, %39
br i1 %40, label %41, label %45
41: ; preds = %33
%42 = load float*, float** %4
%43 = getelementptr inbounds float, float* %42, i64 1
%44 = load float, float* %43
br label %49
45: ; preds = %33
%46 = load float*, float** %4
%47 = getelementptr inbounds float, float* %46, i64 2
%48 = load float, float* %47
br label %49
49: ; preds = %45, %41
%50 = phi float [ %44, %41 ], [ %48, %45 ]
br label %51
51: ; preds = %49, %29
%52 = phi float [ %32, %29 ], [ %50, %49 ]
store float %52, float* %6
%53 = load float, float* %6
%54 = fpext float %53 to double
%55 = fcmp olt double %54, 1.000000e-32
br i1 %55, label %56, label %65
56: ; preds = %51
%57 = load i8*, i8** %3
%58 = getelementptr inbounds i8, i8* %57, i64 3
store i8 0, i8* %58
%59 = load i8*, i8** %3
%60 = getelementptr inbounds i8, i8* %59, i64 2
store i8 0, i8* %60
%61 = load i8*, i8** %3
%62 = getelementptr inbounds i8, i8* %61, i64 1
store i8 0, i8* %62
%63 = load i8*, i8** %3
%64 = getelementptr inbounds i8, i8* %63, i64 0
store i8 0, i8* %64
br label %101
65: ; preds = %51
%66 = load float, float* %6
%67 = call i64 @frexp(float %66, i32* %5)
%68 = sitofp i64 %67 to float
%69 = fmul float %68, 2.560000e+02
%70 = load float, float* %6
%71 = fdiv float %69, %70
store float %71, float* %7
%72 = load float*, float** %4
%73 = getelementptr inbounds float, float* %72, i64 0
%74 = load float, float* %73
%75 = load float, float* %7
%76 = fmul float %74, %75
%77 = fptoui float %76 to i8
%78 = load i8*, i8** %3
%79 = getelementptr inbounds i8, i8* %78, i64 0
store i8 %77, i8* %79
%80 = load float*, float** %4
%81 = getelementptr inbounds float, float* %80, i64 1
%82 = load float, float* %81
%83 = load float, float* %7
%84 = fmul float %82, %83
%85 = fptoui float %84 to i8
%86 = load i8*, i8** %3
%87 = getelementptr inbounds i8, i8* %86, i64 1
store i8 %85, i8* %87
%88 = load float*, float** %4
%89 = getelementptr inbounds float, float* %88, i64 2
%90 = load float, float* %89
%91 = load float, float* %7
%92 = fmul float %90, %91
%93 = fptoui float %92 to i8
%94 = load i8*, i8** %3
%95 = getelementptr inbounds i8, i8* %94, i64 2
store i8 %93, i8* %95
%96 = load i32, i32* %5
%97 = add nsw i32 %96, 128
%98 = trunc i32 %97 to i8
%99 = load i8*, i8** %3
%100 = getelementptr inbounds i8, i8* %99, i64 3
store i8 %98, i8* %100
br label %101
101: ; preds = %65, %56
ret void
}

LLVM IR (-O2):
define void @stbiw__linear_to_rgbe(i8* nocapture writeonly %0, float* nocapture %1) {
%3 = alloca i32
%4 = bitcast i32* %3 to i8*
%5 = load float, float* %1
%6 = getelementptr inbounds float, float* %1, i64 1
%7 = load float, float* %6
%8 = getelementptr inbounds float, float* %1, i64 2
%9 = load float, float* %8
%10 = fcmp ogt float %7, %9
%11 = select i1 %10, float %7, float %9
%12 = fcmp ogt float %5, %11
%13 = select i1 %12, float %5, float %11
%14 = fpext float %13 to double
%15 = fcmp olt double %14, 1.000000e-32
br i1 %15, label %16, label %18
16: ; preds = %2
%17 = getelementptr inbounds i8, i8* %0, i64 1
store i8 0, i8* %17
store i8 0, i8* %0
br label %36
18: ; preds = %2
%19 = call i64 @frexp(float %13, i32* %3)
%20 = sitofp i64 %19 to float
%21 = fmul float %20, 2.560000e+02
%22 = fdiv float %21, %13
%23 = load float, float* %1
%24 = fmul float %23, %22
%25 = fptoui float %24 to i8
store i8 %25, i8* %0
%26 = load float, float* %6
%27 = fmul float %26, %22
%28 = fptoui float %27 to i8
%29 = getelementptr inbounds i8, i8* %0, i64 1
store i8 %28, i8* %29
%30 = load float, float* %8
%31 = fmul float %22, %30
%32 = fptoui float %31 to i8
%33 = load i32, i32* %3
%34 = trunc i32 %33 to i8
%35 = xor i8 %34, -128
br label %36
36: ; preds = %18, %16
%37 = phi i8 [ 0, %16 ], [ %35, %18 ]
%38 = phi i8 [ 0, %16 ], [ %32, %18 ]
%39 = getelementptr inbounds i8, i8* %0, i64 3
store i8 %37, i8* %39
%40 = getelementptr inbounds i8, i8* %0, i64 2
store i8 %38, i8* %40
ret void
}