28 |
|
|
29 |
/**************************************************************/ |
/**************************************************************/ |
30 |
|
|
|
/* |
|
|
#include "mpi_C.h" |
|
|
#include "Paso.h" |
|
|
#include "Pattern.h" |
|
|
*/ |
|
31 |
#include "PasoUtil.h" |
#include "PasoUtil.h" |
32 |
#include "Pattern_coupling.h" |
#include "Pattern_coupling.h" |
33 |
|
|
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"); |