1 files changed, 91 insertions, 0 deletions
diff --git a/src/EOS_GIF.c b/src/EOS_GIF.c
new file mode 100755
index 0000000..4947aa9
--- /dev/null
+++ b/src/EOS_GIF.c
@@ -0,0 +1,91 @@
+#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;
+}