similar to: mdf no sound issue

Did a bit of tuning in SVN. Could you try that? What frame size, filter length, sampling rate are you using? All the process you describe is actually the tip of the iceberg I think. There is something making the adaptation diverge, which causes Eh and Yh to be large and so on. The part I changed is probably the one causing divergence. If it doesn't work, reducing the coefficients even more

Hi I have a question about the AEC and preprocessor. I have seen that when near end is talking for a longer time (about 10 s) without being interrupted by the far end the residual echo (or rather leak_estimate) increases making the preprocessor suppress the near end talk when it shouldn't. Is there a way to make the leakage estimate only update when far-end is present or similar in order to

Hi Ted, Thanks a lot for this analysis. > In FLOAT_DIVU() it hangs at the following: > while (a.m >= b.m) > { > e++; > a.m >>= 1; > } > for the case where a and b are both zero (yes, division by zero). > This happens from mdf.c: True, that needs to be fixed even after I fix the rest. > leak_estimate =

Hi I just upgraded to http://svn.xiph.org/trunk/speex r10357 and got this build warning: alfredh@io:$ make -s mdf.o libspeex/mdf.c: In function 'speex_echo_cancel': libspeex/mdf.c:321: warning: statement with no effect libspeex/mdf.c:317: warning: `adapt_rate' might be used uninitialized in this function Is this intentional? In any case here is a simple fix: Index:

> I am writing to gain a better understanding of the limitations of speex echo > cancellation, esp. with respect to the fixed point implementation. > If these limitations have been documented elsewhere already, please let me > know! Nothing officially documented, sorry. > I observe experimentally that when one or both of the echo or ref data for > speex_echo_cancel() have

Jean Marc, I ran across the silence issue that results from the echo canceller producing "bad" floating point data. Unfortunately this hasn't been done in a way i can reproduce, but I found that calling reset didn't fix the silence, unless I also reset: st->PHI[i] st->Eh[i] st->Yh[i] So it may be good to add these to the reset function. I was using these in

Jean Marc, Ok so I tested with the new code, same result- frame size= 160, filter length=160 ms/1280 samples, 8000 Hz diverged after 9564 calls/packets or 191 seconds. My system stays relatively in synch but is not perfect- difficult to measure if there is any clock drift, but it probably is. In order to make it break faster I use an open air usb microphone that is part of a logitech notebook

Thanks for your prompt and helpful reply. In FLOAT_DIVU() it hangs at the following: while (a.m >= b.m) { e++; a.m >>= 1; } for the case where a and b are both zero (yes, division by zero). This happens from mdf.c: leak_estimate = FLOAT_EXTRACT16(FLOAT_SHL(FLOAT_DIVU(st->Pey, st->Pyy),14)); where st->Pey and st->Pyy are both zero, which happens for the

Jean-Marc, I have traced the second infinite loop further. When st->adapted becomes true (mdf.c line 623), the first Yf[i] value is 4, the leak_estimate is 0xd4e, the resulting r is 3. The first value in st->Rf is 0, so e is 1, and r is set to e>>1, or 0. A little later there is a divide by r, and there is the hang. It seems that the 0 in Rf[0] is the problem, but I am not

Dear All, I would like to plot text with a box around it. I used strwidth and strheight to compute the size of the box which is plotted with rect: z <- rnorm(10) # horizontal text works plot(rnorm(10)) x1 <- 5 y1 <- 0 label <- "Label" cha <- paste(" ", label, " ", sep = "") xh <- strwidth(cha, cex = par("cex")) yh <-

Current SVN compilation fails on Win32 (Intel C++ compiler 9.0): ..\..\libspeex\mdf.c(317): error: declaration may not appear after executable statement in block float adapt_rate; Just moved it up to the variable declaration block to solve it locally, but as it's not valid C, I thought I'd mention it ;) PS: With the same compiler, AGC fails horribly (Zlast goes to several

Dear all, I noticed the following behaviour of plot.lm: > fm1 <- lm(time~dist, data=hills, weights=c(0,0,rep(1,33))) > par(mfrow=c(2,2)) > plot(fm1) Warning messages: 1: longer object length is not a multiple of shorter object length in: res/(sd * (1 - hat)) 2: longer object length is not a multiple of shorter object length in: (res/(sd * (1 - hat)))^2 * hat which seems to be

Skipped content of type multipart/alternative-------------- next part -------------- A non-text attachment was scrubbed... Name: preproc_patch_dth_10_26_05.patch Type: application/octet-stream Size: 8774 bytes Desc: not available Url : http://lists.xiph.org/pipermail/speex-dev/attachments/20051026/c27a3ed6/preproc_patch_dth_10_26_05.obj

I append below a suggested update for p.adjust(). 1. A new method "yh" for control of FDR is included which is valid for any dependency structure. Reference is Benjamini, Y., and Yekutieli, D. (2001). The control of the false discovery rate in multiple testing under dependency. Annals of Statistics 29, 1165-1188. 2. I've re-named the "fdr" method to "bh" but

Hello Jean-Marc, When you count the adaptation step size in case where the filter has already had minimal adaptation, you use the following expression: r = (0.7*r + 0.3*15*RER*e)/e*(power[i] + 10) = (0.7*leak_estimate*Yf[i] + 0.3*15*RER*(Rf[i] + 1))/((Rf[i] + 1)*(power[i] + 10)). Why do we need this weighted sum and the component 0.3*15*RER*(Rf[i] + 1)? Why use the correlation-based RER if

Hi, Is it possible to program ssh tunneling in C or C++ rather than to use shell command? Are there any APIs for ssh tunneling available? I was told by friends that the openssh is working on it, please point me which version of openssh support API ssh tunneling. Thank you. Kind regards, J. yh

Dear R-users, One of the main reasons I moved from GAUSS to R (as an econometrician) was because of the existence of the library LOCFIT for local polynomial regression. While doing some checking between my former `GAUSS code' and my new `R code', I came to realize LOCFIT is not quite doing what I want. I wrote the following example script:

> On Tue, Jul 15, 2008 at 10:33 AM, Suresh Siddha > <suresh.b.siddha at intel.com> wrote: > > On Sun, Jul 13, 2008 at 10:19:35PM -0700, Yinghai Lu wrote: > >> > >> fix for pv. > >> > >> Signed-off-by: Yinghai Lu <yhlu.kernel at gmail.com> > >> > >> --- > >> arch/x86/kernel/paravirt.c | 5 ---- > >>

> On Tue, Jul 15, 2008 at 10:33 AM, Suresh Siddha > <suresh.b.siddha at intel.com> wrote: > > On Sun, Jul 13, 2008 at 10:19:35PM -0700, Yinghai Lu wrote: > >> > >> fix for pv. > >> > >> Signed-off-by: Yinghai Lu <yhlu.kernel at gmail.com> > >> > >> --- > >> arch/x86/kernel/paravirt.c | 5 ---- > >>

thanks to all for the 2d scatter plot. i have one more. how do i plot 'hist(y1, col="red") and hist(y2,col="blue") in the same window? thanks again.