search for: cscalfn

Dear all, Somebody kindly pointed out a problem in my WaveThresh3 code. I can't figure out what is wrong. I have whittled down the code quite a bit into an example case which repeats the problem. There is one R function called "ScalingFunction" and one C function called "CScalFn.c". The idea is that the R function calls the C routine repeatedly to compute a wavelet scaling function. It *seems* that after some iterations of the C routine later on in the code some NAs mysteriously enter into some of the vectors (e.g. v) which then gets rejected by the next call to...