Отправлено GM 23 февраля 2005 г. 14:09 В ответ на: У Вас не определен массив Y[ ]. Это то же самое что Y2[ ]? отправлено GM 23 февраля 2005 г. 13:27

1) Зачем Вы повторяете вычисление частных произведений во внутреннем цикле
a0=0;
for(i=0;i<=k;i++)
a0+=(long)H[i]*(long)Y[k1-i];
a0>>=14;
temp[k]=(short)a0;

Не проще ли сделать так

a0=0;
for(i=0;i {
a0+=(long)H[i]*(long)Y[k1-i];
temp[k]=(short)(a0>>14);
}

Или даже
a0=0;
for(i=0;i {
a0+=(long)H[i]*(long)Y[k1-i];
temp[k]=(short)(a0>>14)*(short)(a0>>14);
}
Тогда останется всего два цикла, время выполнения существенно сократится. Кстати, за сколько циклов выполняется Ваш вариант кода?
Насколько вы хотите его сократить? И что за процессор, конкретно?

2) Вместо фрагмента в конце внешнего цикла
a0=0;
for(k=0;k a0+=(long)temp[k]*(long)temp[k];

можно написать эквивалентный ассемблерный код
rptz #_IDIM-1
squra *ar4,a

Удивительно, но здесь количество строк ассемблера меньше строк Си-кода!

3) Настало время описать словами, что вы хотите сделать...

4) Вот еще, напоследок. Учтите, что Вы пишете long, а компилятор все равно работает с word, а не с double word.

• Ответ на все вопросы — nyavro   (24.02.2005 09:01, 1207 байт)
  • Весь проект меня не интересует, вот когда вы опубликуете весь код(:-)... — GM   (24.02.2005 14:40, 2982 байт)
    • Разворачивание циклов — nyavro   (25.02.2005 11:36, 651 байт)
      • Я рад, что смог дистанционно помочь. С Вас стакан(:-)...Можно посмотреть на окончательный вариант кода этой подпрограммы? — GM   (25.02.2005 13:21, 1844 байт)
        • Окончательный вариант — nyavro   (02.03.2005 14:38, 2463 байт)
          • Ну вот, более-менее пригладили... — GM   (02.03.2005 20:43, 1484 байт)
• Опять где-то глюкануло, повторю, что не прошло — GM   (23.02.2005 14:46, 660 байт)