/[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 2699 by artak, Wed Sep 30 05:43:20 2009 UTC revision 2701 by artak, Wed Sep 30 23:21:08 2009 UTC
# Line 391  void Paso_Solver_solveAMG(Paso_Solver_AM Line 391  void Paso_Solver_solveAMG(Paso_Solver_AM
391               amg->solver=(void*) ptr;               amg->solver=(void*) ptr;
392         #else               #else      
393           #ifdef MKL           #ifdef MKL
           ptr=(Paso_MKL_Handler *)(amg->solver);  
394            temp=Paso_SparseMatrix_alloc(MATRIX_FORMAT_BLK1 + MATRIX_FORMAT_OFFSET1, amg->A->pattern,1,1, FALSE);            temp=Paso_SparseMatrix_alloc(MATRIX_FORMAT_BLK1 + MATRIX_FORMAT_OFFSET1, amg->A->pattern,1,1, FALSE);
395            #pragma omp parallel for private(i) schedule(static)            #pragma omp parallel for private(i) schedule(static)
396            for (i=0;i<amg->A->len;++i) {            for (i=0;i<amg->A->len;++i) {
397                 temp->val[i]=amg->A->val[i];                 temp->val[i]=amg->A->val[i];
398            }            }
           Paso_MKL1(ptr,temp,x,b,verbose);  
399            amg->solver=(void*) ptr;            amg->solver=(void*) ptr;
400            Paso_SparseMatrix_free(temp);            Paso_SparseMatrix_free(temp);
401           #else           #else

Legend:
Removed from v.2699  
changed lines
  Added in v.2701

  ViewVC Help
Powered by ViewVC 1.1.26