/*@@
   @file      ApplyCartoon.c
   @date      6 Jan 2003
   @author    David Rideout
   @desc 
              Applies the Cartoon boundary condition to all variables 
	      selected for a bc.
   @enddesc 
   @version   $Header$
 @@*/

#include "cctk.h"
#include "cctk_Arguments.h"
#include "util_Table.h"
#include "Cartoon2D.h"
#include "Cartoon2D_tensors.h"
#include <stdlib.h>

static const char *rcsid = "$Header$";

CCTK_FILEVERSION(BetaThorns_Cartoon2D_ApplyCartoon_c);

/********************************************************************
 *********************     Local Data Types   ***********************
 ********************************************************************/
/* #define DEBUG 1 */
#define TENSORTYPE_BUFF_SIZE 7
#define PROLONG_BUFF_SIZE 1000

/********************************************************************
 ********************* Local Routine Prototypes *********************
 ********************************************************************/

/********************************************************************
 ***************** Aliased Routine Prototypes ***********************
 ********************************************************************/

/********************************************************************
 ***************** Scheduled Routine Prototypes *********************
 ********************************************************************/

void Cartoon_ApplyBoundaries(CCTK_ARGUMENTS);

/********************************************************************
 ********************* Other Routine Prototypes *********************
 ********************************************************************/

/********************************************************************
 *********************     Local Data   *****************************
 ********************************************************************/

/********************************************************************
 *********************     Aliased Routines   **********************
 ********************************************************************/
 
/********************************************************************
 *********************     Scheduled Routines   **********************
 ********************************************************************/

 /*@@
   @routine    Cartoon_ApplyBoundaries
   @date       6 Jan 2003
   @author     David Rideout
   @desc  
               This will apply the Cartoon boundary condition to all
               variables selected for any (physical) boundary condition.
   @enddesc 
   @calls      
   @history 
   @endhistory 
   @var        CCTK_ARGUMENTS
   @vdesc      Cactus argument list
   @vtype      CCTK_*
   @vio        in
   @endvar
   @returntype void
   @returndesc
   @endreturndesc
@@*/

void Cartoon_ApplyBoundaries(CCTK_ARGUMENTS) 
{

  DECLARE_CCTK_ARGUMENTS;
  int num_vars, err, i, gi, group_tags_table;
  CCTK_INT * vars;
  char tensortype[TENSORTYPE_BUFF_SIZE];
  char prolongtype[PROLONG_BUFF_SIZE];
  int prolongmethod;

  /* Check grid size */
  if(cctk_gsh[1] != 2*cctk_nghostzones[1]+1)
  {
    CCTK_WARN(0, "Grid size in y direction inappropriate.");
  }

  /* Allocate memory to hold selected bcs */
  num_vars = Boundary_SelectedGVs(cctkGH, 0, NULL, NULL, NULL, NULL, NULL);
#ifdef DEBUG
  printf("Cartoon_ApplyBoundaries: num_vars is %d\n",num_vars);
#endif
  vars = malloc(num_vars*sizeof *vars);

  /* get all selected vars */
  err = Boundary_SelectedGVs(cctkGH, num_vars, vars, NULL, NULL, NULL, NULL);
  if (err != num_vars)
  {
    CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING,
	       "Boundary_SelectedGVs returned %d selected variables, but %d "
	       "expected", err, num_vars);
  }

  /* Apply CartoonBC to each of them. */
  /* One should probably check to see that the entire group has been
   * selected, since Cartoon operates on an entire group at a time.
   * For now I'll just skip a variable if it is not a 'group leader'.
   */
  for (i=0; i<num_vars; ++i) {
#ifdef DEBUG
    printf("Cartoon_ApplyBoundaries: i=%d applying cartoon to vi %d\n",i,
	   vars[i]);
#endif
    gi = CCTK_GroupIndexFromVarI(vars[i]);
    if (gi<0) {
      CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING,
		 "Invalid variable index %d selected for a boundary "
		 "condition", gi);
    }
    if (vars[i] != CCTK_FirstVarIndexI(gi)) 
    {
      /* not a group leader -- skip to next var index */
      continue;
    }

    /* Here one should check that the entire group is registered,
     * using CCTK_NumVarsInGroupI.
     */

    /* Get table handle for group tags table */
    group_tags_table = CCTK_GroupTagsTableI(gi);
    if (group_tags_table<0) {
      CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING,
		 "Tags table for variable group %s is %d", CCTK_GroupName(gi),
		 group_tags_table);
    } 
    
    /* Get tensor type from group tags table */
    err = Util_TableGetString(group_tags_table, TENSORTYPE_BUFF_SIZE,
			      tensortype, "tensortypealias");
    if (err<0) 
    {
      CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING,
		 "Error in TAGS table for variable group %s",
		 CCTK_GroupName(gi));
    }
#ifdef DEBUG
    printf("Cartoon_ApplyBoundaries: tensor type is %s\n",tensortype);
#endif

    /* Get prolongation type from group tags table */

    err = Util_TableGetString(group_tags_table, PROLONG_BUFF_SIZE,
                              prolongtype, "Prolongation");

    if (err == UTIL_ERROR_TABLE_NO_SUCH_KEY) 
    {
      /* Use the default */

      prolongmethod = PROLONG_LAGRANGE;
    }
    else if (err < 0)
    {
      CCTK_VWarn(0, __LINE__, __FILE__, CCTK_THORNSTRING,
		 "Error (%d) in TAGS table for variable group %s", err,
		 CCTK_GroupName(gi));
    }
    else
    {
      if (CCTK_Equals(prolongtype, "None")) 
      {
        prolongmethod = PROLONG_NONE; /* But why? */
      } 
      else if (CCTK_Equals(prolongtype, "Lagrange")) 
      {
        prolongmethod = PROLONG_LAGRANGE;
      } 
      else if (CCTK_Equals(prolongtype, "TVD")) 
      {
        prolongmethod = PROLONG_ENO;
      } 
      else if (CCTK_Equals(prolongtype, "ENO")) 
      {
        prolongmethod = PROLONG_ENO;
      } 
      else 
      {
        CCTK_VWarn(0, __LINE__, __FILE__, CCTK_THORNSTRING,
                   "Error in TAGS table for variable group %s",
                   CCTK_GroupName(gi));
      }
    }    
    
    /* Here one should check that the group's size is correct for the
     * specified tensortype.
     */

    /* Call BndCartoon2DVI, passing the appropriate tensor type integer 
       macro */
    if (CCTK_Equals(tensortype, "scalar"))
    {
      BndCartoon2DVI(cctkGH, TENSORTYPE_SCALAR, prolongmethod, vars[i]);
    } else if ((CCTK_Equals(tensortype, "u")) || (CCTK_Equals(tensortype, "d")))
    {
      BndCartoon2DVI(cctkGH, TENSORTYPE_U, prolongmethod, vars[i]);
    } else if (CCTK_Equals(tensortype, "dd_sym"))
    {
      BndCartoon2DVI(cctkGH, TENSORTYPE_DDSYM, prolongmethod, vars[i]);
    } else 
    {
      CCTK_VWarn(0, __LINE__, __FILE__, CCTK_THORNSTRING,
		 "invalid tensor type for group %s", CCTK_GroupName(gi));
    }  
  }

  /* Free data */
  free(vars);
}