BCPeriodic.c File Reference

Periodic boundary conditions. More...

#include <stdlib.h>
#include <basic.h>
#include <arrayfunctions.h>
#include <mpivars.h>
#include <boundaryconditions.h>

Go to the source code of this file.

Functions
int	BCPeriodicU (void *b, void *m, int ndims, int nvars, int *size, int ghosts, double *phi, double waqt)

Detailed Description

Periodic boundary conditions.

Author

Debojyoti Ghosh

Definition in file BCPeriodic.c.

Function Documentation

int BCPeriodicU	(	void *	b,
		void *	m,
		int	ndims,
		int	nvars,
		int *	size,
		int	ghosts,
		double *	phi,
		double	waqt
	)

Applies periodic boundary conditions: Implemented by copying the solution from the other end of the domain into the physical boundary ghost points.

Note**: This function only acts if the the number of processors is 1 along the spatial dimension this boundary corresponds to. If there are more than 1 processors along this dimension, periodicity is handled by MPIExchangeBoundariesnD() to minimize communication.

Parameters

b	Boundary object of type DomainBoundary
m	MPI object of type MPIVariables
ndims	Number of spatial dimensions
nvars	Number of variables/DoFs per grid point
size	Integer array with the number of grid points in each spatial dimension
ghosts	Number of ghost points
phi	The solution array on which to apply the boundary condition
waqt	Current solution time

Definition at line 19 of file BCPeriodic.c.

{
  DomainBoundary *boundary = (DomainBoundary*) b;
  MPIVariables   *mpi      = (MPIVariables*)   m;

  int dim   = boundary->dim;
  int face  = boundary->face;

  if ((boundary->on_this_proc) && (mpi->iproc[dim] == 1)) {
    int bounds[ndims], index1[ndims], index2[ndims];
    _ArraySubtract1D_(bounds,boundary->ie,boundary->is,ndims);
    _ArraySetValue_(index1,ndims,0);
    _ArraySetValue_(index2,ndims,0);
    int done = 0;
    while (!done) {
      int p1 = 0, p2 = 0;
      _ArrayCopy1D_(index1,index2,ndims);
      if (face == 1) {
        index2[dim] = index1[dim] + size[dim]-ghosts;
        _ArrayIndex1DWO_(ndims,size,index1,boundary->is,ghosts,p1);
        _ArrayIndex1D_(ndims,size,index2,ghosts,p2);
      } else if (face == -1) {
        _ArrayIndex1DWO_(ndims,size,index1,boundary->is,ghosts,p1);
        _ArrayIndex1D_(ndims,size,index1,ghosts,p2);
      }
      _ArrayCopy1D_((phi+nvars*p2),(phi+nvars*p1),nvars);
      _ArrayIncrementIndex_(ndims,bounds,index1,done);
    }
  }

  return(0);
}