120121122123124125126127128129130
data[k][0]=1.0+ndev.dev(); data[k][1]=1.0+ndev.dev(); } else if (j == 1) { y[k]=-1.0; data[k][0]=-1.0+ndev.dev(); data[k][1]=1.0+ndev.dev(); } else if (j == 2) { y[k]=1.0; data[k][0]=-1.0+ndev.dev(); data[k][1]=-1.0+ndev.dev(); } else {
123124125126127128129130131132133
y[k]=-1.0; data[k][0]=-1.0+ndev.dev(); data[k][1]=1.0+ndev.dev(); } else if (j == 2) { y[k]=1.0; data[k][0]=-1.0+ndev.dev(); data[k][1]=-1.0+ndev.dev(); } else { y[k]=-1.0; data[k][0]=1.0+ndev.dev(); data[k][1]=-1.0+ndev.dev();
124125126127128129130131132133134
data[k][0]=-1.0+ndev.dev(); data[k][1]=1.0+ndev.dev(); } else if (j == 2) { y[k]=1.0; data[k][0]=-1.0+ndev.dev(); data[k][1]=-1.0+ndev.dev(); } else { y[k]=-1.0; data[k][0]=1.0+ndev.dev(); data[k][1]=-1.0+ndev.dev(); }
127128129130131132133134135136137
y[k]=1.0; data[k][0]=-1.0+ndev.dev(); data[k][1]=-1.0+ndev.dev(); } else { y[k]=-1.0; data[k][0]=1.0+ndev.dev(); data[k][1]=-1.0+ndev.dev(); } } }
128129130131132133134135136137138
data[k][0]=-1.0+ndev.dev(); data[k][1]=-1.0+ndev.dev(); } else { y[k]=-1.0; data[k][0]=1.0+ndev.dev(); data[k][1]=-1.0+ndev.dev(); } } } // Linear kernel
155156157158159160161162163164165
nerror=0; for (j=0;j<4;j++) { // Four quadrants for (i=0;i<M/4;i++) { if (j == 0) { yy=1.0; x[0]=1.0+ndev.dev(); x[1]=1.0+ndev.dev(); } else if (j == 1) { yy=-1.0; x[0]=-1.0+ndev.dev(); x[1]=1.0+ndev.dev();
156157158159160161162163164165166
for (j=0;j<4;j++) { // Four quadrants for (i=0;i<M/4;i++) { if (j == 0) { yy=1.0; x[0]=1.0+ndev.dev(); x[1]=1.0+ndev.dev(); } else if (j == 1) { yy=-1.0; x[0]=-1.0+ndev.dev(); x[1]=1.0+ndev.dev(); } else if (j == 2) {
159160161162163164165166167168169
yy=1.0; x[0]=1.0+ndev.dev(); x[1]=1.0+ndev.dev(); } else if (j == 1) { yy=-1.0; x[0]=-1.0+ndev.dev(); x[1]=1.0+ndev.dev(); } else if (j == 2) { yy=1.0; x[0]=-1.0+ndev.dev(); x[1]=-1.0+ndev.dev();
160161162163164165166167168169170
x[0]=1.0+ndev.dev(); x[1]=1.0+ndev.dev(); } else if (j == 1) { yy=-1.0; x[0]=-1.0+ndev.dev(); x[1]=1.0+ndev.dev(); } else if (j == 2) { yy=1.0; x[0]=-1.0+ndev.dev(); x[1]=-1.0+ndev.dev(); } else {
163164165166167168169170171172173
yy=-1.0; x[0]=-1.0+ndev.dev(); x[1]=1.0+ndev.dev(); } else if (j == 2) { yy=1.0; x[0]=-1.0+ndev.dev(); x[1]=-1.0+ndev.dev(); } else { yy=-1.0; x[0]=1.0+ndev.dev(); x[1]=-1.0+ndev.dev();