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

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

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

revision 3488 by gross, Tue Dec 7 07:36:12 2010 UTC revision 3489 by caltinay, Wed Mar 30 00:46:04 2011 UTC
# Line 78  Paso_Solver_RILU* Paso_Solver_getRILU(Pa Line 78  Paso_Solver_RILU* Paso_Solver_getRILU(Pa
78    index_t iPtr,*index, *where_p;    index_t iPtr,*index, *where_p;
79    dim_t i,k;    dim_t i,k;
80    Paso_SparseMatrix * schur=NULL;    Paso_SparseMatrix * schur=NULL;
81    double A11,A12,A13,A21,A22,A23,A31,A32,A33,D,time0=0,time1=0,time2=0;    double A11,A12,A13,A21,A22,A23,A31,A32,A33,D,time0=0,time1=0;/*,time2=0;*/
82        
83    
84    /* identify independend set of rows/columns */    /* identify independend set of rows/columns */
# Line 105  Paso_Solver_RILU* Paso_Solver_getRILU(Pa Line 105  Paso_Solver_RILU* Paso_Solver_getRILU(Pa
105       #pragma omp parallel for private(i) schedule(static)       #pragma omp parallel for private(i) schedule(static)
106       for (i=0;i<n;++i) mis_marker[i]=-1;       for (i=0;i<n;++i) mis_marker[i]=-1;
107       Paso_Pattern_mis(A_p->pattern,mis_marker);       Paso_Pattern_mis(A_p->pattern,mis_marker);
108       time2=Esys_timer()-time0;       /*time2=Esys_timer()-time0;*/
109       if (Esys_noError()) {       if (Esys_noError()) {
110          #pragma omp parallel for private(i) schedule(static)          #pragma omp parallel for private(i) schedule(static)
111          for (i = 0; i < n; ++i) counter[i]=mis_marker[i];          for (i = 0; i < n; ++i) counter[i]=mis_marker[i];

Legend:
Removed from v.3488  
changed lines
  Added in v.3489

  ViewVC Help
Powered by ViewVC 1.1.26