p61743098
/
sagecal
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

bug fix

This commit is contained in:
Sarod Yatawatta 2017-11-21 21:48:28 +01:00
parent b24279190c
commit 15ec822ec4
1 changed files with 16 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

32
src/lib/predict_model.cu
View File

@ -902,28 +902,28 @@ kernel_correct_residuals(int B, int N, int Nb, int boff, int F, int nchunk, doub
for (int cf=0; cf<F; cf++) { for (int cf=0; cf<F; cf++) {
cuDoubleComplex L1[4],L2[4]; cuDoubleComplex L1[4],L2[4];
L1[0].x=x[cf*8*Nb+n]; L1[0].x=x[cf*8*Nb+8*n];
L1[0].y=x[cf*8*Nb+n+1]; L1[0].y=x[cf*8*Nb+8*n+1];
L1[1].x=x[cf*8*Nb+n+2]; L1[1].x=x[cf*8*Nb+8*n+2];
L1[1].y=x[cf*8*Nb+n+3]; L1[1].y=x[cf*8*Nb+8*n+3];
L1[2].x=x[cf*8*Nb+n+4]; L1[2].x=x[cf*8*Nb+8*n+4];
L1[2].y=x[cf*8*Nb+n+5]; L1[2].y=x[cf*8*Nb+8*n+5];
L1[3].x=x[cf*8*Nb+n+6]; L1[3].x=x[cf*8*Nb+8*n+6];
L1[3].y=x[cf*8*Nb+n+7]; L1[3].y=x[cf*8*Nb+8*n+7];
/* L2=G1*L1 */ /* L2=G1*L1 */
amb(G1,L1,L2); amb(G1,L1,L2);
/* L1=L2*G2^H */ /* L1=L2*G2^H */
ambt(L2,G2,L1); ambt(L2,G2,L1);
x[cf*8*Nb+n]=L1[0].x; x[cf*8*Nb+8*n]=L1[0].x;
x[cf*8*Nb+n+1]=L1[0].y; x[cf*8*Nb+8*n+1]=L1[0].y;
x[cf*8*Nb+n+2]=L1[1].x; x[cf*8*Nb+8*n+2]=L1[1].x;
x[cf*8*Nb+n+3]=L1[1].y; x[cf*8*Nb+8*n+3]=L1[1].y;
x[cf*8*Nb+n+4]=L1[2].x; x[cf*8*Nb+8*n+4]=L1[2].x;
x[cf*8*Nb+n+5]=L1[2].y; x[cf*8*Nb+8*n+5]=L1[2].y;
x[cf*8*Nb+n+6]=L1[3].x; x[cf*8*Nb+8*n+6]=L1[3].x;
x[cf*8*Nb+n+7]=L1[3].y; x[cf*8*Nb+8*n+7]=L1[3].y;
} }
} }