claunia/flac
1
0
Fork 0
mirror of https://github.com/claunia/flac.git synced 2025-12-16 18:54:26 +00:00
Code Issues Packages Projects Releases Wiki Activity

add more debug info, clip bps computation to 0.0

Browse Source
This commit is contained in:
Josh Coalson
2001-02-08 00:26:45 +00:00
parent c4b7c84676
commit 9f77a19607
1 changed files with 29 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

41
src/libFLAC/lpc.c
View File

@@ -149,19 +149,19 @@ void FLAC__lpc_compute_residual_from_qlp_coefficients(const int32 data[], unsign
sum += qlp_coeff[j] * (*(--history)); sum += qlp_coeff[j] * (*(--history));
#ifdef FLAC_OVERFLOW_DETECT #ifdef FLAC_OVERFLOW_DETECT
sumo += (int64)qlp_coeff[j] * (int64)(*history); sumo += (int64)qlp_coeff[j] * (int64)(*history);
if(sumo > 2147483647ll || sumo < -2147483648ll) if(sumo > 2147483647ll || sumo < -2147483648ll) {
fprintf(stderr,"FLAC__lpc_compute_residual_from_qlp_coefficients: OVERFLOW, sumo=%lld\n",sumo); fprintf(stderr,"FLAC__lpc_compute_residual_from_qlp_coefficients: OVERFLOW, i=%u, j=%u, c=%d, d=%d, sumo=%lld\n",i,j,qlp_coeff[j],*history,sumo);
}
#endif #endif
} }
*(residual++) = *(data++) - (sum >> lp_quantization); *(residual++) = *(data++) - (sum >> lp_quantization);
} }
/* Here's a slightly slower but clearer version: /* Here's a slower but clearer version:
for(i = 0; i < data_len; i++) { for(i = 0; i < data_len; i++) {
sum = 0; sum = 0;
history = &(data[i]);
for(j = 0; j < order; j++) for(j = 0; j < order; j++)
sum += qlp_coeff[j] * (*(--history)); sum += qlp_coeff[j] * data[i-j-1];
residual[i] = data[i] - (sum >> lp_quantization); residual[i] = data[i] - (sum >> lp_quantization);
} }
*/ */
@@ -173,7 +173,8 @@ void FLAC__lpc_restore_signal(const int32 residual[], unsigned data_len, const i
int64 sumo; int64 sumo;
#endif #endif
unsigned i, j; unsigned i, j;
int32 sum, *history; int32 sum;
const int32 *history;
#ifdef FLAC_OVERFLOW_DETECT_VERBOSE #ifdef FLAC_OVERFLOW_DETECT_VERBOSE
fprintf(stderr,"FLAC__lpc_restore_signal: data_len=%d, order=%u, lpq=%d",data_len,order,lp_quantization); fprintf(stderr,"FLAC__lpc_restore_signal: data_len=%d, order=%u, lpq=%d",data_len,order,lp_quantization);
@@ -188,17 +189,27 @@ void FLAC__lpc_restore_signal(const int32 residual[], unsigned data_len, const i
sumo = 0; sumo = 0;
#endif #endif
sum = 0; sum = 0;
history = data+i; history = data;
for(j = 0; j < order; j++) { for(j = 0; j < order; j++) {
sum += qlp_coeff[j] * (*(--history)); sum += qlp_coeff[j] * (*(--history));
#ifdef FLAC_OVERFLOW_DETECT #ifdef FLAC_OVERFLOW_DETECT
sumo += (int64)qlp_coeff[j] * (int64)(*history); sumo += (int64)qlp_coeff[j] * (int64)(*history);
if(sumo > 2147483647ll || sumo < -2147483648ll) if(sumo > 2147483647ll || sumo < -2147483648ll) {
fprintf(stderr,"FLAC__lpc_restore_signal: OVERFLOW, sumo=%lld\n",sumo); fprintf(stderr,"FLAC__lpc_restore_signal: OVERFLOW, i=%u, j=%u, c=%d, d=%d, sumo=%lld\n",i,j,qlp_coeff[j],*history,sumo);
}
#endif #endif
} }
*(data++) = *(residual++) + (sum >> lp_quantization);
}
/* Here's a slower but clearer version:
for(i = 0; i < data_len; i++) {
sum = 0;
for(j = 0; j < order; j++)
sum += qlp_coeff[j] * data[i-j-1];
data[i] = residual[i] + (sum >> lp_quantization); data[i] = residual[i] + (sum >> lp_quantization);
} }
*/
} }
real FLAC__lpc_compute_expected_bits_per_residual_sample(real lpc_error, unsigned total_samples) real FLAC__lpc_compute_expected_bits_per_residual_sample(real lpc_error, unsigned total_samples)
@@ -210,10 +221,16 @@ real FLAC__lpc_compute_expected_bits_per_residual_sample(real lpc_error, unsigne
escale = 0.5 * M_LN2 * M_LN2 / (real)total_samples; escale = 0.5 * M_LN2 * M_LN2 / (real)total_samples;
if(lpc_error > 0.0) if(lpc_error > 0.0) {
return 0.5 * log(escale * lpc_error) / M_LN2; real bps = 0.5 * log(escale * lpc_error) / M_LN2;
else if(bps >= 0.0)
return bps;
else
return 0.0;
}
else {
return 0.0; return 0.0;
}
} }
unsigned FLAC__lpc_compute_best_order(const real lpc_error[], unsigned max_order, unsigned total_samples, unsigned bits_per_signal_sample) unsigned FLAC__lpc_compute_best_order(const real lpc_error[], unsigned max_order, unsigned total_samples, unsigned bits_per_signal_sample)