flux_calc_kernel_c.c

/*Crown Copyright 2012 AWE.
*
* This file is part of CloverLeaf.
*
* CloverLeaf is free software: you can redistribute it and/or modify it under
* the terms of the GNU General Public License as published by the
* Free Software Foundation, either version 3 of the License, or (at your option)
* any later version.
*
* CloverLeaf is distributed in the hope that it will be useful, but
* WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
* FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
* details.
*
* You should have received a copy of the GNU General Public License along with
* CloverLeaf. If not, see http://www.gnu.org/licenses/. */

/**
 *  @brief C flux kernel.
 *  @author Wayne Gaudin
 *  @details The edge volume fluxes are calculated based on the velocity fields.
 */

#include <stdio.h>
#include <stdlib.h>
#include "ftocmacros.h"
#include <math.h>

void flux_calc_kernel_c_(int *xmin,int *xmax,int *ymin,int *ymax,
                          double *dbyt,
                          double *xarea,
                          double *yarea,
                          double *xvel0,
                          double *yvel0,
                          double *xvel1,
                          double *yvel1,
                          double *vol_flux_x,
                          double *vol_flux_y)
{
  int x_min=*xmin;
  int x_max=*xmax;
  int y_min=*ymin;
  int y_max=*ymax;
  double dt=*dbyt;

  int j,k;
  
 {

  for (k=y_min;k<=y_max;k++) {
#pragma ivdep
    for (j=x_min;j<=x_max+1;j++) {
      vol_flux_x[FTNREF2D(j  ,k  ,x_max+5,x_min-2,y_min-2)]=0.25*dt*xarea[FTNREF2D(j  ,k  ,x_max+5,x_min-2,y_min-2)]
                            *(xvel0[FTNREF2D(j  ,k  ,x_max+5,x_min-2,y_min-2)]
                             +xvel0[FTNREF2D(j  ,k+1,x_max+5,x_min-2,y_min-2)]
                             +xvel1[FTNREF2D(j  ,k  ,x_max+5,x_min-2,y_min-2)]
                             +xvel1[FTNREF2D(j  ,k+1,x_max+5,x_min-2,y_min-2)]);
    }
  }

  for (k=y_min;k<=y_max+1;k++) {
#pragma ivdep
    for (j=x_min;j<=x_max;j++) {
      vol_flux_y[FTNREF2D(j  ,k  ,x_max+4,x_min-2,y_min-2)]=0.25*dt*yarea[FTNREF2D(j  ,k  ,x_max+4,x_min-2,y_min-2)]
                            *(yvel0[FTNREF2D(j  ,k  ,x_max+5,x_min-2,y_min-2)]
                             +yvel0[FTNREF2D(j+1,k  ,x_max+5,x_min-2,y_min-2)]
                             +yvel1[FTNREF2D(j  ,k  ,x_max+5,x_min-2,y_min-2)]
                             +yvel1[FTNREF2D(j+1,k  ,x_max+5,x_min-2,y_min-2)]);
    }
  }

 }

}