void mant_exp(log_base, x, scientific, m, p) double log_base, x;
int scientific;
double *m;
int *p;
{
  int sign = 1;
  double l10;
  int power;
  if (x == 0) {
    if (m)
      *m = 0;
    if (p)
      *p = 0;
    return;
  }
  if (x < 0) {
    sign = (-1);
    x = (-x);
  }
  l10 = log10(x) / log_base;
  power = floor(l10);
  if (scientific) {
    power = 3 * floor(power / 3.0);
  }
  if (m)
    *m = sign * pow(10.0, (l10 - power) * log_base);
  if (p)
    *p = power;
}


x86 assembly (clang, -O0):
<mant_exp>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
mov    %edi,-0x14(%rbp)
mov    %rsi,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
movl   $0x1,-0x2c(%rbp)
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    119c <mant_exp+0x6c>
jp     119c <mant_exp+0x6c>
cmpq   $0x0,-0x20(%rbp)
je     1182 <mant_exp+0x52>
mov    -0x20(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
cmpq   $0x0,-0x28(%rbp)
je     1197 <mant_exp+0x67>
mov    -0x28(%rbp),%rax
movl   $0x0,(%rax)
jmp    128a <mant_exp+0x15a>
xorps  %xmm0,%xmm0
ucomisd -0x10(%rbp),%xmm0
jbe    11d2 <mant_exp+0xa2>
movl   $0xffffffff,-0x2c(%rbp)
movsd  -0x10(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x10(%rbp)
movsd  -0x10(%rbp),%xmm0
call   1040 <log10@plt>
divsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
call   1050 <floor@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x3c(%rbp)
cmpl   $0x0,-0x14(%rbp)
je     122d <mant_exp+0xfd>
cvtsi2sdl -0x3c(%rbp),%xmm0
movsd  0xdf2(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1050 <floor@plt>
movaps %xmm0,%xmm1
movsd  0xdde(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x3c(%rbp)
cmpq   $0x0,-0x20(%rbp)
je     1276 <mant_exp+0x146>
cvtsi2sdl -0x2c(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm1
cvtsi2sdl -0x3c(%rbp),%xmm0
subsd  %xmm0,%xmm1
mulsd  -0x8(%rbp),%xmm1
movsd  0xdab(%rip),%xmm0
call   1030 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0x48(%rbp),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
cmpq   $0x0,-0x28(%rbp)
je     128a <mant_exp+0x15a>
mov    -0x3c(%rbp),%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,(%rax)
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<mant_exp>:
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x28,%rsp
mov    %rdx,%r14
mov    %rsi,%r15
xorpd  %xmm3,%xmm3
ucomisd %xmm3,%xmm1
jne    115f <mant_exp+0x2f>
jp     115f <mant_exp+0x2f>
test   %r15,%r15
je     1158 <mant_exp+0x28>
movq   $0x0,(%r15)
xor    %ebx,%ebx
jmp    1222 <mant_exp+0xf2>
mov    %edi,%ebp
movapd 0xe97(%rip),%xmm2
xorpd  %xmm1,%xmm2
movsd  %xmm0,(%rsp)
movapd %xmm1,%xmm0
cmpltsd %xmm3,%xmm0
andpd  %xmm0,%xmm2
movapd %xmm1,0x10(%rsp)
andnpd %xmm1,%xmm0
orpd   %xmm2,%xmm0
call   1040 <log10@plt>
divsd  (%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
call   1050 <floor@plt>
cvttsd2si %xmm0,%ebx
test   %ebp,%ebp
je     11ca <mant_exp+0x9a>
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
divsd  0xe57(%rip),%xmm0
call   1050 <floor@plt>
mulsd  0xe4a(%rip),%xmm0
cvttsd2si %xmm0,%ebx
test   %r15,%r15
je     1222 <mant_exp+0xf2>
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
movsd  0x8(%rsp),%xmm1
subsd  %xmm0,%xmm1
mulsd  (%rsp),%xmm1
movsd  0xe2b(%rip),%xmm0
call   1030 <pow@plt>
movapd 0xe06(%rip),%xmm1
xorpd  %xmm0,%xmm1
xorpd  %xmm2,%xmm2
movapd 0x10(%rsp),%xmm3
cmpltsd %xmm2,%xmm3
movapd %xmm3,%xmm2
andnpd %xmm0,%xmm2
andpd  %xmm1,%xmm3
orpd   %xmm2,%xmm3
movlpd %xmm3,(%r15)
test   %r14,%r14
je     122a <mant_exp+0xfa>
mov    %ebx,(%r14)
add    $0x28,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @mant_exp(double %0, double %1, i32 %2, double* %3, i32* %4) {
%6 = alloca double
%7 = alloca double
%8 = alloca i32
%9 = alloca double*
%10 = alloca i32*
%11 = alloca i32
%12 = alloca double
%13 = alloca i32
store double %0, double* %6
store double %1, double* %7
store i32 %2, i32* %8
store double* %3, double** %9
store i32* %4, i32** %10
store i32 1, i32* %11
%14 = load double, double* %7
%15 = fcmp oeq double %14, 0.000000e+00
br i1 %15, label %16, label %27
16: ; preds = %5
%17 = load double*, double** %9
%18 = icmp ne double* %17, null
br i1 %18, label %19, label %21
19: ; preds = %16
%20 = load double*, double** %9
store double 0.000000e+00, double* %20
br label %21
21: ; preds = %19, %16
%22 = load i32*, i32** %10
%23 = icmp ne i32* %22, null
br i1 %23, label %24, label %26
24: ; preds = %21
%25 = load i32*, i32** %10
store i32 0, i32* %25
br label %26
26: ; preds = %24, %21
br label %71
27: ; preds = %5
%28 = load double, double* %7
%29 = fcmp olt double %28, 0.000000e+00
br i1 %29, label %30, label %33
30: ; preds = %27
store i32 -1, i32* %11
%31 = load double, double* %7
%32 = fneg double %31
store double %32, double* %7
br label %33
33: ; preds = %30, %27
%34 = load double, double* %7
%35 = call double @log10(double %34)
%36 = load double, double* %6
%37 = fdiv double %35, %36
store double %37, double* %12
%38 = load double, double* %12
%39 = call double @llvm.floor.f64(double %38)
%40 = fptosi double %39 to i32
store i32 %40, i32* %13
%41 = load i32, i32* %8
%42 = icmp ne i32 %41, 0
br i1 %42, label %43, label %50
43: ; preds = %33
%44 = load i32, i32* %13
%45 = sitofp i32 %44 to double
%46 = fdiv double %45, 3.000000e+00
%47 = call double @llvm.floor.f64(double %46)
%48 = fmul double 3.000000e+00, %47
%49 = fptosi double %48 to i32
store i32 %49, i32* %13
br label %50
50: ; preds = %43, %33
%51 = load double*, double** %9
%52 = icmp ne double* %51, null
br i1 %52, label %53, label %65
53: ; preds = %50
%54 = load i32, i32* %11
%55 = sitofp i32 %54 to double
%56 = load double, double* %12
%57 = load i32, i32* %13
%58 = sitofp i32 %57 to double
%59 = fsub double %56, %58
%60 = load double, double* %6
%61 = fmul double %59, %60
%62 = call double @pow(double 1.000000e+01, double %61)
%63 = fmul double %55, %62
%64 = load double*, double** %9
store double %63, double* %64
br label %65
65: ; preds = %53, %50
%66 = load i32*, i32** %10
%67 = icmp ne i32* %66, null
br i1 %67, label %68, label %71
68: ; preds = %65
%69 = load i32, i32* %13
%70 = load i32*, i32** %10
store i32 %69, i32* %70
br label %71
71: ; preds = %26, %68, %65
ret void
}

LLVM IR (-O2):
define void @mant_exp(double %0, double %1, i32 %2, double* writeonly %3, i32* writeonly %4) {
%6 = fcmp oeq double %1, 0.000000e+00
br i1 %6, label %7, label %12
7: ; preds = %5
%8 = icmp eq double* %3, null
br i1 %8, label %10, label %9
9: ; preds = %7
store double 0.000000e+00, double* %3
br label %10
10: ; preds = %9, %7
%11 = icmp eq i32* %4, null
br i1 %11, label %41, label %39
12: ; preds = %5
%13 = fcmp olt double %1, 0.000000e+00
%14 = fneg double %1
%15 = select i1 %13, double %14, double %1
%16 = tail call double @log10(double %15)
%17 = fdiv double %16, %0
%18 = tail call double @llvm.floor.f64(double %17)
%19 = fptosi double %18 to i32
%20 = icmp eq i32 %2, 0
br i1 %20, label %27, label %21
21: ; preds = %12
%22 = sitofp i32 %19 to double
%23 = fdiv double %22, 3.000000e+00
%24 = tail call double @llvm.floor.f64(double %23)
%25 = fmul double %24, 3.000000e+00
%26 = fptosi double %25 to i32
br label %27
27: ; preds = %21, %12
%28 = phi i32 [ %26, %21 ], [ %19, %12 ]
%29 = icmp eq double* %3, null
br i1 %29, label %37, label %30
30: ; preds = %27
%31 = sitofp i32 %28 to double
%32 = fsub double %17, %31
%33 = fmul double %32, %0
%34 = tail call double @pow(double 1.000000e+01, double %33)
%35 = fneg double %34
%36 = select i1 %13, double %35, double %34
store double %36, double* %3
br label %37
37: ; preds = %30, %27
%38 = icmp eq i32* %4, null
br i1 %38, label %41, label %39
39: ; preds = %37, %10
%40 = phi i32 [ 0, %10 ], [ %28, %37 ]
store i32 %40, i32* %4
br label %41
41: ; preds = %39, %37, %10
ret void
}