Исправление bormant, (текущая версия) :
почему?
https://godbolt.org/z/Cvbn5O
Попробуйте, достаточно ли будет такого изменения:
+ double s = 0.0;
for (k = 0; k < N; ++k)
- *C += *a++ * *t++;
+ s += *a++ * *t++;
- ++C;
+ *C++ = s;
Исходная версия bormant, :