132 |
void Paso_Pattern_RS(Paso_SparseMatrix* A, index_t* mis_marker, double theta) |
void Paso_Pattern_RS(Paso_SparseMatrix* A, index_t* mis_marker, double theta) |
133 |
{ |
{ |
134 |
index_t index_offset=(A->pattern->type & PATTERN_FORMAT_OFFSET1 ? 1:0); |
index_t index_offset=(A->pattern->type & PATTERN_FORMAT_OFFSET1 ? 1:0); |
135 |
dim_t i,j; |
dim_t i; |
136 |
index_t iptr; |
index_t iptr; |
137 |
double threshold,min_offdiagonal; |
double threshold,min_offdiagonal; |
|
bool_t flag; |
|
138 |
dim_t n=A->pattern->numOutput; |
dim_t n=A->pattern->numOutput; |
139 |
if (A->pattern->type & PATTERN_FORMAT_SYM) { |
if (A->pattern->type & PATTERN_FORMAT_SYM) { |
140 |
Paso_setError(TYPE_ERROR,"Paso_Pattern_RS: symmetric matrix pattern is not supported yet"); |
Paso_setError(TYPE_ERROR,"Paso_Pattern_RS: symmetric matrix pattern is not supported yet"); |