ApplyBoundaryConditions.c File Reference

Apply physical boundary conditions to domain. More...

#include <stdio.h>
#include <basic.h>
#include <arrayfunctions.h>
#include <mpivars.h>
#include <boundaryconditions.h>
#include <hypar.h>

Go to the source code of this file.

Functions
int	ApplyBoundaryConditions (void *s, void *m, double *x, double *xref, double waqt)
	Applies the boundary conditions specified for each boundary zone. More...

Detailed Description

Apply physical boundary conditions to domain.

Author

Debojyoti Ghosh Contains the function that applies the physical boundary conditions to each boundary zone.

Definition in file ApplyBoundaryConditions.c.

Function Documentation

int ApplyBoundaryConditions	(	void *	s,
		void *	m,
		double *	x,
		double *	xref,
		double	waqt
	)

Applies the boundary conditions specified for each boundary zone.

The solver object (of type HyPar) contains an oject of type DomainBoundary that contains all the boundary information (dimension, extent, face, type, etc). This function iterates through each of the boundary zones (HyPar::boundary[HyPar::nBoundaryZones]) and calls the corresponding boundary condition function.

The variable flag indicates if the array x is the solution, or a delta-solution (from implicit time-integration methods).

Parameters

s	Object of type HyPar containing solver-related variables
m	Object of type MPIVariables containing MPI-related variables
x	The solution vector on which the boundary conditions are to be applied
xref	Reference solution vector, if needed
waqt	Current simulation time

Definition at line 28 of file ApplyBoundaryConditions.c.

{
  HyPar           *solver   = (HyPar*)          s;
  DomainBoundary  *boundary = (DomainBoundary*) solver->boundary;
  MPIVariables    *mpi      = (MPIVariables*)   m;
  int             nb        = solver->nBoundaryZones;

  int* dim_local;
#if defined(HAVE_CUDA)
  if (solver->use_gpu) {
    dim_local = solver->gpu_dim_local;
  } else {
#endif
    dim_local = solver->dim_local;
#if defined(HAVE_CUDA)
  }
#endif

  /* Apply domain boundary conditions to x */
  int n;
  for (n = 0; n < nb; n++) {
    boundary[n].BCFunctionU(&boundary[n],mpi,solver->ndims,solver->nvars,
                            dim_local,solver->ghosts,x,waqt);
  }

  return(0);
}