IBIdentifyMode.c

Go to the documentation of this file.

#include <string.h>
#include <immersedboundaries.h>

int IBIdentifyMode(
                    double *X,    
                    int    *dim,  
                    void   *ib    
                  )
{
  ImmersedBoundary  *IB = (ImmersedBoundary*) ib;
  Body3D            *body   = IB->body;

  double *x = X, *y = x + dim[0], *z = y + dim[1];

  double grid_xmin = x[0],
         grid_xmax = x[dim[0]-1],
         grid_ymin = y[0],
         grid_ymax = y[dim[1]-1],
         grid_zmin = z[0],
         grid_zmax = z[dim[2]-1];

  if      ( (grid_xmin > body->xmin) && (grid_xmax < body->xmax) ) strcpy(IB->mode,_IB_YZ_);
  else if ( (grid_ymin > body->ymin) && (grid_ymax < body->ymax) ) strcpy(IB->mode,_IB_XZ_);
  else if ( (grid_zmin > body->zmin) && (grid_zmax < body->zmax) ) strcpy(IB->mode,_IB_XY_);
  else                                                             strcpy(IB->mode,_IB_3D_);

  return(0);
}