47 |
dim_t i,j; |
dim_t i,j; |
48 |
/*double threshold=0.05;*/ |
/*double threshold=0.05;*/ |
49 |
index_t iptr,*index,*where_p,diagptr; |
index_t iptr,*index,*where_p,diagptr; |
50 |
bool_t flag,fail; |
bool_t fail; |
51 |
dim_t n=A->pattern->numOutput; |
dim_t n=A->pattern->numOutput; |
52 |
if (A->pattern->type & PATTERN_FORMAT_SYM) { |
if (A->pattern->type & PATTERN_FORMAT_SYM) { |
53 |
Paso_setError(TYPE_ERROR,"Paso_Pattern_mis: symmetric matrix pattern is not supported yet"); |
Paso_setError(TYPE_ERROR,"Paso_Pattern_mis: symmetric matrix pattern is not supported yet"); |
162 |
for (iptr=A->pattern->ptr[i]-index_offset;iptr<A->pattern->ptr[i+1]-index_offset; ++iptr) { |
for (iptr=A->pattern->ptr[i]-index_offset;iptr<A->pattern->ptr[i+1]-index_offset; ++iptr) { |
163 |
if(A->val[iptr-index_offset] < threshold){ |
if(A->val[iptr-index_offset] < threshold){ |
164 |
mis_marker[i]=IS_CONNECTED_TO_MIS; |
mis_marker[i]=IS_CONNECTED_TO_MIS; |
165 |
break; |
#ifndef _OPENMP |
166 |
|
break; |
167 |
|
#endif |
168 |
} |
} |
169 |
} |
} |
170 |
} |
} |