test_grid.c File Reference

C unit tests for grid bounding-box helper routines. More...

#include "test_support.h"
#include "grid.h"
#include <stdlib.h>

Include dependency graph for test_grid.c:

Go to the source code of this file.

Functions
static PetscErrorCode	TestComputeLocalBoundingBoxUniformGrid (void)
	Tests local bounding-box construction on a uniform grid.
static PetscErrorCode	TestGatherAndBroadcastBoundingBoxes (void)
	Tests collective gather and broadcast of rank-local bounding boxes.
int	main (int argc, char **argv)
	Runs the unit-grid PETSc test binary.

Detailed Description

C unit tests for grid bounding-box helper routines.

Definition in file test_grid.c.

Function Documentation

TestComputeLocalBoundingBoxUniformGrid()

static PetscErrorCode TestComputeLocalBoundingBoxUniformGrid ( void  )

static

Tests local bounding-box construction on a uniform grid.

Definition at line 15 of file test_grid.c.

{
    SimCtx *simCtx = NULL;
    UserCtx *user = NULL;
    BoundingBox bbox;
    PetscReal expected_max = 0.0;
 
    PetscFunctionBeginUser;
    PetscCall(PicurvCreateMinimalContexts(&simCtx, &user, 4, 4, 4));
    PetscCall(ComputeLocalBoundingBox(user, &bbox));
    expected_max = ((PetscReal)(user->IM - 1) / (PetscReal)user->IM) + 1.0e-6;
 
    PetscCall(PicurvAssertRealNear(-1.0e-6, bbox.min_coords.x, 1.0e-10, "bbox min x"));
    PetscCall(PicurvAssertRealNear(-1.0e-6, bbox.min_coords.y, 1.0e-10, "bbox min y"));
    PetscCall(PicurvAssertRealNear(-1.0e-6, bbox.min_coords.z, 1.0e-10, "bbox min z"));
    PetscCall(PicurvAssertRealNear(expected_max, bbox.max_coords.x, 1.0e-10, "bbox max x"));
    PetscCall(PicurvAssertRealNear(expected_max, bbox.max_coords.y, 1.0e-10, "bbox max y"));
    PetscCall(PicurvAssertRealNear(expected_max, bbox.max_coords.z, 1.0e-10, "bbox max z"));
 
    PetscCall(PicurvDestroyMinimalContexts(&simCtx, &user));
    PetscFunctionReturn(0);
}

Here is the call graph for this function:

Here is the caller graph for this function:

TestGatherAndBroadcastBoundingBoxes()

static PetscErrorCode TestGatherAndBroadcastBoundingBoxes ( void  )

static

Tests collective gather and broadcast of rank-local bounding boxes.

Definition at line 41 of file test_grid.c.

{
    SimCtx *simCtx = NULL;
    UserCtx *user = NULL;
    BoundingBox local_bbox;
    BoundingBox *boxes = NULL;
    PetscMPIInt rank = 0, size = 1;
 
    PetscFunctionBeginUser;
    PetscCallMPI(MPI_Comm_rank(PETSC_COMM_WORLD, &rank));
    PetscCallMPI(MPI_Comm_size(PETSC_COMM_WORLD, &size));
 
    PetscCall(PicurvCreateMinimalContexts(&simCtx, &user, 4, 4, 4));
    PetscCall(ComputeLocalBoundingBox(user, &local_bbox));
    PetscCall(GatherAllBoundingBoxes(user, &boxes));
 
    if (rank == 0) {
        PetscCall(PicurvAssertBool((PetscBool)(boxes != NULL), "root must receive gathered bounding boxes"));
    } else {
        PetscCall(PicurvAssertBool((PetscBool)(boxes == NULL), "non-root should have NULL pre-broadcast bbox pointer"));
    }
 
    PetscCall(BroadcastAllBoundingBoxes(user, &boxes));
    PetscCall(PicurvAssertBool((PetscBool)(boxes != NULL), "all ranks must have bbox array after broadcast"));
 
    for (PetscMPIInt r = 0; r < size; ++r) {
        PetscCall(PicurvAssertBool((PetscBool)(boxes[r].min_coords.x <= boxes[r].max_coords.x), "bbox x range is valid"));
        PetscCall(PicurvAssertBool((PetscBool)(boxes[r].min_coords.y <= boxes[r].max_coords.y), "bbox y range is valid"));
        PetscCall(PicurvAssertBool((PetscBool)(boxes[r].min_coords.z <= boxes[r].max_coords.z), "bbox z range is valid"));
    }
 
    PetscCall(PicurvAssertRealNear(local_bbox.min_coords.x, boxes[rank].min_coords.x, 1.0e-10, "rank-local bbox min x"));
    PetscCall(PicurvAssertRealNear(local_bbox.max_coords.x, boxes[rank].max_coords.x, 1.0e-10, "rank-local bbox max x"));
 
    free(boxes);
    PetscCall(PicurvDestroyMinimalContexts(&simCtx, &user));
    PetscFunctionReturn(0);
}

Here is the call graph for this function:

Here is the caller graph for this function:

main()

int main	(	int	argc,
		char **	argv
	)

Runs the unit-grid PETSc test binary.

Definition at line 83 of file test_grid.c.

{
    PetscErrorCode ierr;
    const PicurvTestCase cases[] = {
        {"compute-local-bbox-uniform-grid", TestComputeLocalBoundingBoxUniformGrid},
        {"gather-and-broadcast-bboxes", TestGatherAndBroadcastBoundingBoxes},
    };
 
    ierr = PetscInitialize(&argc, &argv, NULL, "PICurv grid tests");
    if (ierr) {
        return (int)ierr;
    }
 
    ierr = PicurvRunTests("unit-grid", cases, sizeof(cases) / sizeof(cases[0]));
    if (ierr) {
        PetscFinalize();
        return (int)ierr;
    }
 
    ierr = PetscFinalize();
    return (int)ierr;
}

Here is the call graph for this function: