#! perl -slw
use Config; package F128;
use Inline C => Config => BUILD_NOISY => 1;
use Inline C => <<END_C,  NAME => 'F128', CLEAN_AFTER_BUILD =>0, TYPEMAPS => './F128.typemap';
#include "../C/mytypes.h"
#define CLASS "F128"

typedef struct { U64 scale; U64 precision; } F128;

F128 *new( I64 scale, U64 precision ) {
    F128 *n = malloc( sizeof( F128 ) );
    n->scale = scale;
    n->precision = precision;
    return n;
}

F128 *add( F128 *a, F128 *b ) {
    if( a->scale != b->scale ) {
        while( a->scale > b->scale ) {
            a->scale >>= 1;
            a->precision <<= 1;
        }
        while( a->scale < b->scale ) {
            a->scale <<= 1;
            a->precision >>= 1;
        }
    }
    a->precision += b->precision;
    return a;
}

F128 *multiply( F128 *a, F128 *b ) {
    a->precision *= b->precision;
    a->scale += b->scale;
    return a;
}

SV *toString( F128 *a ) {
    char *out = malloc( a->scale );
    U32 c = sprintf( out, "%I64ue%I64i", a->precision, a->scale );
    return newSVpv( out, c );
}

void DESTROY( F128 *a ) {
    free( a );
    return;
}

END_C
package main;
my $Fa = F128::new( -10, 45 );
my $Fb = F128::new( -10, 45 );
my $Fc = F128::new( -10, 45 );

print $Fa->toString;

$Fa->add( $Fb );
print $Fa->toString;

$Fc->multiply( $Fb );
print $Fc->toString;