#include <ctype.h>
#include <string.h>
#include <stdlib.h>
#include <stdio.h>

#include "cctk.h"
#include "cctk_Parameters.h"
#include "cctk_Arguments.h"

#include "EOS_GIF.h"

void CCTK_FCALL CCTK_FNAME(EOS_GIF_Pressure) (const CCTK_INT* nelems,
                                              const CCTK_REAL* const rho,
                                              const CCTK_REAL* const eps,
                                              const CCTK_REAL* press);
void CCTK_FCALL CCTK_FNAME(EOS_GIF_IntEn)    (const CCTK_INT* nelems,
                                              const CCTK_REAL* const rho,
                                              const CCTK_REAL* const eps,
                                              const CCTK_REAL* inten);
void CCTK_FCALL CCTK_FNAME(EOS_GIF_DPDRho)   (const CCTK_INT* nelems,
                                              const CCTK_REAL* const rho,
                                              const CCTK_REAL* const eps,
                                              const CCTK_REAL* dpdrho);
void CCTK_FCALL CCTK_FNAME(EOS_GIF_DPDIE)    (const CCTK_INT* nelems,
                                              const CCTK_REAL* const rho,
                                              const CCTK_REAL* const eps,
                                              const CCTK_REAL* dpdie);
void CCTK_FCALL CCTK_FNAME(EOS_GIF_cs2)      (const CCTK_INT* nelems,
                                              const CCTK_REAL* const rho,
                                              const CCTK_REAL* const eps,
                                              const CCTK_REAL* cs2);
void CCTK_FCALL CCTK_FNAME(EOS_GIF_Temp)     (const CCTK_INT* nelems,
		                              const CCTK_REAL* const rho,
					      const CCTK_REAL* const eps,
					      const CCTK_REAL* temperature);
CCTK_INT EOS_GIF_SetArray(const CCTK_INT param_table,
                          const CCTK_INT n_elems,
                          const CCTK_POINTER* indep_vars,
                          const CCTK_INT* which_deps_to_set,
                          CCTK_POINTER* dep_vars)
{

  DECLARE_CCTK_PARAMETERS;
  
  CCTK_INT ierr = 0, var, dim;

/*   printf("A call to the ideal fluid EOS.\n" */
/*          "n_elems %d.\n", n_elems); */

  for (var = 0; var < N_DEPS; ++var)
  {
    if (which_deps_to_set[var])
    {
      switch (var)
      {
        case (0):
          CCTK_FNAME(EOS_GIF_Pressure) (&n_elems,
                                        (const CCTK_REAL*)indep_vars[0], 
                                        (const CCTK_REAL*)indep_vars[1], 
                                        (const CCTK_REAL*)dep_vars[0]);
          break;
        case (1):
          CCTK_FNAME(EOS_GIF_DPDRho)   (&n_elems,
                                        (const CCTK_REAL*)indep_vars[0], 
                                        (const CCTK_REAL*)indep_vars[1], 
                                        (const CCTK_REAL*)dep_vars[1]);
          break;
        case (2):
          CCTK_FNAME(EOS_GIF_DPDIE)    (&n_elems,
                                        (const CCTK_REAL*)indep_vars[0], 
                                        (const CCTK_REAL*)indep_vars[1], 
                                        (const CCTK_REAL*)dep_vars[2]);
          break;
        case (3):
          CCTK_FNAME(EOS_GIF_cs2)      (&n_elems,
                                        (const CCTK_REAL*)indep_vars[0],  
                                        (const CCTK_REAL*)indep_vars[1], 
                                        (const CCTK_REAL*)dep_vars[3]);
          break;
	case (4):
	  CCTK_FNAME(EOS_GIF_Temp)    (&n_elems,
			               (const CCTK_REAL*)indep_vars[0],
				       (const CCTK_REAL*)indep_vars[1],
				       (const CCTK_REAL*)dep_vars[4]);
	  break;
      }
    }
  }
  
  return ierr;
}