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);
|
