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

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

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

revision 3834 by gross, Wed Feb 15 07:09:09 2012 UTC revision 3835 by gross, Wed Feb 15 23:19:01 2012 UTC
# Line 1118  void Paso_Preconditioner_AMG_mergeSolve( Line 1118  void Paso_Preconditioner_AMG_mergeSolve(
1118    index_t size = A->mpi_info->size;    index_t size = A->mpi_info->size;
1119    index_t i, n, p, n_block;    index_t i, n, p, n_block;
1120    index_t *counts, *offset, *dist;    index_t *counts, *offset, *dist;
1121      #ifdef ESYS_MPI
1122      index_t count;
1123      #endif
1124    n_block = amg->n_block;    n_block = amg->n_block;
   
1125    A_D = Paso_Preconditioner_AMG_mergeSystemMatrix(A);    A_D = Paso_Preconditioner_AMG_mergeSystemMatrix(A);
1126    
1127    /* First, gather x and b into rank 0 */    /* First, gather x and b into rank 0 */
# Line 1138  void Paso_Preconditioner_AMG_mergeSolve( Line 1139  void Paso_Preconditioner_AMG_mergeSolve(
1139      offset[i] = p*n_block;      offset[i] = p*n_block;
1140    }    }
1141    #ifdef ESYS_MPI    #ifdef ESYS_MPI
1142    {    count = counts[rank];
1143       index_t count = counts[rank];    MPI_Gatherv(amg->b_C, count, MPI_DOUBLE, b, counts, offset, MPI_DOUBLE, 0, A->mpi_info->comm);
1144       MPI_Gatherv(amg->b_C, count, MPI_DOUBLE, b, counts, offset, MPI_DOUBLE, 0, A->mpi_info->comm);    MPI_Gatherv(amg->x_C, count, MPI_DOUBLE, x, counts, offset, MPI_DOUBLE, 0, A->mpi_info->comm);
      MPI_Gatherv(amg->x_C, count, MPI_DOUBLE, x, counts, offset, MPI_DOUBLE, 0, A->mpi_info->comm);  
   }  
1145    #endif    #endif
1146    
1147    if (rank == 0) {    if (rank == 0) {

Legend:
Removed from v.3834  
changed lines
  Added in v.3835

  ViewVC Help
Powered by ViewVC 1.1.26