461 |
register double sum_row=0; |
register double sum_row=0; |
462 |
register double main_row=0; |
register double main_row=0; |
463 |
register dim_t kdeg=0; |
register dim_t kdeg=0; |
464 |
const register index_t koffset=A->mainBlock->pattern->ptr[i]+A->col_coupleBlock->pattern->ptr[i]; |
register const index_t koffset=A->mainBlock->pattern->ptr[i]+A->col_coupleBlock->pattern->ptr[i]; |
465 |
|
|
466 |
/* collect information for row i: */ |
/* collect information for row i: */ |
467 |
#pragma ivdep |
#pragma ivdep |
532 |
|
|
533 |
const double threshold = remote_threshold[2*i+1]; |
const double threshold = remote_threshold[2*i+1]; |
534 |
register dim_t kdeg=0; |
register dim_t kdeg=0; |
535 |
const register index_t koffset=koffset_0+A->row_coupleBlock->pattern->ptr[i]; |
register const index_t koffset=koffset_0+A->row_coupleBlock->pattern->ptr[i]; |
536 |
if (remote_threshold[2*i]>0) { |
if (remote_threshold[2*i]>0) { |
537 |
#pragma ivdep |
#pragma ivdep |
538 |
for (iptr=A->row_coupleBlock->pattern->ptr[i];iptr<A->row_coupleBlock->pattern->ptr[i+1]; ++iptr) { |
for (iptr=A->row_coupleBlock->pattern->ptr[i];iptr<A->row_coupleBlock->pattern->ptr[i+1]; ++iptr) { |
593 |
register double main_row=0; |
register double main_row=0; |
594 |
register index_t rtmp_offset=-A->mainBlock->pattern->ptr[i]; |
register index_t rtmp_offset=-A->mainBlock->pattern->ptr[i]; |
595 |
register dim_t kdeg=0; |
register dim_t kdeg=0; |
596 |
const register index_t koffset=A->mainBlock->pattern->ptr[i]+A->col_coupleBlock->pattern->ptr[i]; |
register const index_t koffset=A->mainBlock->pattern->ptr[i]+A->col_coupleBlock->pattern->ptr[i]; |
597 |
|
|
598 |
/* collect information for row i: */ |
/* collect information for row i: */ |
599 |
for (iptr=A->mainBlock->pattern->ptr[i];iptr<A->mainBlock->pattern->ptr[i+1]; ++iptr) { |
for (iptr=A->mainBlock->pattern->ptr[i];iptr<A->mainBlock->pattern->ptr[i+1]; ++iptr) { |
685 |
|
|
686 |
const double threshold2 = remote_threshold[2*i+1]*remote_threshold[2*i+1]; |
const double threshold2 = remote_threshold[2*i+1]*remote_threshold[2*i+1]; |
687 |
register dim_t kdeg=0; |
register dim_t kdeg=0; |
688 |
const register index_t koffset=koffset_0+A->row_coupleBlock->pattern->ptr[i]; |
register const index_t koffset=koffset_0+A->row_coupleBlock->pattern->ptr[i]; |
689 |
if (remote_threshold[2*i]>0) { |
if (remote_threshold[2*i]>0) { |
690 |
#pragma ivdep |
#pragma ivdep |
691 |
for (iptr=A->row_coupleBlock->pattern->ptr[i];iptr<A->row_coupleBlock->pattern->ptr[i+1]; ++iptr) { |
for (iptr=A->row_coupleBlock->pattern->ptr[i];iptr<A->row_coupleBlock->pattern->ptr[i+1]; ++iptr) { |