/[escript]/trunk/paso/src/Pattern_coupling.c
ViewVC logotype

Diff of /trunk/paso/src/Pattern_coupling.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1902 by artak, Wed Oct 22 03:54:14 2008 UTC revision 1913 by phornby, Thu Oct 23 08:27:33 2008 UTC
# Line 28  Line 28 
28    
29  /**************************************************************/  /**************************************************************/
30    
31    /*
32  #include "mpi_C.h"  #include "mpi_C.h"
33  #include "Paso.h"  #include "Paso.h"
34  #include "PasoUtil.h"  #include "PasoUtil.h"
35  #include "Pattern.h"  #include "Pattern.h"
36  #include "Solver.h"  */
37    #include "Pattern_coupling.h"
38    
39    
40  /***************************************************************/  /***************************************************************/
# Line 60  void Paso_Pattern_coup(Paso_SparseMatrix Line 62  void Paso_Pattern_coup(Paso_SparseMatrix
62       /* is there any vertex available ?*/       /* is there any vertex available ?*/
63       while (Paso_Util_isAny(n,mis_marker,IS_AVAILABLE)) {       while (Paso_Util_isAny(n,mis_marker,IS_AVAILABLE)) {
64    
65             #pragma omp parallel for private(naib,i,iptr,flag) schedule(static)             #pragma omp parallel for private(i,iptr,flag) schedule(static)
66             for (i=0;i<n;++i) {             for (i=0;i<n;++i) {
67                if (mis_marker[i]==IS_AVAILABLE) {                if (mis_marker[i]==IS_AVAILABLE) {
68                   flag=IS_IN_MIS;                   flag=IS_IN_MIS;
# Line 87  void Paso_Pattern_coup(Paso_SparseMatrix Line 89  void Paso_Pattern_coup(Paso_SparseMatrix
89                  }                  }
90              }              }
91                        
92                #pragma omp parallel for private(naib,i,iptr) schedule(static)                #pragma omp parallel for private(i,iptr) schedule(static)
93                for (i=0;i<n;i++) {                for (i=0;i<n;i++) {
94                 if (mis_marker[i]==IS_AVAILABLE) {                 if (mis_marker[i]==IS_AVAILABLE) {
95                   diagptr=A->pattern->ptr[i];                   diagptr=A->pattern->ptr[i];
# Line 136  void Paso_Pattern_RS(Paso_SparseMatrix* Line 138  void Paso_Pattern_RS(Paso_SparseMatrix*
138  {  {
139    index_t index_offset=(A->pattern->type & PATTERN_FORMAT_OFFSET1 ? 1:0);    index_t index_offset=(A->pattern->type & PATTERN_FORMAT_OFFSET1 ? 1:0);
140    dim_t i,j;    dim_t i,j;
141    index_t naib,iptr;    index_t iptr;
142    double threshold,min_offdiagonal;    double threshold,min_offdiagonal;
143    bool_t flag;    bool_t flag;
144    dim_t n=A->pattern->numOutput;    dim_t n=A->pattern->numOutput;

Legend:
Removed from v.1902  
changed lines
  Added in v.1913

  ViewVC Help
Powered by ViewVC 1.1.26