datatypes/fx/scfx_rep.cpp

100 : m_mant( min_mant ), m_wp(), m_sign(), m_state(), m_msw(), m_lsw(),
107 : m_mant( min_mant ), m_wp(), m_sign(), m_state(), m_msw(), m_lsw(),
113 	m_wp = m_msw = m_lsw = 2;
131 : m_mant( min_mant ), m_wp(), m_sign(), m_state(), m_msw(), m_lsw(),
137 	m_wp = m_msw = m_lsw = 2;
147 : m_mant( min_mant ), m_wp(), m_sign(), m_state(), m_msw(), m_lsw(),
170             m_msw = 2;
181 : m_mant( min_mant ), m_wp(), m_sign(), m_state(), m_msw(), m_lsw(),
187 	m_wp = m_msw = m_lsw = 2;
196 	    m_msw = 2;
207 : m_mant( min_mant ), m_wp( 0 ), m_sign(), m_state( normal ), m_msw( 0 ),
235 : m_mant( min_mant ), m_wp(), m_sign(), m_state(), m_msw(), m_lsw(),
262 : m_mant( min_mant ), m_wp(), m_sign(), m_state(), m_msw(), m_lsw(),
280 : m_mant( min_mant ), m_wp(), m_sign(), m_state(), m_msw(), m_lsw(),
323 : m_mant( min_mant ), m_wp(), m_sign(), m_state(), m_msw(), m_lsw(),
354   m_msw( a.m_msw ), m_lsw( a.m_lsw ), m_r_flag( false )
664 	    m_msw = m_lsw = 0;
762 	m_mant[m_msw] |=  -1 << scfx_find_msb( m_mant[m_msw] );
763 	for( int i = m_msw + 1; i < m_mant.size(); ++ i )
809     int msb = scfx_find_msb( m_mant[m_msw] );
811     int exp = (m_msw - m_wp) * bits_in_word + msb;
835         m0 = m_mant[m_msw] & ~( 1 << mantissa0_size );
836 	if( m_msw > m_lsw )
838 	    m1 = m_mant[m_msw - 1];
839 	    if( m_msw - 1 > m_lsw )
840 	        guard = m_mant[m_msw - 2] >> ( bits_in_word - 1 );
845 	m0 = ( m_mant[m_msw] >> -shift ) & ~( 1 << mantissa0_size );
846 	m1 = m_mant[m_msw] << ( bits_in_word + shift );
847 	if( m_msw > m_lsw )
849 	    m1 |= m_mant[m_msw - 1] >> -shift;
850 	    guard = ( m_mant[m_msw - 1] >> ( -shift - 1 ) ) & 1;
855 	m0 = ( m_mant[m_msw] << shift ) & ~( 1 << mantissa0_size );
856 	if( m_msw > m_lsw )
858 	    m0 |= m_mant[m_msw - 1] >> ( bits_in_word - shift );
859 	    m1 = m_mant[m_msw - 1] << shift;
860 	    if( m_msw - 1 > m_lsw )
862 	        m1 |= m_mant[m_msw - 2] >> ( bits_in_word - shift );
863 		guard = ( m_mant[m_msw - 2] >> (bits_in_word - shift - 1) )
935     for( i = int_part.m_lsw; i <= int_part.m_msw && i < int_part.m_wp; i ++ )
941     for( i = frac_part.m_msw;
946     if( frac_part.m_msw == frac_part.size() - 1 )
956 	double int_wl = ( int_part.m_msw - int_part.m_wp ) * bits_in_word
957 	              + scfx_find_msb( int_part.m_mant[int_part.m_msw] ) + 1;
1001 	double frac_wl = ( frac_part.m_wp - frac_part.m_msw ) * bits_in_word
1002 	               - scfx_find_msb( frac_part.m_mant[frac_part.m_msw] )
1009 	if( frac_part.m_msw == frac_part.size() - 1 )
1023 	    int n = frac_part.m_mant[frac_part.m_msw + 1];
1036 	    frac_part.m_mant[frac_part.m_msw + 1] = 0;
1123 	    msb = ( b.m_msw - b.m_wp ) * bits_in_word
1124 		+ scfx_find_msb( b.m_mant[ b.m_msw ] ) + 1;
1572     int len_lhs = lhs.m_msw - lhs.m_lsw + 1;
1573     int len_rhs = rhs.m_msw - rhs.m_lsw + 1;
1661     int msb_lhs = scfx_find_msb( lhs.m_mant[lhs.m_msw] )
1662 	        + ( lhs.m_msw - lhs.m_wp ) * bits_in_word;
1663     int msb_rhs = scfx_find_msb( rhs.m_mant[rhs.m_msw] )
1664 	        + ( rhs.m_msw - rhs.m_wp ) * bits_in_word;
1681     if( remainder.m_msw == remainder.size() - 1 )
1685     int msw_diff = rhs.m_msw - remainder.m_msw;
1696 	    sub_with_index( remainder.m_mant, remainder.m_msw, remainder.m_lsw,
1697 			    rhs.m_mant, rhs.m_msw, rhs.m_lsw );
1746 	if( m_msw == size() - 1 &&
1747 	    scfx_find_msb( m_mant[m_msw] ) >= bits_in_word - shift_bits )
1802     word a_word = a.m_mant[a.m_msw];
1803     word b_word = b.m_mant[b.m_msw];
1816     int a_msw = a.m_msw - a.m_wp;
1817     int b_msw = b.m_msw - b.m_wp;
1827     int a_i = a.m_msw;
1828     int b_i = b.m_msw;
2267 	int upper_bound_lhs = lhs.m_msw - lhs.m_wp;
2269 	int upper_bound_rhs = rhs.m_msw - rhs.m_wp;
2309     // special case: rhs.m_mant[rhs.m_msw + 1] == 1
2310     if( rhs.m_msw < rhs.size() - 1 && rhs.m_mant[rhs.m_msw + 1 ] != 0 )
2315     int lhs_size = lhs.m_msw - lhs.m_lsw + 1;
2316     int rhs_size = rhs.m_msw - rhs.m_lsw + 1;
2320     int lhs_index = lhs.m_msw;
2321     int rhs_index = rhs.m_msw;
2367     half_word* hw = (half_word*) &m_mant[m_msw];
2369     half_word* hw = ( (half_word*) &m_mant[m_msw] ) + 1;
2378     for( int i = 0, end = ( m_msw - m_wp + 1 ) * 2; i < end; i ++ )
2380     for( int i = 0, end = -( m_msw - m_wp + 1 ) * 2; i > end; i -- )
2474     for( int j = m_lsw; j <= m_msw; j ++ )
2814     os << "msw   = " << m_msw << ::std::endl;
2844     int msb = ( m_msw - m_wp ) * bits_in_word
2845             + scfx_find_msb( m_mant[ m_msw ] ) + 1;
2887     wl_effective = ( m_msw - m_lsw + 1 ) * bits_in_word;
2893     int msb = scfx_find_msb( m_mant[m_msw] );
2896     wl_effective = ( m_msw * bits_in_word + msb ) -
2903     int wi = m_msw - ( wl - 1 ) / bits_in_word;