/[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 3833 by lgao, Wed Feb 15 03:27:58 2012 UTC revision 3834 by gross, Wed Feb 15 07:09:09 2012 UTC
# Line 1116  void Paso_Preconditioner_AMG_mergeSolve( Line 1116  void Paso_Preconditioner_AMG_mergeSolve(
1116    double* b=NULL;    double* b=NULL;
1117    index_t rank = A->mpi_info->rank;    index_t rank = A->mpi_info->rank;
1118    index_t size = A->mpi_info->size;    index_t size = A->mpi_info->size;
1119    index_t i, n, p, count, n_block;    index_t i, n, p, n_block;
1120    index_t *counts, *offset, *dist;    index_t *counts, *offset, *dist;
1121    
1122    n_block = amg->n_block;    n_block = amg->n_block;
# Line 1137  void Paso_Preconditioner_AMG_mergeSolve( Line 1137  void Paso_Preconditioner_AMG_mergeSolve(
1137      counts[i] = (dist[i+1] - p)*n_block;      counts[i] = (dist[i+1] - p)*n_block;
1138      offset[i] = p*n_block;      offset[i] = p*n_block;
1139    }    }
   count = counts[rank];  
1140    #ifdef ESYS_MPI    #ifdef ESYS_MPI
1141    MPI_Gatherv(amg->b_C, count, MPI_DOUBLE, b, counts, offset, MPI_DOUBLE, 0, A->mpi_info->comm);    {
1142    MPI_Gatherv(amg->x_C, count, MPI_DOUBLE, x, counts, offset, MPI_DOUBLE, 0, A->mpi_info->comm);       index_t count = counts[rank];
1143         MPI_Gatherv(amg->b_C, count, MPI_DOUBLE, b, counts, offset, MPI_DOUBLE, 0, A->mpi_info->comm);
1144         MPI_Gatherv(amg->x_C, count, MPI_DOUBLE, x, counts, offset, MPI_DOUBLE, 0, A->mpi_info->comm);
1145      }
1146    #endif    #endif
1147    
1148    if (rank == 0) {    if (rank == 0) {

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

  ViewVC Help
Powered by ViewVC 1.1.26