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

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

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

revision 2122 by artak, Wed Dec 3 02:52:28 2008 UTC revision 2131 by artak, Thu Dec 4 06:09:50 2008 UTC
# Line 140  Paso_Solver_AMG* Paso_Solver_getAMG(Paso Line 140  Paso_Solver_AMG* Paso_Solver_getAMG(Paso
140                }                }
141                                
142                /* Compute row-sum for getting rs(A_FF)^-1*/                /* Compute row-sum for getting rs(A_FF)^-1*/
143                #pragma omp parallel for private(i,iPtr,j) schedule(static)                #pragma omp parallel for private(i,iPtr,j,S) schedule(static)
144                for (i = 0; i < out->n_F; ++i) {                for (i = 0; i < out->n_F; ++i) {
145                  S=0;                  S=0;
146                  for (iPtr=A_p->pattern->ptr[out->rows_in_F[i]];iPtr<A_p->pattern->ptr[out->rows_in_F[i] + 1]; ++iPtr) {                  for (iPtr=A_p->pattern->ptr[out->rows_in_F[i]];iPtr<A_p->pattern->ptr[out->rows_in_F[i] + 1]; ++iPtr) {
# Line 186  Paso_Solver_AMG* Paso_Solver_getAMG(Paso Line 186  Paso_Solver_AMG* Paso_Solver_getAMG(Paso
186    
187                              schur_withFillIn=Paso_SparseMatrix_alloc(A_p->type,Paso_Pattern_binop(PATTERN_FORMAT_DEFAULT, schur->pattern, Paso_Pattern_multiply(PATTERN_FORMAT_DEFAULT,out->A_CF->pattern,out->A_FC->pattern)),1,1);                              schur_withFillIn=Paso_SparseMatrix_alloc(A_p->type,Paso_Pattern_binop(PATTERN_FORMAT_DEFAULT, schur->pattern, Paso_Pattern_multiply(PATTERN_FORMAT_DEFAULT,out->A_CF->pattern,out->A_FC->pattern)),1,1);
188                                                            
189                              fprintf(stderr,"Sparsity of Schure: %dx%d LEN %d Percentage %f\n",schur_withFillIn->pattern->numOutput,schur_withFillIn->pattern->numInput,schur_withFillIn->len,schur_withFillIn->len/(1.*schur_withFillIn->pattern->numOutput*schur_withFillIn->pattern->numInput));                              /*fprintf(stderr,"Sparsity of Schure: %dx%d LEN %d Percentage %f\n",schur_withFillIn->pattern->numOutput,schur_withFillIn->pattern->numInput,schur_withFillIn->len,schur_withFillIn->len/(1.*schur_withFillIn->pattern->numOutput*schur_withFillIn->pattern->numInput));*/
190                                                            
191                              /* copy values over*/                              /* copy values over*/
192                              #pragma omp parallel for private(i,iPtr,j,iPtr_s,index,where_p) schedule(static)                              #pragma omp parallel for private(i,iPtr,j,iPtr_s,index,where_p) schedule(static)

Legend:
Removed from v.2122  
changed lines
  Added in v.2131

  ViewVC Help
Powered by ViewVC 1.1.26