diff options
Diffstat (limited to '1.4.23-rc4/codecs/lpc10/pitsyn.c')
| -rw-r--r-- | 1.4.23-rc4/codecs/lpc10/pitsyn.c | 583 |
1 files changed, 0 insertions, 583 deletions
diff --git a/1.4.23-rc4/codecs/lpc10/pitsyn.c b/1.4.23-rc4/codecs/lpc10/pitsyn.c deleted file mode 100644 index ea8177e22..000000000 --- a/1.4.23-rc4/codecs/lpc10/pitsyn.c +++ /dev/null @@ -1,583 +0,0 @@ -/* - -$Log$ -Revision 1.16 2004/06/26 03:50:14 markster -Merge source cleanups (bug #1911) - -Revision 1.15 2003/11/23 22:14:32 markster -Various warning cleanups - -Revision 1.14 2003/02/12 13:59:15 matteo -mer feb 12 14:56:57 CET 2003 - -Revision 1.1.1.1 2003/02/12 13:59:15 matteo -mer feb 12 14:56:57 CET 2003 - -Revision 1.2 2000/01/05 08:20:39 markster -Some OSS fixes and a few lpc changes to make it actually work - - * Revision 1.2 1996/08/20 20:40:12 jaf - * Removed all static local variables that were SAVE'd in the Fortran - * code, and put them in struct lpc10_decoder_state that is passed as an - * argument. - * - * Removed init function, since all initialization is now done in - * init_lpc10_decoder_state(). - * - * Revision 1.1 1996/08/19 22:31:12 jaf - * Initial revision - * - -*/ - -/* -- translated by f2c (version 19951025). - You must link the resulting object file with the libraries: - -lf2c -lm (in that order) -*/ - -#include "f2c.h" - -#ifdef P_R_O_T_O_T_Y_P_E_S -extern int pitsyn_(integer *order, integer *voice, integer *pitch, real *rms, real *rc, integer *lframe, integer *ivuv, integer *ipiti, real *rmsi, real *rci, integer *nout, real *ratio, struct lpc10_decoder_state *st); -#endif - -/* ***************************************************************** */ - -/* PITSYN Version 53 */ - -/* $Log$ - * Revision 1.16 2004/06/26 03:50:14 markster - * Merge source cleanups (bug #1911) - * - * Revision 1.15 2003/11/23 22:14:32 markster - * Various warning cleanups - * - * Revision 1.14 2003/02/12 13:59:15 matteo - * mer feb 12 14:56:57 CET 2003 - * - * Revision 1.1.1.1 2003/02/12 13:59:15 matteo - * mer feb 12 14:56:57 CET 2003 - * - * Revision 1.2 2000/01/05 08:20:39 markster - * Some OSS fixes and a few lpc changes to make it actually work - * - * Revision 1.2 1996/08/20 20:40:12 jaf - * Removed all static local variables that were SAVE'd in the Fortran - * code, and put them in struct lpc10_decoder_state that is passed as an - * argument. - * - * Removed init function, since all initialization is now done in - * init_lpc10_decoder_state(). - * - * Revision 1.1 1996/08/19 22:31:12 jaf - * Initial revision - * */ -/* Revision 1.2 1996/03/25 18:49:07 jaf */ -/* Added commments about which indices of array arguments are read or */ -/* written. */ - -/* Rearranged local variable declarations to indicate which need to be */ -/* saved from one invocation to the next. Added entry INITPITSYN to */ -/* reinitialize local state variables, if desired. */ - -/* Added lots of comments about proving that the maximum number of pitch */ -/* periods (NOUT) that can be returned is 16. The call to STOP that */ -/* could happen if NOUT got too large was removed as a result. */ - -/* Also proved that the total number of samples returned from N calls, */ -/* each with identical values of LFRAME, will always be in the range */ -/* N*LFRAME-MAXPIT+1 to N*LFRAME. */ - -/* Revision 1.1 1996/02/07 14:48:18 jaf */ -/* Initial revision */ - - -/* ***************************************************************** */ - -/* Synthesize a single pitch epoch */ - -/* Input: */ -/* ORDER - Synthesis order (number of RC's) */ -/* VOICE - Half frame voicing decisions */ -/* Indices 1 through 2 read. */ -/* LFRAME - Length of speech buffer */ -/* Input/Output: */ -/* PITCH - Pitch */ -/* This value should be in the range MINPIT (20) to MAXPIT */ -/* (156), inclusive. */ -/* PITCH can be modified under some conditions. */ -/* RMS - Energy (can be modified) */ -/* RMS is changed to 1 if the value passed in is less than 1. */ -/* RC - Reflection coefficients */ -/* Indices 1 through ORDER can be temporarily overwritten with */ -/* RCO, and then replaced with original values, under some */ -/* conditions. */ -/* Output: */ -/* IVUV - Pitch epoch voicing decisions */ -/* Indices (I) of IVUV, IPITI, and RMSI are written, */ -/* and indices (J,I) of RCI are written, */ -/* where I ranges from 1 to NOUT, and J ranges from 1 to ORDER. */ -/* IPITI - Pitch epoch length */ -/* RMSI - Pitch epoch energy */ -/* RCI - Pitch epoch RC's */ -/* NOUT - Number of pitch periods in this frame */ -/* This is at least 0, at least 1 if MAXPIT .LT. LFRAME (this */ -/* is currently true on every call), and can never be more than */ -/* (LFRAME+MAXPIT-1)/PITCH, which is currently 16 with */ -/* LFRAME=180, MAXPIT=156, and PITCH .GE. 20, as SYNTHS */ -/* guarantees when it calls this subroutine. */ -/* RATIO - Previous to present energy ratio */ -/* Always assigned a value. */ - -/* Subroutine */ int pitsyn_(integer *order, integer *voice, - integer *pitch, real *rms, real *rc, integer *lframe, integer *ivuv, - integer *ipiti, real *rmsi, real *rci, integer *nout, real *ratio, - struct lpc10_decoder_state *st) -{ - /* Initialized data */ - - real *rmso; - logical *first; - - /* System generated locals */ - integer rci_dim1 = 0, rci_offset, i__1, i__2; - real r__1; - - /* Builtin functions */ - double log(doublereal), exp(doublereal); - - /* Local variables */ - real alrn, alro, yarc[10], prop; - integer i__, j, vflag, jused, lsamp; - integer *jsamp; - real slope; - integer *ipito; - real uvpit; - integer ip, nl, ivoice; - integer *ivoico; - integer istart; - real *rco; - real xxy; - -/* Arguments */ -/* $Log$ - * Revision 1.16 2004/06/26 03:50:14 markster - * Merge source cleanups (bug #1911) - * - * Revision 1.15 2003/11/23 22:14:32 markster - * Various warning cleanups - * - * Revision 1.14 2003/02/12 13:59:15 matteo - * mer feb 12 14:56:57 CET 2003 - * - * Revision 1.1.1.1 2003/02/12 13:59:15 matteo - * mer feb 12 14:56:57 CET 2003 - * - * Revision 1.2 2000/01/05 08:20:39 markster - * Some OSS fixes and a few lpc changes to make it actually work - * - * Revision 1.2 1996/08/20 20:40:12 jaf - * Removed all static local variables that were SAVE'd in the Fortran - * code, and put them in struct lpc10_decoder_state that is passed as an - * argument. - * - * Removed init function, since all initialization is now done in - * init_lpc10_decoder_state(). - * - * Revision 1.1 1996/08/19 22:31:12 jaf - * Initial revision - * */ -/* Revision 1.3 1996/03/29 22:03:47 jaf */ -/* Removed definitions for any constants that were no longer used. */ - -/* Revision 1.2 1996/03/26 19:34:33 jaf */ -/* Added comments indicating which constants are not needed in an */ -/* application that uses the LPC-10 coder. */ - -/* Revision 1.1 1996/02/07 14:43:51 jaf */ -/* Initial revision */ - -/* LPC Configuration parameters: */ -/* Frame size, Prediction order, Pitch period */ -/* Local variables that need not be saved */ -/* LSAMP is initialized in the IF (FIRST) THEN clause, but it is */ -/* not used the first time through, and it is given a value before -*/ -/* use whenever FIRST is .FALSE., so it appears unnecessary to */ -/* assign it a value when FIRST is .TRUE. */ -/* Local state */ -/* FIRST - .TRUE. only on first call to PITSYN. */ -/* IVOICO - Previous VOICE(2) value. */ -/* IPITO - Previous PITCH value. */ -/* RMSO - Previous RMS value. */ -/* RCO - Previous RC values. */ - -/* JSAMP - If this routine is called N times with identical values of */ -/* LFRAME, then the total length of all pitch periods returned */ -/* is always N*LFRAME-JSAMP, and JSAMP is always in the range 0 -*/ -/* to MAXPIT-1 (see below for why this is so). Thus JSAMP is */ -/* the number of samples "left over" from the previous call to */ -/* PITSYN, that haven't been "used" in a pitch period returned */ -/* from this subroutine. Every time this subroutine is called, -*/ -/* it returns pitch periods with a total length of at most */ -/* LFRAME+JSAMP. */ - -/* IVOICO, IPITO, RCO, and JSAMP need not be assigned an initial value */ -/* with a DATA statement, because they are always initialized on the */ -/* first call to PITSYN. */ - -/* FIRST and RMSO should be initialized with DATA statements, because */ -/* even on the first call, they are used before being initialized. */ - /* Parameter adjustments */ - if (rc) { - --rc; - } - if (rci) { - rci_dim1 = *order; - rci_offset = rci_dim1 + 1; - rci -= rci_offset; - } - if (voice) { - --voice; - } - if (ivuv) { - --ivuv; - } - if (ipiti) { - --ipiti; - } - if (rmsi) { - --rmsi; - } - - /* Function Body */ - ivoico = &(st->ivoico); - ipito = &(st->ipito); - rmso = &(st->rmso); - rco = &(st->rco[0]); - jsamp = &(st->jsamp); - first = &(st->first_pitsyn); - - if (*rms < 1.f) { - *rms = 1.f; - } - if (*rmso < 1.f) { - *rmso = 1.f; - } - uvpit = 0.f; - *ratio = *rms / (*rmso + 8.f); - if (*first) { - lsamp = 0; - ivoice = voice[2]; - if (ivoice == 0) { - *pitch = *lframe / 4; - } - *nout = *lframe / *pitch; - *jsamp = *lframe - *nout * *pitch; - -/* SYNTHS only calls this subroutine with PITCH in the range -20 */ -/* to 156. LFRAME = MAXFRM = 180, so NOUT is somewhere in th -e */ -/* range 1 to 9. */ - -/* JSAMP is "LFRAME mod PITCH", so it is in the range 0 to */ -/* (PITCH-1), or 0 to MAXPIT-1=155, after the first call. */ - - i__1 = *nout; - for (i__ = 1; i__ <= i__1; ++i__) { - i__2 = *order; - for (j = 1; j <= i__2; ++j) { - rci[j + i__ * rci_dim1] = rc[j]; - } - ivuv[i__] = ivoice; - ipiti[i__] = *pitch; - rmsi[i__] = *rms; - } - *first = FALSE_; - } else { - vflag = 0; - lsamp = *lframe + *jsamp; - slope = (*pitch - *ipito) / (real) lsamp; - *nout = 0; - jused = 0; - istart = 1; - if (voice[1] == *ivoico && voice[2] == voice[1]) { - if (voice[2] == 0) { -/* SSUV - - 0 , 0 , 0 */ - *pitch = *lframe / 4; - *ipito = *pitch; - if (*ratio > 8.f) { - *rmso = *rms; - } - } -/* SSVC - - 1 , 1 , 1 */ - slope = (*pitch - *ipito) / (real) lsamp; - ivoice = voice[2]; - } else { - if (*ivoico != 1) { - if (*ivoico == voice[1]) { -/* UV2VC2 - - 0 , 0 , 1 */ - nl = lsamp - *lframe / 4; - } else { -/* UV2VC1 - - 0 , 1 , 1 */ - nl = lsamp - *lframe * 3 / 4; - } - ipiti[1] = nl / 2; - ipiti[2] = nl - ipiti[1]; - ivuv[1] = 0; - ivuv[2] = 0; - rmsi[1] = *rmso; - rmsi[2] = *rmso; - i__1 = *order; - for (i__ = 1; i__ <= i__1; ++i__) { - rci[i__ + rci_dim1] = rco[i__ - 1]; - rci[i__ + (rci_dim1 << 1)] = rco[i__ - 1]; - rco[i__ - 1] = rc[i__]; - } - slope = 0.f; - *nout = 2; - *ipito = *pitch; - jused = nl; - istart = nl + 1; - ivoice = 1; - } else { - if (*ivoico != voice[1]) { -/* VC2UV1 - - 1 , 0 , 0 */ - lsamp = *lframe / 4 + *jsamp; - } else { -/* VC2UV2 - - 1 , 1 , 0 */ - lsamp = *lframe * 3 / 4 + *jsamp; - } - i__1 = *order; - for (i__ = 1; i__ <= i__1; ++i__) { - yarc[i__ - 1] = rc[i__]; - rc[i__] = rco[i__ - 1]; - } - ivoice = 1; - slope = 0.f; - vflag = 1; - } - } -/* Here is the value of most variables that are used below, depending -on */ -/* the values of IVOICO, VOICE(1), and VOICE(2). VOICE(1) and VOICE(2 -) */ -/* are input arguments, and IVOICO is the value of VOICE(2) on the */ -/* previous call (see notes for the IF (NOUT .NE. 0) statement near th -e */ -/* end). Each of these three values is either 0 or 1. These three */ -/* values below are given as 3-bit long strings, in the order IVOICO, -*/ -/* VOICE(1), and VOICE(2). It appears that the code above assumes tha -t */ -/* the bit sequences 010 and 101 never occur, but I wonder whether a -*/ -/* large enough number of bit errors in the channel could cause such a - */ -/* thing to happen, and if so, could that cause NOUT to ever go over 1 -1? */ - -/* Note that all of the 180 values in the table are really LFRAME, but - */ -/* 180 has fewer characters, and it makes the table a little more */ -/* concrete. If LFRAME is ever changed, keep this in mind. Similarly -, */ -/* 135's are 3*LFRAME/4, and 45's are LFRAME/4. If LFRAME is not a */ -/* multiple of 4, then the 135 for NL-JSAMP is actually LFRAME-LFRAME/ -4, */ -/* and the 45 for NL-JSAMP is actually LFRAME-3*LFRAME/4. */ - -/* Note that LSAMP-JSAMP is given as the variable. This was just for -*/ -/* brevity, to avoid adding "+JSAMP" to all of the column entries. */ -/* Similarly for NL-JSAMP. */ - -/* Variable | 000 001 011,010 111 110 100,101 */ -/* ------------+-------------------------------------------------- */ -/* ISTART | 1 NL+1 NL+1 1 1 1 */ -/* LSAMP-JSAMP | 180 180 180 180 135 45 */ -/* IPITO | 45 PITCH PITCH oldPITCH oldPITCH oldPITCH */ -/* SLOPE | 0 0 0 seebelow 0 0 */ -/* JUSED | 0 NL NL 0 0 0 */ -/* PITCH | 45 PITCH PITCH PITCH PITCH PITCH */ -/* NL-JSAMP | -- 135 45 -- -- -- */ -/* VFLAG | 0 0 0 0 1 1 */ -/* NOUT | 0 2 2 0 0 0 */ -/* IVOICE | 0 1 1 1 1 1 */ - -/* while_loop | once once once once twice twice */ - -/* ISTART | -- -- -- -- JUSED+1 JUSED+1 */ -/* LSAMP-JSAMP | -- -- -- -- 180 180 */ -/* IPITO | -- -- -- -- oldPITCH oldPITCH */ -/* SLOPE | -- -- -- -- 0 0 */ -/* JUSED | -- -- -- -- ?? ?? */ -/* PITCH | -- -- -- -- PITCH PITCH */ -/* NL-JSAMP | -- -- -- -- -- -- */ -/* VFLAG | -- -- -- -- 0 0 */ -/* NOUT | -- -- -- -- ?? ?? */ -/* IVOICE | -- -- -- -- 0 0 */ - - -/* UVPIT is always 0.0 on the first pass through the DO WHILE (.TRUE.) - */ -/* loop below. */ - -/* The only possible non-0 value of SLOPE (in column 111) is */ -/* (PITCH-IPITO)/FLOAT(LSAMP) */ - -/* Column 101 is identical to 100. Any good properties we can prove -*/ -/* for 100 will also hold for 101. Similarly for 010 and 011. */ - -/* SYNTHS calls this subroutine with PITCH restricted to the range 20 -to */ -/* 156. IPITO is similarly restricted to this range, after the first -*/ -/* call. IP below is also restricted to this range, given the */ -/* definitions of IPITO, SLOPE, UVPIT, and that I is in the range ISTA -RT */ -/* to LSAMP. */ - - while(TRUE_) { - -/* JUSED is the total length of all pitch periods curr -ently */ -/* in the output arrays, in samples. */ - -/* An invariant of the DO I = ISTART,LSAMP loop below, - under */ -/* the condition that IP is always in the range 1 thro -ugh */ -/* MAXPIT, is: */ - -/* (I - MAXPIT) .LE. JUSED .LE. (I-1) */ - -/* Note that the final value of I is LSAMP+1, so that -after */ -/* the DO loop is complete, we know: */ - -/* (LSAMP - MAXPIT + 1) .LE. JUSED .LE. LSAMP */ - - i__1 = lsamp; - for (i__ = istart; i__ <= i__1; ++i__) { - r__1 = *ipito + slope * i__; - ip = (integer)(r__1 + .5f); - if (uvpit != 0.f) { - ip = (integer)uvpit; - } - if (ip <= i__ - jused) { - ++(*nout); - -/* The following check is no longer nece -ssary, now that */ -/* we can prove that NOUT will never go -over 16. */ - -/* IF (NOUT .GT. 16) STOP 'PITSYN: too many epochs' -*/ - - ipiti[*nout] = ip; - *pitch = ip; - ivuv[*nout] = ivoice; - jused += ip; - prop = (jused - ip / 2) / (real) lsamp; - i__2 = *order; - for (j = 1; j <= i__2; ++j) { - alro = (real)log((rco[j - 1] + 1) / (1 - rco[j - 1])); - alrn = (real)log((rc[j] + 1) / (1 - rc[j])); - xxy = alro + prop * (alrn - alro); - xxy = (real)exp(xxy); - rci[j + *nout * rci_dim1] = (xxy - 1) / (xxy + 1); - } - rmsi[*nout] = (real)(log(*rmso) + prop * (log(*rms) - log(*rmso))); - rmsi[*nout] = (real)exp(rmsi[*nout]); - } - } - if (vflag != 1) { - goto L100; - } - -/* I want to prove what range UVPIT must lie in after -the */ -/* assignments to it below. To do this, I must determ -ine */ -/* what range (LSAMP-ISTART) must lie in, after the */ -/* assignments to ISTART and LSAMP below. */ - -/* Let oldLSAMP be the value of LSAMP at this point in - the */ -/* execution. This is 135+JSAMP in state 110, or 45+J -SAMP in */ -/* states 100 or 101. */ - -/* Given the loop invariant on JUSED above, we know th -at: */ - -/* (oldLSAMP - MAXPIT + 1) .LE. JUSED .LE. oldLSAMP */ - -/* ISTART is one more than this. */ - -/* Let newLSAMP be the value assigned to LSAMP below. - This */ -/* is 180+JSAMP. Thus (newLSAMP-oldLSAMP) is either 4 -5 or */ -/* 135, depending on the state. */ - -/* Thus, the range of newLSAMP-ISTART is: */ - -/* (newLSAMP-(oldLSAMP+1)) .LE. newLSAMP-ISTART */ -/* .LE. (newLSAMP-(oldLSAMP - MAXPIT + 2)) */ - -/* or: */ - -/* 46 .LE. newLSAMP-ISTART .LE. 133+MAXPIT .EQ. 289 */ - -/* Therefore, UVPIT is in the range 23 to 144 after th -e first */ -/* assignment to UVPIT below, and after the conditiona -l */ -/* assignment, it is in the range 23 to 90. */ - -/* The important thing is that it is in the range 20 t -o 156, */ -/* so that in the loop above, IP is always in this ran -ge. */ - - vflag = 0; - istart = jused + 1; - lsamp = *lframe + *jsamp; - slope = 0.f; - ivoice = 0; - uvpit = (real) ((lsamp - istart) / 2); - if (uvpit > 90.f) { - uvpit /= 2; - } - *rmso = *rms; - i__1 = *order; - for (i__ = 1; i__ <= i__1; ++i__) { - rc[i__] = yarc[i__ - 1]; - rco[i__ - 1] = yarc[i__ - 1]; - } - } -L100: - *jsamp = lsamp - jused; - } -/* Given that the maximum pitch period MAXPIT .LT. LFRAME (this is -*/ -/* currently true on every call, since SYNTHS always sets */ -/* LFRAME=180), NOUT will always be .GE. 1 at this point. */ - if (*nout != 0) { - *ivoico = voice[2]; - *ipito = *pitch; - *rmso = *rms; - i__1 = *order; - for (i__ = 1; i__ <= i__1; ++i__) { - rco[i__ - 1] = rc[i__]; - } - } - return 0; -} /* pitsyn_ */ |