Commit b6ba7b7a authored by Martin Reinecke's avatar Martin Reinecke
Browse files

fix

parent 1d4deeb6
...@@ -969,7 +969,7 @@ template<typename T, typename Serv> void x2grid_c ...@@ -969,7 +969,7 @@ template<typename T, typename Serv> void x2grid_c
{ {
complex<T> tmp(v*ku[cu]); complex<T> tmp(v*ku[cu]);
size_t cv=0; size_t cv=0;
for (; cv<supp-3; cv+=4) for (; cv+3<supp; cv+=4)
{ {
ptr[cv ] += tmp*kv[cv ]; ptr[cv ] += tmp*kv[cv ];
ptr[cv+1] += tmp*kv[cv+1]; ptr[cv+1] += tmp*kv[cv+1];
...@@ -1027,7 +1027,7 @@ template<typename T, typename Serv> void grid2x_c ...@@ -1027,7 +1027,7 @@ template<typename T, typename Serv> void grid2x_c
{ {
complex<T> tmp(0); complex<T> tmp(0);
size_t cv=0; size_t cv=0;
for (; cv<supp-3; cv+=4) for (; cv+3<supp; cv+=4)
tmp += ptr[cv ]*kv[cv ] tmp += ptr[cv ]*kv[cv ]
+ ptr[cv+1]*kv[cv+1] + ptr[cv+1]*kv[cv+1]
+ ptr[cv+2]*kv[cv+2] + ptr[cv+2]*kv[cv+2]
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment