- for(k = 1; k <= K; k++){
- for(i = 1; i <= pow(2,k-1); i++){
- matno=get_y(k,2*i-1);
- CDEBUG2(surf_sdp_out,"%d %d 1 1 1\n", matno , block_num);
- addentry(constraints, &C, matno, block_num, 1, 1, 1.0, C.blocks[block_num].blocksize);
-
- matno=get_y(k,2*i);
- CDEBUG2(surf_sdp_out,"%d %d 2 2 1\n", matno , block_num);
- addentry(constraints, &C, matno, block_num, 2, 2, 1.0, C.blocks[block_num].blocksize);
-
- matno=get_y(k-1,i);
- CDEBUG2(surf_sdp_out,"%d %d 1 2 1\n", matno , block_num);
- addentry(constraints, &C, matno, block_num, 1, 2, 1.0, C.blocks[block_num].blocksize);
-
- matno=get_y(k-1,i);
- CDEBUG2(surf_sdp_out,"%d %d 2 1 1\n", matno , block_num);
- addentry(constraints, &C, matno, block_num, 2, 1, 1.0, C.blocks[block_num].blocksize);
-
+ for (k = 1; k <= K; k++) {
+ for (i = 1; i <= pow(2, k - 1); i++) {
+ matno = get_y(k, 2 * i - 1);
+ CDEBUG2(surf_sdp_out, "%d %d 1 1 1", matno, block_num);
+ addentry(constraints, &C, matno, block_num, 1, 1, 1.0,
+ C.blocks[block_num].blocksize);
+
+ matno = get_y(k, 2 * i);
+ CDEBUG2(surf_sdp_out, "%d %d 2 2 1", matno, block_num);
+ addentry(constraints, &C, matno, block_num, 2, 2, 1.0,
+ C.blocks[block_num].blocksize);
+
+ matno = get_y(k - 1, i);
+ CDEBUG2(surf_sdp_out, "%d %d 1 2 1", matno, block_num);
+ addentry(constraints, &C, matno, block_num, 1, 2, 1.0,
+ C.blocks[block_num].blocksize);
+
+ matno = get_y(k - 1, i);
+ CDEBUG2(surf_sdp_out, "%d %d 2 1 1", matno, block_num);
+ addentry(constraints, &C, matno, block_num, 2, 1, 1.0,
+ C.blocks[block_num].blocksize);
+