|
From: Martin R. <ru...@us...> - 2004-08-27 17:03:03
|
Update of /cvsroot/foo/foo/libfoo/modules/orthodox In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv17755 Modified Files: FOOMLookup.m FOOMReadTranspSnd.m FOOMTransposeSnd.m Log Message: corrected sinc lookup loops Index: FOOMLookup.m =================================================================== RCS file: /cvsroot/foo/foo/libfoo/modules/orthodox/FOOMLookup.m,v retrieving revision 1.6 retrieving revision 1.7 diff -C2 -d -r1.6 -r1.7 *** FOOMLookup.m 27 Aug 2004 07:00:46 -0000 1.6 --- FOOMLookup.m 27 Aug 2004 17:02:47 -0000 1.7 *************** *** 169,173 **** } tdelta = 1 - tdelta; ! for (si = 1; si < tot; si++, ti2 += tinc) { sum += s[si] * (base[ti2] + delta[ti2] * tdelta); --- 169,173 ---- } tdelta = 1 - tdelta; ! for (si = 1; si <= tot; si++, ti2 += tinc) { sum += s[si] * (base[ti2] + delta[ti2] * tdelta); Index: FOOMTransposeSnd.m =================================================================== RCS file: /cvsroot/foo/foo/libfoo/modules/orthodox/FOOMTransposeSnd.m,v retrieving revision 1.5 retrieving revision 1.6 diff -C2 -d -r1.5 -r1.6 *** FOOMTransposeSnd.m 15 Aug 2004 17:42:37 -0000 1.5 --- FOOMTransposeSnd.m 27 Aug 2004 17:02:47 -0000 1.6 *************** *** 199,203 **** P = 1. / lp - P; ! for (fti = P * tinc, si = 1; si < ptot; si++, fti += ptinc) { ti = (int)fti; --- 199,203 ---- P = 1. / lp - P; ! for (fti = P * tinc, si = 1; si <= ptot; si++, fti += ptinc) { ti = (int)fti; *************** *** 235,239 **** tdelta = 1 - tdelta; ! for (si = 1; si < tot; si++, ti2 += tinc) { sum += s[si] * (base[ti2] + delta[ti2] * tdelta); --- 235,239 ---- tdelta = 1 - tdelta; ! for (si = 1; si <= tot; si++, ti2 += tinc) { sum += s[si] * (base[ti2] + delta[ti2] * tdelta); Index: FOOMReadTranspSnd.m =================================================================== RCS file: /cvsroot/foo/foo/libfoo/modules/orthodox/FOOMReadTranspSnd.m,v retrieving revision 1.5 retrieving revision 1.6 diff -C2 -d -r1.5 -r1.6 *** FOOMReadTranspSnd.m 15 Aug 2004 17:42:37 -0000 1.5 --- FOOMReadTranspSnd.m 27 Aug 2004 17:02:47 -0000 1.6 *************** *** 288,292 **** } P = 1. / lp - P; ! for (fti = P * tinc, si = 1; si < ptot; si++, fti += ptinc) { ti = (int) fti; --- 288,292 ---- } P = 1. / lp - P; ! for (fti = P * tinc, si = 1; si <= ptot; si++, fti += ptinc) { ti = (int) fti; *************** *** 315,319 **** } tdelta = 1 - tdelta; ! for (si = 1; si < tot; si++, ti2 += tinc) { sum += s[si] * (base[ti2] + delta[ti2] * tdelta); --- 315,319 ---- } tdelta = 1 - tdelta; ! for (si = 1; si <= tot; si++, ti2 += tinc) { sum += s[si] * (base[ti2] + delta[ti2] * tdelta); |