#if !DECODER_ONLY
/* Lzip - Data compressor based on the LZMA algorithm
Copyright (C) 2008, 2009, 2010, 2011 Antonio Diaz Diaz.
This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with this program. If not, see .
*/
#define _FILE_OFFSET_BITS 64
#include
#include
#include
#include
#include "lzip.h"
#include "encoder.h"
#include "fast_encoder.h"
bool Fmatchfinder::read_block()
{
if( !at_stream_end && stream_pos < buffer_size )
{
const int size = buffer_size - stream_pos;
const int rd = readblock( infd, buffer + stream_pos, size );
stream_pos += rd;
if( rd != size && errno ) throw Error( "Read error" );
at_stream_end = ( rd < size );
}
return pos < stream_pos;
}
Fmatchfinder::Fmatchfinder( const int ifd )
:
partial_data_pos( 0 ),
prev_positions( new int32_t[num_prev_positions] ),
pos( 0 ),
cyclic_pos( 0 ),
key4( 0 ),
stream_pos( 0 ),
match_len_limit_( 16 ),
infd( ifd ),
at_stream_end( false )
{
const int dict_size = 65536;
const int buffer_size_limit = ( 16 * dict_size ) + before_size + after_size;
buffer_size = dict_size;
buffer = (uint8_t *)malloc( buffer_size );
if( !buffer ) exit(-1);
if( read_block() && !at_stream_end && buffer_size < buffer_size_limit )
{
buffer_size = buffer_size_limit;
buffer = (uint8_t *)realloc( buffer, buffer_size );
if( !buffer ) exit(-1);
read_block();
}
if( at_stream_end && stream_pos < dict_size )
dictionary_size_ = max( (int)min_dictionary_size, stream_pos );
else dictionary_size_ = dict_size;
pos_limit = buffer_size;
if( !at_stream_end ) pos_limit -= after_size;
prev_pos_chain = new int32_t[dictionary_size_];
for( int i = 0; i < num_prev_positions; ++i ) prev_positions[i] = -1;
}
void Fmatchfinder::reset()
{
const int size = stream_pos - pos;
if( size > 0 ) memmove( buffer, buffer + pos, size );
partial_data_pos = 0;
stream_pos -= pos;
pos = 0;
cyclic_pos = 0;
key4 = 0;
for( int i = 0; i < num_prev_positions; ++i ) prev_positions[i] = -1;
read_block();
}
void Fmatchfinder::move_pos()
{
if( ++cyclic_pos >= dictionary_size_ ) cyclic_pos = 0;
if( ++pos >= pos_limit )
{
if( pos > stream_pos )
internal_error( "pos > stream_pos in Fmatchfinder::move_pos" );
if( !at_stream_end )
{
const int offset = pos - dictionary_size_ - before_size;
const int size = stream_pos - offset;
memmove( buffer, buffer + offset, size );
partial_data_pos += offset;
pos -= offset;
stream_pos -= offset;
for( int i = 0; i < num_prev_positions; ++i )
if( prev_positions[i] >= 0 ) prev_positions[i] -= offset;
for( int i = 0; i < dictionary_size_; ++i )
if( prev_pos_chain[i] >= 0 ) prev_pos_chain[i] -= offset;
read_block();
}
}
}
int Fmatchfinder::longest_match_len( int * const distance )
{
int len_limit = match_len_limit_;
if( len_limit > available_bytes() )
{
len_limit = available_bytes();
if( len_limit < 4 ) return 0;
}
const uint8_t * const data = buffer + pos;
key4 = ( ( key4 << 4 ) ^ data[3] ) & ( num_prev_positions - 1 );
int newpos = prev_positions[key4];
prev_positions[key4] = pos;
int32_t * ptr0 = prev_pos_chain + cyclic_pos;
int maxlen = 0;
for( int count = 4; ; )
{
if( newpos < (pos - dictionary_size_ + 1) || newpos < 0 || --count < 0 )
{ *ptr0 = -1; break; }
const uint8_t * const newdata = buffer + newpos;
int len = 0;
if( newdata[maxlen] == data[maxlen] )
while( len < len_limit && newdata[len] == data[len] ) ++len;
const int delta = pos - newpos;
if( maxlen < len ) { maxlen = len; *distance = delta - 1; }
int32_t * const newptr = prev_pos_chain +
( cyclic_pos - delta +
( ( cyclic_pos >= delta ) ? 0 : dictionary_size_ ) );
if( len < len_limit )
{
*ptr0 = newpos;
ptr0 = newptr;
newpos = *ptr0;
}
else
{
*ptr0 = *newptr;
break;
}
}
return maxlen;
}
void Fmatchfinder::longest_match_len()
{
int len_limit = match_len_limit_;
if( len_limit > available_bytes() )
{
len_limit = available_bytes();
if( len_limit < 4 ) return;
}
const uint8_t * const data = buffer + pos;
key4 = ( ( key4 << 4 ) ^ data[3] ) & ( num_prev_positions - 1 );
const int newpos = prev_positions[key4];
prev_positions[key4] = pos;
int32_t * const ptr0 = prev_pos_chain + cyclic_pos;
if( newpos < (pos - dictionary_size_ + 1) || newpos < 0 ) *ptr0 = -1;
else
{
const uint8_t * const newdata = buffer + newpos;
if( newdata[len_limit-1] != data[len_limit-1] ||
memcmp( newdata, data, len_limit - 1 ) ) *ptr0 = newpos;
else
{
int idx = cyclic_pos - pos + newpos;
if( idx < 0 ) idx += dictionary_size_;
*ptr0 = prev_pos_chain[idx];
}
}
}
// Return value == number of bytes advanced (len).
// *disp returns the distance to encode.
// ( *disp == -1 && len == 1 ) means literal.
int FLZ_encoder::sequence_optimizer( const int reps[num_rep_distances],
int * const disp, const State & state )
{
const int main_len = read_match_distances();
int replen = 0;
int rep_index = 0;
for( int i = 0; i < num_rep_distances; ++i )
{
const int len = fmatchfinder.true_match_len( 0, reps[i] + 1, max_match_len );
if( len > replen ) { replen = len; rep_index = i; }
}
if( replen > min_match_len && replen + 4 > main_len )
{
*disp = rep_index;
move_pos( replen, true );
return replen;
}
if( main_len > min_match_len ||
( main_len == min_match_len && match_distance < modeled_distances ) )
{
*disp = num_rep_distances + match_distance;
move_pos( main_len, true );
return main_len;
}
const uint8_t cur_byte = fmatchfinder[0];
const uint8_t match_byte = fmatchfinder[-reps[0]-1];
*disp = -1;
if( match_byte == cur_byte )
{
const uint8_t prev_byte = fmatchfinder[-1];
const int pos_state = fmatchfinder.data_position() & pos_state_mask;
int price = price0( bm_match[state()][pos_state] );
if( state.is_char() )
price += literal_encoder.price_symbol( prev_byte, cur_byte );
else
price += literal_encoder.price_matched( prev_byte, cur_byte, match_byte );
const int short_rep_price = price1( bm_match[state()][pos_state] ) +
price1( bm_rep[state()] ) +
price0( bm_rep0[state()] ) +
price0( bm_len[state()][pos_state] );
if( short_rep_price < price ) *disp = 0;
}
fmatchfinder.move_pos();
return 1;
}
// End Of Stream mark => (dis == 0xFFFFFFFFU, len == min_match_len)
void FLZ_encoder::full_flush( const State & state )
{
const int pos_state = fmatchfinder.data_position() & pos_state_mask;
range_encoder.encode_bit( bm_match[state()][pos_state], 1 );
range_encoder.encode_bit( bm_rep[state()], 0 );
encode_pair( 0xFFFFFFFFU, min_match_len, pos_state );
range_encoder.flush();
File_trailer trailer;
trailer.data_crc( crc() );
trailer.data_size( fmatchfinder.data_position() );
trailer.member_size( range_encoder.member_position() + File_trailer::size() );
for( int i = 0; i < File_trailer::size(); ++i )
range_encoder.put_byte( trailer.data[i] );
range_encoder.flush_data();
}
FLZ_encoder::FLZ_encoder( Fmatchfinder & mf, const File_header & header,
const int outfd )
:
crc_( 0xFFFFFFFFU ),
fmatchfinder( mf ),
range_encoder( outfd ),
len_encoder( fmatchfinder.match_len_limit() ),
rep_match_len_encoder( fmatchfinder.match_len_limit() ),
num_dis_slots( 2 * real_bits( fmatchfinder.dictionary_size() - 1 ) )
{
for( int i = 0; i < File_header::size; ++i )
range_encoder.put_byte( header.data[i] );
}
bool FLZ_encoder::encode_member( const long long member_size )
{
const long long member_size_limit =
member_size - File_trailer::size() - max_marker_size;
int rep_distances[num_rep_distances];
State state;
for( int i = 0; i < num_rep_distances; ++i ) rep_distances[i] = 0;
if( fmatchfinder.data_position() != 0 ||
range_encoder.member_position() != File_header::size )
return false; // can be called only once
if( !fmatchfinder.finished() ) // encode first byte
{
const uint8_t prev_byte = 0;
const uint8_t cur_byte = fmatchfinder[0];
range_encoder.encode_bit( bm_match[state()][0], 0 );
literal_encoder.encode( range_encoder, prev_byte, cur_byte );
crc32.update( crc_, cur_byte );
move_pos( 1 );
}
while( true )
{
if( fmatchfinder.finished() ) { full_flush( state ); return true; }
const int pos_state = fmatchfinder.data_position() & pos_state_mask;
int dis;
const int len = sequence_optimizer( rep_distances, &dis, state );
if( len <= 0 ) return false;
bool bit = ( dis < 0 && len == 1 );
range_encoder.encode_bit( bm_match[state()][pos_state], !bit );
if( bit ) // literal byte
{
const uint8_t prev_byte = fmatchfinder[-len-1];
const uint8_t cur_byte = fmatchfinder[-len];
crc32.update( crc_, cur_byte );
if( state.is_char() )
literal_encoder.encode( range_encoder, prev_byte, cur_byte );
else
{
const uint8_t match_byte = fmatchfinder[-len-rep_distances[0]-1];
literal_encoder.encode_matched( range_encoder,
prev_byte, cur_byte, match_byte );
}
state.set_char();
}
else // match or repeated match
{
crc32.update( crc_, fmatchfinder.ptr_to_current_pos() - len, len );
mtf_reps( dis, rep_distances );
bit = ( dis < num_rep_distances );
range_encoder.encode_bit( bm_rep[state()], bit );
if( bit )
{
bit = ( dis == 0 );
range_encoder.encode_bit( bm_rep0[state()], !bit );
if( bit )
range_encoder.encode_bit( bm_len[state()][pos_state], len > 1 );
else
{
range_encoder.encode_bit( bm_rep1[state()], dis > 1 );
if( dis > 1 )
range_encoder.encode_bit( bm_rep2[state()], dis > 2 );
}
if( len == 1 ) state.set_short_rep();
else
{
rep_match_len_encoder.encode( range_encoder, len, pos_state );
state.set_rep();
}
}
else
{
encode_pair( dis - num_rep_distances, len, pos_state );
state.set_match();
}
}
if( range_encoder.member_position() >= member_size_limit )
{
full_flush( state );
return true;
}
}
}
#endif