CFD/hypre
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fixed problem when a grid is distributed accross only a subset of

Browse Source 
of the processors.  This happens often in the context of "level solves"
in an AMR code.
This commit is contained in:
falgout 1999-06-30 23:04:53 +00:00
parent 2eea77ce6e
commit 44c159b60f
2 changed files with 20 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
struct_matrix_vector/struct_grid.c
View File

@ -177,7 +177,7 @@ hypre_AssembleStructGrid( hypre_StructGrid *grid )
int ierr = 0;
MPI_Comm comm = hypre_StructGridComm(grid);
hypre_BoxArray *boxes = hypre_StructGridBoxes(grid);
hypre_BoxArray *boxes;
int global_size;
int local_size;
hypre_BoxNeighbors *neighbors;
@ -192,6 +192,15 @@ hypre_AssembleStructGrid( hypre_StructGrid *grid )
hypre_BoxArray *periodic_all_boxes;
int *periodic_processes;
/* If I don't own any boxes, need to add a "zero-sized" */
if (hypre_StructGridBoxes(grid) == NULL)
{
hypre_Index ilower, iupper;
hypre_SetIndex(ilower, 99, 99, 99);
hypre_SetIndex(iupper, -99, -99, -99);
hypre_SetStructGridExtents(grid, ilower, iupper);
}
boxes = hypre_StructGridBoxes(grid);
if (hypre_StructGridAllBoxes(grid) == NULL)

11
struct_mv/struct_grid.c
View File

@ -177,7 +177,7 @@ hypre_AssembleStructGrid( hypre_StructGrid *grid )
int ierr = 0;
MPI_Comm comm = hypre_StructGridComm(grid);
hypre_BoxArray *boxes = hypre_StructGridBoxes(grid);
hypre_BoxArray *boxes;
int global_size;
int local_size;
hypre_BoxNeighbors *neighbors;
@ -192,6 +192,15 @@ hypre_AssembleStructGrid( hypre_StructGrid *grid )
hypre_BoxArray *periodic_all_boxes;
int *periodic_processes;
/* If I don't own any boxes, need to add a "zero-sized" */
if (hypre_StructGridBoxes(grid) == NULL)
{
hypre_Index ilower, iupper;
hypre_SetIndex(ilower, 99, 99, 99);
hypre_SetIndex(iupper, -99, -99, -99);
hypre_SetStructGridExtents(grid, ilower, iupper);
}
boxes = hypre_StructGridBoxes(grid);
if (hypre_StructGridAllBoxes(grid) == NULL)